算法-- 回溯算法 及 示例

回溯法介绍

回溯法(英语:backtracking)也称试探法,回溯法有“通用的解题方法”之称。它可以系统的搜索一个问题的所有解或者任意解。

回溯法是一个既带有系统性又带有跳跃性的的搜索算法。它在包含问题的所有解的解空间树中,按照深度优先的策略,从根结点

出发搜索解空间树。算法搜索至解空间树的任一结点时,总是先判断该结点是否肯定不包含问题的解。如果肯定不包含,则跳过

对以该结点为根的子树的系统搜索,逐层向其祖先结点回溯。否则,进入该子树,继续按深度优先的策略进行搜索。回溯法在用来

求问题的所有解时,要回溯到根,且根结点的所有子树都已被搜索遍才结束。而回溯法在用来求问题的任一解时,只要搜索到问题

的一个解就可以结束。这种以深度优先的方式系统地搜索问题的解的算法称为回溯法,它适用于解一些组合数较大的问题.

回溯法通常用最简单的递归方法来实现,在反复重复上述的步骤后可能出现两种情况:

1 .找到一个可能存在的正确的答案。

2. 在尝试了所有可能的分步方法后宣告该问题没有答案。

在最坏的情况下,回溯法会导致一次复杂度为指数时间的计算。

适用范围

在包含问题的所有解的解空间树中,按照深度优先搜索的策略,从根结点出发深度探索解空间树。

 适用范围:

1,问题的解用向量表示

  X = (x1, x2, …, xn)

2,需要搜索一个或一组解

3,满足约束条件的最优解

回溯法的基本思想

对于用回溯法求解的问题,首先要将问题进行适当的转化,得出状态空间树。这棵树的每条完整路径都代表了一种解的可能。通过深度优先搜索

这棵树,枚举每种可能的解的情况;从而得出结果。但是,回溯法中通过构造约束函数,可以大大提升程序效率,因为在深度优先搜索的过程中,

不断的将每个解(并不一定是完整的,事实上这也就是构造约束函数的意义所在)与约束函数进行对照从而删除一些不可能的解,这样就不必继续

把解的剩余部分列出从而节省部分时间。

回溯法中,首先需要明确下面三个概念

 1,约束函数:约束函数是根据题意定出的。通过描述合法解的一般特征用于去除不合法的解,从而避免继续搜索出这个不合法解的剩余部分。因此,

       约束函数是对于任何状态空间树上的节点都有效、等价的。

2,状态空间树:刚刚已经提到,状态空间树是一个对所有解的图形描述。树上的每个子节点的解都只有一个部分与父节点不同。

3,扩展节点、活结点、死结点:所谓扩展节点,就是当前正在求出它的子节点的节点,在DFS中,只允许有一个扩展节点。活结点就是通过与约束函数

      的对照,节点本身和其父节点均满足约束函数要求的节点;死结点反之。由此很容易知道死结点是不必求出其子节点的(没有意义)。

利用回溯法解题的具体步骤

首先,要通过读题完成下面三个步骤

(1)描述解的形式,定义一个解空间,它包含问题的所有解,这一步主要明确问题的解空间树。

(2)构造状态空间树。

(3)构造约束函数(用于杀死节点)。

然后就要通过DFS思想完成回溯,具体流程如下:

(1)设置初始化的方案(给变量赋初值,读入已知数据等)。

(2)变换方式去试探,若全部试完则转(7)。

(3)判断此法是否成功(通过约束函数),不成功则转(2)。

(4)试探成功则前进一步再试探。

(5)正确方案还未找到则转(2)。

(6)已找到一种方案则记录并打印。

(7)退回一步(回溯),若未退到头则转(2)。

(8)已退到头则结束或打印无解。

总结起来就是:

针对所给问题,确定问题的解空间 –> 确定结点的扩展搜索规则–> 以DFS方式搜索解空间,并在搜索过程中用剪枝函数避免无效搜索。

回溯法基本框架

设问题的解是一个n维向量(a1,a2,………,an),约束条件是bi(i=1,2,3,…..,n)之间满足某种条件,记为f(bi)。

递归实现

int a[n];
初始化数组a[]等操作  
 ...
void dfs(int cur)
{
	int i;
     if(cur>n)
       //统计、输出结果等;
      else
     {
        for(i = 下界; i <= 上界; ++i)  // 枚举i所有可能的路径
        {
           if(fun(i))                 // 满足限界函数和约束条件
              {
                 a[cur] = i;
                  ...                 // 其他操作,设置标志等
                 dfs(cur+1);
                 //回溯前的清理工作(如a[cur]置空值,标志置0等);
              }
          }
      }
}

典型问题应用

素数环问题,hdoj1016:http://acm.hdu.edu.cn/showproblem.php?pid=1016

题目大意:

将从1到n这n个整数围成一个圆环,若其中任意2个相邻的数字相加,结果均为素数,那么这个环就成为素数环。

要求输出:从整数1开始。

分析问题,可以构造解空间树,比较顺利的想到 DFS 、回溯。

代码如下:

#include<stdio.h>
#include <string.h>

#define M 40

int isPrime[M];////素数表,下标为素数的置为1,否则0
int vis[M>>1];// vis 标识 1-n,是否被选
int res[M>>1];// 存储解向量

int cnt;// 测试样例个数

void prime()//求出1-40的所有素数
{
	int i, j;
	for(i=1; i<M; ++i)
	{
		int ok = 1;
		for(j=2; j*j<=i; ++j)
		{
			if(i%j == 0)
			{
				ok = 0;
				break;
			}
		}
		if(ok)
			isPrime[i]=1;

	}
}

void dfs(int cur, int n)
{
	int i;
	if(cur == n && isPrime[res[n-1] + res[0]])//别忘了测试边界,最后一个和第一个数 构成的环
	{
		for(i=0; i<n-1; ++i)
			printf("%d ", res[i]);
		printf("%d\n", res[i]);
	}
	else
	{
		for(i=2; i<=n; i++)// 尝试每个i, 1始终在排头,因此从2开始计算
		{
			if(!vis[i] && isPrime[res[cur-1] + i])// i未用过且和前一个数和为素数
			{
				res[cur] = i;
				vis[i] = 1; // 设置标志
				dfs(cur+1, n);
				vis[i] = 0; // 回溯, 清除标识
			}
		}
	}
}


int main()
{
	int n;
	//freopen("in.txt", "r", stdin);
	prime();
	cnt = 0;
	while(scanf("%d", &n) != EOF)
	{
		++cnt;
		printf("Case %d:\n", cnt);
		memset(vis, 0, sizeof(vis));
		res[0] = 1;
		dfs(1, n);
		printf("\n");
	}
	return 0;
}


代码中 prime()函数 求出1-40的所有素数,因为只要测试 1-19 因此 可以事前 把1-38的素数存储到一个 素数表里。这样计算时间更快。

下面是最近更新的回溯法的相关例子,http://blog.csdn.net/daniel_ustc/article/details/17041933

参考:

 http://zh.wikipedia.org/zh-cn/%E5%9B%9E%E6%BA%AF%E6%B3%95

http://www.newsjz.com/wxqgr/xxas/as2/sfjc/201103/337.html

    原文作者:回溯法
    原文地址: https://blog.csdn.net/daniel_ustc/article/details/17040315
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