結果
| 問題 |
No.1127 変形パスカルの三角形
|
| ユーザー |
|
| 提出日時 | 2022-05-15 09:24:29 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 62 ms / 1,500 ms |
| コード長 | 740 bytes |
| コンパイル時間 | 195 ms |
| コンパイル使用メモリ | 82,208 KB |
| 実行使用メモリ | 67,252 KB |
| 最終ジャッジ日時 | 2024-07-23 21:27:07 |
| 合計ジャッジ時間 | 3,214 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 30 |
ソースコード
P = 10 ** 9 + 7
a,b = map(int,input().split())
a %= P
b %= P
n,k = map(int,input().split())
C = n + 5
fact = [1] * C
fact_inv = [1] * C
for i in range(2,C):
fact[i] = fact[i-1] * i % P
fact_inv[-1] = pow(fact[-1],P - 2,P)
for i in range(C - 2,1,-1):
fact_inv[i] = fact_inv[i + 1] * (i + 1) % P
def calc(n,k):
return fact[n] * fact_inv[k] % P * fact_inv[n - k] % P
if k == 1:
ans = a
elif k == n + 1:
ans = b
else:
ans = calc(n - 1,k - 1) * a + calc(n - 1,n - k + 1) * b
print(ans % P)
ans = 0
for k in range(1,n + 1 + 1):
if k == 1:
tmp = a
elif k == n + 1:
tmp = b
else:
tmp = calc(n - 1,k - 1) * a + calc(n - 1,n - k + 1) * b
ans += (tmp % P) ** 2
ans %= P
print(ans)