//Let's join Kaede Takagaki Fan Club !! #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define int long long //#define L __int128 typedef long long ll; typedef pair P; typedef pair P1; typedef pair P2; #define pu push #define pb push_back #define eb emplace_back #define mp make_pair #define eps 1e-7 #define INF 1000000000 #define a first #define b second #define fi first #define sc second //#define rng(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define rep(i,x) for(int i=0;i bool chmax(t&a,u b){if(a bool chmin(t&a,u b){if(b using vc=vector; template ostream& operator<<(ostream& os,const pair& p){ return os<<"{"< ostream& operator<<(ostream& os,const vc& v){ os<<"{"; for(auto e:v)os<> 30)) * 0xbf58476d1ce4e5b9; x = (x ^ (x >> 27)) * 0x94d049bb133111eb; return x ^ (x >> 31); } size_t operator()(uint64_t x) const { static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count(); return splitmix64(x + FIXED_RANDOM); } //don't make x negative! size_t operator()(pair x)const{ return operator()(uint64_t(x.first)<<32|x.second); } }; //unordered_set -> dtype, null_type //unordered_map -> dtype(key), dtype(value) using namespace __gnu_pbds; template using hash_table=gp_hash_table; template void g(T &a){ cin >> a; } template void o(const T &a,bool space=false){ cout << a << (space?' ':'\n'); } //ios::sync_with_stdio(false); const ll mod = 998244353; //const ll mod = 1000000007; mt19937_64 mt(chrono::steady_clock::now().time_since_epoch().count()); template void add(T&a,T b){ a+=b; if(a >= mod) a-=mod; } ll modpow(ll x,ll n){ ll res=1; while(n>0){ if(n&1) res=res*x%mod; x=x*x%mod; n>>=1; } return res; } #define _sz 400005 ll F[_sz],R[_sz]; void make(){ F[0] = 1; for(int i=1;i<_sz;i++) F[i] = F[i-1]*i%mod; R[_sz-1] = modpow(F[_sz-1], mod-2); for(int i=_sz-2;i>=0;i--) R[i] = R[i+1] * (i+1) % mod; } ll C(int a,int b){ if(b < 0 || a < b) return 0; return F[a]*R[b]%mod*R[a-b]%mod; } template struct ntt{ inline void add(int &a, int b) { a += b; if(a >= md) a -= md; } inline void sub(int &a, int b) { a -= b; if(a < 0) a += md; } inline int mul(int a, int b) { return (int)((ll)a*b%md); } inline int power(int a, long long b) { int res = 1; while (b > 0) { if (b & 1) res = mul(res, a); a = mul(a, a); b >>= 1; } return res; } inline int inv(int a) { a %= md; if (a < 0) a += md; int b = md, u = 0, v = 1; while (a) { int t = b / a; b -= t * a; swap(a, b); u -= t * v; swap(u, v); } assert(b == 1); if (u < 0) u += md; return u; } int max_base, root; vector dw, idw; ntt() { int tmp = md - 1; max_base = 0; while (tmp % 2 == 0) { tmp /= 2; max_base++; } root = 2; while (power(root, (md-1)>>1) == 1) root++; dw.resize(max_base); idw.resize(max_base); rep(i, max_base){ sub(dw[i], power(root, (md-1) >> (i+2))); idw[i] = inv(dw[i]); } } void fft(vector &a, bool inv) { const int n = a.size(); assert((n & (n - 1)) == 0); assert(__builtin_ctz(n) <= max_base); if(!inv){ for(int m=n;m>>=1;){ int w = 1; for(int s=0,k=0; s=y?x-y:x+md-y); a[i] = (x+y>=md?x+y-md:x+y); } w = mul(w, dw[__builtin_ctz(++k)]); } } } else{ for(int m=1;m=y?x-y:x+md-y); a[j] = mul(a[j], w); a[i] = (x+y>=md?x+y-md:x+y); } w = mul(w, idw[__builtin_ctz(++k)]); } } const int inv_sz = this->inv(n); for(auto &&e:a) e = mul(e, inv_sz); } } vector multiply(vector a, vector b, int eq = 0) { int need = a.size() + b.size() - 1; int nbase = 0; while ((1 << nbase) < need) nbase++; int sz = 1 << nbase; a.resize(sz); b.resize(sz); fft(a, 0); if (eq) b = a; else fft(b, 0); rep(i, sz) a[i] = mul(a[i], b[i]); fft(a, 1); a.resize(need); return a; } vector square(vector a) { return multiply(a, a, 1); } }; ntt<998244353>nt; #define vi vector vi init(vi a, int x){ a.resize(x,0); return a; } vi add(vi a, vi b,int M=-1){ if(a.size() < b.size()) swap(a,b); for(int i=0;i= mod) a[i] -= mod; } if(M >= 0 && a.size() > M) a.resize(M); return a; } vi sub(vi a, vi b,int M=-1){ if(a.size() < b.size()) a.resize(b.size(),0); for(int i=0;i= mod) a[i] -= mod; } if(M >= 0 && a.size() > M) a.resize(M); return a; } //微分 vi diff(vi a,int M=-1){ reverse(a.begin(),a.end()); a.pop_back(); reverse(a.begin(),a.end()); for(int i=0;i= 0 && a.size() > M) a.resize(M); if(a.empty()) a.pb(0); return a; } //積分、定数項は指定可能 vi intg(vi a,int coef=0,int M=-1){ reverse(a.begin(),a.end()); a.push_back(coef); reverse(a.begin(),a.end()); for(int i=1;i= 0 && a.size() > M) a.resize(M); return a; } vi mul_int(vi a,int x,int M=-1){ for(int i=0;i= 0 && a.size() > M) a.resize(M); return a; } vi mul(vi a,vi b,int M=-1){ //modが99...のときはnttできる //任意modの場合は、適当に区切ってもいいしgarnerしてもいいし任せます a = nt.multiply(a,b); if(M >= 0 && a.size() > M) a.resize(M); return a; } vi lw(vi a, int x){ if(a.size() > x) a.resize(x); return a; } vi inv(vi a,int M){ if(a.empty() || a[0] == 0) return vi(); //1の場合をやっておく vi ret(1); ret[0] = modpow(a[0],mod-2); int cur = 1; int nxt = 1; while(cur < M){ auto g = lw(ret, cur); auto f = lw(a, cur*2); f.resize(cur << 1); g.resize(cur << 1); nt.fft(f, 0); nt.fft(g, 0); rep(i, f.size()) f[i] = nt.mul(f[i], g[i]); nt.fft(f, 1); f.erase(begin(f), begin(f) + cur); f.resize(cur << 1), nt.fft(f, 0); rep(i, f.size()) { f[i] = mod - nt.mul(f[i], g[i]); } nt.fft(f, 1); ret.insert(end(ret), begin(f), begin(f)+ cur); nxt++; cur <<= 1; } assert(ret.size() >= M); ret.resize(M); return ret; } //vecで表される多項式の平方根、x^(M-1)の項までを求めてM長のvectorを返す //ただし、vec[0] = 0の場合は無理 //整数の平方根求めるのかったるいのでそれはまた後日 //s1(z)^2 = f(z) (mod z) //sn(z)^2 = f(z) (mod z^n) とする //s(2n)(z) = 1/2 * (sn(z)+f(z)*sn(z)^-1) vi sqrt(vi a,int M){ if(a.empty() || a[0] == 0) return vi(); //1の場合をやっておく vi s[30]; //今日はvec[0] = 1のことしか考えません s[0].pb(1); assert(1LL*s[0][0]*s[0][0]%mod == a[0]); int cur = 1; int nxt = 1; vi vec3; while(cur < M){ int nw = min(cur*2,(int)(a.size())); int sz = vec3.size(); vec3.resize(nw,0); for(int t=sz;t= a.size()) aa.pb(0); else aa.pb(a[aa.size()]); } if(m>1){ vi w = add(q,mul(g,sub(diff(f,2*m-1),mul(f,q,2*m-1),2*m-1),2*m-1),2*m-1); f = add(f,mul(f,sub(aa,intg(w,0,2*m),2*m),2*m),2*m); } else{ vi w = mul(g,diff(f,2*m-1),2*m-1); f = add(f,mul(f,sub(aa,intg(w,0,2*m),2*m),2*m),2*m); } m*=2; while(q.size() < m-1){ if(q.size() >= h.size()) q.pb(0); else q.pb(h[q.size()]); } } f.resize(M); assert(f.size() == M); return f; } //o(ans?"Yes":"No"); void solve(){ int n,k; cin >> n >> k; make(); ll ans = C(2*n, n); if(k == 1); else if(k == 2) ans = (ans+mod-2) % mod; else{ //k-3, k-2, k-1 vcpl[3]; rep(i, 3){ pl[i].resize(k/2+2, 0); for(int x=0;x> t; while(t--) solve(); }