[转帖]常用算法设计方法

simon21

要使计算机能完成人们预定的工作，首先必须为如何完成预定的工作设计一个算法，然后再根据算法编写程序。计算机程序要对问题的每个对象和处理规则给出正确详尽的描述，其中程序的数据结构和变量用来描述问题的对象，程序结构、函数和语句用来描述问题的算法。算法数据结构是程序的两个重要方面。

算法是问题求解过程的精确描述，一个算法由有限条可完全机械地执行的、有确定结果的指令组成。指令正确地描述了要完成的任务和它们被执行的顺序。计算机按算法指令所描述的顺序执行算法的指令能在有限的步骤内终止，或终止于给出问题的解，或终止于指出问题对此输入数据无解。

通常求解一个问题可能会有多种算法可供选择，选择的主要标准是算法的正确性和可靠性，简单性和易理解性。其次是算法所需要的存储空间少和执行更快等。

算法设计是一件非常困难的工作，经常采用的算法设计技术主要有迭代法、穷举搜索法、递推法、贪婪法、回溯法、分治法、动态规划法等等。另外，为了更简洁的形式设计和藐视算法，在算法设计时又常常采用递归技术，用递归描述算法。

一、迭代法

迭代法是用于求方程或方程组近似根的一种常用的算法设计方法。设方程为f(x)=0，用某种数学方法导出等价的形式x=g(x)，然后按以下步骤执行：

（1） 选一个方程的近似根，赋给变量x0；

（2） 将x0的值保存于变量x1，然后计算g(x1)，并将结果存于变量x0；

（3） 当x0与x1的差的绝对值还小于指定的精度要求时，重复步骤（2）的计算。

若方程有根，并且用上述方法计算出来的近似根序列收敛，则按上述方法求得的x0就认为是方程的根。上述算法用C程序的形式表示为：

【算法】迭代法求方程的根

1. { x0=初始近似根；
   
2. do {
   
3. x1=x0；
   
4. x0=g(x1)； /*按特定的方程计算新的近似根*/
   
5. } while ( fabs(x0-x1)>Epsilon)；
   
6. printf(“方程的近似根是%f\n”，x0)；
   
7. }

复制代码

迭代算法也常用于求方程组的根，令

X=（x0，x1，…，xn-1）

设方程组为：

xi=gi(X) (I=0，1，…，n-1)

则求方程组根的迭代算法可描述如下：

【算法】迭代法求方程组的根

1. {for (i=0;i<n;i++)
   
2. x=初始近似根;
   
3. do {
   
4. for (i=0;i<n;i++)
   
5. y=x;
   
6. for (i=0;i<n;i++)
   
7. x=gi(X);
   
8. for (delta=0.0,i=0;i<n;i++)
   
9. if (fabs(y-x)>delta) delta=fabs(y-x)；
   
10. } while (delta>Epsilon)；
    
11. for (i=0;i<n;i++)
    
12. printf(“变量x[%d]的近似根是 %f”，I，x)；
    
13. printf(“\n”)；}

复制代码

具体使用迭代法求根时应注意以下两种可能发生的情况：

（1） 如果方程无解，算法求出的近似根序列就不会收敛，迭代过程会变成死循环，因此在使用迭代算法前应先考察方程是否有解，并在程序中对迭代的次数给予限制；

（2） 方程虽然有解，但迭代公式选择不当，或迭代的初始近似根选择不合理，也会导致迭代失败。

二、穷举搜索法

穷举搜索法是对可能是解的众多候选解按某种顺序进行逐一枚举和检验，并从众找出那些符合要求的候选解作为问题的解。

【问题】 将A、B、C、D、E、F这六个变量排成如图所示的三角形，这六个变量分别取[1，6]上的整数，且均不相同。求使三角形三条边上的变量之和相等的全部解。如图就是一个解。

程序引入变量a、b、c、d、e、f，并让它们分别顺序取1至6的证书，在它们互不相同的条件下，测试由它们排成的如图所示的三角形三条边上的变量之和是否相等，如相等即为一种满足要求的排列，把它们输出。当这些变量取尽所有的组合后，程序就可得到全部可能的解。细节见下面的程序。

【程序1】

1. # include <stdio.h>
   
2. void main() {
   
3.         int a,b,c,d,e,f;
   
4.         for (a=1;a<=6;a++)｛
   
5.                 for (b=1;b<=6;b++) {
   
6.                         if (b==a) continue;
   
7.                         for (c=1;c<=6;c++) {
   
8.                                 if (c==a)||(c==b) continue;
   
9.                                 for (d=1;d<=6;d++) {
   
10.                                         if (d==a)||(d==b)||(d==c) continue;
    
11.                                         for (e=1;e<=6;e++) {
    
12.                                                 if (e==a)||(e==b)||(e==c)||(e==d) continue;
    
13.                                                 f=21-(a+b+c+d+e);
    
14.                                                 if ((a+b+c==c+d+e))&&(a+b+c==e+f+a)) {
    
15.                                                         printf(“%6d,a);
    
16.                                                         printf(“%4d%4d”,b,f);
    
17.                                                         printf(“%2d%4d%4d”,c,d,e);
    
18.                                                         scanf(“%*c”);
    
19.                                                 }
    
20.                                         }
    
21.                                 }
    
22.                         }
    
23.                 }
    
24.         }
    
25. ｝

复制代码

按穷举法编写的程序通常不能适应变化的情况。如问题改成有9个变量排成三角形，每条边有4个变量的情况，程序的循环重数就要相应改变。

simon21

对一组数穷尽所有排列，还有更直接的方法。将一个排列看作一个长整数，则所有排列对应着一组整数。将这组整数按从小到大的顺序排列排成一个整数，从对应最小的整数开始。按数列的递增顺序逐一列举每个排列对应的每个整数，这能更有效地完成排列的穷举。从一个排列找出对应数列的下一个排列可在当前排列的基础上作部分调整来实现。倘若当前排列为1，2，4，6，5，3，并令其对应的长整数为124653。要寻找比长整数124653更大的排列，可从该排列的最后一个数字顺序向前逐位考察，当发现排列中的某个数字比它前一个数字大时，如本例中的6比它的前一位数字4大，这说明还有对应更大整数的排列。但为了顺序从小到大列举出所有的排列，不能立即调整得太大，如本例中将数字6与数字4交换得到的排列126453就不是排列124653的下一个排列。为了得到排列124653的下一个排列，应从已经考察过的那部分数字中选出比数字大，但又是它们中最小的那一个数字，比如数字5，与数字4交换。该数字也是从后向前考察过程中第一个比4大的数字。5与4交换后，得到排列125643。在前面数字1，2，5固定的情况下，还应选择对应最小整数的那个排列，为此还需将后面那部分数字的排列顺序颠倒，如将数字6，4，3的排列顺序颠倒，得到排列1，2，5，3，4，6，这才是排列1，2，4，6，5，3的下一个排列。按以上想法编写的程序如下。

【程序2】

1. # include <stdio.h>
   
2. # define SIDE_N 3
   
3. # define LENGTH 3
   
4. # define VARIABLES 6
   
5. int A,B,C,D,E,F;
   
6. int *pt[]={&A,&B,&C,&D,&E,&F};
   
7. int *side[SIDE_N][LENGTH]={&A,&B,&C,&C,&D,&E,&E,&F,&A};
   
8. int side_total[SIDE_N];
   
9. main{}
   
