bzoj 3527: [Zjoi2014]力

都说这是个FFT模板题 可我这种蒟蒻还是看了大半天。。。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<queue>
#include<algorithm>
#include<vector>
#include<complex>
#define M 1000009
#define EPS 1e-10
#define MO 10000
#define ll long long
#define E complex<double>
#define pi acos(-1)
using namespace std;
ll read()
{
    char ch=getchar();
    ll x=0,f=1;
    for(;ch<'0'||ch>'9';ch=getchar())
        if(ch=='-')
          f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar())
        x=x*10+ch-'0';
    return x*f;
}
E g[M],f[M],_f[M];
int n,L,R[M];
void fft(E a[M],int f)
{
  for(int i=0;i<n;i++)
    if(i>R[i])
      swap(a[i],a[R[i]]);
  for(int i=1;i<n;i<<=1)
    {
      E wn(cos(pi/i),f*sin(pi/i));
      for(int j=0;j<n;j+=(i<<1))
        {
          E W(1,0);
          for(int k=0;k<i;k++,W*=wn)
            {
              E x=a[j+k],y=W*a[i+j+k];
              a[j+k]=x+y;
              a[i+j+k]=x-y;
            }
        }
    }
  if(f==-1)
    for(int i=0;i<n;i++)
      a[i]/=n;
}
int main()
{
   n=read();
   n--;
   for(int i=0;i<=n;i++)
     {
       double a1;
       scanf("%lf",&a1);
       f[i]=a1;
       _f[n-i]=a1;
     }
   for(int i=1;i<=n;i++)
     g[i]=(1.0/i/i);
   int m=2*n;
   for(n=1;n<=m;n<<=1)
     L++;
   for(int i=0;i<n;i++)
     R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
   fft(f,1);
   fft(_f,1);
   fft(g,1);
   for(int i=0;i<n;i++)
     {
       f[i]*=g[i];
       _f[i]*=g[i];
     }
   fft(f,-1);
   fft(_f,-1);
   for(int i=0;i<=m/2;i++)
     printf("%.3lf
",f[i].real()-_f[m/2-i].real());
   return 0;
}