結果
| 問題 |
No.1191 数え上げを愛したい(数列編)
|
| コンテスト | |
| ユーザー |
umezo
|
| 提出日時 | 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;
}
umezo