基础算法系列总结:回溯算法(解火力网问题)
理论辅助:
回溯算法也叫试探法,它是一种系统地搜索问题的解的方法。回溯算法的基本思想是:从一条路往前走,能进则进,不能进则退回来,换一条路再试。用回溯算法解决问题的一般步骤为:
1、定义一个解空间,它包含问题的解。
2、利用适于搜索的方法组织解空间。
3、利用深度优先法搜索解空间。
4、利用限界函数避免移动到不可能产生解的子空间。
问题的解空间通常是在搜索问题的解的过程中动态产生的,这是回溯算法的一个重要特性。
还是那个基调,不喜欢纯理论的东西,喜欢使用例子来讲诉理论,在算法系列总结:动态规划(解公司外包成本问题) 的那一节里面 我们举得是经典的0-1背包问题,在回溯算法里面也有一些很经典的问题,当然,动态规划的0-1背包问题其实也可以使用回溯算法来解。在诸如此类似的求最优解的问题中,大部分其实都可以用回溯法来解决,可以认为回溯算法一个”通用解题法“,这是由他试探性的行为决定的,就好比求一个最优解,我可能没有很好的概念知道怎么做会更快的求出这个最优解,但是我可以尝试所有的方法,先试探性的尝试每一个组合,看看到底通不通,如果不通,则折回去,由最近的一个节点继续向前尝试其他的组合,如此反复。这样所有解都出来了,在做一下比较,能求不出最优解吗?
例子先行,现在我们来看看经典的N后问题
问题描述:在n*n格的棋盘上放置彼此不受攻击的n个皇后。按照国际象棋的规矩,皇后可以攻击与之处在同一行或同一列或同一斜线上的棋子。n后问题等价于在n*n格的棋盘上方置n个皇后,任何2个皇后不放在同一行或同一列或同一斜线上。我们需要求的是可放置的总数。
基本思路: 用n元组x[1;n]表示n后问题的解。其中,x[i]表示皇后i放置在棋盘的第i行的第x[i]列。由于不容许将2个皇后放在同一列上,所以解向量中的x[i]互不相同。2个皇后不能放在同一斜线上是问题的隐约束。对于一般的n后问题,这一隐约束条件可以化成显约束的形式。如果将n*n 格的棋盘看做二维方阵,其行号从上到下,列号从左到右依次编号为1,2,...n。从棋盘左上角到右下角的主对角线及其平行线(即斜率为-1的各斜线)上,2个下标值的差(行号-列号)值相等。同理,斜率为+1的每条斜线上,2个下标值的和(行号+列号)值相等。因此,若2个皇后放置的位置分别是(i,j)和(k,l),且 i-j = k -l 或 i+j = k+l,则说明这2个皇后处于同一斜线上。以上2个方程分别等价于i-k = j-l 和 i-k =l-j。由此可知,只要|i-k|=|l-j|成立,就表明2个皇后位于同一条斜线上。
代码:
#include <math.h>
#include<stdlib.h>
static int n,x[1000 ];
static long sum;
int Place(int k)
{
for(int j=1;j <k; j++ )
if((abs(k-j) == abs(x[j]-x[k]))||(x[j]==x[k])) return 0 ;
return 1 ;
}
void Backtrak(int t)
{
if(t>n) sum++ ;
else
for(int i=1; i <= n; i++ )
{
x[t] = i;
if(Place(t))Backtrak(t+1 );
}
}
int main()
{
int nn;
while(scanf("%d",&nn)!= EOF)
{
n= nn;
sum=0 ;
for(int i=0;i<=n;i++ )
x[i]=0 ;
Backtrak(1 );
printf("%d\n" ,sum);
}
}
这段代码有必要解释一下,Place(int)即尝试看是否可以,如果不可以则回退到t+1层,再尝试其他的组合。
这里也道出了回溯算法的核心思想:但当探索到某一步时,发现原先选择并不优或达不到目标,就退回一步重新选择
算法实践:
问题描述:在一个n*n的网格里,每个网格可能为“墙壁”(用‘X’表示)和“街道”(用‘.’表示)。现在在街道放置碉堡,每个碉堡可以向上下左右四个方向开火,子弹射程无限远。墙壁可以阻挡子弹。问最多能放置多少个碉堡,使它们彼此不会互相摧毁。
如下面四张图,墙壁用黑正方形表示,街道用空白正方形表示,圆球就代表碉堡。1,2,3是正确的,4,5是错误的。以为4,5里面在某一行或者某一列有两个碉堡,这样他们就会互相攻击了。意思明白了吗?可能我的表达很不清晰,呵呵….
4 ——————输入的n值
.X..
....
XX..
....
2
XX
.X
3
.X.
X.X
.X.
4
....
....
....
....
Sample output:
5
1
5
4
初拿到这个问题,你会不会想到回溯算法呢?有人说遍历墙的位置,然后再墙的上下左右四个格子放置碉堡会得到最优解,这个我没有验证过,细细的用笔画了画,好像是这么回事,但是很多时候要知道最优解用什么方法是很难发现的,利用通用解题方法回溯法,我们可以在一片茫然的时候开始我们的编程
首先我们来分析一下这个问题:使用回溯法,我们尝试每一种可能放置的情况,然后进行判断是否满足要求,若不满足,尝试放到下一个单元格,如此反复,最终,我们将所有可能放置的情况全部遍历出来了,连所有情况都出来了,难不成还找不到最优解吗?哈哈。。说做就做…
char map[ 4 ][ 4 ];
int best,n;
int canput( int row, int col)
{
int i;
for (i = row - 1 ; i >= 0 ; i -- )
{
if (map[i][col] == ' o ' ) return 0 ;
if (map[i][col] == ' x ' ) break ;
}
for (i = col - 1 ; i >= 0 ; i -- )
{
if (map[row][i] == ' o ' ) return 0 ;
if (map[row][i] == ' x ' ) break ;
}
return 1 ;
}
void solve( int k, int tot)
{
int x,y;
if (k == n * n)
{
if (tot > best)
{
best = tot; return ;
}
}
else
{
x = k / n;
y = k % n;
if ((map[x][y] == ' . ' ) && (canput(x,y) ) )
{
map[x][y] = ' o ' ;
solve(k + 1 ,tot + 1 );
map[x][y] = ' . ' ;
}
solve(k + 1 ,tot);
}
}
int main()
{
int i,j;
scanf( " %d " , & n);
while (n > 0 )
{
for (i = 0 ;i < n;i ++ )
for (j = 0 ;j < n;j ++ )
scanf( " %1s " , & map[i][j]);
best = 0 ;
solve( 0 , 0 );
printf( " %d\n " ,best);
n = 0 ;
scanf( " %d " , & n);
}
return 0 ;
}
对上面的代码做一下点解释,canput是做检验的,检验放在某个地点到底行不行得通,solve才是真正进行递归回溯的函数。。
算法系列目录: