結果
問題 |
No.1275 綺麗な式
|
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
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)