十五、多项式乘法与快速傅里叶变换

最新推荐文章于 2024-09-06 16:54:11 发布
最新推荐文章于 2024-09-06 16:54:11 发布
阅读量3.8w 收藏 52
点赞数 9
分类专栏： 十五大经典算法研究 文章标签： fft 算法 matlab winapi delete 图像处理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/v_july_v/article/details/6684636
版权

十五、多项式乘法与快速傅里叶变换

前言

    经典算法研究系列，已经写到第十五章了，本章，咱们来介绍多项式的乘法以及快速傅里叶变换算法。本博客之前也已详细介绍过离散傅里叶变换（请参考：十、从头到尾彻底理解傅里叶变换算法、上，及十、从头到尾彻底理解傅里叶变换算法、下），这次咱们从多项式乘法开始，然后介绍FFT算法的原理与实现。同时，本文虽涉及到不少数学公式和定理（当然，我会尽量舍去一些与本文咱们要介绍的中心内容无关的定理或证明，只为保证能让读者易于接受或理解），但尽量保证通俗易懂，以让读者能看个明白。

    有朋友建议，算法专一种，就ok，没必要各个都学习。但个人实在抑制不住自己的兴趣，就是想写，当没法做到强迫自己不写时，这个经典算法研究系列便一直这么写下来了。

第一节、多项式乘法

    我们知道，有两种表示多项式的方法，即系数表示法和点值表示法。什么是系数表示法?所谓的系数表示法，举个例子如下图所示，A（x）=6x^3 + 7x^2 - 10x + 9，B（x）=-2x^3+4x-5，则C（x）=A（x）*B（x）就是普通的多项式相乘的算法，系数与系数相乘，这就是所谓的系数表示法。

    那么，何谓点值表示法?。一个次数为n次的多项式A（x）的点值表示就是n个点值所形成的集合：{(x0,y0), (x1,y1),..., (xn-1,yn-1)}。其中所有xk各不相同，且当k=0，1,……，n-1时，有y（k）=A（xk）。

    一个多项式可以由很多不同的点值表示，这是由于任意n个相异点x0，x1，...，xn-1组成的集合，都可以看做是这种点值法的表示方法。对于一个用系数形式表示的多项式来说，在原则上计算其点值表示是简单易行的，我们只需要选取n个相异点x0，x1，...，xn-1，然后对k=0，1，...，n-1，求出A（xk），然后根据霍纳法则，求出这n个点的值所需要的时间为O（n^2）。

    然，稍后，你将看到，如果巧妙的选取xk的话，适当的利用点值表示可以使多项式的乘法可以在线性时间内完成。所以，如果我们能恰到好处的利用系数表示法与点值表示法的相互转化，那么我们可以加速多项式乘法（下面，将详细阐述这个过程），达到O（n*logn）的最佳时间复杂度。

    前面说了，当用系数表示法，即用一般的算法表示多项式时，两个 n次多项式的乘法需要用 O（n^2）的时间才能完成。但采用点值表示法时，多项式乘法的运行时间仅为O（n）。接下来，咱们要做的是，如何利用系数表示法与点值表示法的相互转化来实现多项式系数表示法的O（n*logn）的快速乘法。

