結果
問題 |
No.3004 ヤング図形
|
ユーザー |
![]() |
提出日時 | 2025-01-18 02:27:50 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,092 bytes |
コンパイル時間 | 3,510 ms |
コンパイル使用メモリ | 282,848 KB |
実行使用メモリ | 13,636 KB |
最終ジャッジ日時 | 2025-01-18 02:28:32 |
合計ジャッジ時間 | 38,714 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 TLE * 4 |
ソースコード
#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 int MOD=998244353; 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; } ll nCr(ll n,ll r){ if(n<r) return 0; if(n<0 || r<0) return 0; ll x=1,y=1; for(ll i=1;i<=r;i++){ x=x*(n+1-i)%MOD; y=y*i%MOD; } return x*modpow(y,MOD-2)%MOD; } map<ll,ll> m; ll kai(ll n){ if(m.contains(n)) return m[n]; if(n<=1) return m[n]=1; ll t=kai(n/2); if(n%2==0) return m[n]=t*t%MOD*nCr(n,n/2)%MOD; return m[n]=t*t%MOD*nCr(n,n/2)%MOD*((n+1)/2)%MOD; } int main(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); int k; cin>>k; vector<ll> L(k),M(k); ll sum=0; rep(i,k){ cin>>L[i]>>M[i]; sum+=L[i]*M[i]; } ll ans=1; rep(i,k){ ll tmp=L[i]*M[i]; rep(j,M[i]) ans=ans*nCr(tmp-j*L[i],L[i])%MOD; ans=ans*modpow(kai(M[i]),MOD-2)%MOD; ans=ans*nCr(sum,tmp)%MOD; sum-=tmp; } cout<<ans<<endl; return 0; }