結果
| 問題 |
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;
}