【快速傅里叶变换】FFT高精乘法

比较优美的算法，虽然我至今还推倒不出那个定理（多项式差值的唯一性还有逆运算的公式），嘛，记住算了。

 

可以用来优化高精度乘法，做到NlogN,不过这个NlogN真的很蛋疼，我四位一压才可以做到300000一秒出，因为涉及复数和实数运算，常数太大了。

 

基本上是照着算导写的，常数小优化就是那些复数根可以预先处理出来，代码小优化就是DFT和DFT-1可以写在一起，只要把单位复根的y取反即可。

 

看了一下h8oj发现rank1的神牛比我快3倍，难道是底层优化？

rank4。。。

 

program fft;{$inline+} const s=4; type cp=record x,y:real end; arr=array[0..1 shl 17] of cp; var c:array[0..100000] of int64; a,b,w,tt:arr; n,i,p:longint; st:ansistring; wt:cp; operator *(var a,b:cp)c:cp;inline; begin c.x:=a.x*b.x-a.y*b.y;c.y:=a.x*b.y+a.y*b.x end; operator +(var a,b:cp)c:cp;inline;begin c.x:=a.x+b.x;c.y:=a.y+b.y end; operator -(var a,b:cp)c:cp;inline;begin c.x:=a.x-b.x;c.y:=a.y-b.y end; procedure dft(var a:arr;s,t:longint);inline; begin if n&