深度优先搜索

深度优先搜索

        搜索算法是全国信息学竞赛必考的算法之一，也是最常见的算法之一。很多复杂的算法都在搜索算法的基础上发展起来的。在竞赛中，有一些问题一般难以找到公式或规律，特别是信息学竞赛试题难度较大，根本就没有公式可循，这时一般采用搜索法利用计算机高速运行的特点来编程求解，如何分析此题是否需要用搜索算法来解，这就显得关键。我们知道，所谓搜索，就是在数据集合中寻找满足某种条件的数据对象。
         搜索算法是利用计算机的高性能来有目的的穷举一个问题的部分或所有的可能情况，从而求出问题的解的一种方法。
         搜索算法被称为“通用解题法”，在算法和人工智能中占有重要地位。很多算法都是搜索算法的扩展，如图论中的某些算法、动态规划、贪心等，只是他们找到了某些规律来进行变相剪枝而已。 我们常用的搜索算法有深度优先搜索depth first search（dfs）、广度优先搜索breadth first search（bfs）和A*算法。其中A*算法属于启发式搜索。

深度优先搜索的概念

        按照深度优先的顺序遍历状态空间，通常用递归或者栈来实现。说白了就是往深里搜，当某个节点的所有子节点都搜索过了就返回父节点继续进行。就像走迷宫，你没有办法用分身术来站在每个走过的位置。不撞南山不回头。

    假设按照以下的顺序来搜索：
         （1） v0 -> v1 ->v4，此时已经到尽头，仍然到不了v6，于是原路返回到v1去搜索其它路径；
         （2） 返回到v1后继续搜索v2，于是搜索路径是v0 -> v1 -> v2 -> v6，找到目标节点，返回有解。
         这样，搜索2步就到达了。这就是dfs的过程。搜索的过程，实际上就是建立一棵树的过程。如果图1全部搜索完，构建的树如下：
    
    图1的dfs伪代码如下：

/*dfs核心伪代码 */

void dfs(k,…){    //k代表目前dfs的深度
	if 找到解 or 走不下去了  { …. return; }
	for i=1 to n do     //解空间
	……
	枚举下一种情况，dfs(k+1,…)
}

【程序1】特殊的质数肋骨

         农民约翰母牛总是产生最好的肋骨。
         你能通过农民约翰和美国农业部标记在每根肋骨上的数字认出它们。
         农民约翰确定他卖给买方的是真正的质数肋骨,是因为从右边开始切下肋骨,每次还剩下的肋骨上的数字都组成一个质数,举例来说:
         7 3 3 1
         全部肋骨上的数字 7331是质数;三根肋骨 733是质数;二根肋骨 73 是质数;当然,最后一根肋骨 7 也是质数。
         7331 被叫做长度 4 的特殊质数。
         写一个程序对给定的肋骨的数目 N (1<=N<=8),求出所有的特殊质数。
         数字1不被看作一个质数。
输入格式
         单独的一行包含N。
输出格式
         按顺序输出长度为 N 的特殊质数,每行一个。
样例输入
         4
样例输出
         2333
         2339
         2393
         2399
         2939
         3119
         3137
         3733
         3739
         3793
         3797
         5939
         7193
         7331
         7333
         7393

问题分析：
         样例中n=4，即四位数中具有这种性质的数。如果直接搜索枚举，n稍稍大一点，就会超时。那么我们可以直接构造，第一位数字只能是2，3，5，7中的1个，因为其它的都不是素数。我们进行适当的剪枝，缩小解空间，除第一位外，其它位的数字不可能是偶数，也不可能是5，只能是1,3,7,9。构造的搜索树如下（n=3为例）：

         深度优先搜索的过程就是建立一棵搜索树的过程，上图中的根节点0是不存在的，添加以方便表示。
子程序dfs(k)框架如下：

void dfs(x, k){  //搜索第k位上可能的数字，构造前k位数x，此处k为目前dfs的深度
	if  x不是素数 then 返回;
	if  x=n then 输出;
	dfs(x*10+1,k+1);
	dfs(x*10+3,k+1);
	dfs(x*10+7,k+1);
	dfs(x*10+9,k+1);
}

参考代码：
#include<iostream>
using namespace std;
int n;
bool prime(int m) {
	for (int i=2;i*i<=m;++i) 
		if (m%i==0) return false;
	return true;
}