bzoj2194

http://www.lydsy.com/JudgeOnline/problem.php?id=2194

卷积。。。

卷积并不高深，其实卷积就是两个多项式相乘的系数，但是得满足一点条件，就是f[n]=a[i]*b[n-i],就是下标和固定。。。然后这道题下标和不固定，但是我们把b反过来,就是一个卷积了。每次和是固定的

但是输出的时候得输出从n-2n,因为c[n+k]=a[i]*b[n+k-i],n<=n+k<=2*n

#include<bits/stdc++.h>
using namespace std;
#define pi acos(-1)
const int N = ;
complex<double> a[N], b[N];
int n, m, l;
int r[N];
void fft(complex<double> *a, int f)
{
    for(int i = ; i <= n; ++i) if(i < r[i]) swap(a[i], a[r[i]]);
    for(int i = ; i < n; i <<= )
    {
        complex<double> w(cos(pi / i), f * sin(pi / i));
        for(int p = i << , j = ; j <= n; j += p)
        {
            complex<double> t(, );
            for(int k = ; k < i; ++k, t = t * w)
            {
                complex<double> x = a[j + k], y = t * a[j + k + i];
                a[j + k] = x + y; a[j + k + i] = x - y;
            }
        }
    }
}
int main()
{
    scanf("%d", &n); --n;
    for(int i = ; i <= n; ++i)
    {
        int x, y; scanf("%d%d", &x, &y); a[i] = x; b[n - i] = y;
    }
    m =  * n; for(n = ; n <= m; n <<= ) ++l;
    for(int i = ; i <= n; ++i) r[i] = (r[i >> ] >> ) | ((i & ) << (l - ));
    fft(a, ); fft(b, );
    for(int i = ; i <= n; ++i) a[i] = a[i] * b[i];
    fft(a, -);
    for(int i = m / ; i <= m; ++i) printf("%d\n", (int)(a[i].real() / n + 0.5));
    return ;
}