結果
| 問題 |
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 |
ソースコード
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)
ygd.