結果
| 問題 |
No.2528 pop_(backfront or not)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-11-03 22:00:50 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 65 ms / 2,000 ms |
| コード長 | 875 bytes |
| コンパイル時間 | 4,222 ms |
| コンパイル使用メモリ | 255,904 KB |
| 最終ジャッジ日時 | 2025-02-17 18:15:48 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 19 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll=long long;
using mint=modint998244353;
int main()
{
int n;
cin>>n;
vector<vector<mint>> dp(n+1,vector<mint>(2*n+1));
dp[0][0]=1;
vector<mint> fac(4001),invfac(4001);
fac[0]=1;
for(int i=1;i<=4000;i++)fac[i]=fac[i-1]*i;
for(int i=0;i<=4000;i++)invfac[i]=fac[i].inv();
for(int i=1;i<=n;i++){
for(int j=1;j<=i;j++){
if(2*i-j>=3)dp[i][j]+=dp[i-1][j]*fac[2*i-j-1]*invfac[2*i-j-3]*invfac[2];
if(j>=3)dp[i][j]+=dp[i-1][j-2]*fac[j-1]*invfac[j-3]*invfac[2];
if(j>=2&&2*i-j>=2)dp[i][j]+=dp[i-1][j-1]*(j-1)*(2*i-j-1);
dp[i][j]+=dp[i-1][j-1];
dp[i][2*i-j]=dp[i][j];
}
}
for(int i=0;i<=2*n;i++){
cout<<dp[n][min(i,2*n-i)].val()<<endl;
}
return 0;
}