10. { int i,j,t,equal;
    
11. for (j=0;j<VARIABLES;j++)
    
12. *pt[j]=j+1;
    
13. while(1)
    
14. { for (i=0;i<SIDE_N;i++)
    
15. { for (t=j=0;j<LENGTH;j++)
    
16. t+=*side[j];
    
17. side_total=t;
    
18. }
    
19. for (equal=1,i=0;equal&&i<SIDE_N-1;i++)
    
20. if (side_total!=side_total[i+1] equal=0;
    
21. if (equal)
    
22. { for (i=1;i<VARIABLES;i++)
    
23. printf(“%4d”,*pt);
    
24. printf(“\n”);
    
25. scanf(“%*c”);
    
26. }
    
27. for (j=VARIABLES-1;j>0;j--)
    
28. if (*pt[j]>*pt[j-1]) break;
    
29. if (j==0) break;
    
30. for (i=VARIABLES-1;i>=j;i--)
    
31. if (*pt>*pt[i-1]) break;
    
32. t=*pt[j-1];* pt[j-1] =* pt; *pt=t;
    
33. for (i=VARIABLES-1;i>j;i--,j++)
    
34. { t=*pt[j]; *pt[j] =* pt; *pt=t; }
    
35. }
    
36. }

复制代码

从上述问题解决的方法中，最重要的因素就是确定某种方法来确定所有的候选解。下面再用一个示例来加以说明。

【问题】 背包问题

问题描述：有不同价值、不同重量的物品n件，求从这n件物品中选取一部分物品的选择方案，使选中物品的总重量不超过指定的限制重量，但选中物品的价值之和最大。

设n个物品的重量和价值分别存储于数组w[ ]和v[ ]中，限制重量为tw。考虑一个n元组（x0，x1，…，xn-1），其中xi=0 表示第i个物品没有选取，而xi=1则表示第i个物品被选取。显然这个n元组等价于一个选择方案。用枚举法解决背包问题，需要枚举所有的选取方案，而根据上述方法，我们只要枚举所有的n元组，就可以得到问题的解。

显然，每个分量取值为0或1的n元组的个数共为2n个。而每个n元组其实对应了一个长度为n的二进制数，且这些二进制数的取值范围为0～2n-1。因此，如果把0～2n-1分别转化为相应的二进制数，则可以得到我们所需要的2n个n元组。

【算法】

1. maxv=0;
   
2. for (i=0;i<2n;i++)
   
3. { B[0..n-1]=0;
   
4. 把i转化为二进制数，存储于数组B中;
   
5. temp_w=0;
   
6. temp_v=0;
   
7. for (j=0;j<n;j++)
   
8. { if (B[j]==1)
   
9. { temp_w=temp_w+w[j];
   
10. temp_v=temp_v+v[j];
    
11. }
    
12. if ((temp_w<=tw)&&(temp_v>maxv))
    
13. { maxv=temp_v;
    
14. 保存该B数组；
    
15. }
    
16. }
    
17. }

复制代码

三、递推法

递推法是利用问题本身所具有的一种递推关系求问题解的一种方法。设要求问题规模为N的解，当N=1时，解或为已知，或能非常方便地得到解。能采用递推法构造算法的问题有重要的递推性质，即当得到问题规模为i-1的解后，由问题的递推性质，能从已求得的规模为1，2，…，i-1的一系列解，构造出问题规模为I的解。这样，程序可从i=0或i=1出发，重复地，由已知至i-1规模的解，通过递推，获得规模为i的解，直至得到规模为N的解。

simon21

【问题】 阶乘计算

问题描述：编写程序，对给定的n（n≦100），计算并输出k的阶乘k！（k=1，2，…，n）的全部有效数字。

由于要求的整数可能大大超出一般整数的位数，程序用一维数组存储长整数，存储长整数数组的每个元素只存储长整数的一位数字。如有m位成整数N用数组a[ ]存储：

N=a[m]×10m-1+a[m-1]×10m-2+ … +a[2]×101+a[1]×100

并用a[0]存储长整数N的位数m，即a[0]=m。按上述约定，数组的每个元素存储k的阶乘k！的一位数字，并从低位到高位依次存于数组的第二个元素、第三个元素……。例如，5！=120，在数组中的存储形式为：

3
0
2
1
……

首元素3表示长整数是一个3位数，接着是低位到高位依次是0、2、1，表示成整数120。

计算阶乘k！可采用对已求得的阶乘(k-1)！连续累加k-1次后求得。例如，已知4！=24，计算5！，可对原来的24累加4次24后得到120。细节见以下程序。

1. # include <stdio.h>
   
2. # include <malloc.h>
   
3. # define MAXN 1000
   
4. 
   
5. void pnext(int a[ ],int k)
   
6. { int *b,m=a[0],i,j,r,carry;
   
7. b=(int * ) malloc(sizeof(int)* (m+1));
   
8. for ( i=1;i<=m;i++) b=a;
   
9. for ( j=1;j<=k;j++)
   
10. { for ( carry=0,i=1;i<=m;i++)
    
11. { r=(i<a[0]?a+b:a)+carry;
    
12. a=r%10;
    
13. carry=r/10;
    
14. }
    
15. if (carry) a[++m]=carry;
    
16. }
    
17. free(b);
    
18. a[0]=m;
    
19. }
    
20. 
    
21. void write(int *a,int k)
    
22. { int i;
    
23. printf(“%4d！=”,k);
    
24. for (i=a[0];i>0;i--)
    
25. printf(“%d”,a);
    
26. printf(“\n\n”);
    
27. }
    
28. 
    
29. void main()
    
30. { int a[MAXN],n,k;
    
31. printf(“Enter the number n: “);
    
32. scanf(“%d”,&n);
    
33. a[0]=1;
    
34. a[1]=1;
    
35. write(a,1);
    
36. for (k=2;k<=n;k++)
    
37. { pnext(a,k);
    
38. write(a,k);
    
39. getchar();
    
40. }
    
41. }

复制代码

四、递归

递归是设计和描述算法的一种有力的工具，由于它在复杂算法的描述中被经常采用，为此在进一步介绍其他算法设计方法之前先讨论它。

能采用递归描述的算法通常有这样的特征：为求解规模为N的问题，设法将它分解成规模较小的问题，然后从这些小问题的解方便地构造出大问题的解，并且这些规模较小的问题也能采用同样的分解和综合方法，分解成规模更小的问题，并从这些更小问题的解构造出规模较大问题的解。特别地，当规模N=1时，能直接得解。

【问题】 编写计算斐波那契（Fibonacci）数列的第n项函数fib（n）。

斐波那契数列为：0、1、1、2、3、……，即：

fib(0)=0;
fib(1)=1;
fib(n)=fib(n-1)+fib(n-2) （当n>1时）。
写成递归函数有：

1. int fib(int n)
   
2. { if (n==0) return 0;
   
3. if (n==1) return 1;
   
4. if (n>1) return fib(n-1)+fib(n-2);
   
5. }

复制代码

[世纪青年资料站 FW88.COM]

递归算法的执行过程分递推和回归两个阶段。在递推阶段，把较复杂的问题（规模为n）的求解推到比原问题简单一些的问题（规模小于n）的求解。例如上例中，求解fib(n)，把它推到求解fib(n-1)和fib(n-2)。也就是说，为计算fib(n)，必须先计算fib(n-1)和fib(n-2)，而计算fib(n-1)和fib(n-2)，又必须先计算fib(n-3)和fib(n-4)。依次类推，直至计算fib(1)和fib(0)，分别能立即得到结果1和0。在递推阶段，必须要有终止递归的情况。例如在函数fib中，当n为1和0的情况。

在回归阶段，当获得最简单情况的解后，逐级返回，依次得到稍复杂问题的解，例如得到fib(1)和fib(0)后，返回得到fib(2)的结果，……，在得到了fib(n-1)和fib(n-2)的结果后，返回得到fib(n)的结果。

在编写递归函数时要注意，函数中的局部变量和参数知识局限于当前调用层，当递推进入“简单问题”层时，原来层次上的参数和局部变量便被隐蔽起来。在一系列“简单问题”层，它们各有自己的参数和局部变量。

由于递归引起一系列的函数调用，并且可能会有一系列的重复计算，递归算法的执行效率相对较低。当某个递归算法能较方便地转换成递推算法时，通常按递推算法编写程序。例如上例计算斐波那契数列的第n项的函数fib(n)应采用递推算法，即从斐波那契数列的前两项出发，逐次由前两项计算出下一项，直至计算出要求的第n项。

simon21

【问题】 组合问题

问题描述：找出从自然数1、2、……、n中任取r个数的所有组合。例如n=5，r=3的所有组合为： （1）5、4、3 （2）5、4、2 （3）5、4、1

（4）5、3、2 （5）5、3、1 （6）5、2、1

（7）4、3、2 （8）4、3、1 （9）4、2、1

（10）3、2、1

分析所列的10个组合，可以采用这样的递归思想来考虑求组合函数的算法。设函数为void comb(int m,int k)为找出从自然数1、2、……、m中任取k个数的所有组合。当组合的第一个数字选定时，其后的数字是从余下的m-1个数中取k-1数的组合。这就将求m个数中取k个数的组合问题转化成求m-1个数中取k-1个数的组合问题。设函数引入工作数组a[ ]存放求出的组合的数字，约定函数将确定的k个数字组合的第一个数字放在a[k]中，当一个组合求出后，才将a[ ]中的一个组合输出。第一个数可以是m、m-1、……、k，函数将确定组合的第一个数字放入数组后，有两种可能的选择，因还未去顶组合的其余元素，继续递归去确定；或因已确定了组合的全部元素，输出这个组合。细节见以下程序中的函数comb。

【程序】

1. # include <stdio.h>
   
2. # define MAXN 100
   
3. 
   
4. int a[MAXN];
   
5. void comb(int m,int k)
   
6. { int i,j;
   
7. for (i=m;i>=k;i--)
   
8. { a[k]=i;
   
9. if (k>1)
   
10. comb(i-1,k-1);
    
11. else
    
12. { for (j=a[0];j>0;j--)
    
13. printf(“%4d”,a[j]);
    
14. printf(“\n”);
    
15. }
    
16. }
    
17. }
    
18. 
    
19. void main()
    
20. { a[0]=3;
    
21. comb(5,3);
    
22. }

复制代码

【问题】 背包问题

问题描述：有不同价值、不同重量的物品n件，求从这n件物品中选取一部分物品的选择方案，使选中物品的总重量不超过指定的限制重量，但选中物品的价值之和最大。

设n件物品的重量分别为w0、w1、…、wn-1，物品的价值分别为v0、v1、…、vn-1。采用递归寻找物品的选择方案。设前面已有了多种选择的方案，并保留了其中总价值最大的方案于数组option[ ]，该方案的总价值存于变量maxv。当前正在考察新方案，其物品选择情况保存于数组cop[ ]。假定当前方案已考虑了前i-1件物品，现在要考虑第i件物品；当前方案已包含的物品的重量之和为tw；至此，若其余物品都选择是可能的话，本方案能达到的总价值的期望值为tv。算法引入tv是当一旦当前方案的总价值的期望值也小于前面方案的总价值maxv时，继续考察当前方案变成无意义的工作，应终止当前方案，立即去考察下一个方案。因为当方案的总价值不比maxv大时，该方案不会被再考察，这同时保证函数后找到的方案一定会比前面的方案更好。

对于第i件物品的选择考虑有两种可能：

（1） 考虑物品i被选择，这种可能性仅当包含它不会超过方案总重量限制时才是可行的。选中后，继续递归去考虑其余物品的选择。

（2） 考虑物品i不被选择，这种可能性仅当不包含物品i也有可能会找到价值更大的方案的情况。

按以上思想写出递归算法如下：

try(物品i，当前选择已达到的重量和，本方案可能达到的总价值tv)

{ /*考虑物品i包含在当前方案中的可能性*/

if(包含物品i是可以接受的)

{ 将物品i包含在当前方案中；

if (i<n-1)

try(i+1,tw+物品i的重量,tv);

else

/*又一个完整方案，因为它比前面的方案好，以它作为最佳方案*/

以当前方案作为临时最佳方案保存;

恢复物品i不包含状态；

}

/*考虑物品i不包含在当前方案中的可能性*/

if (不包含物品i仅是可男考虑的)

if (i<n-1)

try(i+1,tw,tv-物品i的价值)；

else

/*又一个完整方案，因它比前面的方案好，以它作为最佳方案*/

以当前方案作为临时最佳方案保存;

}

为了理解上述算法，特举以下实例。设有4件物品，它们的重量和价值见表：

物品	重量	价值
0	5	4
1	3	4
2	2	3
3	1	1

并设限制重量为7。则按以上算法，下图表示找解过程。由图知，一旦找到一个解，算法就进一步找更好的佳。如能判定某个查找分支不会找到更好的解，算法不会在该分支继续查找，而是立即终止该分支，并去考察下一个分支。

Try(0,0,12)

Try(1,5,12)

Try(1,0,8)

Try(2,5,8)

Try(3,7,8)

Try(2,3,8)

Try(3,5,8)

不能得到更好的解

不能得到更好的解

超重

不能得到更好的解

得到解：(1,0,1,0)

maxv=7

得到解：(0,1,1,1)

maxv=8

不能得到更好的解

超重

simon21

按上述算法编写函数和程序如下：

【程序】

1. # include <stdio.h>
   
2. # define N 100
   
3. 
   
4. double limitW,totV,maxV;
   
5. int option[N],cop[N];
   
6. struct { double weight;
   
7. double value;
   
8. }a[N];
   
9. int n;
   
10. 
    
11. void find(int i,double tw,double tv)
    
12. { int k;
    
13. /*考虑物品i包含在当前方案中的可能性*/
    
14. if (tw+a.weight<=limitW)
    
15. { cop=1;
    
16. if (i<n-1) find(i+1,tw+a.weight,tv);
    
17. else
    
18. { for (k=0;k<n;k++)
    
19. option[k]=cop[k];
    
20. maxv=tv;
    
21. }
    
22. cop=0;
    
23. }
    
24. /*考虑物品i不包含在当前方案中的可能性*/
    
25. if (tv-a.value>maxV)
    
26. if (i<n-1) find(i+1,tw,tv-a.value);
    
27. else
    
28. { for (k=0;k<n;k++)
    
29. option[k]=cop[k];
    
30. maxv=tv-a.value;
    
31. }
    
32. }
    
33. 
    
34. void main()
    
35. { int k;
    
36. double w,v;
    
37. printf(“输入物品种数\n”);
    
38. scanf((“%d”,&n);
    
39. printf(“输入各物品的重量和价值\n”);
    
40. for (totv=0.0,k=0;k<n;k++)
    
41. { scanf(“%1f%1f”,&w,&v);
    
42. a[k].weight=w;
    
43. a[k].value=v;
    
44. totV+=V;
    
45. }
    
46. printf(“输入限制重量\n”);
    
47. scanf(“%1f”,&limitV);
    
48. maxv=0.0;
    
49. for (k=0;k<n;k++) cop[k]=0;
    
50. find(0,0.0,totV);
    
51. for (k=0;k<n;k++)
    
52. if (option[k]) printf(“%4d”,k+1);
    
53. printf(“\n总价值为%.2f\n”,maxv);
    
54. }

复制代码

作为对比，下面以同样的解题思想，考虑非递归的程序解。为了提高找解速度，程序不是简单地逐一生成所有候选解，而是从每个物品对候选解的影响来形成值得进一步考虑的候选解，一个候选解是通过依次考察每个物品形成的。对物品i的考察有这样几种情况：当该物品被包含在候选解中依旧满足解的总重量的限制，该物品被包含在候选解中是应该继续考虑的；反之，该物品不应该包括在当前正在形成的候选解中。同样地，仅当物品不被包括在候选解中，还是有可能找到比目前临时最佳解更好的候选解时，才去考虑该物品不被包括在候选解中；反之，该物品不包括在当前候选解中的方案也不应继续考虑。对于任一值得继续考虑的方案，程序就去进一步考虑下一个物品。

【程序】

1. # include <stdio.h>
   
2. # define N 100
   
3. 
   
4. double limitW;
   
5. int cop[N];
   
6. struct ele { double weight;
   
7. double value;
   
8. } a[N];
   
9. int k,n;
   
10. struct { int flg;
    
11. double tw;
    
12. double tv;
    
13. }twv[N];
    
14. void next(int i,double tw,double tv)
    
15. { twv.flg=1;
    
16. twv.tw=tw;
    
17. twv.tv=tv;
    
18. }
    
19. double find(struct ele *a,int n)
    
20. { int i,k,f;
    
21. double maxv,tw,tv,totv;
    
22. maxv=0;
    
23. for (totv=0.0,k=0;k<n;k++)
    
24. totv+=a[k].value;
    
25. next(0,0.0,totv);
    
26. i=0;
    
27. While (i>=0)
    
28. { f=twv.flg;
    
29. tw=twv.tw;
    
30. tv=twv.tv;
    
31. switch(f)
    
32. { case 1: twv.flg++;
    
33. if (tw+a.weight<=limitW)
    
34. if (i<n-1)
    
35. { next(i+1,tw+a.weight,tv);
    
36. i++;
    
37. }
    
38. else
    
39. { maxv=tv;
    
40. for (k=0;k<n;k++)
    
41. cop[k]=twv[k].flg!=0;
    
42. }
    
43. break;
    
44. case 0: i--;
    
45. break;
    
46. default: twv.flg=0;
    
47. if (tv-a.value>maxv)
    
48. if (i<n-1)
    
49. { next(i+1,tw,tv-a.value);
    
50. i++;
    
51. }
    
52. else
    
53. { maxv=tv-a.value;
    
54. for (k=0;k<n;k++)
    
55. cop[k]=twv[k].flg!=0;
    
56. }
    
57. break;
    
58. }
    
59. }
    
60. return maxv;
    
61. }
    
62. 
    
63. void main()
    
64. { double maxv;
    
65. printf(“输入物品种数\n”);
    
66. scanf((“%d”,&n);
    
67. printf(“输入限制重量\n”);
    
68. scanf(“%1f”,&limitW);
    
69. printf(“输入各物品的重量和价值\n”);
    
70. for (k=0;k<n;k++)
    
71. scanf(“%1f%1f”,&a[k].weight,&a[k].value);
    
72. maxv=find(a,n);
    
73. printf(“\n选中的物品为\n”);
    
74. for (k=0;k<n;k++)
    
75. if (option[k]) printf(“%4d”,k+1);
    
76. printf(“\n总价值为%.2f\n”,maxv);
    
77. }

复制代码

五、回溯法

回溯法也称为试探法，该方法首先暂时放弃关于问题规模大小的限制，并将问题的候选解按某种顺序逐一枚举和检验。当发现当前候选解不可能是解时，就选择下一个候选解；倘若当前候选解除了还不满足问题规模要求外，满足所有其他要求时，继续扩大当前候选解的规模，并继续试探。如果当前候选解满足包括问题规模在内的所有要求时，该候选解就是问题的一个解。在回溯法中，放弃当前候选解，寻找下一个候选解的过程称为回溯。扩大当前候选解的规模，以继续试探的过程称为向前试探。

simon21

1、回溯法的一般描述
可用回溯法求解的问题P，通常要能表达为：对于已知的由n元组（x1，x2，…，xn）组成的一个状态空间E={（x1，x2，…，xn）∣xi∈Si ，i=1，2，…，n}，给定关于n元组中的一个分量的一个约束集D，要求E中满足D的全部约束条件的所有n元组。其中Si是分量xi的定义域，且 |Si| 有限，i=1，2，…，n。我们称E中满足D的全部约束条件的任一n元组为问题P的一个解。

解问题P的最朴素的方法就是枚举法，即对E中的所有n元组逐一地检测其是否满足D的全部约束，若满足，则为问题P的一个解。但显然，其计算量是相当大的。

我们发现，对于许多问题，所给定的约束集D具有完备性，即i元组（x1，x2，…，xi）满足D中仅涉及到x1，x2，…，xi的所有约束意味着j（j<i）元组（x1，x2，…，xj）一定也满足D中仅涉及到x1，x2，…，xj的所有约束，i=1，2，…，n。换句话说，只要存在0≤j≤n-1，使得（x1，x2，…，xj）违反D中仅涉及到x1，x2，…，xj的约束之一，则以（x1，x2，…，xj）为前缀的任何n元组（x1，x2，…，xj，xj+1，…，xn）一定也违反D中仅涉及到x1，x2，…，xi的一个约束，n≥i>j。因此，对于约束集D具有完备性的问题P，一旦检测断定某个j元组（x1，x2，…，xj）违反D中仅涉及x1，x2，…，xj的一个约束，就可以肯定，以（x1，x2，…，xj）为前缀的任何n元组（x1，x2，…，xj，xj+1，…，xn）都不会是问题P的解，因而就不必去搜索它们、检测它们。回溯法正是针对这类问题，利用这类问题的上述性质而提出来的比枚举法效率更高的算法。

回溯法首先将问题P的n元组的状态空间E表示成一棵高为n的带权有序树T，把在E中求问题P的所有解转化为在T中搜索问题P的所有解。树T类似于检索树，它可以这样构造：

设Si中的元素可排成xi(1) ，xi(2) ，…，xi(mi-1) ，|Si| =mi，i=1，2，…，n。从根开始，让T的第I层的每一个结点都有mi个儿子。这mi个儿子到它们的双亲的边，按从左到右的次序，分别带权xi+1(1) ，xi+1(2) ，…，xi+1(mi) ，i=0，1，2，…，n-1。照这种构造方式，E中的一个n元组（x1，x2，…，xn）对应于T中的一个叶子结点，T的根到这个叶子结点的路径上依次的n条边的权分别为x1，x2，…，xn，反之亦然。另外，对于任意的0≤i≤n-1，E中n元组（x1，x2，…，xn）的一个前缀I元组（x1，x2，…，xi）对应于T中的一个非叶子结点，T的根到这个非叶子结点的路径上依次的I条边的权分别为x1，x2，…，xi，反之亦然。特别，E中的任意一个n元组的空前缀（），对应于T的根。

因而，在E中寻找问题P的一个解等价于在T中搜索一个叶子结点，要求从T的根到该叶子结点的路径上依次的n条边相应带的n个权x1，x2，…，xn满足约束集D的全部约束。在T中搜索所要求的叶子结点，很自然的一种方式是从根出发，按深度优先的策略逐步深入，即依次搜索满足约束条件的前缀1元组（x1i）、前缀2元组（x1，x2）、…，前缀I元组（x1，x2，…，xi），…，直到i=n为止。

在回溯法中，上述引入的树被称为问题P的状态空间树；树T上任意一个结点被称为问题P的状态结点；树T上的任意一个叶子结点被称为问题P的一个解状态结点；树T上满足约束集D的全部约束的任意一个叶子结点被称为问题P的一个回答状态结点，它对应于问题P的一个解。

【问题】 组合问题

问题描述：找出从自然数1、2、……、n中任取r个数的所有组合。

例如n=5，r=3的所有组合为：

（1）1、2、3 （2）1、2、4 （3）1、2、5

（4）1、3、4 （5）1、3、5 （6）1、4、5

（7）2、3、4 （8）2、3、5 （9）2、4、5

（10）3、4、5

则该问题的状态空间为：

E={（x1，x2，x3）∣xi∈S ，i=1，2，3 } 其中：S={1，2，3，4，5}

约束集为： x1<x2<x3

显然该约束集具有完备性。

问题的状态空间树T：

⊙ ⊙ ⊙ ⊙ ⊙

⊙ ⊙ ⊙ ⊙ ⊙ ⊙ ⊙ ⊙ ⊙ ⊙

simon21

2、回溯法的方法

对于具有完备约束集D的一般问题P及其相应的状态空间树T，利用T的层次结构和D的完备性，在T中搜索问题P的所有解的回溯法可以形象地描述为：

从T的根出发，按深度优先的策略，系统地搜索以其为根的子树中可能包含着回答结点的所有状态结点，而跳过对肯定不含回答结点的所有子树的搜索，以提高搜索效率。具体地说，当搜索按深度优先策略到达一个满足D中所有有关约束的状态结点时，即“激活”该状态结点，以便继续往深层搜索；否则跳过对以该状态结点为根的子树的搜索，而一边逐层地向该状态结点的祖先结点回溯，一边“杀死”其儿子结点已被搜索遍的祖先结点，直到遇到其儿子结点未被搜索遍的祖先结点，即转向其未被搜索的一个儿子结点继续搜索。

在搜索过程中，只要所激活的状态结点又满足终结条件，那么它就是回答结点，应该把它输出或保存。由于在回溯法求解问题时，一般要求出问题的所有解，因此在得到回答结点后，同时也要进行回溯，以便得到问题的其他解，直至回溯到T的根且根的所有儿子结点均已被搜索过为止。

例如在组合问题中，从T的根出发深度优先遍历该树。当遍历到结点（1，2）时，虽然它满足约束条件，但还不是回答结点，则应继续深度遍历；当遍历到叶子结点（1，2，5）时，由于它已是一个回答结点，则保存（或输出）该结点，并回溯到其双亲结点，继续深度遍历；当遍历到结点（1，5）时，由于它已是叶子结点，但不满足约束条件，故也需回溯。

3、回溯法的一般流程和技术

在用回溯法求解有关问题的过程中，一般是一边建树，一边遍历该树。在回溯法中我们一般采用非递归方法。下面，我们给出回溯法的非递归算法的一般流程：
NYNYNNYY
建立根结点root
建立root的第一个孩子结点node
建树完毕？
Node是叶子？
Node是解？
处理解
回溯node=parent(node)
Node还有孩子？
建立node的孩子结点node=parent(node)
建立node的孩子结点node=parent(node)
结束
开始

在用回溯法求解问题，也即在遍历状态空间树的过程中，如果采用非递归方法，则我们一般要用到栈的数据结构。这时，不仅可以用栈来表示正在遍历的树的结点，而且可以很方便地表示建立孩子结点和回溯过程。

例如在组合问题中，我们用一个一维数组Stack[ ]表示栈。开始栈空，则表示了树的根结点。如果元素1进栈，则表示建立并遍历（1）结点；这时如果元素2进栈，则表示建立并遍历（1，2）结点；元素3再进栈，则表示建立并遍历（1，2，3）结点。这时可以判断它满足所有约束条件，是问题的一个解，输出（或保存）。这时只要栈顶元素（3）出栈，即表示从结点（1，2，3）回溯到结点（1，2）。

【问题】 组合问题

问题描述：找出从自然数1，2，…，n中任取r个数的所有组合。

采用回溯法找问题的解，将找到的组合以从小到大顺序存于a[0]，a[1]，…，a[r-1]中，组合的元素满足以下性质：

（1） a[i+1]>a，后一个数字比前一个大；

（2） a-i<=n-r+1。

按回溯法的思想，找解过程可以叙述如下：

首先放弃组合数个数为r的条件，候选组合从只有一个数字1开始。因该候选解满足除问题规模之外的全部条件，扩大其规模，并使其满足上述条件（1），候选组合改为1，2。继续这一过程，得到候选组合1，2，3。该候选解满足包括问题规模在内的全部条件，因而是一个解。在该解的基础上，选下一个候选解，因a[2]上的3调整为4，以及以后调整为5都满足问题的全部要求，得到解1，2，4和1，2，5。由于对5不能再作调整，就要从a[2]回溯到a[1]，这时，a[1]=2，可以调整为3，并向前试探，得到解1，3，4。重复上述向前试探和向后回溯，直至要从a[0]再回溯时，说明已经找完问题的全部解。按上述思想写成程序如下：

【程序】

1. # define MAXN 100
   
2. int a[MAXN];
   
3. void comb(int m,int r)
   
4. { int i,j;
   
5. i=0;
   
6. a=1;
   
7. do {
   
8. if (a-i<=m-r+1
   
9. { if (i==r-1)
   
10. { for (j=0;j<r;j++)
    
11. printf(“%4d”,a[j]);
    
12. printf(“\n”);
    
13. }
    
14. a++;
    
15. continue;
    
16. }
    
17. else
    
18. { if (i==0)
    
19. return;
    
20. a[--i]++;
    
21. }
    
22. } while (1)
    
23. }
    
24. 
    
25. main()
    
26. { comb(5,3);
    
27. }

复制代码

simon21

【问题】 填字游戏

问题描述：在3×3个方格的方阵中要填入数字1到N（N≥10）内的某9个数字，每个方格填一个整数，似的所有相邻两个方格内的两个整数之和为质数。试求出所有满足这个要求的各种数字填法。

可用试探发找到问题的解，即从第一个方格开始，为当前方格寻找一个合理的整数填入，并在当前位置正确填入后，为下一方格寻找可填入的合理整数。如不能为当前方格找到一个合理的可填证书，就要回退到前一方格，调整前一方格的填入数。当第九个方格也填入合理的整数后，就找到了一个解，将该解输出，并调整第九个的填入的整数，寻找下一个解。

为找到一个满足要求的9个数的填法，从还未填一个数开始，按某种顺序（如从小到大的顺序）每次在当前位置填入一个整数，然后检查当前填入的整数是否能满足要求。在满足要求的情况下，继续用同样的方法为下一方格填入整数。如果最近填入的整数不能满足要求，就改变填入的整数。如对当前方格试尽所有可能的整数，都不能满足要求，就得回退到前一方格，并调整前一方格填入的整数。如此重复执行扩展、检查或调整、检查，直到找到一个满足问题要求的解，将解输出。

回溯法找一个解的算法：

1. { int m=0,ok=1;
   
2. int n=8;
   
3. do{
   
4. if (ok) 扩展;
   
5. else 调整;
   
6. ok=检查前m个整数填放的合理性;
   
7. } while ((!ok||m!=n)&&(m!=0))
   
8. if (m!=0) 输出解;
   
9. else 输出无解报告；
   
10. }

复制代码

如果程序要找全部解，则在将找到的解输出后，应继续调整最后位置上填放的整数，试图去找下一个解。相应的算法如下：

回溯法找全部解的算法：

1. { int m=0,ok=1;
   
2. int n=8;
   
3. do{
   
4. if (ok)
   
5. { if (m==n)
   
6. { 输出解；
   
7. 调整；
   
8. }
   
9. else 扩展;
   
10. }
    
11. else 调整;
    
12. ok=检查前m个整数填放的合理性;
    
13. } while (m!=0);
    
14. }

复制代码

为了确保程序能够终止，调整时必须保证曾被放弃过的填数序列不会再次实验，即要求按某种有许模型生成填数序列。给解的候选者设定一个被检验的顺序，按这个顺序逐一形成候选者并检验。从小到大或从大到小，都是可以采用的方法。如扩展时，先在新位置填入整数1，调整时，找当前候选解中下一个还未被使用过的整数。将上述扩展、调整、检验都编写成程序，细节见以下找全部解的程序。

【程序】

1. # include <stdio.h>
   
2. # define N 12
   
3. 
   
4. void write(int a[ ])
   
5. { int i,j;
   
6. for (i=0;i<3;i++)
   
7. { for (j=0;j<3;j++)
   
8. printf(“%3d”,a[3*i+j]);
   
9. printf(“\n”);
   
10. }
    
11. scanf(“%*c”);
    
12. }
    
13. int b[N+1];
    
14. int a[10];
    
15. int isprime(int m)
    
16. { int i;
    
17. int primes[ ]={2,3,5,7,11,17,19,23,29,-1};
    
18. if (m==1||m%2=0) return 0;
    
19. for (i=0;primes>0;i++)
    
20. if (m==primes) return 1;
    
21. for (i=3;i*i<=m;)
    
22. { if (m%i==0) return 0;
    
23. i+=2;
    
24. }
    
25. return 1;
    
26. }
    
27. int checkmatrix[ ][3]={ {-1},{0,-1},{1,-1},{0,-1},{1,3,-1},{2,4,-1},{3,-1},{4,6,-1},{5,7,-1}};
    
28. int selectnum(int start)
    
29. { int j;
    
30. for (j=start;j<=N;j++)
    
31. if (b[j]) return j
    
32. return 0;
    
33. }
    
34. 
    
35. int check(int pos)
    
36. { int i,j;
    
37. if (pos<0) return 0;
    
38. for (i=0;(j=checkmatrix[pos])>=0;i++)
    
39. if (!isprime(a[pos]+a[j])
    
40. return 0;
    
41. return 1;
    
42. }
    
43. 
    
44. int extend(int pos)
    
45. { a[++pos]=selectnum(1);
    
46. b[a][pos]]=0;
    
47. return pos;
    
48. }
    
49. 
    
50. int change(int pos)
    
51. { int j;
    
52. while (pos>=0&&(j=selectnum(a[pos]+1))==0)
    
53. b[a[pos--]]=1;
    
54. if (pos<0) return –1
    
55. b[a[pos]]=1;
    
56. a[pos]=j;
    
57. b[j]=0;
    
58. return pos;
    
59. }
    
60. 
    
61. void find()
    
62. { int ok=0,pos=0;
    
63. a[pos]=1;
    
64. b[a[pos]]=0;
    
65. do {
    
66. if (ok)
    
67. if (pos==8)
    
68. { write(a);
    
69. pos=change(pos);
    
70. }
    
71. else pos=extend(pos);
    
72. else pos=change(pos);
    
73. ok=check(pos);
    
74. } while (pos>=0)
    
75. }
    
76. 
    
77. void main()
    
78. { int i;
    
79. for (i=1;i<=N;i++)
    
80. b=1;
    
81. find();
    
82. }

复制代码

simon21

【问题】 n皇后问题

问题描述：求出在一个n×n的棋盘上，放置n个不能互相捕捉的国际象棋“皇后”的所有布局。

这是来源于国际象棋的一个问题。皇后可以沿着纵横和两条斜线4个方向相互捕捉。如图所示，一个皇后放在棋盘的第4行第3列位置上，则棋盘上凡打“×”的位置上的皇后就能与这个皇后相互捕捉。

1	×	×	×	×
2	×	×	×	×
3	×	×	×	×
4	×	×	Q	×
5	×	×	×	×
6	×	×	×	×
7	×	×	×	×
8	×	×	×	×

从图中可以得到以下启示：一个合适的解应是在每列、每行上只有一个皇后，且一条斜线上也只有一个皇后。

求解过程从空配置开始。在第1列至第m列为合理配置的基础上，再配置第m+1列，直至第n列配置也是合理时，就找到了一个解。接着改变第n列配置，希望获得下一个解。另外，在任一列上，可能有n种配置。开始时配置在第1行，以后改变时，顺次选择第2行、第3行、…、直到第n行。当第n行配置也找不到一个合理的配置时，就要回溯，去改变前一列的配置。得到求解皇后问题的算法如下：

1. { 输入棋盘大小值n；
   
2. m=0;
   
3. good=1;
   
4. do {
   
5. if (good)
   
6. if (m==n)
   
7. { 输出解；
   
8. 改变之，形成下一个候选解;
   
9. }
   
10. else 扩展当前候选接至下一列；
    
11. else 改变之，形成下一个候选解；
    
12. good=检查当前候选解的合理性；
    
13. } while (m!=0);
    
14. }

复制代码

simon21

在编写程序之前，先确定边式棋盘的数据结构。比较直观的方法是采用一个二维数组，但仔细观察就会发现，这种表示方法给调整候选解及检查其合理性带来困难。更好的方法乃是尽可能直接表示那些常用的信息。对于本题来说，“常用信息”并不是皇后的具体位置，而是“一个皇后是否已经在某行和某条斜线合理地安置好了”。因在某一列上恰好放一个皇后，引入一个一维数组（col[ ]），值col表示在棋盘第i列、col行有一个皇后。例如：col[3]=4，就表示在棋盘的第3列、第4行上有一个皇后。另外，为了使程序在找完了全部解后回溯到最初位置，设定col[0]的初值为0当回溯到第0列时，说明程序已求得全部解，结束程序运行。
为使程序在检查皇后配置的合理性方面简易方便，引入以下三个工作数组：

（1） 数组a[ ]，a[k]表示第k行上还没有皇后；

（2） 数组b[ ]，b[k]表示第k列右高左低斜线上没有皇后；

（3） 数组 c[ ]，c[k]表示第k列左高右低斜线上没有皇后；

棋盘中同一右高左低斜线上的方格，他们的行号与列号之和相同；同一左高右低斜线上的方格，他们的行号与列号之差均相同。

初始时，所有行和斜线上均没有皇后，从第1列的第1行配置第一个皇后开始，在第m列col[m]行放置了一个合理的皇后后，准备考察第m+1列时，在数组a[ ]、b[ ]和c[ ]中为第m列，col[m]行的位置设定有皇后标志；当从第m列回溯到第m-1列，并准备调整第m-1列的皇后配置时，清除在数组a[ ]、b[ ]和c[ ]中设置的关于第m-1列，col[m-1]行有皇后的标志。一个皇后在m列，col[m]行方格内配置是合理的，由数组a[ ]、b[ ]和c[ ]对应位置的值都为1来确定。细节见以下程序：

【程序】

1. # include <stdio.h>
   
2. # include <stdlib.h>
   
3. # define MAXN 20
   
4. 
   
5. int n,m,good;
   
6. int col[MAXN+1],a[MAXN+1],b[2*MAXN+1],c[2*MAXN+1];
   
7. 
   
8. void main()
   
9. { int j;
   
10. char awn;
    
11. printf(“Enter n: “); scanf(“%d”,&n);
    
12. for (j=0;j<=n;j++) a[j]=1;
    
13. for (j=0;j<=2*n;j++) cb[j]=c[j]=1;
    
14. m=1; col[1]=1; good=1; col[0]=0;
    
15. do {
    
16. if (good)
    
17. if (m==n)
    
18. { printf(“列\t行”);
    
19. for (j=1;j<=n;j++)
    
20. printf(“%3d\t%d\n”,j,col[j]);
    
21. printf(“Enter a character (Q/q for exit)!\n”);
    
22. scanf(“%c”,&awn);
    
23. if (awn==’Q’||awn==’q’) exit(0);
    
24. while (col[m]==n)
    
25. { m--;
    
26. a[col[m]]=b[m+col[m]]=c[n+m-col[m]]=1;
    
27. }
    
28. col[m]++;
    
29. }
    
30. else
    
31. { a[col[m]]=b[m+col[m]]=c[n+m-col[m]]=0;
    
32. col[++m]=1;
    
33. }
    
34. else
    
35. { while (col[m]==n)
    
36. { m--;
    
37. a[col[m]]=b[m+col[m]]=c[n+m-col[m]]=1;
    
38. }
    
39. col[m]++;
    
40. }
    
41. good=a[col[m]]&&b[m+col[m]]&&c[n+m-col[m]];
    
42. } while (m!=0);
    
43. }

复制代码

试探法找解算法也常常被编写成递归函数，下面两程序中的函数queen_all()和函数queen_one()能分别用来解皇后问题的全部解和一个解。

【程序】

1. # include <stdio.h>
   
2. # include <stdlib.h>
   
3. # define MAXN 20
   
4. 
   
5. int n;
   
6. int col[MAXN+1],a[MAXN+1],b[2*MAXN+1],c[2*MAXN+1];
   
7. 
   
8. void main()
   
9. { int j;
   
10. printf(“Enter n: “); scanf(“%d”,&n);
    
11. for (j=0;j<=n;j++) a[j]=1;
    
12. for (j=0;j<=2*n;j++) cb[j]=c[j]=1;
    
13. queen_all(1,n);
    
14. }
    
15. 
    
16. void queen_all(int k,int n)
    
17. { int i,j;
    
18. char awn;
    
19. for (i=1;i<=n;i++)
    
20. if (a&&b[k+i]&&c[n+k-i])
    
21. { col[k]=i;
    
22. a=b[k+i]=c[n+k-i]=0;
    
23. if (k==n)
    
24. { printf(“列\t行”);
    
25. for (j=1;j<=n;j++)
    
26. printf(“%3d\t%d\n”,j,col[j]);
    
27. printf(“Enter a character (Q/q for exit)!\n”);
    
28. scanf(“%c”,&awn);
    
29. if (awn==’Q’||awn==’q’) exit(0);
    
30. }
    
31. queen_all(k+1,n);
    
32. a=b[k+i]=c[n+k-i];
    
33. }
    
34. }

复制代码

采用递归方法找一个解与找全部解稍有不同，在找一个解的算法中，递归算法要对当前候选解最终是否能成为解要有回答。当它成为最终解时，递归函数就不再递归试探，立即返回；若不能成为解，就得继续试探。设函数queen_one()返回1表示找到解，返回0表示当前候选解不能成为解。细节见以下函数。

【程序】

1. # define MAXN 20
   
2. int n;
   
3. int col[MAXN+1],a[MAXN+1],b[2*MAXN+1],c[2*MAXN+1];
   
4. int queen_one(int k,int n)
   
5. { int i,found;
   
6. i=found=0;
   
7. While (!found&&i<n)
   
8. { i++;
   
9. if (a&&b[k+i]&&c[n+k-i])
   
10. { col[k]=i;
    
11. a=b[k+i]=c[n+k-i]=0;
    
12. if (k==n) return 1;
    
13. else
    
14. found=queen_one(k+1,n);
    
15. a=b[k+i]=c[n+k-i]=1;
    
16. }
    
17. }
    
18. return found;
    
19. }

复制代码

simon21

六、贪婪法

贪婪法是一种不追求最优解，只希望得到较为满意解的方法。贪婪法一般可以快速得到满意的解，因为它省去了为找最优解要穷尽所有可能而必须耗费的大量时间。贪婪法常以当前情况为基础作最优选择，而不考虑各种可能的整体情况，所以贪婪法不要回溯。

例如平时购物找钱时，为使找回的零钱的硬币数最少，不考虑找零钱的所有各种发表方案，而是从最大面值的币种开始，按递减的顺序考虑各币种，先尽量用大面值的币种，当不足大面值币种的金额时才去考虑下一种较小面值的币种。这就是在使用贪婪法。这种方法在这里总是最优，是因为银行对其发行的硬币种类和硬币面值的巧妙安排。如只有面值分别为1、5和11单位的硬币，而希望找回总额为15单位的硬币。按贪婪算法，应找1个11单位面值的硬币和4个1单位面值的硬币，共找回5个硬币。但最优的解应是3个5单位面值的硬币。

【问题】 装箱问题

问题描述：装箱问题可简述如下：设有编号为0、1、…、n-1的n种物品，体积分别为v0、v1、…、vn-1。将这n种物品装到容量都为V的若干箱子里。约定这n种物品的体积均不超过V，即对于0≤i＜n，有0＜vi≤V。不同的装箱方案所需要的箱子数目可能不同。装箱问题要求使装尽这n种物品的箱子数要少。

若考察将n种物品的集合分划成n个或小于n个物品的所有子集，最优解就可以找到。但所有可能划分的总数太大。对适当大的n，找出所有可能的划分要花费的时间是无法承受的。为此，对装箱问题采用非常简单的近似算法，即贪婪法。该算法依次将物品放到它第一个能放进去的箱子中，该算法虽不能保证找到最优解，但还是能找到非常好的解。不失一般性，设n件物品的体积是按从大到小排好序的，即有v0≥v1≥…≥vn-1。如不满足上述要求，只要先对这n件物品按它们的体积从大到小排序，然后按排序结果对物品重新编号即可。装箱算法简单描述如下：

{ 输入箱子的容积；

输入物品种数n；

按体积从大到小顺序，输入各物品的体积；

预置已用箱子链为空；

预置已用箱子计数器box_count为0；

for (i=0;i<n;i++)

{ 从已用的第一只箱子开始顺序寻找能放入物品i 的箱子j；

if （已用箱子都不能再放物品i）

{ 另用一个箱子，并将物品i放入该箱子；

box_count++；

}

else

将物品i放入箱子j；

}

}

上述算法能求出需要的箱子数box_count，并能求出各箱子所装物品。下面的例子说明该算法不一定能找到最优解，设有6种物品，它们的体积分别为：60、45、35、20、20和20单位体积，箱子的容积为100个单位体积。按上述算法计算，需三只箱子，各箱子所装物品分别为：第一只箱子装物品1、3；第二只箱子装物品2、4、5；第三只箱子装物品6。而最优解为两只箱子，分别装物品1、4、5和2、3、6。

若每只箱子所装物品用链表来表示，链表首结点指针存于一个结构中，结构记录尚剩余的空间量和该箱子所装物品链表的首指针。另将全部箱子的信息也构成链表。以下是按以上算法编写的程序。

【程序】

# include <stdio.h>

# include <stdlib.h>

typedef struct ele

{ int vno;

struct ele *link;

} ELE;

typedef struct hnode

{ int remainder;

ELE *head;

Struct hnode *next;

} HNODE;


void main()

{ int n, i, box_count, box_volume, *a;

HNODE *box_h, *box_t, *j;

ELE *p, *q;

Printf(“输入箱子容积\n”);

Scanf(“%d”,&box_volume);

Printf(“输入物品种数\n”);

Scanf(“%d”,&n);

A=(int *)malloc(sizeof(int)*n);

Printf(“请按体积从大到小顺序输入各物品的体积：”);

For (i=0;i<n;i++) scanf(“%d”,a+i);

Box_h=box_t=NULL;

Box_count=0;

For (i=0;i<n;i++)

{ p=(ELE *)malloc(sizeof(ELE));

p->vno=i;

for (j=box_h;j!=NULL;j=j->next)

if (j->remainder>=a) break;

if (j==NULL)

{ j=(HNODE *)malloc(sizeof(HNODE));

j->remainder=box_volume-a;

j->head=NULL;

if (box_h==NULL) box_h=box_t=j;

else box_t=boix_t->next=j;

j->next=NULL;

box_count++;

}

else j->remainder-=a;

for (q=j->next;q!=NULL&&q->link!=NULL;q=q->link);

if (q==NULL)

{ p->link=j->head;

j->head=p;

}

else

{ p->link=NULL;

q->link=p;

}

}

printf(“共使用了%d只箱子”，box_count);

printf(“各箱子装物品情况如下：”);

for (j=box_h,i=1;j!=NULL;j=j->next,i++)

{ printf(“第%2d只箱子，还剩余容积%4d，所装物品有；\n”,I,j->remainder);

for (p=j->head;p!=NULL;p=p->link)

printf(“%4d”,p->vno+1);

printf(“\n”);

}

}

simon21

【问题】 马的遍历

问题描述：在8×8方格的棋盘上，从任意指定的方格出发，为马寻找一条走遍棋盘每一格并且只经过一次的一条路径。

马在某个方格，可以在一步内到达的不同位置最多有8个，如图所示。如用二维数组board[ ][ ]表示棋盘，其元素记录马经过该位置时的步骤号。另对马的8种可能走法（称为着法）设定一个顺序，如当前位置在棋盘的（i，j）方格，下一个可能的位置依次为（i+2，j+1）、（i+1，j+2）、（i-1，j+2）、（i-2，j+1）、（i-2，j-1）、（i-1，j-2）、（i+1，j-2）、（i+2，j-1），实际可以走的位置尽限于还未走过的和不越出边界的那些位置。为便于程序的同意处理，可以引入两个数组，分别存储各种可能走法对当前位置的纵横增量。

对于本题，一般可以采用回溯法，这里采用Warnsdoff策略求解，这也是一种贪婪法，其选择下一出口的贪婪标准是在那些允许走的位置中，选择出口最少的那个位置。如马的当前位置（i，j）只有三个出口，他们是位置（i+2，j+1）、（i-2，j+1）和（i-1，j-2），如分别走到这些位置，这三个位置又分别会有不同的出口，假定这三个位置的出口个数分别为4、2、3，则程序就选择让马走向（i-2，j+1）位置。

由于程序采用的是一种贪婪法，整个找解过程是一直向前，没有回溯，所以能非常快地找到解。但是，对于某些开始位置，实际上有解，而该算法不能找到解。对于找不到解的情况，程序只要改变8种可能出口的选择顺序，就能找到解。改变出口选择顺序，就是改变有相同出口时的选择标准。以下程序考虑到这种情况，引入变量start，用于控制8种可能着法的选择顺序。开始时为0，当不能找到解时，就让start增1，重新找解。细节以下程序。

【程序】

# include <stdio.h>

int delta_i[ ]={2,1,-1,-2,-2,-1,1,2};

int delta_j[ ]={1,2,2,1,-1,-2,-2,-1};

int board[8][8];

int exitn(int i,int j,int s,int a[ ])

{ int i1,j1,k,count;

for (count=k=0;k<8;k++)

{ i1=i+delta_i[(s+k)%8];

j1=i+delta_j[(s+k)%8];

if (i1>=0&&i1<8&&j1>=0&&j1<8&&board[I1][j1]==0)

a[count++]=(s+k)%8;

}

return count;

}


int next(int i,int j,int s)

{ int m,k,mm,min,a[8],b[8],temp;

m=exitn(i,j,s,a);

if (m==0) return –1;

for (min=9,k=0;k<m;k++)

{ temp=exitn(I+delta_i[a[k]],j+delta_j[a[k]],s,b);

if (temp<min)

{ min=temp;

kk=a[k];

}

}

return kk;

}


void main()

{ int sx,sy,i,j,step,no,start;

for (sx=0;sx<8;sx++)

for (sy=0;sy<8;sy++)

{ start=0;

do {

for (i=0;i<8;i++)

for (j=0;j<8;j++)

board[j]=0;

board[sx][sy]=1;

I=sx; j=sy;

For (step=2;step<64;step++)

{ if ((no=next(i,j,start))==-1) break;

I+=delta_i[no];

j+=delta_j[no];

board[j]=step;

}

if (step>64) break;

start++;

} while(step<=64)

for (i=0;i<8;i++)

{ for (j=0;j<8;j++)

printf(“%4d”,board[j]);

printf(“\n\n”);

}

scanf(“%*c”);

}

}

simon21

1、分治法的基本思想

任何一个可以用计算机求解的问题所需的计算时间都与其规模N有关。问题的规模越小，越容易直接求解，解题所需的计算时间也越少。例如，对于n个元素的排序问题，当n=1时，不需任何计算；n=2时，只要作一次比较即可排好序；n=3时只要作3次比较即可，…。而当n较大时，问题就不那么容易处理了。要想直接解决一个规模较大的问题，有时是相当困难的。

分治法的设计思想是，将一个难以直接解决的大问题，分割成一些规模较小的相同问题，以便各个击破，分而治之。

如果原问题可分割成k个子问题（1<k≤n），且这些子问题都可解，并可利用这些子问题的解求出原问题的解，那么这种分治法就是可行的。由分治法产生的子问题往往是原问题的较小模式，这就为使用递归技术提供了方便。在这种情况下，反复应用分治手段，可以使子问题与原问题类型一致而其规模却不断缩小，最终使子问题缩小到很容易直接求出其解。这自然导致递归过程的产生。分治与递归像一对孪生兄弟，经常同时应用在算法设计之中，并由此产生许多高效算法。

2、分治法的适用条件

分治法所能解决的问题一般具有以下几个特征：

（1）该问题的规模缩小到一定的程度就可以容易地解决；

（2）该问题可以分解为若干个规模较小的相同问题，即该问题具有最优子结构性质；

（3）利用该问题分解出的子问题的解可以合并为该问题的解；

（4）该问题所分解出的各个子问题是相互独立的，即子问题之间不包含公共的子子问题。

上述的第一条特征是绝大多数问题都可以满足的，因为问题的计算复杂性一般是随着问题规模的增加而增加；第二条特征是应用分治法的前提，它也是大多数问题可以满足的，此特征反映了递归思想的应用；第三条特征是关键，能否利用分治法完全取决于问题是否具有第三条特征，如果具备了第一条和第二条特征，而不具备第三条特征，则可以考虑贪心法或动态规划法。第四条特征涉及到分治法的效率，如果各子问题是不独立的，则分治法要做许多不必要的工作，重复地解公共的子问题，此时虽然可用分治法，但一般用动态规划法较好。

3、分治法的基本步骤

分治法在每一层递归上都有三个步骤：

（1）分解：将原问题分解为若干个规模较小，相互独立，与原问题形式相同的子问题；

（2）解决：若子问题规模较小而容易被解决则直接解，否则递归地解各个子问题；

（3）合并：将各个子问题的解合并为原问题的解。

它的一般的算法设计模式如下：

Divide_and_Conquer（P）

if |P|≤n0

then return（ADHOC（P））

将P分解为较小的子问题P1、P2、…、Pk

for i←1 to k

do

yi ← Divide-and-Conquer（Pi） △ 递归解决Pi

T ← MERGE（y1，y2，…，yk） △ 合并子问题

Return（T）

其中 |P| 表示问题P的规模；n0为一阈值，表示当问题P的规模不超过n0时，问题已容易直接解出，不必再继续分解。ADHOC（P）是该分治法中的基本子算法，用于直接解小规模的问题P。因此，当P的规模不超过n0时，直接用算法ADHOC（P）求解。

算法MERGE（y1，y2，…，yk）是该分治法中的合并子算法，用于将P的子问题P1、P2、…、Pk的相应的解y1、y2、…、yk合并为P的解。

根据分治法的分割原则，原问题应该分为多少个子问题才较适宜？各个子问题的规模应该怎样才为适当？这些问题很难予以肯定的回答。但人们从大量实践中发现，在用分治法设计算法时，最好使子问题的规模大致相同。换句话说，将一个问题分成大小相等的k个子问题的处理方法是行之有效的。许多问题可以取k=2。这种使子问题规模大致相等的做法是出自一种平衡子问题的思想，它几乎总是比子问题规模不等的做法要好。

分治法的合并步骤是算法的关键所在。有些问题的合并方法比较明显，有些问题合并方法比较复杂，或者是有多种合并方案；或者是合并方案不明显。究竟应该怎样合并，没有统一的模式，需要具体问题具体分析。

simon21

【问题】 大整数乘法

问题描述：

通常，在分析一个算法的计算复杂性时，都将加法和乘法运算当作是基本运算来处理，即将执行一次加法或乘法运算所需的计算时间当作一个仅取决于计算机硬件处理速度的常数。

这个假定仅在计算机硬件能对参加运算的整数直接表示和处理时才是合理的。然而，在某些情况下，我们要处理很大的整数，它无法在计算机硬件能直接表示的范围内进行处理。若用浮点数来表示它，则只能近似地表示它的大小，计算结果中的有效数字也受到限制。若要精确地表示大整数并在计算结果中要求精确地得到所有位数上的数字，就必须用软件的方法来实现大整数的算术运算。

请设计一个有效的算法，可以进行两个n位大整数的乘法运算。

设X和Y都是n位的二进制整数，现在要计算它们的乘积XY。我们可以用小学所学的方法来设计一个计算乘积XY的算法，但是这样做计算步骤太多，显得效率较低。如果将每2个1位数的乘法或加法看作一步运算，那么这种方法要作O(n2)步运算才能求出乘积XY。下面我们用分治法来设计一个更有效的大整数乘积算法。


图6-3 大整数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)。利用式（3），并考虑到X和Y的符号对结果的影响，我们给出大整数相乘的完整算法MULT如下：

function MULT(X，Y，n); {X和Y为2个小于2n的整数，返回结果为X和Y的乘积XY}

begin

S=SIGN(X)*SIGN(Y); {S为X和Y的符号乘积}

X=ABS(X);

Y=ABS(Y); {X和Y分别取绝对值}

if n=1 then

if (X=1)and(Y=1) then return(S)

else return(0)

else begin

A=X的左边n/2位;

B=X的右边n/2位;

C=Y的左边n/2位;

D=Y的右边n/2位;

ml=MULT(A,C,n/2);

m2=MULT(A-B,D-C,n/2);

m3=MULT(B,D,n/2);

S=S*(m1*2n+(m1+m2+m3)*2n/2+m3);

return(S);

end;

end;

上述二进制大整数乘法同样可应用于十进制大整数的乘法以提高乘法的效率减少乘法次数。

simon21

【问题】 最接近点对问题

问题描述：

在应用中，常用诸如点、圆等简单的几何对象代表现实世界中的实体。在涉及这些几何对象的问题中，常需要了解其邻域中其他几何对象的信息。例如，在空中交通控制问题中，若将飞机作为空间中移动的一个点来看待，则具有最大碰撞危险的2架飞机，就是这个空间中最接近的一对点。这类问题是计算几何学中研究的基本问题之一。下面我们着重考虑平面上的最接近点对问题。

最接近点对问题的提法是：给定平面上n个点，找其中的一对点，使得在n个点的所有点对中，该点对的距离最小。

严格地说，最接近点对可能多于1对。为了简单起见，这里只限于找其中的一对。

这个问题很容易理解，似乎也不难解决。我们只要将每一点与其他n-1个点的距离算出，找出达到最小距离的两个点即可。然而，这样做效率太低，需要O(n2)的计算时间。我们能否找到问题的一个O (nlogn)算法。

这个问题显然满足分治法的第一个和第二个适用条件，我们考虑将所给的平面上n个点的集合S分成2个子集S1和S2，每个子集中约有n/2个点，然后在每个子集中递归地求其最接近的点对。在这里，一个关键的问题是如何实现分治法中的合并步骤，即由S1和S2的最接近点对，如何求得原集合S中的最接近点对，因为S1和S2的最接近点对未必就是S的最接近点对。如果组成S的最接近点对的2个点都在S1中或都在S2中，则问题很容易解决。但是，如果这2个点分别在S1和S2中，则对于S1中任一点p，S2中最多只有n/2个点与它构成最接近点对的候选者，仍需做n2/4次计算和比较才能确定S的最接近点对。因此，依此思路，合并步骤耗时为O(n2)。整个算法所需计算时间T(n)应满足：

T(n)=2T(n/2)+O(n2)


它的解为T(n)=O(n2)，即与合并步骤的耗时同阶，显示不出比用穷举的方法好。从解递归方程的套用公式法，我们看到问题出在合并步骤耗时太多。这启发我们把注意力放在合并步骤上。

为了使问题易于理解和分析，我们先来考虑一维的情形。此时S中的n个点退化为x轴上的n个实数x1、x2、…、xn。最接近点对即为这n个实数中相差最小的2个实数。我们显然可以先将x1、x2、…、xn排好序，然后，用一次线性扫描就可以找出最接近点对。这种方法主要计算时间花在排序上，因此如在排序算法中所证明的，耗时为O(nlogn)。然而这种方法无法直接推广到二维的情形。因此，对这种一维的简单情形，我们还是尝试用分治法来求解，并希望能推广到二维的情形。

假设我们用x轴上某个点m将S划分为2个子集S1和S2，使得S1={x∈S | x≤m}；S2={x∈S | x>m}。这样一来，对于所有p∈S1和q∈S2有p<q。

递归地在S1和S2上找出其最接近点对{p1，p2}和{q1，q2}，并设δ=min{|p1-p2|，|q1-q2|}，S中的最接近点对或者是{p1，p2}，或者是{q1，q2}，或者是某个{p3，q3}，其中p3∈S1且q3∈S2。如图1所示。


图1 一维情形的分治法

我们注意到，如果S的最接近点对是{p3，q3}，即 | p3-q3 | < δ，则p3和q3两者与m的距离不超过δ，即 | p3-m | < δ，| q3-m | < δ，也就是说，p3∈(m-δ，m)，q3∈(m，m+δ)。由于在S1中，每个长度为δ的半闭区间至多包含一个点（否则必有两点距离小于δ），并且m是S1和S2的分割点，因此(m-δ，m)中至多包含S中的一个点。同理，(m，m+δ)中也至多包含S中的一个点。由图1可以看出，如果(m-δ，m)中有S中的点，则此点就是S1中最大点。同理，如果(m，m+δ)中有S中的点，则此点就是S2中最小点。因此，我们用线性时间就能找到区间(m-δ，m)和(m，m+δ)中所有点，即p3和q3。从而我们用线性时间就可以将S1的解和S2的解合并成为S的解。也就是说，按这种分治策略，合并步可在O(n)时间内完成。这样是否就可以得到一个有效的算法了呢？
还有一个问题需要认真考虑，即分割点m的选取，及S1和S2的划分。选取分割点m的一个基本要求是由此导出集合S的一个线性分割，即S=S1∪S2 ，S1∩S2=Φ，且S1 {x | x≤m}；S2 {x | x>m}。容易看出，如果选取m=[max（S）+min（S）]/2，可以满足线性分割的要求。选取分割点后，再用O(n)时间即可将S划分成S1={x∈S | x≤m}和S2={x∈S | x>m}。然而，这样选取分割点m，有可能造成划分出的子集S1和S2的不平衡。例如在最坏情况下，|S1|=1，|S2|=n-1，由此产生的分治法在最坏情况下所需的计算时间T（n）应满足递归方程：

T（n）=T（n-1）+O（n）

它的解是T（n）=O（n2）。这种效率降低的现象可以通过分治法中“平衡子问题”的方法加以解决。也就是说，我们可以通过适当选择分割点m，使S1和S2中有大致相等个数的点。自然地，我们会想到用S的n个点的坐标的中位数来作分割点。在选择算法中介绍的选取中位数的线性时间算法使我们可以在O（n）时间内确定一个平衡的分割点m。

至此，我们可以设计出一个求一维点集S中最接近点对的距离的算法pair如下。

Float pair（S）;

{ if | S | =2 δ= | x[2]－x[1] | /*x[1..n]存放的是S中n个点的坐标*/

else

{ if ( | S | =1) δ=∞

else

{ m=S中各点的坐标值的中位数;

构造S1和S2，使S1={x∈S | x≤m}，S2={x∈S | x＞m};

δ1=pair(S1);

δ2=pair(S2);

p=max(S1);

q=min(S2);

δ=min(δ1，δ2，q-p);

}

return(δ);

}