#include #include #include #include #include using namespace std; using u32 = uint32_t; constexpr u32 O = 3, M = 998'244'353, M2 = M - 2; inline constexpr u32 mod_mul(uint64_t x, u32 y) { return static_cast(x * y % M); } map, u32> cache; inline u32 mod_pow(u32 a, u32 n) { auto key = make_pair(a, n); if(cache.count(key)) { return cache[key]; } u32 res = 1; for(; n; n>>=1) { if(n & 1) { res = mod_mul(res, a); } a = mod_mul(a, a); } return cache[key] = res; } void myfmt(vector &a, bool inv) { u32 n = static_cast(a.size()); if(n == 1) { return; } size_t m = n / 2; vector a0(m), a1(m); for(size_t i=0, j=0; i= M) { a[k] -= M; } pz = mod_mul(pz, z); } for(size_t i=0, k=m; k= M) { a[k] -= M; } pz = mod_mul(pz, z); } } void fmt(vector &a) { myfmt(a, false); } void ifmt(vector &a) { myfmt(a, true); u32 n = static_cast(a.size()); u32 inv = mod_pow(n, M2); for(size_t i=0; i convol(vector &a, vector &b) { size_t n = 1; while(n < a.size() + b.size()) { n <<= 1; } a.resize(n); b.resize(n); fmt(a); fmt(b); vector c(n); for(size_t i=0; i> que; for(int i=0; i(--a % M), 1}); } while(que.size() > 1) { vector a = que.front(); que.pop(); vector b = que.front(); que.pop(); que.emplace(convol(a, b)); } vector v = que.front(); for(int q=0; q