結果
問題 |
No.3004 ヤング図形
|
ユーザー |
|
提出日時 | 2025-01-31 19:22:52 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 634 ms / 4,000 ms |
コード長 | 1,354 bytes |
コンパイル時間 | 3,919 ms |
コンパイル使用メモリ | 208,652 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2025-01-31 19:23:09 |
合計ジャッジ時間 | 15,752 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll = long long; #define rep(i,n) for(int i=0;i<n;i++) #define all(A) A.begin(),A.end() #include<atcoder/modint> using namespace atcoder; using mint = modint998244353; vector<mint> fact, factinv, inv; const ll mod = 998244353; void prenCkModp(ll n) { fact.resize(n + 5); factinv.resize(n + 5); inv.resize(n + 5); fact[0] = fact[1] = 1; factinv[0] = factinv[1] = 1; inv[1] = 1; for (ll i = 2; i < n + 5; i++) { fact[i] = (fact[i - 1] * i); inv[i] = mod - (inv[mod % i] * (mod / i)); factinv[i] = (factinv[i - 1] * inv[i]); } } mint nCk(ll n, ll k) { if (n < k||k<0) return 0; return fact[n] * (factinv[k] * factinv[n - k]); } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); ll K; cin>>K; vector<ll> L(K),M(K); ll S=0; rep(i,K){ cin>>L[i]>>M[i]; S+=L[i]*M[i]; } int l=0; mint an=1,bn=1; map<ll,ll> U; ll m=0; rep(k,K){ U[M[k]]++; m=max(m,M[k]); } for(int i=1;i<=S;i++){ an*=i; if(i<=m){ if(U.count(i))bn*=mint(an).pow(U[i]); } if(l<K&&L[l]==i){ mint p=mint(an).pow(M[l]); l++; bn*=p; } } an*=bn.inv(); cout<<an.val()<<endl; }