题意;求两个数相乘。

  第一次写非递归的fft,因为一个数组开小了调了两天TAT。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;
#define PI 3.1415926535897932384
#define MAXN 1200000
#pragma optimize("O2")
struct Complex
{
double x,y;
Complex(){};
Complex(double x,double y):x(x),y(y){};
void init(double x,double y)
{
this->x=x;this->y=y;
}
Complex operator +(Complex a)
{
Complex ret(x + a.x, y + a.y);
return ret;
}
Complex operator -(Complex a)
{
Complex ret(x - a.x, y - a.y);
return ret;
}
Complex operator *(Complex a)
{
Complex ret(x*a.x - y*a.y, x*a.y + y*a.x);
return ret;
}
Complex operator /(double v)
{
Complex ret(x/v,y/v);
return ret;
}
Complex operator =(double x)
{
this->x=x;
return *this;
}
void operator -=(Complex a)
{
this->x-=a.x;
this->y-=a.y;
}
void operator +=(Complex a)
{
this->x+=a.x;
this->y+=a.y;
}
};
char ss[MAXN];
int t1[MAXN],t2[MAXN];
int pp[MAXN];
Complex ww[MAXN][];
Complex g1[MAXN],g2[MAXN];
long long res[MAXN];
int n;
void dft(Complex g[],int len,int sign)
{
Complex t;
for (int i=;i<len;i<<=)
{
for (int j=;j<len;j+=(i<<))
{
for (int k=;k<i;k++)
{
t=g[j+k];
g[j+k]=g[j+k]+g[j+k+i]*ww[k*((n>>)/i)][sign];
g[j+k+i]=t-g[j+k+i]*ww[k*((n>>)/i)][sign];
}
}
}
} int main()
{
freopen("input.txt","r",stdin);
// freopen("output.txt","w",stdout);
int i,j,k,x,y,z,m;
int n1,n2;
int nn;
scanf("%d",&nn);
while (nn--)
{ scanf("%s",ss);
m=strlen(ss);
for (i=;i<m;i++)
t1[(m-i-)/]=t1[(m-i-)/]*+ss[i]-'';
n1=(m-)/+;
scanf("%s",ss);
m=strlen(ss);
for (i=;i<m;i++)
t2[(m-i-)/]=t2[(m-i-)/]*+ss[i]-'';
n2=(m-)/+;
n=max(n1,n2);
while (n!=(n&(-n)))n-=n&(-n);
n<<=;
for (i=;i<n;i++)
for (j=n>>,x=;j;j>>=,x<<=)
if (i&j)pp[i]+=x;
for (i=;i<n1;i++)g1[pp[i]].x=t1[i];
for (i=;i<n2;i++)g2[pp[i]].x=t2[i];
for (i=;i<=n;i++)
{
ww[i][].init(cos(*PI*i/n), -sin(*PI*i/n));
ww[i][].x=ww[i][].x;
ww[i][].y=-ww[i][].y;
}
dft(g1,n,);
dft(g2,n,);
for (i=;i<n;i++)g2[i]=g1[i]*g2[i];
for (i=;i<n;i++)g1[pp[i]]=g2[i];
dft(g1,n,);
for (i=;i<n;i++)res[i]=(long long)(g1[i].x/n+0.5);
for (i=;i<n;i++)res[i+]+=res[i]/,res[i]%=;
m=n;
while (m && !res[m-])m--;
printf("%d",(int)res[m-]);
for (i=m-;i>=;i--)
printf("%05d",(int)res[i]);
printf("\n");
memset(t1,,sizeof(t1[])*n);
memset(t2,,sizeof(t2[])*n);
memset(pp,,sizeof(pp[])*n);
memset(res,,sizeof(res[])*n);
memset(g1,,sizeof(g1[])*n);
memset(g2,,sizeof(g2[])*n);
}
}

