#include #include #include #include #include using namespace std; using u32 = uint32_t; constexpr u32 O = 3, D = 119, S = 23, M = D << S | 1, M2 = M - 2; inline constexpr u32 mod_mul(u32 x, uint64_t y) { return static_cast(x * y % M); } inline constexpr u32 mod_pow(u32 a, u32 n) { u32 res = 1; for(; n; n>>=1) { if(n & 1) { res = mod_mul(res, a); } a = mod_mul(a, a); } return res; } template struct Seq { u32 arr[N+N]; constexpr Seq() : arr() { arr[0] = mod_pow(O, D); arr[N] = mod_pow(arr[0], M-2); for(size_t i=0; i(); void myfmt(vector &a, bool inv, size_t ln) { 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, size_t ln) { myfmt(a, false, ln); } void ifmt(vector &a, size_t ln) { myfmt(a, true, ln); 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, ln = 0; while(n < a.size() + b.size() - 1) { n <<= 1; ++ln; } a.resize(n); b.resize(n); fmt(a, ln); fmt(b, ln); 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