IDA*搜索算法学习笔记

IDA*

紫书第207页；

紫书在介绍迭代加深搜索时提到了一种搜索算法IDA*，个人认为IDA算法是A算法和迭代加深算法的结合体；

总的来说，IDA*算法的基本思路就是：

1.主函数枚举最大搜索树的层数，如果再当前的层数中可以出现答案，则这个层数就为答案的最小层数；

2.dfs函数一般有两个参数，一个是当前层数，一个是最大层数；dfs函数的前面部分一般是判断当前剩余层数是否可以得出解（这个就是俗称的“乐观估价函数”）；如果不能得出解，则直接减枝；

总结：

IDA*算法弥补了dfs只能找出解，不能找出最优解的缺陷，同时又比单向的bfs快许多；唯一的难点就是“乐观估价函数”并不好求；

---

紫书第210页的一道题就用到了IDA*算法，比bfs更加简洁，更快；

The Rotation Game

这道题的难点在于转向的实现，运用先预处理棋盘，再进行操作，更加简洁；

还有个关键点是“乐观估价函数”，因为一次转，最多只能改变中间8个数字的1个，当1，2，3的数量加上剩余的转向次数都不大于8时，说明无法得出解，剪枝；

代码：

#include<bits/stdc++.h>
#define LL long long
#define pa pair<int,int>
#define lson k<<1
#define rson k<<1|1
#define inf 0x3f3f3f3f
//ios::sync_with_stdio(false);
using namespace std;
const int N=200100;
const int M=1000100;
const LL mod=998244353;
int num[30];
int ans[8]={6,7,8,11,12,15,16,17};
int ma[8][7]={ {0,2,6,11,15,20,22},
                {1,3,8,12,17,21,23},
                {10,9,8,7,6,5,4},
                {19,18,17,16,15,14,13},
                {23,21,17,12,8,3,1},
                {22,20,15,11,6,2,0},
                {13,14,15,16,17,18,19},
                {4,5,6,7,8,9,10}};
int ok;
int b;
char path[N];
int s1,s2,s3;
bool judge(){
	s1=s2=s3=0;
	for(int i=0;i<8;i++){
		if(num[ans[i]]==1) s1++;
		if(num[ans[i]]==2) s2++;
		if(num[ans[i]]==3) s3++;
	}
	if(s1==8){
		b=1;
		return true;
	}
	if(s2==8){
		b=2;
		return true;
	}
	if(s3==8){
		b=3;
		return true;
	}
	return false;
}
void move(int p){
	int t=num[ma[p][0]];
	for(int i=0;i<6;i++) num[ma[p][i]]=num[ma[p][i+1]];
	num[ma[p][6]]=t;
}
void dfs(int p,int q){//p是当前层数，q是最大层数
	if(ok) return;
	if(judge()){//符合 
		ok=1;
		return;
	}
	if(q-p+s1<8&&q-p+s2<8&&q-p+s3<8) return;//直接退出，找不到符合条件的方法 
	int num1[30];//保留num
	memcpy(num1,num,sizeof(num));
	for(int i=0;i<8;i++){
		if(ok) return;
		move(i);
		path[p]=i+'A';
		dfs(p+1,q);
		memcpy(num,num1,sizeof(num1));
	}
	return;
}
int main(){
	ios::sync_with_stdio(false);
	while(cin>>num[0]){
		if(num[0]==0) break;
		ok=0;
		for(int i=1;i<=23;i++) cin>>num[i];
		int a=0;
		for(int i=0;;i++){
			dfs(0,i);
			if(ok){
				a=i;
				break;
			}
		}
		if(!a){
			cout<<"No moves needed"<<endl;
			cout<<b<<endl;
		}
		else{
			for(int i=0;i<a;i++) cout<<path[i];
			cout<<endl;
			cout<<b<<endl;
		}
	}
	return 0;
}