FFT的一种迭代实现

struct Complex
{
    double x,y;
    Complex(double x1=0.0 ,double y1=0.0)
    {
        x=x1;
        y=y1;
    }
    Complex operator -(const Complex &b)const
    {
        return Complex(x-b.x,y-b.y);
    }
    Complex operator +(const Complex &b)const
    {
        return Complex(x+b.x,y+b.y);
    }
    Complex operator *(const Complex &b)const
    {
        return Complex(x*b.x-y*b.y,x*b.y+y*b.x);
    }
};

void change(Complex y[],int len)
{
    int i,j,k;
    ,j=len/;i<len-;i++)
    {
        if(i<j)swap(y[i],y[j]);
        k=len/;
        while(j>=k)
        {
            j-=k;
            k/=;
        }
        j+=k;
    }
}

void FFT(Complex y[],int len,int on)
{
    change(y,len);
    ;h<=len;h<<=)
    {
        Complex wn(cos(on**pi/h),sin(on**pi/h));
        ;j<=len;j+=h)
        {
            Complex w(,);
            ;k++)
            {
                Complex u=y[k];
                Complex t=w*y[k+h/];
                y[k]=u+t;
                y[k+h/]=u-t;
                w=w*wn;
            }
        }
    }
    )
        ;i<len;i++)
        y[i].x/=len;
}