結果
| 問題 |
No.3043 括弧列の数え上げ
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2025-02-28 22:33:40 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 641 ms / 2,000 ms |
| コード長 | 1,610 bytes |
| コンパイル時間 | 6,202 ms |
| コンパイル使用メモリ | 334,236 KB |
| 実行使用メモリ | 120,924 KB |
| 最終ジャッジ日時 | 2025-02-28 22:34:09 |
| 合計ジャッジ時間 | 25,701 ms |
|
ジャッジサーバーID (参考情報) |
judge6 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 45 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include<atcoder/all>
using namespace atcoder;
using mint=atcoder::modint998244353;
template<long long mod,long long MAX_N>
struct factional_prime{
long long inv_[MAX_N+1];
long long fac_[MAX_N+1];
long long fac_inv_[MAX_N+1];
factional_prime(){
inv_[0]=0;inv_[1]=fac_[0]=fac_[1]=fac_inv_[0]=fac_inv_[1]=1;
for(long long i=2;i<=MAX_N;i++){
inv_[i]=((mod-mod/i)*inv_[mod%i])%mod;
fac_[i]=(fac_[i-1]*i)%mod;
fac_inv_[i]=(fac_inv_[i-1]*inv_[i])%mod;
}
}
long long inv(long long n){
if(n<0)return 0;
return inv_[n];
}
long long fac(long long n){
if(n<0)return 0;
return fac_[n];
}
long long finv(long long n){
if(n<0)return 0;
return fac_inv_[n];
}
long long nCr(long long n,long long r){
if(n<r||n<0||r<0)return 0;
return ((fac_[n]*fac_inv_[n-r])%mod*fac_inv_[r])%mod;
}
long long nPr(long long n,long long r){
if(n<r||n<0||r<0)return 0;
return (fac_[n]*fac_inv_[n-r])%mod;
}
};
factional_prime<998244353,5000000> fp;
int main(){
int N;cin>>N;
if(N%2){
cout<<0<<endl;
return 0;
}
mint ans=0;
for(int i=0;i<N;i++){
for(int j=i+1;j<N;j++){
if((j-i-1)%2)continue;
int m1=(j-i-1)/2;
mint tmp=mint(fp.nCr(2*m1,m1))/(m1+1);
int m2=N-(j-i+1);m2/=2;
tmp*=mint(fp.nCr(2*m2,m2))/(m2+1);
ans+=m1*tmp;
// cout<<i<<" "<<j<<" : "<<ans.val()<<endl;
}
}
cout<<ans.val()<<endl;
}