結果
問題 | No.2383 Naphthol |
ユーザー |
|
提出日時 | 2023-07-14 21:54:12 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,789 bytes |
コンパイル時間 | 1,820 ms |
コンパイル使用メモリ | 195,020 KB |
最終ジャッジ日時 | 2025-02-15 13:58:25 |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 WA * 6 |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll=long long; using pii=pair<int,int>; #define all(a) a.begin(),a.end() #define pb push_back #define sz(a) ((int)a.size()) const int N=210005,mod=998244353; int add(int x, int y){x+=y; if(x>=mod) x-=mod; return x;} int sub(int x, int y){x-=y; if(x<0) x+=mod; return x;} int mul(int x, int y){return ((ll)x)*y%mod;} int Pow(int x, ll y=mod-2){int res=1; for(; y; x=mul(x,x),y>>=1) if(y&1) res=mul(res,x); return res;} int fac[N],inv[N],ifac[N]; inline int C(int n, int m){if(m<0||m>n) return 0; return mul(fac[n],mul(ifac[m],ifac[n-m]));} void init_comb(){ fac[0]=inv[1]=ifac[0]=1; for(int i=1; i<N; ++i) fac[i]=mul(fac[i-1],i); for(int i=2; i<N; ++i) inv[i]=mul(inv[mod%i],mod-mod/i); for(int i=1; i<N; ++i) ifac[i]=mul(ifac[i-1],inv[i]); } int n,k,cnt[4]; signed main(){ ios_base::sync_with_stdio(0),cin.tie(0); init_comb(); cin >> n >> k; if(2*n+4-k<k) k=2*n+4-k; if(n==1){ if(k<=1) cout << "1\n"; else if(k==2) cout << "3\n"; else cout << "3\n"; return 0; } cnt[0]=C(2*n+4,k); cnt[1]=k&1?0:C(n+2,k/2); if(n&1){ for(int i=0; i<=2; ++i) if((k+i)%2==0) cnt[2]=add(cnt[2],mul(C(2,i),C(n+1,(k-i)/2))); } else{ if(k%2==0) cnt[2]=C(n+2,k/2); } if(n&1){ if(k%4==0) cnt[3]=C((n+1)/2,k/4); if(k%4==2) cnt[3]=C((n+1)/2,k/4); } else{ if(k%4==0) cnt[3]=C((n+2)/2,k/4); } cnt[1]=sub(cnt[1],cnt[3]),cnt[2]=sub(cnt[2],cnt[3]); for(int i=1; i<=3; ++i) cnt[0]=sub(cnt[0],cnt[i]); //cout << cnt[0] << ' ' << cnt[1] << ' ' << cnt[2] << ' ' << cnt[3] << "\n"; int res=0; for(int i=0; i<4; ++i) res=add(res,mul(cnt[i],Pow(1<<(2-__builtin_popcount(i))))); cout << res << "\n"; }