题意;求两个数相乘。

  第一次写非递归的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. BS开发平台,一小时搞定功能强大的统计分析页面

    BS开发平台,一小时搞定功能强大的统计分析页面,拥有强大的功能和详细的权限控制. 1.组织数据,分析需求(实际耗时大约20分钟)  2.建立需要的业务数据表(大致10分钟)3. 运行代码工具,生产需要 ...

  2. unity3d最新面试题与参考答案汇总

    1.在类的构造函数前加上static会报什么错?为什么? 构造函数格式为 public+类名,如果加上static会报错(静态构造函数不能有访问修饰符)原因:静态构造函数不允许访问修饰符,也不接受任何 ...

  3. 20151221jquery学习笔记--验证插件

    验证插件(validate.js),是一款验证常规表单数据合法性的插件.使用它,极大的解放了在表单上繁杂的验证过程,并且错误提示显示的完善也增加了用户体验.一. 使用 validate.js 插件官网 ...

  4. 在Eclipse的DDMS中查看手机data文件夹中的内容

    在模拟器状态下,在Eclipse的DDMS中的File Explorer下查看手机data/data中的内容是件非常轻松的事情,特别是查看databases 但是在真机模式下就全然不是那么一回事了,在 ...

  5. AudioStreamer使用之快速点击下/上一首按钮,音频会重复的问题解决。

    针对此问题,我的做法是在AudioStreamer.h文件新加入一个属性 @property (nonatomic,strong) NSURL *musicUrl; 和一个单例方法: +(id)sha ...

  6. linux 列出安装包内的文件

    rpm系列: 1. 文件类型  rpm -qpl xxx.rpm2. 已安装文件  rpm -ql packagename debain系列: dpkg -c xxxx.deb

  7. 转JSONObject put,accumulate,element的区别

        public Object put (Object key, Object value) 将value映射到key下.如果此JSONObject对象之前存在一个value在这个key下,当前的 ...

  8. C++ GUI Programming with Qt4 笔记 -- chap1

    1. Hello Qt #include <QApplication> #include <QLabel> int main(int argc, char *argv[]){ ...

  9. 页面mask css

    <html> <head> <style type="text/css"> .share_mask { position: fixed; top ...

  10. SQL优化(2)

    建表时候数据库引擎的选择也可以达到优化的效果 InnoDB: 基于磁盘的资源是InnoDB表空间数据文件和它的日志文件,InnoDB 表的大小只受限于操作系统文件的大小,一般为 2GB MyISAM: ...