#include #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 fact; std::vector rev; std::vector 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=x); if (x>t; rep(i,0,t) solve(); } void solve(){ int N; ll B,C; cin>>N>>B>>C; combination table(N,mod); vector> f(2); rep(i,0,N+1) f[i&1].push_back(table.Comb(N,i)); vector ans={1}; auto convolution=[&](vector l,vector r)->vector{ int n=l.size(),m=r.size(); vector res(n+m-1); rep(i,0,n) rep(j,0,m){ res[i+j]=(res[i+j]+l[i]*r[j])%mod; } return res; }; rep(d,0,63){ int ind=(B&1)-(C&1); vector 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<