結果
問題 | No.802 だいたい等差数列 |
ユーザー |
|
提出日時 | 2023-01-25 23:51:31 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 125 ms / 2,000 ms |
コード長 | 1,224 bytes |
コンパイル時間 | 2,061 ms |
コンパイル使用メモリ | 194,248 KB |
最終ジャッジ日時 | 2025-02-10 06:51:49 |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define fi first#define se second#define pb push_backusing vi = vector <int>;using ll = long long;using pii = pair <int, int>;//~ const ll mod = 998244353;const ll mod = 1e9 + 7;ll qpow(ll a, ll b, ll m = mod) { ll r = 1, t = a;for(; b; b /= 2) { if(b & 1) r = r * t % m; t = t * t % m; } return r; }const int N = 2e6 + 11;ll inv[N], fac[N], ifac[N];ll binom(int n, int m) {if(n < m || m < 0) return 0; // n >= 0return fac[n] * ifac[m] % mod * ifac[n - m] % mod;}ll f(int m, int n) {return binom(m + n, n);}int main() {inv[1] = fac[0] = ifac[0] = 1;for(int i = 2; i < N; i ++)inv[i] = (mod - mod / i) * inv[mod % i] % mod;for(int i = 1; i < N; i ++) {fac[i] = i * fac[i - 1] % mod;ifac[i] = inv[i] * ifac[i - 1] % mod;}ios :: sync_with_stdio(false);int n; cin >> n;ll m, d1, d2; cin >> m >> d1 >> d2;m -= (n - 1) * d1; m --;int d = d2 - d1;ll ans = 0;if(m >= 0) {ans = f(m, n);for(int i = 1; i < n && m - i * (d + 1) >= 0; i ++) {int sgn = ((i & 1) ? -1 : 1);ans += sgn * binom(n - 1, i) * f(m - i * (d + 1), n);ans %= mod;}}ans %= mod;if(ans < 0) ans += mod;cout << ans << '\n';}