穷举搜索法是对可能是解的众多候选解按某种顺序进行逐一枚举和检验,并从众找出那些符合要求的候选解作为问题的解。
【问题】
将A、B、C、D、E、F这六个变量排成如图所示的三角形,这六个变量分别取[1,6]上的整数,且均不相同。求使三角形三条边上的变量之和相等的全部解。如图就是一个解。
程序引入变量a、b、c、d、e、f,并让它们分别顺序取1至6的证书,在它们互不相同的条件下,测试由它们排成的如图所示的三角形三条边上的变量之和是否相等,如相等即为一种满足要求的排列,把它们输出。当这些变量取尽所有的组合后,程序就可得到全部可能的解。细节见下面的程序。
【程序1】
# include <stdio.h> |
按穷举法编写的程序通常不能适应变化的情况。如问题改成有9个变量排成三角形,每条边有4个变量的情况,程序的循环重数就要相应改变。
对一组数穷尽所有排列,还有更直接的方法。将一个排列看作一个长整数,则所有排列对应着一组整数。将这组整数按从小到大的顺序排列排成一个整数,从对应最小的整数开始。按数列的递增顺序逐一列举每个排列对应的每个整数,这能更有效地完成排列的穷举。从一个排列找出对应数列的下一个排列可在当前排列的基础上作部分调整来实现。倘若当前排列为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】
# include <stdio.h> # define SIDE_N 3 # define LENGTH 3 # define VARIABLES 6 int A,B,C,D,E,F; int *pt[]={&A,&B,&C,&D,&E,&F}; int *side[SIDE_N][LENGTH]={&A,&B,&C,&C,&D,&E,&E,&F,&A}; int side_total[SIDE_N]; main{} { int i,j,t,equal; for (j=0;j<VARIABLES;j++) *pt[j]=j+1; while(1) { for (i=0;i<SIDE_N;i++) { for (t=j=0;j<LENGTH;j++) t+=*side[j]; side_total=t; } for (equal=1,i=0;equal&&i<SIDE_N-1;i++) if (side_total!=side_total[i+1] equal=0; if (equal) { for (i=1;i<VARIABLES;i++) printf(“%4d”,*pt); printf(“\n”); scanf(“%*c”); } for (j=VARIABLES-1;j>0;j--) if (*pt[j]>*pt[j-1]) break; if (j==0) break; for (i=VARIABLES-1;i>=j;i--) if (*pt>*pt[i-1]) break; t=*pt[j-1];* pt[j-1] =* pt; *pt=t; for (i=VARIABLES-1;i>j;i--,j++) { t=*pt[j]; *pt[j] =* pt; *pt=t; } } } |
从上述问题解决的方法中,最重要的因素就是确定某种方法来确定所有的候选解。下面再用一个示例来加以说明。
【问题】 背包问题
问题描述:有不同价值、不同重量的物品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元组。
【算法】
maxv=0; for (i=0;i<2n;i++) { B[0..n-1]=0; 把i转化为二进制数,存储于数组B中; temp_w=0; temp_v=0; for (j=0;j<n;j++) { if (B[j]==1) { temp_w=temp_w+w[j]; temp_v=temp_v+v[j]; } if ((temp_w<=tw)&&(temp_v>maxv)) { maxv=temp_v; 保存该B数组; } } } |