spoj VFMUL FFT快速傅立叶变换模板题的更多相关文章

  1. FFT(快速傅立叶变换):HDU 1402 A * B Problem Plus

    Calculate A * B. Input Each line will contain two integers A and B. Process to end of file. Note: th ...

  2. FFT快速傅立叶变换的工作原理

    实数DFT,复数DFT,FFTFFT是计算DFT的快速算法,但是它是基于复数的,所以计算实数DFT的时候需要将其转换为复数的格式,下图展示了实数DFT和虚数DFT的情况,实数DFT将时域中N点信号转换 ...

  3. FFT快速傅立叶变换

    //最近突然发现博客园支持\(\rm\LaTeX\),非常高兴啊! 话说离省选只有不到五天了还在学新东西确实有点逗…… 切到正题,FFT还是非常神奇的一个东西,能够反直觉地把两个多项式相乘的时间复杂度 ...

  4. FFT快速傅立叶变换:解析wav波频图、Time Domain、Frequency Domain

    您好,此教程将教大家使用scipy.fft分析wav文件的波频图.Time Domain.Frequency Domain. 实际案例:声音降噪,去除高频. 结果: 波频图: Time Domain:

  5. 【bzoj2179】FFT快速傅立叶 FFT模板

    2016-06-01  09:34:54 很久很久很久以前写的了... 今天又比较了一下效率,貌似手写复数要快很多. 贴一下模板: #include<iostream> #include& ...

  6. 离散傅立叶变换与快速傅立叶变换(DFT与FFT)

    自从去年下半年接触三维重构以来,听得最多的词就是傅立叶变换,后来了解到这个变换在图像处理里面也是重点中的重点. 本身自己基于高数知识的理解是傅立叶变换是将一个函数变为一堆正余弦函数的和的变换.而图像处 ...

  7. 快速傅立叶变换(FFT)算法

    已知多项式f(x)=a0+a1x+a2x2+...+am-1xm-1, g(x)=b0+b1x+b2x2+...+bn-1xn-1.利用卷积的蛮力算法,得到h(x)=f(x)g(x),这一过程的时间复 ...

  8. $\mathcal{FFT}$·$\mathcal{Fast \ \ Fourier \ \ Transformation}$快速傅立叶变换

    \(2019.2.18upd:\) \(LINK\) 之前写的比较适合未接触FFT的人阅读--但是有几个地方出了错,大家可以找一下233 啊-本来觉得这是个比较良心的算法没想到这么抽搐这个算法真是将一 ...

  9. BZOJ 2194 快速傅立叶变换之二 | FFT

    BZOJ 2194 快速傅立叶变换之二 题意 给出两个长为\(n\)的数组\(a\)和\(b\),\(c_k = \sum_{i = k}^{n - 1} a[i] * b[i - k]\). 题解 ...

随机推荐

  1. ant有什么用

    内容摘要: ANT是一个基于Java的自动化脚本引擎,脚本格式为XML.除了做Java编译相关任务外,ANT还可以通过插件实现很多应用的调用. 1)ANT的基本概念: 2)ANT的安装:解包,设置路径 ...

  2. iOS开发经验总结(下)

    四十.AFNetworking 传送 form-data 将JSON的数据,转化为NSData, 放入Request的body中. 发送到服务器就是form-data格式. 四十一.非空判断注意 BO ...

  3. Pyhon之类学习1

    #!/usr/bin/python # Filename: class.py __metaclass__=type class Person: def set_name(self,name): sel ...

  4. tcl/tk实例详解——返回一个文件夹下所有文件的绝对路径

    http://blog.csdn.net/dulixin/article/details/2133840 #所有代码如下,使用注释的方式讲解脚本#修改好文件夹和保存结果路径,可以把本文件直接拷贝进tc ...

  5. PHP计算2点经纬度之间的距离

    hp] view plaincopy function getDistanceBetweenPointsNew($latitude1, $longitude1, $latitude2, $longit ...

  6. php生成验证码图片

    0:效果图 1:index.php用来生成验证码图片 <?php session_start(); header ('Content-Type: image/png'); $image=imag ...

  7. eclipse上传显示svn上传者名

    这里来记录下eclipse上传显示svn上传者名. 如图所示,修改即可.

  8. datejs lib

    // Get today's date Date.today(); // Add 5 days to today Date.today().add(5).days(); // Get Friday o ...

  9. 响应者链 hittest:withEvent: 方法的使用

    关于响应者链部分的基础内容 参考http://www.cnblogs.com/wendingding/p/3795171.html 这里我要说明的是 关于- (UIView *)hitTest:(CG ...

  10. ZIP压缩文件夹中上个月的文件,并将备份文件拷贝到服务器

    遍历文件夹的子文件夹下的所有文件,将上个月的文件集中到一起,然互压缩,并copy到服务器的映射磁盘. static void Main(string[] args) { //原始文件存放的位置 Dir ...