結果
| 問題 |
No.2528 pop_(backfront or not)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-11-03 22:50:18 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 280 ms / 2,000 ms |
| コード長 | 996 bytes |
| コンパイル時間 | 3,797 ms |
| コンパイル使用メモリ | 255,132 KB |
| 最終ジャッジ日時 | 2025-02-17 18:40:01 |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 19 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using ll = long long;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
using mint = atcoder::modint998244353;
const int MOD=998244353;
const int MAX_N=2e6;
vector<ll> fac(MAX_N+1,1);
vector<ll> finv(MAX_N+1,1);
vector<ll> inv(MAX_N+1,1);
void comb_setup(){
for(int i=2;i<=MAX_N;i++){
fac.at(i)=(fac.at(i-1)*i)%MOD;
inv.at(i)=MOD-(inv.at(MOD%i)*(MOD/i))%MOD;
finv.at(i)=(finv.at(i-1)*inv.at(i))%MOD;
}
}
mint comb(int n,int k){
if(n<k) return 0;
return fac.at(n)*(finv.at(n-k)*finv.at(k)%MOD)%MOD;
}
vector<vector<mint>> mm(5000,vector<mint>(5000,-1));
mint dv(int l,int r){
if(l==0&&r==0) return 1;
if(l==0||r==0) return 0;
if(mm.at(l).at(r)!=-1) return mm.at(l).at(r);
mint rt=0;
rt+=dv(l-1,r-1)*((l-1)*(r-1)+1);
if(l>2) rt+=dv(l-2,r)*comb(l-1,2);
if(r>2) rt+=dv(l,r-2)*comb(r-1,2);
return mm.at(l).at(r)=rt;
}
int main(){
comb_setup();
int n;
cin>>n;
rep(i,2*n+1){
cout<<dv(i,2*n-i).val()<<"\n";
}
}