結果
問題 |
No.1191 数え上げを愛したい(数列編)
|
ユーザー |
![]() |
提出日時 | 2020-08-28 20:02:32 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 17 ms / 2,000 ms |
コード長 | 899 bytes |
コンパイル時間 | 2,335 ms |
コンパイル使用メモリ | 192,652 KB |
最終ジャッジ日時 | 2025-01-13 16:30:28 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 26 |
ソースコード
#define rep(i, n) for (int i = 0; i < (int)(n); i++) #define ALL(v) v.begin(), v.end() typedef long long ll; #include <bits/stdc++.h> using namespace std; const ll MAX = 310000; const ll MOD = 998244353; ll fac[MAX], finv[MAX], inv[MAX]; void init(){ 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; } } ll nCr(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; } ll nPr(int n,int k){ if(n<k) return 0; if(n<0 || k<0) return 0; return fac[n]*finv[n-k]%MOD; } int main(){ ll n,m,a,b; cin>>n>>m>>a>>b; init(); ll ans=0; for(int i=1;i<=m;i++){ if(min(b,m-i)-(n-1)*(a-1)<0) break; ans=(ans+nCr(min(b,m-i)-(n-1)*(a-1),n-1))%MOD; } cout<<(ans*nPr(n,n))%MOD<<endl; return 0; }