結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0