結果
| 問題 |
No.802 だいたい等差数列
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-10-26 17:32:12 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 148 ms / 2,000 ms |
| コード長 | 746 bytes |
| コンパイル時間 | 527 ms |
| コンパイル使用メモリ | 82,200 KB |
| 実行使用メモリ | 110,260 KB |
| 最終ジャッジ日時 | 2024-10-26 17:32:18 |
| 合計ジャッジ時間 | 5,947 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 30 |
ソースコード
N, M, D1, D2 = map(int, input().split())
MOD = 10**9 + 7
num = 2*10**6
frac = [1] * num
inv_frac = [1] * num
for i in range(2, num):
frac[i] = (frac[i-1] * i) % MOD
inv_frac[-1] = pow(frac[-1], MOD-2, MOD)
for i in range(num-1, 0, -1):
inv_frac[i-1] = (inv_frac[i] * i) % MOD
def comb(n, r):
if r < 0 or r > n:
return 0
return (frac[n] * (inv_frac[n-r] * inv_frac[r])%MOD) % MOD
def calc(n, m):
# 合計 m個を n人で分ける
# n人それぞれの最小個数は0個
# {m+1}_H_{n-1} = {m+n-1}_C_{n-1}
return comb(n+m-1, n-1)
M = M - 1 - D1*(N-1)
ans = calc(N+1, M)
for k in range(1, N):
m = M - (D2-D1+1)*k
ans += (-1)**(k%2) * comb(N-1, k) * calc(N+1, m)
ans %= MOD
print(ans)