結果
問題 | No.802 だいたい等差数列 |
ユーザー |
![]() |
提出日時 | 2019-03-18 18:13:35 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 130 ms / 2,000 ms |
コード長 | 1,271 bytes |
コンパイル時間 | 889 ms |
コンパイル使用メモリ | 101,584 KB |
実行使用メモリ | 50,432 KB |
最終ジャッジ日時 | 2024-07-18 13:01:16 |
合計ジャッジ時間 | 4,937 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include <algorithm>#include <iostream>#include <iomanip>#include <cstring>#include <string>#include <vector>#include <random>#include <queue>#include <cmath>#include <stack>#include <set>#include <map>typedef long long ll;using namespace std;const int MAX_N = 2000000;const ll MOD = 1000000007LL;ll inv[MAX_N], fac[MAX_N], fiv[MAX_N];void COMinit() {inv[1] = fac[1] = fiv[1] = inv[0] = fac[0] = fiv[0] = 1;for (ll i = 2; i < MAX_N; i++) {fac[i] = fac[i - 1] * i % MOD; // n!inv[i] = inv[MOD % i] * (MOD - MOD / i) % MOD; // n^-1fiv[i] = fiv[i - 1] * inv[i] % MOD; // (n!)^-1}}ll com(ll n, ll k) {if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * fiv[k] % MOD * fiv[n-k] % MOD;}int main() {cin.sync_with_stdio(false);cin.tie(0);ll n, m, d1, d2;cin >> n >> m >> d1 >> d2;ll d = d2 - d1;ll x = m - 1 - d1 * (n - 1);if (x < 0) {cout << 0 << '\n';return 0;}COMinit();ll ans = 0;for (ll i = 0; i < n; i++) {ll tmp = com(n-1, i) * com(x - i*(d+1) + n, n) % MOD;if (i % 2 == 0) ans += tmp;else ans += MOD - tmp;ans %= MOD;}cout << ans << "\n";return 0;}