# AtCoder Beginner Contest 145

URL: https://atcoder.jp/contests/abc145 (opens new window)

# D

x,y = gets.chomp.split(" ").map(&:to_i)

# パスカルの三角形を書くとわかりやすい
#        - (1, 2) - (2, 4)
# (0, 0)          - (3, 3)
#        - (2, 1) - (4, 3)

# n_C_k において, 下のような式が成り立つ
n = (x + y) / 3
k = y - n

if (x + y) % 3 != 0
  puts 0
  exit
end

# これより下は組み合わせを高速に解くための関数群 -----

MOD = 10 ** 9 + 7
# 繰り返し二乗法(x.pow(n, mod)でもOK)
def pow(x, n, mod = MOD)
  ans = 1
  # 下位ビットから繰り返し計算を行う
  while n > 0 do
    ans = ans * x % mod if n.odd? # 下位ビットが1かどうか
    x = x * x % mod
    n >>= 1 # 右ビットシフト
  end
  return ans
end

# (mod p)下での逆元
# フェルマーの小定理より a^(p-1) = 1 (mod p) (aは整数, pは素数, aとpは互いに素)
# a * a^(p-2) = 1 なので, aの逆元が a^(p-2) となる
def inv(a, mod = MOD)
  pow(a, mod - 2, mod)
end

# 階乗
# 再利用できるようグローバル変数化
$fact = [1] # 0! = 1
def fact(n, mod = MOD)
  # 計算済ならそれを返す
  return $fact[n] unless $fact[n].nil?

  # 階乗を求めながら$factへ格納
  f = $fact[-1]
  ($fact.size).upto(n) do |i|
    f = (f * i) % mod
    $fact << f
  end
  return $fact[n]
end

# 順列
def perm(n, r, mod = MOD)
  # n*(n-1)*(n-2)*...*(n-r+1)
  return n.downto(n - r + 1).reduce(1) {|ans, x| (ans * x) % mod }
end

# 組み合わせ
def comb(n, r, mod = MOD)
  return 0 unless 0 <= r && r <= n
  return 1 if r == 0 || r == n
  # nCk = n!/(n-k)! * 1/k!
  return perm(n, r, mod) * inv(fact(r, mod) % mod, mod) % mod
end

# 関数群ここまで -----

puts comb(n, k)
Last Updated: 2020/08/09 18:33