結果
問題 |
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; }