結果

問題 No.1275 綺麗な式
ユーザー ygd.
提出日時 2020-10-30 22:37:31
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 739 bytes
コンパイル時間 155 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 310,280 KB
最終ジャッジ日時 2024-07-22 02:06:04
合計ジャッジ時間 61,221 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 12 RE * 48
権限があれば一括ダウンロードができます

ソースコード

diff #

def cmb(n, r, p):
  if (r < 0) or (n < r):
    return 0
  r = min(r, n - r)
  return fac[n]*finv[r]*finv[n-r]%p

def perm(n,r,p):
  if (r < 0) or (n < r):
    return 0
  return fac[n]*finv[n-r]%p

N = pow(10,7)+10
MOD = pow(10,9)+7

fac = [-1]*(N+1); fac[0] = 1; fac[1] = 1 #階乗
finv = [-1]*(N+1); finv[0] = 1; finv[1] = 1 #階乗の逆元
inv = [-1]*(N+1); inv[0] = 0; inv[1] = 1 #逆元
for i in range(2,N+1):
  fac[i] = fac[i-1]*i%MOD
  inv[i] = MOD - inv[MOD%i]*(MOD//i)%MOD
  finv[i] = finv[i-1]*inv[i]%MOD

a,b = map(int,input().split())
n = int(input())
ans = 0
m = n//2
#print(m)
for j in range(m+1):
  #print(n-2*j)
  temp = 2*cmb(n,2*j,MOD)*pow(a,n-2*j,MOD)*pow(b,j,MOD)%MOD
  #print(temp)
  ans = (ans+temp)%MOD
print(ans)
  
0