結果
問題 | No.2468 Mercurialist |
ユーザー |
|
提出日時 | 2023-09-23 23:59:23 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 40 ms / 2,000 ms |
コード長 | 999 bytes |
コンパイル時間 | 2,181 ms |
コンパイル使用メモリ | 166,940 KB |
実行使用メモリ | 8,192 KB |
最終ジャッジ日時 | 2024-07-16 23:49:40 |
合計ジャッジ時間 | 3,616 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include<bits/stdc++.h> using namespace std; typedef long long ll; #define fi first #define se second const ll mod=998244353; const int N=3e5+5; const int iu=3e5; ll x,y,z,k; ll f[N],inf[N]; ll pw(ll x,ll y){ if(y==0) return 1; if(y%2) return x*pw(x,y-1)%mod; ll res=pw(x,y/2); return res*res%mod; } ll C(ll x,ll y){ return f[x]*inf[x-y]%mod*inf[y]%mod; } int main(){ ios::sync_with_stdio(false);cin.tie(0); cin >> x >> y >> z >> k; f[0]=1; for(int i=1; i<=iu ;i++) f[i]=f[i-1]*i%mod; inf[iu]=pw(f[iu],mod-2); for(int i=iu; i>=1 ;i--) inf[i-1]=inf[i]*i%mod; ll n=x+y+z; ll tot=C(x+y+z,x)*C(y+z,y)%mod*f[y]%mod; ll cum=0; for(int i=1; i<=n-x+1 ;i++){ ll frog=C(n-i,x-1); ll ways=0; int st=i-k; if(st<=0) ways=f[y+z]*inf[z]%mod; else if(st<=y) ways=pw(n-i-x+k+1,st)*f[y+z-st]%mod*inf[z]%mod; else ways=pw(n-i-x+k+1,y); cum=(cum+frog*ways)%mod; //cout << i << ' ' << frog << ' ' << ways << ' ' << n-i-x+k+1 << endl;; } cum=cum*pw(tot,mod-2)%mod; cout << cum << '\n'; }