結果
| 問題 |
No.802 だいたい等差数列
|
| コンテスト | |
| ユーザー |
ayaoni
|
| 提出日時 | 2020-11-26 10:36:34 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 131 ms / 2,000 ms |
| コード長 | 617 bytes |
| コンパイル時間 | 314 ms |
| コンパイル使用メモリ | 82,108 KB |
| 実行使用メモリ | 97,440 KB |
| 最終ジャッジ日時 | 2024-07-23 20:44:03 |
| 合計ジャッジ時間 | 4,485 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 30 |
ソースコード
import sys
def MI(): return map(int,sys.stdin.readline().rstrip().split())
N,M,D1,D2 = MI()
mod = 10**9+7
D = D2-D1
MM = M-1-(N-1)*D1
fac = [0]*(N+M)
fac[0] = 1
for i in range(1,N+M):
fac[i] = (i*fac[i-1]) % mod
fac_inv = [0]*(N+M)
a = pow(fac[-1],mod-2,mod)
for i in range(N+M-1,-1,-1):
fac_inv[i] = a
a *= i
a %= mod
def nCr(n,r):
if n < r:
return 0
return (fac[n]*fac_inv[r]*fac_inv[n-r]) % mod
ans = 0
for i in range(N):
X = MM-i*(D+1)
a = nCr(N-1,i)*nCr(X+N,N)
a %= mod
if i % 2 == 0:
ans += a
else:
ans -= a
ans %= mod
print(ans)
ayaoni