算法笔记：高级搜索

一.折半搜索

这类题目标志都一般很明显，物品的种类n>=40，物品特性（比如价格、重量）m>=10^9。因此，使用背包dp、状态dp等dp都因为空间太大而不行，直接暴搜时间复杂度有不允许。

将搜索物品分为两半，搜索前面一半，将搜索结果记录并排好序；再搜索剩下一半，使用二分查找排好序的记录的搜索结果。

例题：

洛谷：P4799 [CEOI2015 Day2] 世界冰球锦标赛

给出 Bobek 的预算和每场比赛的票价，试求：如果总票价不超过预算，他有多少种观赛方案。如果存在以其中一种方案观看某场比赛而另一种方案不观看，则认为这两种方案不同。

输入格式

第一行，两个正整数 N 和 M(1≤N≤40,1≤M≤10^18)，表示比赛的个数和 Bobek 那家徒四壁的财产。

第二行，N 个以空格分隔的正整数，均不超过 10^16，代表每场比赛门票的价格。

输出格式

输出一行，表示方案的个数。由于 N 十分大，注意：答案 ≤2^40。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
int n;
ll m,ans = 0;
ll p[41];
int cnt = 0;
ll w[1<<20];
void dfs1(int cur,ll s){
	if(cur>=n/2){
		w[cnt++] = s;
		return ;
	}
	dfs1(cur+1,s);
	if(s+p[cur]<=m){
		dfs1(cur+1,s+p[cur]);
	}
	
}
void dfs2(int cur,ll s){
	if(cur>=n){
		ans += (upper_bound(w,w+cnt,m-s)-w);
		return ;
	}
	dfs2(cur+1,s);
	if(s+p[cur]<=m){
		dfs2(cur+1,s+p[cur]);
	}
	
}
int main(){
	cin>>n>>m;
	for(int i = 0;i<n;i++){
		cin>>p[i];
	}
	//sort(p,p+n);
	dfs1(0,0);
	sort(w,w+cnt);
	dfs2(n/2,0);
	cout<<ans;
	
}

二.双向搜索

对于无法用动态规划、搜索剪枝，同时暴搜超时的题目，可以尝试双向搜索，即开始状态和终止状态同时开始搜索，如果达到一个相同的状态说明搜索完成。

相当将一颗搜索树变成两颗，如果两颗搜索树的某个节点（状态）相遇（状态相同），说明搜索结束。这样可以有效减少搜索树的大小。

洛谷P1379八数码难题

题目描述

在 3×3的棋盘上，摆有八个棋子，每个棋子上标有 1至 8 的某一数字。棋盘中留有一个空格，空格用 00 来表示。空格周围的棋子可以移到空格中。要求解的问题是：给出一种初始布局（初始状态）和目标布局（为了使题目简单,设目标状态为 123804765），找到一种最少步骤的移动方法，实现从初始布局到目标布局的转变。

输入格式

输入初始状态，一行九个数字，空格用 0 表示。

输出格式

只有一行，该行只有一个数字，表示从初始状态到目标状态需要的最少移动次数。保证测试数据中无特殊无法到达目标状态数据。

#include<bits/stdc++.h>
using namespace std;
int tar = 123804765;
int sta;
unordered_map<int, int>um1, um2;
int ans = 0;
int temp1[9];//num -> pos
int temp2[9];//pos ->num
int nums[9] = {1, 10, 100, 1000, 10000, 100000, 1000000, 10000000, 100000000};
vector<int> um[9] = {{1, 3}, {0, 2, 4}, {5, 1}, {0, 4, 6}, {1, 3, 5, 7}, {2, 4, 8}, {3, 7}, {4, 6, 8}, {5, 7}};//记录八数码中每个格字挨着哪些格子
queue<int>q1, q2;
int bfs() {
	while (true) {
        //从开始状态搜索
		int t = q1.front(), cur = t;
		q1.pop();
		for (int i = 0; i < 9; i++) {
			if (t == 0) {
				temp1[0] = 8;
				temp2[8] = 0;
			} else {
				temp1[t % 10] = i;
				temp2[i] = t % 10;
				t /= 10;
			}

		}
		for (int i = um[temp1[0]].size() - 1; i >= 0; i--) {
			int num = cur + (nums[temp1[0]] - nums[um[temp1[0]][i]]) * temp2[um[temp1[0]][i]];
			if (um2.count(num)) {
                //遇到了相同状态，搜索结束
				return um2[num] + um1[cur] + 1;
			}
			if (!um1.count(num)) {
				q1.push(num);
				um1[num] = um1[cur] + 1;

			}
		}

        //从终止状态搜索
		t = q2.front();
		cur = t;
		q2.pop();
		for (int i = 0; i < 9; i++) {
			if (t == 0) {
				temp1[0] = 8;
				temp2[8] = 0;
			} else {
				temp1[t % 10] = i;
				temp2[i] = t % 10;
				t /= 10;
			}

		}
		for (int i = um[temp1[0]].size() - 1; i >= 0; i--) {
			int num = cur + (nums[temp1[0]] - nums[um[temp1[0]][i]]) * temp2[um[temp1[0]][i]];
			if (um1.count(num)) {
                //遇到了相同状态，搜索结束
				return um1[num] + um2[cur] + 1;
			}
			if (!um2.count(num)) {
				q2.push(num);
				um2[num] = um2[cur] + 1;

			}
		}

	}
	return 0;
}
int main() {
	cin >> sta;
	if (sta == tar) {
		cout << 0;
		return 0;
	}
	um2[tar] = 0;
	um1[sta] = 0;
	q1.push(sta);
	q2.push(tar);

	cout << bfs();
}

三.启发式搜索（A*）

未完待续