結果
問題 | No.802 だいたい等差数列 |
ユーザー |
![]() |
提出日時 | 2020-08-06 17:53:12 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,492 bytes |
コンパイル時間 | 1,238 ms |
コンパイル使用メモリ | 115,932 KB |
最終ジャッジ日時 | 2025-01-12 15:25:29 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 14 WA * 5 TLE * 11 |
ソースコード
#include<iostream>#include<string>#include<algorithm>#include<vector>#include<iomanip>#include<math.h>#include<complex>#include<queue>#include<deque>#include<stack>#include<map>#include<set>#include<bitset>#include<functional>#include<assert.h>#include<numeric>using namespace std;#define REP(i,m,n) for(int i=(int)(m) ; i < (int) (n) ; ++i )#define rep(i,n) REP(i,0,n)using ll = long long;const int inf=1e9+7;const ll longinf=1LL<<60 ;const ll mod=1e9+7 ;#define pai 3.141592653589793238462643383279const int MAX = 10000000;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(int n, int 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(){int n, m, d1, d2; cin >> n >> m >> d1 >> d2;int zenbu = m-1-(n-1)*d1;if(zenbu<0){cout << 0 << endl; return 0;}COMinit(); ll ans = 0;rep(i, n){if(i%2){ans-=COM(zenbu+n-(d2-d1+1)*i, n)*COM(n-1, i)%mod; ans%=mod;}else{ans+=COM(zenbu+n-(d2-d1+1)*i, n)*COM(n-1, i); ans%=mod;}//cout << ans << endl;}cout << ans << endl;}