結果
| 問題 |
No.3004 ヤング図形
|
| コンテスト | |
| ユーザー |
umezo
|
| 提出日時 | 2025-01-18 02:35:45 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 1,101 ms / 4,000 ms |
| コード長 | 1,058 bytes |
| コンパイル時間 | 5,252 ms |
| コンパイル使用メモリ | 275,928 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2025-01-18 02:36:22 |
| 合計ジャッジ時間 | 23,964 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 25 |
ソースコード
#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;
}
int kai(int n){
if(n<=1) return 1;
ll t=kai(n/2);
if(n%2==0) return t*t%MOD*nCr(n,n/2)%MOD;
return 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];
if(L[i]>1){
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;
}
umezo