#include #include #define PI2 6.28318530717958647692 const int k = 17; const int m = (1 << 17); typedef double complex cmplx; cmplx C[1<<17]; cmplx w[1<<16]; int read_uint(){ int x=0; char c=0; while(c<'0') c = getchar_unlocked(); do { x *= 10; x += c-'0'; c = getchar_unlocked(); } while(c>='0'); return x; } void write_int(int x){ int i; static char buf[12]; if(x==0){ puts("0"); return; } for(i=10; x; i--){ buf[i] = '0' + x % 10; x /= 10; } puts(buf+1+i); } unsigned minus(unsigned x) { int y = 1 << (31-__builtin_clz(x)); return x^(y-1); } void fft(int k, cmplx *A, const cmplx *w){ const int m = 1 << k; int u = 1; int v = m/4; int i, j; if(k&1){ for(j=0; j>= 1; } for(i=k&~1;i>0;i-=2){ int jh; for(jh=0;jh>= 2; } } void ifft(int k, cmplx *A, const cmplx *w){ const int m = 1 << k; int u = m/4; int v = 1; int i, j; for(i=2;i<=k;i+=2){ int jh; for(jh=0;jh>= 2; v <<= 2; } if(k&1){ for(j = 0;j>1, z); genw(i|b, b>>1, z*w[b]); } } void convolver(int k, cmplx *A, const cmplx *w){ int i; const int m = 1 << k; fft(k, A, w); A[0] = 4 * creal(A[0]) * cimag(A[0]) * I; A[1] = 4 * creal(A[1]) * cimag(A[1]) * I; for(i = 2; i < m; i+=2){ int y = 1 << (31-__builtin_clz(i)); int j = i^(y-1); A[i] = (A[i] + conj(A[j]))*(A[i] - conj(A[j])); A[j] = -conj(A[i]); } for(i = 0; i < m; i+=2){ A[i/2] = (A[i]+A[i^1] - (A[i]-A[i^1])*w[i/2]*I)/(4*m); } ifft(k-1, A, w); } int main(){ int l, mm, q; int i, j; const double arg = -PI2/m; const int B = 100009; for(i=1, j=m/4; j; i<<=1, j>>=1){ w[i] = cexp(I * (arg * j)); } genw(0, m/4, 1); l = read_uint(); mm = read_uint(); read_uint(); for(i=0;i