uoj34多项式乘法

这是一道模板题。
给你两个多项式，请输出乘起来后的多项式。

输入格式
第一行两个整数 nn 和 mm，分别表示两个多项式的次数。
第二行 n+1n+1 个整数，分别表示第一个多项式的 00 到 nn 次项前的系数。
第三行 m+1m+1 个整数，分别表示第一个多项式的 00 到 mm 次项前的系数。

输出格式
一行 n+m+1n+m+1 个整数，分别表示乘起来后的多项式的 00 到 n+mn+m 次项前的系数。

样例一
input
1 2
1 2
1 2 1

output
1 4 5 2

explanation
(1+2x)⋅(1+2x+x2)=1+4x+5x2+2x3(1+2x)⋅(1+2x+x2)=1+4x+5x2+2x3

限制与约定
0≤n,m≤1050≤n,m≤105，保证输入中的系数大于等于 00 且小于等于 99。

时间限制：1s
空间限制：256MB

分析：
如题所言，确实是一道FFT模板题，
寻求曲神帮助
虽然曲神神的有点不可思议，但人还是很亲切的
第一次写，代码极为丑陋

这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>

using namespace std;

const int N=301000;  //数组大小避免开成2的次方 
const double pi=acos(-1.0);
struct node{
    double a,b;
    node (double aa=0,double bb=0)
    {
        a=aa;b=bb;
    }
};
node a[N],b[N],omega[N],a_omega[N]; 
int n,m,k,l;

node operator +(const node &x,const node &y){return node(x.a+y.a,x.b+y.b);}
node operator -(const node &x,const node &y){return node(x.a-y.a,x.b-y.b);}
node operator *(const node &x,const node &y){return node(x.a*y.a-x.b*y.b,x.a*y.b+x.b*y.a);}
node operator *(const node &x,const double &y){return node(x.a*y,x.b*y);}

void init(int n)
{
    for (int i=0;i<n;i++)
    {
        omega[i]=node(cos(2.0*pi*i/n),sin(2.0*pi*i/n));
        a_omega[i]=node(cos(2.0*pi*i/n),-sin(2.0*pi*i/n));
    }
}

void FFT(int n,node *a,node *w)
{
    int i,j=0,k;
    for (i=0;i<n;i++)  //很显然这就是优美的rev操作 
    {
        if (i>j) swap(a[i],a[j]);
        for (int l=n>>1;(j^=l)<l;l>>=1);
    }
    for (int i=2;i<=n;i<<=1)  //合并 
    {
        int m=i>>1;  //一半区间的长度 
        for (int j=0;j<n;j+=i)
            for (k=0;k<m;k++)
            {
                node z=a[j+m+k]*w[n/i*k];
                a[j+m+k]=a[j+k]-z;
                a[j+k]=a[j+k]+z;
            } 
    }
}

int main()
{
    int fn;
    scanf("%d%d",&n,&m);
    for (int i=0;i<=n;i++) scanf("%lf",&a[i].a);
    for (int i=0;i<=m;i++) scanf("%lf",&b[i].a);
    fn=1;
    while (fn<=n+m) fn<<=1;   //找到大于n+m的最小2次幂
    init(fn); 
    FFT(fn,a,omega);
    FFT(fn,b,omega);
    for (int i=0;i<=fn;i++)
        a[i]=a[i]*b[i];
    FFT(fn,a,a_omega);
    for (int i=0;i<=n+m;i++) printf("%d ",(int)(a[i].a/fn+0.5));  //四舍五入 
}