算法分析与设计上机实验指导书.docx

第 PAGE \* Arabic 16 页第 PAGE \* Arabic 16 页第 PAGE \* Arabic 16 页 天津商业大学信息工程学院 算法分析与设计 实验指导书 杨亮 2011/3/1 PAGEI / NUMPAGES23 目录 TOC \o 1-3 \h \z \u 实验1：大整数运算（分治策略） 1 程序1-1. 未使用分治策略 2 程序1-2. 使用分治策略 3 实验2：最长公共子序列（动态规划） 5 程序2-1 自底向上方法 6 实验3：图的最小生成树（贪心策略） 9 程序 3-1 Kruskal算法 13 程序 3-2 Prim算法 16 实验4：八皇后问题（解空间搜索方法） 19 程序4-1 回溯法与剪枝 20 算法分析与设计 实验指导书 PAGE 1 实验1：大整数运算（分治策略） 实验目的：通过大整数相乘的算法实现来加深对分治策略的了解。 技术描述：在大整数相乘的算法中，如果不使用分之策略所需的时间复杂度为O(n2)，而使用了动态规划思想进行重新设计后时间复杂度可以降低到O(n1.59)。 时间：100分钟 实验内容：编程完成相应的大整数加法，普通乘法和使用分治策略设计的大整数乘法，比较其相应的性能差别，并完成实验报告。 算法描述：设X和Y都是n位的二进制整数，现在要计算它们的乘积XY。我们可以用小学所学的方法来设计一个计算乘积XY的算法，但是这样做计算步骤太多，显得效率较低。如果将每2个1位数的乘法或加法看作一步运算，那么这种方法要作O(n2)步运算才能求出乘积XY。下面我们用 分治法来设计一个更有效的大整数乘积算法。 图1-1 大整数X和Y的分段　 我们将n位的二进制整数X和Y各分为2段，每段的长为n/2位(为简单起见，假设n是2的幂)，如图6-3所示。 由此，X=A2n/2+B ，Y=C2n/2+D。这样，X和Y的乘积为： XY=(A2n/2+B)(C2n/2+D)=AC2n+(AD+CB)2n/2+BD???? (1) 如果按式(1)计算XY，则我们必须进行4次n/2位整数的乘法(AC，AD，BC和BD)，以及3次不超过n位的整数加法(分别对应于式(1)中的加号)，此外还要做2次移位(分别对应于式(1)中乘2n和乘2n/2)。所有这些加法和移位共用O(n)步运算。设T(n)是2个n位整数相乘所需的运算总数，则由式(1)，我们有: ?????????????????????? (2) 由此可得T(n)=O(n2)。因此，用(1)式来计算X和Y的乘积并不比小学生的方法更有效。要想改进算法的计算复杂性，必须减少乘法次数。为此我们把XY写成另一种形式: XY=AC2n+[(A-B)(D-C)+AC+BD]2n/2+BD??? ???? (3) 虽然，式(3)看起来比式(1)复杂些，但它仅需做3次n/2位整数的乘法(AC，BD和(A-B)(D-C))，6次加、减法和2次移位。由此可得: ???????????????????????????????? (4) 用 解递归方程的 套用公式法马上可得其解为T(n)=O(nlog3)=O(n1.59) 程序1-1. 未使用分治策略 #include stdio.h #include string.h #define MAX_LEN 200 unsigned an1[MAX_LEN+10]; unsigned an2[MAX_LEN+10]; unsigned aResult[MAX_LEN* 2 + 10]; char szLine1[MAX_LEN+10]; char szLine2[MAX_LEN+10]; int main(){ gets( szLine1); //gets函数读取一行 gets( szLine2); int i, j; int nLen1 = strlen( szLine1); memset( an1, 0, sizeof(an1)); memset( an2, 0, sizeof(an2)); memset( aResult, 0, sizeof(aResult)); j = 0; for( i = nLen1 - 1;i = 0 ; i --) an1[j++] = szLine1[i] - 0; int nLen2 = strlen(szLine2); j = 0; for( i = nLen2 - 1;i = 0 ; i --) an2[j++] = szLine2[i] - 0; //每一轮都用an1的一位，去和an2各位相乘，从an1的个位开始 for( i = 0;