洛谷 3803 【模板】多项式乘法（FFT）

题目：https://www.luogu.org/problemnew/show/P3803

第一道FFT！

https://www.cnblogs.com/zwfymqz/p/8244902.html

http://www.cnblogs.com/RabbitHu/p/FFT.html

就是把系数转化为2*n个点值，点值相乘一下，再转化回2*n个系数的过程。

转化为点值的过程就是倍增一样，第一步是w_{1,0}，也就是说x都是1，所以一开始2*n个位置上的点值都是原来的系数；然后变成两个一组取w_{2,0}，w_{2,1}的点值，最后变成2*n个分别取w_{2*n,0}，w_{2*n,1}，......，w_{2*n,2*n-1}的点值。过程就是DFT，证明可见上面博客。

从点值转化回系数的方法和DFT差不多，似乎只要把 x 都变成倒数、做一边刚才的就行。变成倒数的方法就是那个Wn方向变成负的，这样 x^k 就是倒着转的，上面那个角标就一直是和原来相反的了。

最后算答案的时候似乎直接把虚数的部分舍弃了。

对于iDFT的证明还有点不明白。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define db double
using namespace std;
const int N=1e6+;const db pi=acos(-1.0);
int n,m,len,r[N<<];//<<2! for (n+m)<<1
struct cpl{
  db x,y;
}I,a[N<<],b[N<<];
cpl operator+ (cpl a,cpl b){return (cpl){a.x+b.x,a.y+b.y};}
cpl operator- (cpl a,cpl b){return (cpl){a.x-b.x,a.y-b.y};}
cpl operator* (cpl a,cpl b){return (cpl){a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};}
int rdn()
{
  int ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='') ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
  return fx?ret:-ret;
}
void fft(cpl *a,bool fx)
{
  for(int i=;i<len;i++)
    if(i<r[i])swap(a[i],a[r[i]]);
  for(int R=;R<=len;R<<=)//<=
    {
      int m=R>>;
      cpl Wn=(cpl){ cos(pi/m),(fx?-:)*sin(pi/m) };
      for(int i=;i<len;i+=R)
    {
      cpl w=I;
      for(int j=;j<m;j++,w=w*Wn)
        {
          cpl tmp=w*a[i+m+j];
          a[i+m+j]=a[i+j]-tmp;
          a[i+j]=a[i+j]+tmp;
        }
    }
    }
}
int main()
{
  n=rdn(); m=rdn(); I.x=; I.y=;
  for(int i=;i<=n;i++)a[i].x=rdn();
  for(int i=;i<=m;i++)b[i].x=rdn();
  len=;
  while(len<=n+m)len<<=;//<=
  for(int i=;i<len;i++)
    r[i]=(r[i>>]>>)+((i&)?len>>:);
  fft(a,);  fft(b,);
  for(int i=;i<len;i++)
    a[i]=a[i]*b[i];
  fft(a,);
  for(int i=;i<=n+m;i++)
    printf("%d ",int(a[i].x/len+0.5));puts("");
  return ;
}