先上模板 1 #include 2 #include 3 const int P=(1>=1,t=t*t%mod)if(k&1)f=f*t%mod;return f;} 6 int a[1j)tmp=x[i],x[i]=x[j],x[j]=tmp;13 for(l=k>>1;(j^=...
分类：
其他好文 时间：
2014-07-22 08:19:37
阅读次数：
416
Description给出两个n位10进制整数x和y，你需要计算x*y。 Input第一行一个正整数n。 第二行描述一个位数为n的正整数x。 第三行描述一个位数为n的正整数y。 Output输出一行，即x*y的结果。 Sample Input134Sample Output 12HINTn 2 .....
分类：
其他好文 时间：
2014-07-20 00:29:12
阅读次数：
243
/********************************************************************************* 程序名称：快速傅里叶变换(FFT) ** 程序描述：本程序实现快速傅里叶变换 ** 程序作者：宋元瑞 ** 最后修改：2011年4月5...
分类：
编程语言 时间：
2014-07-18 10:12:07
阅读次数：
213
下面的matlab程序分别使用周期图法、相关函数法以及AR谱方法计算信号的功率谱。% power spectrum estimated
clear all;
clc;
close all;
Fs=1000; % 采样频率
nfft = 1024; % fft计算点数
%产生含有噪声的序列
n=0:1/Fs:1;
xn=cos(2*pi*100*n)+3*cos(2*pi*200...
分类：
其他好文 时间：
2014-07-08 16:23:01
阅读次数：
277
题目连接：Codeforces 444B DZY Loves FFT
题目大意：根据题目的算法生成a，b数组，然后对于每个长度的l，求a[i]*b[l-i]的最大值。
解题思路：概率问题，枚举前30大的数，如果有就可以直接输出答案，如果没有，就暴力枚举b数组为1的位置找最大值。
#include
#include
#include
using namespace std;
t...
分类：
其他好文 时间：
2014-07-08 10:41:58
阅读次数：
197
准零基础搞懂FFT快速傅里叶变换及其实现程序（二）上一篇文章我们了解了DFT的原理，FFT是基于DFT的更适合计算机运算的算法，本文我们就正式开始学习FFT的原理。首先我么先来宏观的看一下FFT。如果我们把整个FFT的算法看成一个黑盒子的话，那么它的输入就是时间波形信号，比如声音波形（横轴为时间，纵...
分类：
其他好文 时间：
2014-07-02 21:30:08
阅读次数：
236
题目大意：给出n个数qi，定义 Fj为 令 Ei=Fi/qi，求Ei。看了很久题解，终于有些眉目，因为知道要用FFT，所以思路就很直了其实我们就是要±1/(j-i)^2 ( i-j大于0时为正，小于0时为负 ) 和 qi 的乘积要算到j这个位置上，这个满足卷积，所以用FFT优化，但是j-i有负...
分类：
其他好文 时间：
2014-06-27 12:09:57
阅读次数：
184
为什么要分析算法
分析算法可以预测算法的性能，比较算法之间的优劣，保证算法的正确性，理解算法的理论基础。
成功算法的例子
离散傅立叶变换，如果使用暴力方法，那么算法的复杂度是是N^2，如果使用FFT快速傅立叶变换可以实现O(N logN)复杂度
N-body模拟：使用Barnes-hut算法可以将复杂度降低到N logN
...
分类：
其他好文 时间：
2014-06-03 00:10:12
阅读次数：
364
学了若干天终于学(bei)会了传说中的法法塔感觉也没那么难用嘛fft快速傅里叶变换
在大表课件上写就是解决高精乘的工具
其实很有理有据fft就是用复数的折半引理优化两个多项式相乘的高端东西他能使O(n^2)的多项式相乘优化到O(nlogn)听ak说这也是比较模板的东西
也就不去理解什么证明了（其实是...
分类：
其他好文 时间：
2014-05-19 13:24:27
阅读次数：
296
#include#include#include#include#includeusing
namespace std;typedef double dd;const dd pi=acos(0.0)*2;#define N 400005struct
P{ dd x,y; P(dd A=0...
分类：
其他好文 时间：
2014-05-17 03:38:45
阅读次数：
248
