转载请注明出处,谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents
by---cxlove
继续IDA*搜索,估价函数H仍然是曼哈顿距离,每一次转换会改变4个位置的曼哈顿距离,分别改变1,所以把曼哈顿距离和+3/4便可以作为H函数,表示至少需要多少步,一个DFS的剪枝。
这题最多九步,BFS应该也无压力
可惜没有优化到0MS。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define LD long double
#define LL __int64
#define M 200005
using namespace std;
int T,a[9];
int depth;
char str[10];
bool flag;
int rotation[4][4]={{0,1,4,3},{1,2,5,4},{3,4,7,6},{4,5,8,7}};
int get_h(int *b){
int ans=0;
for(int i=0;i<9;i++)
ans+=abs(i/3-(b[i]-1)/3)+abs(i%3-(b[i]-1)%3);
return (ans+3)/4;
}
void change(int *b,int kind){
if(kind&1){
kind/=2;
int tmp=b[rotation[kind][3]];
for(int i=3;i>0;i--)
b[rotation[kind][i]]=b[rotation[kind][i-1]];
b[rotation[kind][0]]=tmp;
}
else{
kind/=2;
int tmp=b[rotation[kind][0]];
for(int i=1;i<4;i++)
b[rotation[kind][i-1]]=b[rotation[kind][i]];
b[rotation[kind][3]]=tmp;
}
}
void IDAstar(int *b,int tmp_depth,int pre){
if(flag)
return;
if(get_h(b)>tmp_depth)
return;
if(tmp_depth==0&&get_h(b)==0){
flag=true;
return;
}
for(int i=0;i<8;i++){
if(pre>=0&&pre/2==i/2&&(pre%2)^(i%2))
continue;
int tmp[9];
for(int j=0;j<9;j++)
tmp[j]=b[j];
change(tmp,i);
IDAstar(tmp,tmp_depth-1,i);
}
}
int main(){
int cas=0;
while(scanf("%s",str)!=EOF&&strcmp(str,"0000000000")){
T=str[0]-'0';
for(int i=0;i<9;i++)
a[i]=str[i+1]-'0';
flag=false;
for(depth=get_h(a);depth<=T;depth++){
IDAstar(a,depth,-1);
if(flag){
printf("%d. %d\n",++cas,depth);
break;
}
}
if(!flag)
printf("%d. -1\n",++cas);
}
return 0;
}
分享到:
相关推荐
中南大学ACM模板, ACM 模板 JPVision Fighting! To be or not to be , that is a question. alpc48 2008-10-5
hdu acm 教案 搜索入门 hdu acm 教案 搜索入门
HDU的1250,主要是利用高精度加法,但是代码有点繁琐,效率不是很高
杭电ACMhdu1163
HDU1059的代码
hdu1001解题报告
搜索 dfs 解题代码 hdu1241
hdu 1574 passed sorce
HDU的一题........HDU DP动态规
hdu2101AC代码
杭州电子科技大学oj平台上的第1010题,是关于搜索的题目,很不错的题
hdu 5007 Post Robot 字符串枚举。 暴力一下就可以了。
hdu acm 教案 动态规划(1) hdu acm 教案 动态规划(1)
hdu 1166线段树代码
自己做的HDU ACM已经AC的题目
ACM HDU题目分类,我自己总结的大概只有十来个吧
HDU最全ac代码
hdu动态规划算法集锦
hdu题目分类
HDU图论题目分类