五、回溯法 e>\[OwF-x
V(n7hpS
回溯法也称为试探法,该方法首先暂时放弃关于问题规模大小的限制,并将问题的候选解按某种顺序逐一枚举和检验。当发现当前候选解不可能是解时,就选择下一个候选解;倘若当前候选解除了还不满足问题规模要求外,满足所有其他要求时,继续扩大当前候选解的规模,并继续试探。如果当前候选解满足包括问题规模在内的所有要求时,该候选解就是问题的一个解。在回溯法中,放弃当前候选解,寻找下一个候选解的过程称为回溯。扩大当前候选解的规模,以继续试探的过程称为向前试探。 Xl^=&!S>me
1、回溯法的一般描述 fOs}5J
可用回溯法求解的问题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的一个解。 gB,~Y511
解问题P的最朴素的方法就是枚举法,即对E中的所有n元组逐一地检测其是否满足D的全部约束,若满足,则为问题P的一个解。但显然,其计算量是相当大的。 1:5jUUL8
我们发现,对于许多问题,所给定的约束集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的解,因而就不必去搜索它们、检测它们。回溯法正是针对这类问题,利用这类问题的上述性质而提出来的比枚举法效率更高的算法。 B*#lkMr
回溯法首先将问题P的n元组的状态空间E表示成一棵高为n的带权有序树T,把在E中求问题P的所有解转化为在T中搜索问题P的所有解。树T类似于检索树,它可以这样构造: t=\y|Idc
设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的根。 d[ (KgX9
因而,在E中寻找问题P的一个解等价于在T中搜索一个叶子结点,要求从T的根到该叶子结点的路径上依次的n条边相应带的n个权x1,x2,…,xn满足约束集D的全部约束。在T中搜索所要求的叶子结点,很自然的一种方式是从根出发,按深度优先的策略逐步深入,即依次搜索满足约束条件的前缀1元组(x1i)、前缀2元组(x1,x2)、…,前缀I元组(x1,x2,…,xi),…,直到i=n为止。 N0h* |
在回溯法中,上述引入的树被称为问题P的状态空间树;树T上任意一个结点被称为问题P的状态结点;树T上的任意一个叶子结点被称为问题P的一个解状态结点;树T上满足约束集D的全部约束的任意一个叶子结点被称为问题P的一个回答状态结点,它对应于问题P的一个解。 2.6%?E]
【问题】 组合问题 dq[X:3i
问题描述:找出从自然数1、2、……、n中任取r个数的所有组合。 }DiMt4!ZC!
例如n=5,r=3的所有组合为: 7xidBVx
(1)1、2、3 (2)1、2、4 (3)1、2、5 ^6tGj+D9
(4)1、3、4 (5)1、3、5 (6)1、4、5 7:'7EqM
(7)2、3、4 (8)2、3、5 (9)2、4、5 nS'hdeoW
(10)3、4、5 [;yEG$)K
则该问题的状态空间为: dHjJLs_
E={(x1,x2,x3)∣xi∈S ,i=1,2,3 } 其中:S={1,2,3,4,5} *CG-F=
约束集为: x1<x2<x3 luAmq+
显然该约束集具有完备性。 x-0S-1M
cC@B\Q
>.^/Z/[.L
2、回溯法的方法 G`3/${ti
对于具有完备约束集D的一般问题P及其相应的状态空间树T,利用T的层次结构和D的完备性,在T中搜索问题P的所有解的回溯法可以形象地描述为: vqDd][ n
从T的根出发,按深度优先的策略,系统地搜索以其为根的子树中可能包含着回答结点的所有状态结点,而跳过对肯定不含回答结点的所有子树的搜索,以提高搜索效率。具体地说,当搜索按深度优先策略到达一个满足D中所有有关约束的状态结点时,即“激活”该状态结点,以便继续往深层搜索;否则跳过对以该状态结点为根的子树的搜索,而一边逐层地向该状态结点的祖先结点回溯,一边“杀死”其儿子结点已被搜索遍的祖先结点,直到遇到其儿子结点未被搜索遍的祖先结点,即转向其未被搜索的一个儿子结点继续搜索。 Q@zD'G>
在搜索过程中,只要所激活的状态结点又满足终结条件,那么它就是回答结点,应该把它输出或保存。由于在回溯法求解问题时,一般要求出问题的所有解,因此在得到回答结点后,同时也要进行回溯,以便得到问题的其他解,直至回溯到T的根且根的所有儿子结点均已被搜索过为止。 2)oT\m
例如在组合问题中,从T的根出发深度优先遍历该树。当遍历到结点(1,2)时,虽然它满足约束条件,但还不是回答结点,则应继续深度遍历;当遍历到叶子结点(1,2,5)时,由于它已是一个回答结点,则保存(或输出)该结点,并回溯到其双亲结点,继续深度遍历;当遍历到结点(1,5)时,由于它已是叶子结点,但不满足约束条件,故也需回溯。 *V}T}nK7
3、回溯法的一般流程和技术 'o41)p
在用回溯法求解有关问题的过程中,一般是一边建树,一边遍历该树。在回溯法中我们一般采用非递归方法。下面,我们给出回溯法的非递归算法的一般流程: G}ZJ}5h
\!D <u'n
在用回溯法求解问题,也即在遍历状态空间树的过程中,如果采用非递归方法,则我们一般要用到栈的数据结构。这时,不仅可以用栈来表示正在遍历的树的结点,而且可以很方便地表示建立孩子结点和回溯过程。 @AM;58.
例如在组合问题中,我们用一个一维数组Stack[ ]表示栈。开始栈空,则表示了树的根结点。如果元素1进栈,则表示建立并遍历(1)结点;这时如果元素2进栈,则表示建立并遍历(1,2)结点;元素3再进栈,则表示建立并遍历(1,2,3)结点。这时可以判断它满足所有约束条件,是问题的一个解,输出(或保存)。这时只要栈顶元素(3)出栈,即表示从结点(1,2,3)回溯到结点(1,2)。 .;#Wf@V
【问题】 组合问题 iupuhq$]
问题描述:找出从自然数1,2,…,n中任取r个数的所有组合。 k$j4~C'$
采用回溯法找问题的解,将找到的组合以从小到大顺序存于a[0],a[1],…,a[r-1]中,组合的元素满足以下性质: h8# 14?
(1) a[i+1]>a,后一个数字比前一个大; ;la sk4|
(2) a-i<=n-r+1。 "`;$wA
按回溯法的思想,找解过程可以叙述如下: i> }P V
首先放弃组合数个数为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]再回溯时,说明已经找完问题的全部解。按上述思想写成程序如下: !>gi9z,
【程序】 -DWyKR= j"
# define MAXN 100 53>y<
int a[MAXN]; x93t.5E6
void comb(int m,int r) z^;0{q,
{ int i,j; y>4r<YZQ
i=0; Vu6$84>-,
a=1; jx*jYil
do { tkQ#mipAj
if (a-i<=m-r+1 jhU'UAn
{ if (i==r-1) -%R3YU3
{ for (j=0;j<r;j++) 9eG{"0)
printf(“%4d”,a[j]); L;jzDng<
printf(“\n”); X9ua&T2(l
} > J4Tk1//b
a++; t}$WP&XRG<
continue; \9r1JP0
} `XT8}9z!
else Y8$Y]2
{ if (i==0) \OR=+\].9
return; ,0j7qn@tm
a[--i]++; 0^l)9zE
} VTM*=5|c
} while (1) xLZJ[:gr
} ,#
i@jB
e u{
main() 1#]0\Y(
{ comb(5,3); uP;qs8
} }bkQr)us
【问题】 填字游戏 \W%Aeg*c
问题描述:在3×3个方格的方阵中要填入数字1到N(N≥10)内的某9个数字,每个方格填一个整数,似的所有相邻两个方格内的两个整数之和为质数。试求出所有满足这个要求的各种数字填法。 ]'
ck!eG
可用试探发找到问题的解,即从第一个方格开始,为当前方格寻找一个合理的整数填入,并在当前位置正确填入后,为下一方格寻找可填入的合理整数。如不能为当前方格找到一个合理的可填证书,就要回退到前一方格,调整前一方格的填入数。当第九个方格也填入合理的整数后,就找到了一个解,将该解输出,并调整第九个的填入的整数,寻找下一个解。 [9}D+k F
为找到一个满足要求的9个数的填法,从还未填一个数开始,按某种顺序(如从小到大的顺序)每次在当前位置填入一个整数,然后检查当前填入的整数是否能满足要求。在满足要求的情况下,继续用同样的方法为下一方格填入整数。如果最近填入的整数不能满足要求,就改变填入的整数。如对当前方格试尽所有可能的整数,都不能满足要求,就得回退到前一方格,并调整前一方格填入的整数。如此重复执行扩展、检查或调整、检查,直到找到一个满足问题要求的解,将解输出。 nw.,`M,N
回溯法找一个解的算法: hmijp1u
{ int m=0,ok=1; 'r/+za:2
int n=8; PIo8m f/
do{ c;RL<83:
if (ok) 扩展; [Un~]E.'J
else 调整; F},#%_4
ok=检查前m个整数填放的合理性; KM'*+.I
} while ((!ok||m!=n)&&(m!=0)) 7IEG%FY
T
if (m!=0) 输出解; 9V ]{q
else 输出无解报告; yX`5x^wVw
} !Eqp,"ts7
如果程序要找全部解,则在将找到的解输出后,应继续调整最后位置上填放的整数,试图去找下一个解。相应的算法如下: utH/E7^8
回溯法找全部解的算法: dq6|m
}g{
{ int m=0,ok=1; pgp@Zw)r)k
int n=8; }3%L3v&
do{ {hP_"nN#
if (ok) .()|0A B&g
{ if (m==n) +mgm39
{ 输出解; Xyz w.%4c
调整; v@qU<\Y>
} ~4T:v_Q7g
else 扩展; d_ [l{
} AlVBhR`
else 调整; g`n5-D@3
ok=检查前m个整数填放的合理性; YA?46[:
} while (m!=0); i(hL6DLD
} neLQ>WT
L
为了确保程序能够终止,调整时必须保证曾被放弃过的填数序列不会再次实验,即要求按某种有许模型生成填数序列。给解的候选者设定一个被检验的顺序,按这个顺序逐一形成候选者并检验。从小到大或从大到小,都是可以采用的方法。如扩展时,先在新位置填入整数1,调整时,找当前候选解中下一个还未被使用过的整数。将上述扩展、调整、检验都编写成程序,细节见以下找全部解的程序。 6g6BE^o\
【程序】 n7>L&?N#y#
# include <stdio.h> %}`zq8Q;
# define N 12 Z)Zc9SVC
void write(int a[ ]) ]cKxYX)J
{ int i,j; >QA;02
for (i=0;i<3;i++) -wdd'G
{ for (j=0;j<3;j++) L-,C5^
printf(“%3d”,a[3*i+j]); :(?F(Q^
printf(“\n”); VhfMj|
} o`{@':%D`
scanf(“%*c”); grs~<n|o\
} IEP^u
`}
z P`&X:8
int b[N+1]; R?Dc*,
int a[10]; S.! n35
int isprime(int m) O8mmS!
{ int i; pWm==Ds|
int primes[ ]={2,3,5,7,11,17,19,23,29,-1}; 5ltrr(MeD
if (m==1||m%2=0) return 0; ==ZL0 ][
for (i=0;primes>0;i++) phc9esz
if (m==primes) return 1; XwlF[3VbiX
for (i=3;i*i<=m;) qX%oLa
{ if (m%i==0) return 0; Y0?<~Gf
i+=2; k^i\<@v
} "NMSLqO
return 1; Mr}K-C?ge
} r9vO(m~
|z
8Wh
int checkmatrix[ ][3]={ {-1},{0,-1},{1,-1},{0,-1},{1,3,-1}, cvk$ I"q+
{2,4,-1},{3,-1},{4,6,-1},{5,7,-1}}; Ia%S=xU{=
int selectnum(int start) {Y6;/".DM
{ int j; sWGc1jC?.F
for (j=start;j<=N;j++) `JZ`j7f
if (b[j]) return j IwZe2$f
return 0; ] v8 .ym
} j` 5K7~hv
5<RZht$i
int check(int pos) Fu$JI8
{ int i,j; huTWoMU
if (pos<0) return 0; n ]<>$
for (i=0;(j=checkmatrix[pos])>=0;i++) ~6!TMVr
if (!isprime(a[pos]+a[j]) 5f-eWW]!
return 0; tXg>R _\C
return 1; mNDd>4%H_
} z:JQ3D7/we
$OP w$
int extend(int pos) ^)%TQ.
{ a[++pos]=selectnum(1); =K8z8K?
b[a][pos]]=0; t
\;,$i
return pos; {~0r3N4Zl
} }}q_QD_
Xt$o$V
int change(int pos) _@_w6Rh
{ int j; 'g#EBy
while (pos>=0&&(j=selectnum(a[pos]+1))==0) 7|Bg--G1
b[a[pos--]]=1; KwxO%/-}S
if (pos<0) return –1 L^%jR=
b[a[pos]]=1; NU/:jr.W#
a[pos]=j; ,5Nf9z!hk(
b[j]=0; ^,sKj-
return pos; '(-SuaH49
} nm{J
;+NU;f/WM
void find() Qs
#7<NQ
{ int ok=0,pos=0; wxW\L!@
a[pos]=1; Sf"]enwB
b[a[pos]]=0; w\`u|f;Aq
do { 2(s-8E:
if (ok) t`
f.HJe
if (pos==8) fKkH
[
{ write(a); d'UCPg<Y
pos=change(pos); ;%V)lP "o
} E%np-is{1
else pos=extend(pos); l0@+&Xj
else pos=change(pos); d>k"#|
ok=check(pos); X0&[cyP!
} while (pos>=0) D%,AdR"m
} 7
yp}
#Ks2a):8
void main() N~0~1
WQn
{ int i; ;'7gg]
for (i=1;i<=N;i++) $9Pscu bM4
b=1; %DIZgPd\
find(); W, -fnJk
} rhQv,F9
【问题】 n皇后问题 w^N3Ma
问题描述:求出在一个n×n的棋盘上,放置n个不能互相捕捉的国际象棋“皇后”的所有布局。 ;Q8LA",5d
这是来源于国际象棋的一个问题。皇后可以沿着纵横和两条斜线4个方向相互捕捉。如图所示,一个皇后放在棋盘的第4行第3列位置上,则棋盘上凡打“×”的位置上的皇后就能与这个皇后相互捕捉。 Ce-D^9kC
E@N& Y1t
1 2 3 4 5 6 7 8 x17K8De
× × Kq4b`cn{_
× × × K'u66%wAL
× × × xa]e9u%
× × Q × × × × × ['#3GJz-
× × × {6Qd,CX
× × × ! 1wf/C;=
× × I]vCra
× × i-E/#zni
从图中可以得到以下启示:一个合适的解应是在每列、每行上只有一个皇后,且一条斜线上也只有一个皇后。 FAbl5VW'
求解过程从空配置开始。在第1列至第m列为合理配置的基础上,再配置第m+1列,直至第n列配置也是合理时,就找到了一个解。接着改变第n列配置,希望获得下一个解。另外,在任一列上,可能有n种配置。开始时配置在第1行,以后改变时,顺次选择第2行、第3行、…、直到第n行。当第n行配置也找不到一个合理的配置时,就要回溯,去改变前一列的配置。得到求解皇后问题的算法如下: HCOE'24I
{ 输入棋盘大小值n; S{7 R6,B5
m=0; 5FQtlB9F
good=1; DB>.Uf"
do { '~;vp
if (good) e%K
oecq
if (m==n) ZgH(,g,TU
{ 输出解;
~uRL+<.c
改变之,形成下一个候选解; }n[<$*W^
} YfZ96C[a
else 扩展当前候选接至下一列; TL*8h7.(
else 改变之,形成下一个候选解; } Rs@
good=检查当前候选解的合理性; *@SZ0
} while (m!=0); vzPuk|q3
} =*Ru2
在编写程序之前,先确定边式棋盘的数据结构。比较直观的方法是采用一个二维数组,但仔细观察就会发现,这种表示方法给调整候选解及检查其合理性带来困难。更好的方法乃是尽可能直接表示那些常用的信息。对于本题来说,“常用信息”并不是皇后的具体位置,而是“一个皇后是否已经在某行和某条斜线合理地安置好了”。因在某一列上恰好放一个皇后,引入一个一维数组(col[ ]),值col表示在棋盘第i列、col行有一个皇后。例如:col[3]=4,就表示在棋盘的第3列、第4行上有一个皇后。另外,为了使程序在找完了全部解后回溯到最初位置,设定col[0]的初值为0当回溯到第0列时,说明程序已求得全部解,结束程序运行。 VyWPg7}e
为使程序在检查皇后配置的合理性方面简易方便,引入以下三个工作数组: Z`{ZV5
(1) 数组a[ ],a[k]表示第k行上还没有皇后; 4b<>gpQ
(2) 数组b[ ],b[k]表示第k列右高左低斜线上没有皇后; K! e51P
(3) 数组 c[ ],c[k]表示第k列左高右低斜线上没有皇后; Ubf@"B
棋盘中同一右高左低斜线上的方格,他们的行号与列号之和相同;同一左高右低斜线上的方格,他们的行号与列号之差均相同。 ZCuLgCP?Z
初始时,所有行和斜线上均没有皇后,从第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来确定。细节见以下程序: e=#'rDm
【程序】 rUOl+p_47
# include <stdio.h> *CS2ndp
# include <stdlib.h> Y}UVC|Ef
# define MAXN 20 MF]EX
int n,m,good; ^mZ eAW
int col[MAXN+1],a[MAXN+1],b[2*MAXN+1],c[2*MAXN+1]; H(,D5y`k1
V3t;V-Lkt
void main() F* h\ #?
{ int j; 9?L,DThQ
char awn; 9Atnnx]n
printf(“Enter n: “); scanf(“%d”,&n); j$mCU?
for (j=0;j<=n;j++) a[j]=1; lOJ3_8
for (j=0;j<=2*n;j++) cb[j]=c[j]=1; f'28s*n
m=1; col[1]=1; good=1; col[0]=0; {IR-g,B
do { E3P2
if (good) g+ P
if (m==n) Pt^SlX^MM
{ printf(“列\t行”); . ]D7Il
for (j=1;j<=n;j++) #Rx|oSc}
printf(“%3d\t%d\n”,j,col[j]); iwS55o
printf(“Enter a character (Q/q for exit)!\n”); TeXt'G=M
scanf(“%c”,&awn); /lqVMlz\77
if (awn==’Q’||awn==’q’) exit(0); Eu)(@,]we
while (col[m]==n) 3rh@|fg)E
{ m--; [t }\8^y
a[col[m]]=b[m+col[m]]=c[n+m-col[m]]=1; " _{o}8L
} y3{F\K
col[m]++;
##_Jz 5P
} 6L4<c+v_
else lS!uL9t.
{ a[col[m]]=b[m+col[m]]=c[n+m-col[m]]=0; %{*)-_M
col[++m]=1; IqrT@jgN-
} ~&\} qz3
else =/zb$d cz
{ while (col[m]==n) `+?g96
{ m--; H\!p%Y
a[col[m]]=b[m+col[m]]=c[n+m-col[m]]=1; m. EIMuj
} fq,LXQ#G
col[m]++; `%oJa`
} 66 @#V
good=a[col[m]]&&b[m+col[m]]&&c[n+m-col[m]]; v"3($?au0
} while (m!=0); R m^$Dn
} uW4wTAk;qh
试探法找解算法也常常被编写成递归函数,下面两程序中的函数queen_all()和函数queen_one()能分别用来解皇后问题的全部解和一个解。 v62O+{
【程序】 Z36C7 kw
# include <stdio.h> |ryV7VJ8
# include <stdlib.h> <A+n[h
# define MAXN 20 W3aFao>!OZ
int n; jS<(Oo
int col[MAXN+1],a[MAXN+1],b[2*MAXN+1],c[2*MAXN+1]; %f'mW2
void main() bjCO@t
{ int j; >A_:qyGk
printf(“Enter n: “); scanf(“%d”,&n); |18h
p
for (j=0;j<=n;j++) a[j]=1; 9qcA+gz:|
for (j=0;j<=2*n;j++) cb[j]=c[j]=1;
o.|P7{v}
queen_all(1,n); u zgQ_
} JDp{d c
yMVlTO
void queen_all(int k,int n) O~~WP*N
{ int i,j; RF$2p4=[
char awn; |X6/Y@N
for (i=1;i<=n;i++) 9?J
3G,&
if (a&&b[k+i]&&c[n+k-i]) _`- trE.
{ col[k]=i; X" m0||
a=b[k+i]=c[n+k-i]=0; *}<