結果
問題 |
No.2237 Xor Sum Hoge
|
ユーザー |
👑 ![]() |
提出日時 | 2023-02-21 16:29:07 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 238 ms / 10,000 ms |
コード長 | 1,584 bytes |
コンパイル時間 | 4,297 ms |
コンパイル使用メモリ | 268,796 KB |
最終ジャッジ日時 | 2025-02-10 19:43:38 |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 32 |
ソースコード
#include <bits/stdc++.h> #pragma GCC optimize("Ofast") using namespace std; using ll=long long; const int mod=998244353; #define rep(i,a,b) for (ll i=a;i<b;i++) #include<atcoder/all> using namespace atcoder; namespace po167{ struct combination{ int upper; long long MOD; std::vector<long long> fact; std::vector<long long> rev; std::vector<long long> fact_rev; combination(int max,long long mod):upper(max),MOD(mod),fact(max+1),rev(max+1),fact_rev(max+1){ for(long long i=0;i<=max;i++){ if(i<2){ fact[i]=1; fact_rev[i]=1; rev[i]=1; continue; } fact[i]=(fact[i-1]*i)%mod; rev[i]=mod-((mod/i)*rev[mod%i])%mod; fact_rev[i]=(fact_rev[i-1]*rev[i])%mod; } } long long Comb(int x,int y){ assert(upper>=x); if (x<y||y<0||x<0) return 0; return (((fact_rev[y]*fact_rev[x-y])%MOD)*fact[x])%MOD; } long long P(int x,int y){ assert(upper>=x); if (x<y||y<0||x<0) return 0; return (fact_rev[x-y]*fact[x])%MOD; } }; } using po167::combination; void solve(); // oddloop int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int t=1; //cin>>t; rep(i,0,t) solve(); } void solve(){ int N; ll B,C; cin>>N>>B>>C; combination table(N,mod); vector<vector<ll>> f(2); rep(i,0,N+1) f[i&1].push_back(table.Comb(N,i)); vector<ll> ans={1}; rep(d,0,63){ int ind=(B&1)-(C&1); vector<ll> n_ans; while(ind<(int)(ans.size())){ if(ind<0) n_ans.push_back(0); else n_ans.push_back(ans[ind]); ind+=2; } ans=convolution(n_ans,f[C&1]); B>>=1,C>>=1; } cout<<ans[0]<<"\n"; }