結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

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)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0