#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; typedef vector vl; typedef vector vvl; typedef vector vc; typedef vector vs; typedef vector vb; typedef vector vd; typedef pair P; typedef pair pii; typedef vector

vpl; typedef tuple tapu; #define rep(i,n) for(int i=0; i<(n); i++) #define REP(i,a,b) for(int i=(a); i<(b); i++) #define all(x) (x).begin(), (x).end() #define rall(x) (x).rbegin(), (x).rend() const int inf = 1<<30; const ll linf = 1LL<<62; const int MAX = 1020000; ll dy[8] = {1,-1,0,0,1,-1,1,-1}; ll dx[8] = {0,0,1,-1,1,-1,-1,1}; const double pi = acos(-1); const double eps = 1e-7; template inline bool chmin(T1 &a,T2 b){ if(a>b){ a = b; return true; } else return false; } template inline bool chmax(T1 &a,T2 b){ if(a inline void print(T &a){ rep(i,a.size()) cout << a[i] << " "; cout << "\n"; } template inline void print2(T1 a, T2 b){cout << a << " " << b << "\n";} template inline void print3(T1 a, T2 b, T3 c){ cout << a << " " << b << " " << c << "\n"; } ll pcount(ll x) {return __builtin_popcountll(x);} //const int mod = 1e9 + 7; const int mod = 998244353; vector> fft(vector> a, bool inverse = false){ int n = a.size(); int h = 0; while((1<>k & 1) << (h-1-k); if(i < j) swap(a[i], a[j]); } for(int b=1; b w = polar(1.0, -(2*pi) / (2*b) * j * (inverse ? -1 : 1)); for(int k=0; k s = a[j+k]; complex t = a[j+k+b] * w; a[j+k] = s + t; a[j+k+b] = s - t; } } } if(inverse) for(int i=0; i> fft(vector a, bool inverse = false){ int n = a.size(); vector> ret(n); for(int i=0; i(a[i],0); return fft(ret, inverse); } vector convolve(vector a, vector b){ int sz = a.size() + b.size() - 1; int p = 1; while(p < sz) p *= 2; a.resize(p); b.resize(p); vector> A = fft(a); vector> B = fft(b); for(int i=0; i