第二节、多项式的快速乘法 

    在此之前，我们得明白两个概念，求值与插值。通俗的讲，所谓求值（系数->点值），是指系数形式的多项式转换为点值形式的表示方式。而插值（点值->系数）正好是求值的逆过程，即反过来，它是已知点值表示法，而要你转换其为多项式的系数表示法（用n个点值对也可以唯一确定一个不超过n-1次多项式，这个过程称之为插值）。而这个系数表示法与点值表示法的相互转化，即无论是从系数表示法转化到点值表示法所谓的求值，还是从点值表示法转化到系数表示法所谓的插值，求值和插值两种过程的时间复杂度都是O（n*logn）。

    前面，我们已经说了，在多项式乘法中，如果用系数表示法表示多项式，那么多项式乘法的复杂度将为O（n^2），而用点值表示法表示的多项式，那么多项式的乘法的复杂度将是线性复杂度为O（n），即：  适当的利用点值表示可以使多项式的乘法可以在线性时间内完成。此时读者可以发挥你的想象，假设我们以下面这样一种过程来计算多项式的乘法（不过在此之前，你得先把两个要相乘的多项式A（x）和B（x）扩充为次数或者说系数为2n次的多项式），我们将会得到我们想要的结果：

  1. 系数表示法转化为点值表示法。先用系数表示法表示一个多项式，然后对这个多项式进行求值操作，即多项式从系数表示法变成了点值表示法，时间复杂度为O（n*logn）；
  2. 点值表示法计算多项式乘法。用点值表示法表示多项式后，计算多项式的乘法，线性复杂度为O（n）；
  3. 点值表示法转化为系数表示法。最终再次将点值表示法表示的多项式转变为系数表示法表示的多项式，这一过程，为O（n*logn）。

    那么，综上上述三项操作，系数表示法表示的多项式乘法总的时间复杂度已被我们降到了O（n*logn+n+n*logn）=O（N*logN）。

    如下图所示，从左至右，看过去，这个过程即是完成多项式乘法的计算过程。不过，完成这个过程有两种方法，一种就是前面第一节中所说的普通方法，即直接对系数表示法表示的多项式进行乘法运算，复杂度为O（n^2），它体现在下图中得Ordinary multiplication过程。还有一种就是本节上文处所述的三个步骤：1、将多项式由系数表示法转化为点值表示法(点值过程）；2、 利用点值表示法完成多项式乘法；3、将点值表示法再转化为系数表示法（插值过程）。三个步骤下来，总的时间复杂度为O（N*logN）。它体现在下图中的Evaluation + Pointwise multiplication + Interpolation 三个合过程。

    由上图中，你也已经看到了。我们巧妙的利用了关于点值形式的多项式的线性时间乘法算法，来加快了系数形式表示的多项式乘法的运算速度。在此过程中，我们的工作最关键的就在于以O（n*logn）的时间快速把一个多项式从系数形式转换为点值形式（求值，我们即称之为FFT），然后再以O（n*logn）的时间从点值形式转换为系数形式（插值，我们即称之为FFT逆）。最终达到了我们以最终的系数形式表示的多项式的快速乘法为O（N*logN）的时间复杂度。好不令人心生快意。

    对上图，有一点必须说明，项w2n是2n次单位复根。且不容忽视的是，在上述两种表示法即系数表示法和点值表示法相互转换的过程中， 都使用了单位复根，才得以在O（n*logn）的时间内完成求值与插值运算（至于何谓单位复根，请参考相关资料。因为为了保证本文的通俗易懂性，无意引出一大堆公式或定理）。

第三节、FFT

     注：本节有相当部分文字来自算法导论中文版第二版以及维基百科。

     在具体介绍FFT之前，我们得熟悉知道折半定理是怎么一回事儿：如果n>0且n为偶数，那么，n个n次单位复根的平方等于n/2个n/2个单位复根。我们已经知道，通过使用一种称为快速傅里叶变换（FFT）的方法，可以在O（n*logn）的时间内计算出DFTn（a），而若采用直接的方法复杂度为O（n^2）。FFT就是利用了单位复根的特殊性质。

    你将看到，FFT方法运用的策略为分治策略，所谓分治，即分而治之。两个多项式A（x）与B（x）相乘的过程中，FFT用A（x）中偶数下标的系数与奇数下标的系数，分别定义了两个新的次数界为n/2的多项式A[0]（x）和A[1]（x）：

A[0](x) =a0 +a2x +a4x2 +··· +an-2xn/2-1,

A[1](x) =a1 +a3x +a5x2 +··· +an-1xn/2-1.

    注意，A[0]包含了A中所有偶数下标的系数（下标的相应二进制数的最后一位为0），而A[1]包含A中所有奇数下标的系数（下标的相应二进制数的最后一位为1）。有下式：

     这样，求A（x）在处的问题就转换为：

     1）求次数界为n/2的多项式A[0]与A[1]在点的值，然后

     2）将上述结果进行组合。

    下面，我们用N次单位根WN来表示。

WN的性质：

  1. 周期性，WN具有周期N。
  2. 对称性：。

    为了简单起见，我们下面设待变换序列长度n = 2r。 根据上面单位根的对称性，求级数时，可以将求和区间分为两部分：

Fodd(k) 和 Feven(k)是两个分别关于序列奇数号和偶数号序列N/2点变换。由此式只能计算出yk的前N/2个点，对于后N/2个点，注意Fodd(k) 和 Feven(k) 都是周期为N/2的函数，由单位根的对称性，于是有以下变换公式：

  • 。

    这样，一个N点变换就分解成了两个N/2点变换，照这样可继续分解下去。这就是库利-图基快速傅里叶变换算法的基本原理。此时，我们已经不难分析出此时算法的时间复杂度将为O(NlogN)。

     ok，没想到，本文之中还是出现了这么多的公式（没办法，搞这个FFT就是个纯数学活儿。之前买的一本小波与傅里叶分析基础也是如此，就是一本数学公式和定理的聚集书。不过，能看懂更好，实在无法弄懂也只权当做个稍稍了解）。

    好了，下面，咱们来简单理解下FFT中的蝶形算法，本文将告结束。如下图所示：

 

      有人这样解释蝶形算法：对于N（2的x次方）点的离散信号，把它按索引位置分成两个序列，分别为0，2，4，....，2K（记为A）和1，3，5，....，2K-1（记为B），由傅立叶变换可以推出N点的傅立叶变换前半部分的结果为A+B*旋转因子，后半部分的结果为A-B*旋转因子。于是求N点的傅立叶变换就变成分别求两个N/2点序列的傅立叶变换，对每一个N/2点的序列，递归前面的步骤，直到只有两点的序列，就只变成简单的加减关系了。把这些点的加减关系用线连接，看上去就是个蝶形。ok，更多可参考算法导论第30章。

举一个例子，我们知道，4X4的矩阵运算如果按常规算法的话仍要进行64次乘法运算和48次加法，这将耗费较多的时间，于是在H.264中，有一种改进的算法（蝶形算法）可以减少运算的次数。这种矩阵运算算法构造非常巧妙，利用构造的矩阵的整数性质和对称性，可完全将乘法运算转化为加法运算。如下图所示：

     下面的代码来自一本数字图像处理的书上的源代码：

VOID WINAPI FFT(complex<double> * TD, complex<double> * FD, int r)
{
	// 付立叶变换点数
	LONG count;

	// 循环变量
	int i,j,k;

	// 中间变量
	int bfsize,p;

	// 角度
	double angle;

	complex<double> *W,*X1,*X2,*X;

	// 计算付立叶变换点数
	count = 1 << r;

	// 分配运算所需存储器
	W = new complex<double>[count / 2];
	X1 = new complex<double>[count];
	X2 = new complex<double>[count];

	// 计算加权系数
	for(i = 0; i < count / 2; i++)
	{
		angle = -i * PI * 2 / count;
		W[i] = complex<double> (cos(angle), sin(angle));
	}

	// 将时域点写入X1
	memcpy(X1, TD, sizeof(complex<double>) * count);

	// 采用蝶形算法进行快速付立叶变换
	for(k = 0; k < r; k++)
	{
		for(j = 0; j < 1 << k; j++)
		{
			bfsize = 1 << (r-k);
			for(i = 0; i < bfsize / 2; i++)
			{
				p = j * bfsize;
				X2[i + p] = X1[i + p] + X1[i + p + bfsize / 2];
				X2[i + p + bfsize / 2] = (X1[i + p] - X1[i + p + bfsize / 2]) * W[i * (1<<k)];
			}
		}
		X = X1;
		X1 = X2;
		X2 = X;
	}

	// 重新排序
	for(j = 0; j < count; j++)
	{
		p = 0;
		for(i = 0; i < r; i++)
		{
			if (j&(1<<i))
			{
				p+=1<<(r-i-1);
			}
		}
		FD[j]=X1[p];
	}

	// 释放内存
	delete W;
	delete X1;
	delete X2;
}

    待我日后继续斟酌，补充。完。

    updated：关于快速傅立叶变换(FFT)的C++实现与Matlab实验，这里有一篇不错的文章，读者可以看看：http://blog.csdn.net/rappy/article/details/1700829。2012.03.03。

v_JULY_v
关注
专栏目录
多项式算法1：多项式乘法
alpha202的博客
06-02 5265
多项式简介 对于数域F\mathbb FF，若有∀i∈{1,2,3,⋯&ThinSpace;,n}\forall i\in\{1,2,3,\cdots,n \}∀i∈{1,2,3,⋯,n}，则 f(x)=a0+a1x+a2x2+⋯+anx=∑i=1naixif(x)=a_0+a_1x+a_2x^2+\cdots+a_nx=\sum_{i=1}^na_ix^if(x)=a0​+a1​x+a2...
快速傅里叶变换（FFT）求解多项式乘法
To be a better man
05-13 2597
前提：这里是从数学和算法的角度理解FFT，不涉及信号处理背景知识 已知系数表达式求点值表达式 这里n是2的整数幂，不够要补零。 核心：就是我①式计算前面2/n部分的结果可以保留下来用于计算后面2/n部分②式的计算，这样的话就可以把每个区间每一次进行等分成两份，就可以分成logn次，每次需要计算n次（比如A1(x)公式你看里面最大那一项是），最后时间复杂度就是，这样我们就在时间内得到了点值表达式 这样就可以每计算一个值需要递归logn层，计算n个值所以在 时间内得到了点值表达式.
快速傅里叶变换（FFT）及其在多项式乘法中的应用 —— 深入分析与 Python 实现
qq_52010229的博客
09-06 1543
单位根是指在复平面上均匀分布在单位圆上的点，满足 x^n = 1 的复数点。它们可以表示为：这些点的对称性和周期性使得我们可以通过分治法高效计算傅里叶变换。通过本文，我们深入探讨了如何利用快速傅里叶变换FFT进行高效的多项式乘法。我们从最基本的多项式操作出发，介绍了如何通过单位根和分治法将乘法复杂度从 O(n^2) 降低到 O(n log n)。我们还通过 Python 代码展示了 FFT 和 IFFT 的具体实现，并进行了多项式乘法的示例。FFT 在多项式乘法中的应用只是其众多用途之一。
数据结构——多项式的加减乘法以及合并同类项的算法
LJS021013的博客
11-18 1516
数据结构——加减乘法以及合并同类项算法
FFT与多项式乘法
Leo_h
06-08 1万+
鉴于网上关于FFT作用的文章并不少，这里尽量少说废话，直接说如何理解和具体实现FFT计算多项式乘法和大整数乘法
快速傅里叶变换(FFT)
Nerd_的博客
07-25 4726
快速傅里叶变换 问题引入
【FFT实战篇】C++实现：利用快速傅里叶变换快速计算（多项式）乘法
GoWilli的CSDN
03-22 3920
快速傅里叶变换FFT全流程解析，洛谷P3803 FFT多项式乘法
多项式乘法问题的快速傅里叶变换算法简述
Foliciatarier的博客
12-23 738
给定多项式 $A(x) = \Sigma_{i = 0}^n a_i x^i, B(x) = \Sigma_{i = 0}^m b_i x^i$ 的系数 $a_{0 .. n}, b_{0 .. m}$，求其多项式乘积 $C(x) = A(x) B(x) = \Sigma_{i = 0}^{m + n} c_i x^i$ 的各项系数 $c_{0 .. n + m}$
经典算法研究：多项式乘法与快速傅里叶变换
"多项式乘法与快速傅里叶变换-pfc 5.0 manual手册版" 本文档主要涵盖了多项式乘法与快速傅里叶变换（FFT）的算法原理和实现，作者通过一系列的文章深入浅出地介绍了这些经典算法。快速傅里叶变换是信号处理和计算...
多项式乘法与快速傅里叶变换
wenyusuran的专栏
07-07 1844
前言     经典算法研究系列，已经写到第十五章了，本章，咱们来介绍多项式的乘法以及快速傅里叶变换算法。本博客之前也已详细介绍过离散傅里叶变换（请参考：十、从头到尾彻底理解傅里叶变换算法、上，及十、从头到尾彻底理解傅里叶变换算法、下），这次咱们从多项式乘法开始，然后介绍FFT算法的原理与实现。同时，本文虽涉及到不少数学公式和定理（当然，我会尽量舍去一些与本文咱们要介绍的中心内容无关的定理或证明，
一元多项式乘法算法
03-18
NULL 博文链接：https://qingdaoguy.iteye.com/blog/694445
python123求多项式的值_大数乘法—多项式与快速傅里叶变换
weixin_33089993的博客
03-01 951
本章涉及知识点：1、多项式乘法的时间复杂度2、多项式的表示：系数3、多项式的表示：点值4、复数的表示5、单位复数根6、单位复数根的性质—消去引理7、单位复数根的性质—折半引理8、离散傅里叶变换：DFT和IDFT9、快速傅里叶变换：FFT10、FFT求解多项式乘法的步骤11、python编程实战FFT大数乘法12、结果分析一、多项式乘法的时间复杂度数学中，我们可以将任意一个n位的数字写为一个n-1次...
多项式快速幂——快速傅里叶变换
lxz6527的博客
06-11 861
快速傅里叶变换实现
多项式乘法傅立叶变换c++
12-09
多项式乘法傅立叶变换的好处是可以将多项式的乘法运算转化为卷积运算，从而可以通过快速傅立叶变换（FFT）算法来加速计算。在计算机科学领域中，快速傅立叶变换（FFT）是一种高效的算法，可以在较短的时间内完成乘法...
如何将两个多项式相乘
weixin_62675848的博客
10-18 1238
将P1当前项(ci,ei)乘P2多项式，再加到结果多项式中。
多项式乘法（FFT）
热门推荐
zhouyuheng2003的博客
11-21 2万+
多项式乘法（FFT） 前言 作为一名OI选手，至今未写过fft相关的博客，真是一大遗憾，这也导致我并没有真正推过fft的所有式子 这一篇fft的博客我将详细介绍多项式乘法，易于理解，当然，一些公式会有些枯燥，如果是初学者请耐心看完哦 介绍 多项式乘法，顾名思义，首先是讲到多项式，那么什么是多项式呢？ 多项式 首先是多项式的定义，想必大家都知道(你上过初中吧)，而在这里，我们所说的多项式都是单个未知...
【数据结构与算法分析】学习笔记课后答案第三章3.6 3.7多项式加法乘法不同时间复杂度
su_1998的博客
12-12 2005
第三章 表、栈和队列 3.6 编写两个多项式相加的函数 时间复杂度为O(M+N) 【数据结构与算法分析】学习笔记课后答案第三章3.6 3.7多项式加法乘法不同时间复杂度
多项式乘法计算
Lilyの博客 ٩(❛ัᴗ❛ั)
04-23 2345
这是一个多项式的乘法问题。 A(X)=amx^m+am-1x^m-1+.....a1x+a0； B(x)=bnx^n+bn-1x^n-1+.....b1x+b0 这两个多项式的项目分别为m+1和n+1，最高次分别为m和n， 现在需要求解A*B的值。
用C语言实现两个多项式相乘运算
最新发布
10-10
在C语言中，我们可以使用数组来表示多项式的系数，然后通过循环结构来实现两个多项式的乘法。由于多项式乘法涉及到每个系数的逐项相乘并累加，所以通常会创建一个新的数组来存储结果。这里是一个简单的示例： ```c #include <stdio.h> #define MAX_TERMS 50 // 用于存储最大项数 // 函数声明 void multiplyPolynomials(int* poly1, int* poly2, int* result, int n1, int n2); int main() { int poly1[MAX_TERMS], poly2[MAX_TERMS], result[MAX_TERMS * MAX_TERMS]; int n1 = 3, n2 = 4; // 第一个多项式的阶数和第二个多项式的阶数 int numTerms = 0; // 输入两个多项式的系数 printf("Enter coefficients for the first polynomial (separated by spaces): "); for (int i = 0; i < n1 + 1; i++) { scanf("%d", &poly1[i]); } printf("Enter coefficients for the second polynomial (separated by spaces): "); for (int i = 0; i < n2 + 1; i++) { scanf("%d", &poly2[i]); } // 乘法操作 multiplyPolynomials(poly1, poly2, result, n1, n2); // 输出结果多项式的系数 printf("\nProduct of the polynomials:\n"); for (int i = 0; i <= n1 + n2; i++) { if (result[i] != 0) { printf("%dx^%d + ", result[i], i); } numTerms++; if (numTerms % 5 == 0) { // 每五行显示一次 printf("\n"); } } return 0; } // 实现多项式乘法函数 void multiplyPolynomials(int* poly1, int* poly2, int* result, int n1, int n2) { for (int i = 0; i <= n1; i++) { for (int j = 0; j <= n2; j++) { int term = poly1[i] * poly2[j]; // 相乘得到临时结果 result[i+j] += term; // 累加到相应位置 } } } ``` 在这个例子中，`multiplyPolynomials`函数接受两个输入多项式的系数数组、结果数组以及它们的阶数，然后遍历每一个可能的乘积项并将它们加到结果数组对应的位置。
评论 21
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

v_JULY_v

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值