题目大意：求nn个点的无向简单连通图个数，n≤1.3?105n\leq1.3*10^5
递推式：fi=2C2i?∑i?1j=1fj?Cj?1i?1?2C2i?jf_i=2^{C_i^2}-\sum_{j=1}^{i-1}f_j*C_{i-1}^{j-1}*2^{C_{i-j}^2}
推导戳这里
然后两侧同除(i?1)!(i-1)!得到：
fi(i?1)!=2C2i(i?1)!?∑i?1j=1...
分类：
其他好文 时间：
2015-05-27 19:14:58
阅读次数：
160
没什么好说的，今天有考了FFT（虽然不用FFT也能过）但是确实有忘了怎么写FFT了，于是乎只有重新写一遍FFT模板练一下手了。第一部分普通FFT，第二部分数论FFT，记一下模数2^23*7*17+1 #include#include#include#include#includeusing n...
分类：
其他好文 时间：
2015-05-26 23:12:31
阅读次数：
198
利用一维FFT查看正弦波频谱信息，所需的SIMULINK模型如下图1所示。Zero-OrderHold用于采样正弦波，设置方式如图2所示。BUffer用于缓存采样数据，深度可以设置，深度等于FFT的长度，设置方式如图3所示。FFT的点数跟BUFFER的深度相同，比如，buffer深度为1024，那么此处的FFT就采..
分类：
其他好文 时间：
2015-05-26 19:09:29
阅读次数：
284
在simulink中搭建这样的模块还是蛮复杂的，也可能是因为我对Simulink不熟悉，有些高级功能未能用上，所以下图中所示的显得较为复杂。整个过程如下：1，载入一个256*256的图像，如图1所示；2，将图像格式转化为double型；3，将2D矩阵转换为1D矩阵；4，使用buffer，每次读入一行..
分类：
其他好文 时间：
2015-05-25 20:38:20
阅读次数：
159
经常可以在文献上看到2D-FFT能够以两个1D-FFT来实现，今天我用MATLAB证明了，的确如此。MATLAB的代码如下clearall;
clc;
f=ones(256,256);
center_loc=size(f);
rd=2;
f(round(center_loc(1)/2)-rd:round(center_loc(1)/2)+rd,round(center_loc(2)/2)-rd:round(center_loc(2)/2..
分类：
其他好文 时间：
2015-05-25 14:52:54
阅读次数：
333
快速傅里叶变换FFTDFT是信号分析与处理中的一种重要变换。但直接计算DFT的计算量与变换区间长度N的平方成正比，当N较大时，计算量太大，直接用DFT算法进行谱分析和信号的实时处理是不切实际的。
1.直接计算DFT
长度为N的有限长序列x(n)的DFT为：
2.减少运算量的思路和方法
思路：N点DFT的复乘次数等于N2。把N点DFT分解为几个较短的DFT，可使乘法次数大大减少。另外，旋转...
分类：
其他好文 时间：
2015-05-24 08:54:18
阅读次数：
142
上篇讲述了一维FFT的GPU实现（FFT算法实现——基于GPU的基2快速傅里叶变换），后来我又由于需要做了一下二维FFT，大概思路如下。首先看的肯定是公式：如上面公式所描述的，2维FFT只需要拆分成行FFT，和列FFT就行了，其中我在下面的实现是假设原点在F(0,0)，由于我的代码需要原点在中心，所...
分类：
编程语言 时间：
2015-05-22 01:46:24
阅读次数：
369
莫比乌斯反演
莫比乌斯反演在数论中占有重要的地位，许多情况下能大大简化运算。那么我们先来认识莫比乌斯反演公式。
定理：和是定义在非负整数集合上的两个函数，并且满足条件，那么我们得到结论
在上面的公式中有一个函数，它的定义如下：
（1）若，那么
（2）若，均为互异素数，那么
...
分类：
其他好文 时间：
2015-05-17 13:44:53
阅读次数：
279
FFT是离散傅立叶变换的高速算法，能够将一个信号变换到频域。有些信号在时域上是非常难看出什么特征的，可是如果变换到频域之后，就非常easy看出特征了。这就是非常多信号分析採用FFT变换的原因。另外，FFT能够将一个信号的频谱提取出来，这在频谱分析方面也是经经常使用的。 尽管非常多人都知道FFT是什么...
分类：
编程语言 时间：
2015-05-16 17:59:35
阅读次数：
125
FFT基本操作。。。讲解请自己看大学信号转置系列。。。 1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 #define PAU putchar(' ') 8 #define ENT putchar('\n') 9...
分类：
其他好文 时间：
2015-05-16 09:00:08
阅读次数：
151
