結果
問題 | No.802 だいたい等差数列 |
ユーザー |
|
提出日時 | 2024-03-30 15:39:29 |
言語 | C++17(clang) (17.0.6 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 974 bytes |
コンパイル時間 | 7,331 ms |
コンパイル使用メモリ | 162,716 KB |
実行使用メモリ | 261,248 KB |
最終ジャッジ日時 | 2024-09-30 17:31:23 |
合計ジャッジ時間 | 19,585 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 WA * 3 |
other | AC * 5 WA * 25 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/modint>using namespace atcoder;using namespace std;using mint = modint1000000007;using ll = long long;const int MAX = 11000000;const int MOD = 1000000007;long long fac[MAX], finv[MAX], inv[MAX];void COMinit(){fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for(int i = 2; i < MAX; i++){fac[i] = fac[i-1] * i % MOD;inv[i] = MOD - inv[MOD%i] * (MOD/i) % MOD;finv[i] = finv[i-1] * inv[i] % MOD;}}long long COM(long long n, long long k){if(n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n-k] % MOD) % MOD;}int main() {ll N, M, D1, D2;cin >> N >> M >> D1 >> D2;ll D = D2 - D1;ll x =M - (D1 * (N -1)) - 1;if (x < 0) {cout << 0 << endl;return 0;}COMinit();mint ans = 0;for (ll p = 0;p < N;p++) {mint tmp = mint(COM(x + N - (p * (D + 1)), N)) * mint(COM(N - 1, p));}cout << ans.val() << endl;}