結果
| 問題 |
No.3160 Party Game
|
| コンテスト | |
| ユーザー |
umezo
|
| 提出日時 | 2025-05-24 07:59:50 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,125 bytes |
| コンパイル時間 | 2,585 ms |
| コンパイル使用メモリ | 276,456 KB |
| 実行使用メモリ | 50,548 KB |
| 最終ジャッジ日時 | 2025-05-27 22:06:15 |
| 合計ジャッジ時間 | 5,593 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 37 WA * 1 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define ALL(v) v.begin(),v.end()
typedef long long ll;
template <class T> using V=vector<T>;
template <class T> using VV=V<V<T>>;
const int MAX=2000200;
const int 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 modpow(ll x,ll n){
ll ans=1;
while(n){
if(n&1) ans=ans*x%MOD;
x=x*x%MOD;
n/=2;
}
return ans;
}
int main(){
ios::sync_with_stdio(false);
std::cin.tie(nullptr);
init();
ll n,m;
cin>>n>>m;
if(m==1){
cout<<0<<endl;
return 0;
}
ll sum=nCr(m+n,n)-n;
ll a=modpow(sum,MOD-2);
ll ans=0,psum=0;
for(ll i=m;i>0;i--){
if(m-n*i<0) continue;
ll p=nCr(m-n*i+n,n)*a%MOD;
ans=(ans+(p-psum+MOD)%MOD*i%MOD)%MOD;
psum=p;
}
cout<<ans<<endl;
return 0;
}
umezo