結果

問題 No.3004 ヤング図形
ユーザー maeshun
提出日時 2025-01-18 14:48:56
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 3,343 bytes
コンパイル時間 4,355 ms
コンパイル使用メモリ 255,396 KB
実行使用メモリ 814,716 KB
最終ジャッジ日時 2025-01-18 14:49:42
合計ジャッジ時間 36,706 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 20 TLE * 3 MLE * 2
権限があれば一括ダウンロードができます

ソースコード

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

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
#define rep(i, n) for(int i=0;i<(n);++i)
#define rep1(i, n) for(int i=1;i<=(n);i++)
#define ll long long
using mint = modint998244353;
using P = pair<ll,ll>;
using lb = long double;
using T = tuple<ll, ll, ll>;
#ifdef LOCAL
# include <debug_print.hpp>
# define dbg(...) debug_print::multi_print(#__VA_ARGS__, __VA_ARGS__)
#else
# define dbg(...) (static_cast<void>(0))
#endif
// combination mod prime
// https://youtu.be/8uowVvQ_-Mo?t=6002
// https://youtu.be/Tgd_zLfRZOQ?t=9928
struct modinv {
int n; vector<mint> d;
modinv(): n(2), d({0,1}) {}
mint operator()(int i) {
while (n <= i) d.push_back(-d[mint::mod()%n]*(mint::mod()/n)), ++n;
return d[i];
}
mint operator[](int i) const { return d[i];}
} invs;
struct modfact {
int n; vector<mint> d;
modfact(): n(2), d({1,1}) {}
mint operator()(int i) {
while (n <= i) d.push_back(d.back()*n), ++n;
return d[i];
}
mint operator[](int i) const { return d[i];}
} facts;
struct modfactinv {
int n; vector<mint> d;
modfactinv(): n(2), d({1,1}) {}
mint operator()(int i) {
while (n <= i) d.push_back(d.back()*invs(n)), ++n;
return d[i];
}
mint operator[](int i) const { return d[i];}
} ifacts;
mint comb(int n, int k) {
if (n < k || k < 0) return 0;
return facts(n)*ifacts(k)*ifacts(n-k);
}
uint64_t combinations2(uint64_t n, uint64_t k) {
uint64_t r = 1;
for (uint64_t d = 1; d <= k; ++d) {
r *= n--;
r /= d;
}
return r;
}
// vector<vector<mint>> comb(n+1, vector<mint>(n+1));
// rep(i,n+1) comb[i][0] = 1;
// for(int i=1;i<=n;i++){
// for(int j=1;j<=n;j++){
// comb[i][j] = comb[i-1][j] + comb[i-1][j-1];
// }
// }
//Lucasmod
/* ComnCk % p
: O(p^2)
nCk(n,k): nCk % p O(log n)
*/
// struct Comb {
// vector<vector<long long>> com; //
// long long p; // p
// Comb(long long _p) : p(_p) {
// init(p);
// }
// void init(long long p) { //
// com.assign(p, vector<long long>(p));
// com[0][0] = 1;
// for (int i = 1; i < p; i++) {
// com[i][0] = 1;
// for (int j = i; j > 0; j--) {
// com[i][j] = (com[i - 1][j - 1] + com[i - 1][j]) % p;
// }
// }
// }
// long long nCk(long long n, long long k) {
// long long ret = 1;
// while (n > 0) { //
// int ni = n % p;
// int ki = k % p;
// ret *= com[ni][ki];
// ret %= p;
// n /= p;
// k /= p;
// }
// return ret;
// }
// };
int main()
{
int k;
cin >> k;
vector<int> l(k), m(k);
rep(i,k) cin >> l[i] >> m[i];
mint ans = 1;
int sum = 0;
rep(i,k) {
mint now = 1;
int cnt = l[i];
for(int j=1;j<m[i];j++) {
now = comb(cnt+l[i]-1,l[i]-1)*now;
cnt += l[i];
}
ans *= now;
ans *= ifacts(cnt);
sum += cnt;
}
rep(i,sum) ans *= i+1;
cout << ans.val() << endl;
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0