第5讲回溯法专题讲座编程大赛课程Word下载.docx
《第5讲回溯法专题讲座编程大赛课程Word下载.docx》由会员分享,可在线阅读,更多相关《第5讲回溯法专题讲座编程大赛课程Word下载.docx(44页珍藏版)》请在冰豆网上搜索。
图5-14皇后问题回溯实施求解
图(a)为在第1行第1列放置一个皇后的初始状态。
图(b)中,第2个皇后不能放在第1、2列,因而放置在第3列上。
图(c)中,表示第3行的所有各列均不能放置皇后,则返回第2行,第2个皇后需后移。
图(d)中,第2个皇后后移到第4列,第3个皇后放置在第2列。
图(e)中,第4行的所有各列均不能放置皇后,则返回第3行;
第3个皇后后移的所有位置均不能放置皇后,则返回第2行;
第2个皇后已无位可退,则返回第1行;
第1个皇后需后移。
图(f)中,第1个皇后后移至第2格。
图(g)中,第2个皇后不能放在第1,2,3列,因而放置在第4列上。
图(h)中,第3个皇后放在第1列;
第4个皇后不能放置1,2列,于是放置在第3列。
这样经以上回溯,得到4皇后问题的一个解:
2413。
继续以上的回溯探索,可得4皇后问题的另一个解:
3142。
3.回溯算法框架描述
(1)回溯描述
对于一般含参量m,n的搜索问题,回溯法框架描述如下:
输入正整数n,m,(n≥m)
i=1;
a[i]=<
元素初值>
;
while
(1)
{
for(g=1,k=i-1;
k>
=1;
k--)
if(<
约束条件1>
)g=0;
//检测约束条件,不满足则返回
if(g&
&
<
约束条件2>
)
printf(a[1-m]);
//输出一个解
if(i<
n&
g){i++;
取值点>
continue;
}
while(a[i]==<
回溯点>
&
i>
1)i--;
//向前回溯
if(a[i]==n&
i==1)break;
//退出循环,结束
elsea[i]=a[i]+1;
}
具体求解问题的试探搜索范围与要求不同,在应用回溯设计时,需根据问题的具体实际确定数组元素的初值、取值点与回溯点,同时需把问题中的约束条件进行必要的分解,以适应上述回溯流程。
其中实施向前回溯的循环
while(a[i]==<
是向前回溯一步,还是回溯两步或更多步,完全根据a[i]是否达到回溯点来确定。
例如,回溯点是n,i=6,当a[6]=n时回溯到i=5;
若a[5]=n时回溯到i=4;
依此类推,到a[i]达到回溯点则停止。
图5-1所示的4皇后问题迭代回溯过程描述为:
n=4;
i=1;
a[i]=1;
while
(1)
{
g=1;
for(k=i-1;
if(a[i]=a[k]&
abs(a[i]-a[k])=i-k)
g=0;
//检测约束条件,不满足则返回
i==4)
printf(a[1-4]);
if(i<
while(a[i]==n&
//向前回溯
if(a[i]==n&
//退出循环结束
}
以上回溯体现在迭代式i=i-1,因而又称为迭代回溯。
此外,递归也能实现回溯。
(2)递归回溯
intput(intk)
{inti,j,u;
if(k<
=问题规模)
{u=0;
if(<
约束条件>
)
u=1;
//当k时不可操作
if(u==0)//当k时可操作
{if(k==规模)//若已满足规模,则打印出一个解
printf(<
一个解>
);
elseput(k+1);
//调用put(k+1)
}
在调用put(k)时,当检测约束条件知不可操作(记u=1),即再往前不可能得解,此时当然不可能输出解,也不调用put(k+1),而是回溯,返回调用put(k)之处。
这就是递归回溯的机理。
如果是主程序调用put
(1),最后返回到主程序调用put
(1)的后续语句,完成递归。
图5.1所示的4皇后问题迭代回溯过程描述为:
if(k<
=4)
{for(i=1;
i<
=4;
i++)//探索第k行从第1格开始放皇后
{a[k]=i;
for(u=0,j=1;
j<
=k-1;
j++)
if(a[k]==a[j]||abs(a[k]-a[j])==k-j)
//若第k行第i格放不下,则置u=1
if(u==0)//若第k行第i格可放,则检测是否满4行
{if(k==4)//若已放满到4行时,则打印出一个解
{s++;
printf("
"
);
for(j=1;
printf("
%d"
a[j]);
}
//若没放满4行,则放下一行put(k+1)
4.回溯法的效益分析
应用回溯设计求解实际问题,由于解空间的结构差异,很难精确计算与估计回溯产生的结点数,因此回溯法的复杂度是分析回溯法效率时遇到的主要困难。
回溯法产生的结点数通常只有解空间结点数的一小部分,这也是回溯法的计算效率大大高于穷举法的原因所在。
回溯求解过程实质上是一个遍历一棵“状态树”的过程,只是这棵树不是遍历前预先建立的。
回溯算法在搜索过程中,只要所激活的状态结点满足终结条件,应该把它输出或保存。
由于在回溯法求解问题时,一般要求输出问题的所有解,因此在得到结点后,同时也要进行回溯,以便得到问题的其他解,直至回溯到状态树的根且根的所有子结点均已被搜索过为止。
组织解空间便于算法在求解集时更易于搜索,典型的组织方法是图或树。
一旦定义了解空间的组织方法,这个空间即可从开始结点进行搜索。
回溯法的时间通常取决于状态空间树上实际生成的那部分问题状态的数目。
对于元组长度为n的问题,若其状态空间树中结点总数为n!
,则回溯算法的最坏情形的时间复杂度可达O(p(n)n!
);
若其状态空间树中结点总数为2n,则回溯算法的最坏情形的时间复杂度可达O(p(n)2n),其中p(n)为n的多项式。
对于不同的实例,回溯法的计算时间有很大的差异。
对于很多具有大n的求解实例,应用回溯法一般可在很短的时间内求得其解,可见回溯法不失为一种快速有效的算法。
对于某一具体实际问题的回溯求解,常通过计算实际生成结点数的方法即蒙特卡罗方法(Montecarlo)来评估其计算效率。
蒙特卡罗方法的基本思想是在状态空间树上随机选择一条路径(x0,x1,…,xn-1),设X是这一路径上部分向量(x0,x1,…,xk-1)的结点,如果在X处不受限制的子向量数是mk,则认为与X同一层的其他结点不受限制的子向量数也都是mk。
也就是说,若不受限制的x0取值有m0个,则该层上有m0个结点;
若不受限制的x1取值有m1个,则该层上有m0m1个结点;
依此类推。
由于认为在同一层上不受限制的结点数相同,因此,该路径上实际生成的结点数估计为
计算路径上结点数m的蒙特卡罗算法描述如下:
//已知随机路径上取值数据m0,m1,…,mk-1
m=1;
t=1;
for(j=0;
{t=t*m[j];
m=m+t;
printf(“%ld”,m);
把所求得的随机路径上的结点数(或若干条随机路径的结点数的平均值)与状态空间树上的总结点数进行比较,由其比值可以初步看出回溯设计的效益。
在下面的n皇后问题的回溯求解时将具体应用以上蒙特卡罗算法估计回溯设计的效益。
5.2桥本分数式
5.2.1桥本分数式
1.案例提出
日本数学家桥本吉彦教授于1993年10月在我国山东举行的中日美三国数学教育研讨会上向与会者提出以下填数趣题:
把1,2,…,9这9个数字填入下式的9个方格中(数字不得重复),使下面的分数等式成立
.□□□
──+──=──
□□□□□□
桥本教授当即给出了一个解答。
这一分数式填数趣题究竟共有多少个解答?
试求出所有解答。
(等式左边两个分数交换次序只算一个解答)。
这一填数趣题的解是否唯一?
如果不唯一究竟有多少个解?
由人工推算求解难度太大,通过程序设计由计算机来探求更为合适。
2.回溯设计
(1)设计要点
我们采用回溯法逐步调整探求。
把式中9个□规定一个顺序后,先在第一个□中填入一个数字(从1开始递增),然后从小到大选择一个不同于前面□的数字填在第二个□中,依此类推,把九个□都填入没有重复的数字后,检验是否满足等式。
若等式成立,打印所得的解。
然后第九个□中的数字调整增1再试,直到调整为9(不能再增);
返回前一个□中数字调整增1再试;
依此类推,直至第一个□中的数字调整为9时,完成调整探求。
可见,问题的解空间是9位的整数组,其约束条件是9位数中没有相同数字且必须满足分式的要求。
为此,设置a数组,式中每一□位置用一个数组元素来表示:
同时,记式中的3个分母分别为
m1=a
(2)a(3)=a
(2)*10+a(3)
m2=a(5)a(6)=a(5)*10+a(6)
m3=a(8)a(9)=a(8)*10+a(9)
所求分数等式等价于整数等式a
(1)*m2*m3+a(4)*m1*m3=a(7)*m1*m2成立。
这一转化可以把分数的测试转化为整数测试。
注意到等式左侧两分数交换次序只算一个解,为避免解的重复,设a
(1)<
a(4)。
式中9个□各填一个数字,不允许重复。
为判断数字是否重复,设置中间变量g:
先赋值g=1;
若出现某两数字相同(即a(i)=a(k))或a
(1)>a(4),则赋值g=0(重复标记)。
首先从a
(1)=1开始,逐步给a(i)(1≤i≤9)赋值,每一个a(i)赋值从1开始递增至9。
直至a(9)赋值,判断:
若i=9,g=1,a
(1)*m2*m3+a(4)*m1*m3=a(7)*m1*m2同时满足,则为一组解,用n统计解的个数后,格式打印输出这组解。
若i<
9且g=1,表明还不到9个数字,则下一个a(i)从1开始赋值继续。
若a(9)=9,则返回前一个数组元素a(8)增1赋值(此时,a(9)又从1开始)再试。
若a(8)=9,则返回前一个数组元素a(7)增1赋值再试。
依此类推,直到a
(1)=9时,已无法返回,意味着已全部试毕,求解结束。
按以上所描述的回溯的参量:
m=n=9
元素初值:
a[1]=1,数组元素初值取1。
取值点:
a[i]=1,各元素从1开始取值。
回溯点:
a[i]=9,各元素取值至9后回溯。
约束条件1:
a[i]==a[k]||a[1]>a[4],其中(i>k)。
约束条件2:
i=9&
a[1]*m2*m3+a[4]*m1*m3=a[7]*m1*m2
(2)桥本分数式回溯程序设计
//桥本分数式回溯实现
//把1,2,...,9填入□/□□+□/□□=□/□□
#include<
stdio.h>
voidmain()
{intg,i,k,s,a[10];
longm1,m2,m3;
a[1]=1;
s=0;
{g=1;
for(k=i-1;
if(a[i]==a[k]){g=0;
break;
}//两数相同,标记g=0
if(i==9&
g==1&
a[1]<
a[4])
{m1=a[2]*10+a[3];
m2=a[5]*10+a[6];
m3=a[8]*10+a[9];
if(a[1]*m2*m3+a[4]*m1*m3==a[7]*m1*m2)//判断等式
{s++;
printf("
(%2d)"
s);
%d/%ld+%d/"
a[1],m1,a[4]);
%ld=%d/%ld"
m2,a[7],m3);
if(s%2==0)printf("
\n"
9&
g==1)
{i++;
}//不到9个数,往后继续
while(a[i]==9&
//往前回溯
if(a[i]==9&
elsea[i]++;
//至第1个数为9结束
共以上%d个解。
(3)程序运行结果
(1)1/26+5/78=4/39
(2)1/32+5/96=7/84
(3)1/32+7/96=5/48(4)1/78+4/39=6/52
(5)1/96+7/48=5/32(6)2/68+9/34=5/17
(7)2/68+9/51=7/34(8)4/56+7/98=3/21
(9)5/26+9/78=4/13(10)6/34+8/51=9/27
共以上10个解。
3.递归设计
设置桥本分数式递归函数put(k):
当k<
=9时,第k个数字取值a[k]=i(i=1,2,…,9),标记u=0。
a[k]与已取的a[j](j<
k)比较,是否出现重复数字。
若a[k]==a[j],则第k个数字取值不成功,标记u=1;
重新取值。
若保持u=0,第k个数字取值成功:
1)检测k是否到9;
若到9且满足等式,输出一个解。
2)若不到9,或不满足等式要求,则调用put(k+1)。
若a[k]已取到9,返回调用put(k)的k-1状态,即回溯到k-1状态重新取值。
主程序调用put
(1),返回put
(1)时,即输出解的个数s,结束。
(2)递归程序实现
//桥本分数式递归求解
inta[10],s=0;
{intput(intk);
put
(1);
//调用递归函数put
(1)
共有以上%d个解。
//桥本分数式递归函数
{inti,j,u,m1,m2,m3;
=9)
=9;
i++)//探索第k个数字取值i
if(a[k]==a[j])
//出现重复数字,则置u=1
if(u==0)//若第k个数字可为i
{if(k==9&
a[4])//若已9个数字,则检查等式
{m1=a[2]*10+a[3];
if(a[1]*m2*m3+a[4]*m1*m3==a[7]*m1*m2)
%2d:
//输出一个解
printf("
%d/%d+%d/%d"
a[1],m1,a[4],m2);
=%d/%d"
a[7],m3);
if(s%2==0)printf("
//若不到9个数字,则调用put(k+1)
returns;
4.求解说明
以上回溯与递归求解都有回溯功能,所以能快捷地求出所有解。
关于桥本分数式求解,已有应用程序设计得到9个解的报导,遗失了一个解。
可见在程序设计求解时,如果程序中结构欠妥或参量设置不当,都可能造成增解或遗解。
5.2.210数字分数式
1.案例提出
把0,1,2,...,9这10个数字填入下式的10个方格中,要求:
□□□
──+───=──
□□□□□□□
(1)各数字不得重复;
(2)数字“0”不得填在各分数的分子与分母的首位;
(3)式中各分数为最简真分数,即分子分母没有大于1的公因数。
这一分数等式填数趣题究竟共有多少个解答?
试应用回溯求出所有解答。
2.回溯设计
设置a数组表示式中的10个数字,即
m2=a(5)a(6)a(7)=a(5)*100+a(6)*10+a(7)
m3=a(9)a(10)=a(9)*10+a(10)
在上述回溯设计基础上修改若干参数:
数字从9个增加到10个,因而i<
9改为i<
10;
i==9改为i==10;
数组元素取值修改为从“0”开始,即a[1]=0;
a[i]=0;
数字“0”不得在各分数的分子与分母的首位,即“0”只能在a(3),a(6),a(7)与a(10)这4个数字中,因而在输出解的条件中增加a(3)*a(6)*a(7)*a(10)=0。
此外,需增加判断3个分数是否为真分数的测试循环。
3.10数字分数式程序实现
//10数字分数式
{intg,i,k,s,t,u,a[11];
a[1]=0;
}//两数相同,标记g=0
if(i==10&
a[3]*a[6]*a[7]*a[10]==0)
{m1=a[2]*10+a[3];
m2=a[5]*100+a[6]*10+a[7];
m3=a[9]*10+a[10];
if(a[1]*m2*m3+a[4]*m1*m3==a[8]*m1*m2)//判断等式
{t=0;
for(u=2;
u<
u++)//测试3个分数是否为真分数
{if(a[1]%u==0&
m1%u==0){t=1;
if(a[4]%u==0&
m2%u==0){t=1;
if(a[8]%u==0&
m3%u==0){t=1;
}
if(t==0)
{printf("
%d/%ld+%d/"
%ld=%d/%ld\n"
m2,a[8],m3);
10&
a[i]=0;
}//不到10个数,往后继续
//往前回溯
//至第1个数为9结束
4.程序运行结果与说明
4/19+5/608=7/32
以上10数字分数式求解是在9数字分数式设计基础上改动所得,结构完全相同。
请比较以上两个回溯设计的参数变化。
5.3逐位整除数
本节探索一个新颖有趣的案例——逐位整除数,包括高逐位整除数与低逐位整除数。
定义高逐位整除数:
从其高位开始,前1位能被1整除,前2位能被2整除,…,前n位能被n整除。
例如10245就是一个5位高逐位整除数。
定义低逐位整除数:
从其低位(即个位)开始,1位数能被1整除,2位数能被2整除,…,n位数能被n整除。
例如5111120就是一个7位低逐位整除数。
5.3.1高逐位