結果
問題 | No.802 だいたい等差数列 |
ユーザー |
![]() |
提出日時 | 2019-03-18 22:08:33 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,620 ms / 2,000 ms |
コード長 | 1,022 bytes |
コンパイル時間 | 594 ms |
コンパイル使用メモリ | 65,024 KB |
実行使用メモリ | 34,688 KB |
最終ジャッジ日時 | 2024-07-19 06:20:31 |
合計ジャッジ時間 | 57,524 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include <iostream>using namespace std;typedef long long ll;ll inf = 1e9+7;ll func[2000010] = {0};ll inv[2000010] = {0};ll mult(ll n, ll m){if(m==1) return n%inf;else if(m%2==0){ll t = mult(n,m/2);return (t*t)%inf;}else{ll t = mult(n,m-1);return (t*n)%inf;}}void factorial(ll N){for(ll i=0;i<=N;i++){if(i==0){func[i] = 1;inv[i] = 1;}else{func[i] = (i*func[i-1])%inf;inv[i] = mult(func[i],inf-2);}}}ll comb(ll n,ll k){if(n<k) return 0;if(n<0 || k<0) return 1;else return (((func[n]*inv[k])%inf)*inv[n-k])%inf;}ll hcomb(ll n,ll k){if(n==0 && k==0) return 1;return comb(n+k-1,k);}ll N,M,D1,D2;int main(){cin >> N >> M >> D1 >> D2;factorial(2000000);ll ans = 0;ll x = M-1-D1*(N-1);if(x<0){cout << 0 << endl;return 0;}for(ll i=0;i<=N-1;i++){ll r = (comb(N-1,i)*comb(x-i*(D2-D1+1)+N,N))%inf;if(i%2==0) (ans += r)%=inf;else (ans += -r+inf)%=inf;}cout << ans << endl;return 0;}