結果
問題 | No.3004 ヤング図形 |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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 longusing 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=9928struct 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];// }// }//Lucasの定理(素数mod)/* Com:nCk % 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;}