poj1324-Holedox Moving-状态压缩+BFS

朴素暴力方法做的。。。

用 一个1 2 3 4表示 当前蛇身与 上一部分蛇身的位置关系，蛇最长为7，所以可以 用 4进制 表示状态，每一个蛇头 ，其蛇身有 4^7种状态

用状态压缩 标记出现过的状态，，以免重复计算TLE。。。。

别的地方就是朴素的BFS。。。。

写得太挫。。状态数组要开 bool才能过，不然tle

当然本题还可以用双向bfs，A*算法。。。

 在朴素算法的基础上用bitset表示状态 节省了空间 3000kb，时间还是差不多 3.8s 

 

以下是朴素方法  ：  3.9s... 10520kb

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <iostream>
#include <queue>
#include <map>
#include <set>
#include <vector>
using namespace std;
struct node
{
	int hx,hy;
	int body[8]; //1 2 3 4 分别代表上下左右 
	int step;
};
int len;
int n,m,l;
queue<node> qq;
int dirx[]={0,0,1,-1}; //you  zuo xia shang
int diry[]={1,-1,0,0};
int mp[22][22];
bool state[22][22][1<<14];
int trans(int x,int y)			//计算del_x,del_y代表的方位
{
	if (x==-1&&y==0)
		return 1;  //shang 
	if (x==1&&y==0)	
		return 2; //下
	if (x==0&&y==-1)
		return 3 ;		//左
	if (x==0&&y==1)
		return 4;	//右
}
node move_to(node t,int x,int y)		//蛇移动到[X,Y]后整体坐标的变化
{
	node res;
	int i;
	for (i=len;i>1;i--)
	{
		res.body[i]=t.body[i-1];
	}
	int dir=trans(t.hx-x,t.hy-y);
 
		res.body[1]=dir;  
	res.step=t.step+1;
	res.hx=x;
	res.hy=y;
	return res;
	
}
inline void get_xy(int &tx,int &ty,int dir,int x,int y)		//得到坐标值
{
	if (dir==1)//上
	{
		tx=x-1;ty=y;
	}
	if (dir==2)//xia
	{
		tx=x+1;ty=y;
	}
	if (dir==3)//zuo
	{
		tx=x;ty=y-1;
	}
	if (dir==4)//you
	{
		tx=x;ty=y+1;
	}
 return ;

}
void deal_map(node t)		//把蛇身标记为障碍
{
	int i;
	int lastx,lasty;
	for (i=1;i<=len;i++)
	{
		int tx,ty;
		if (i==1)
		{
			get_xy(tx,ty,t.body[i],t.hx,t.hy);
		}
		else
		{
			get_xy(tx,ty,t.body[i],lastx,lasty);
		}
			lastx=tx;
			lasty=ty;
		if (mp[tx][ty]==0)
			mp[tx][ty]=2;	 
	}
		
}
void un_deal_map(node t)		//去掉 蛇身记号
{
int i;
	int lastx,lasty;
	for (i=1;i<=len;i++)
	{
		int tx,ty;
		if (i==1)
		{
			get_xy(tx,ty,t.body[i],t.hx,t.hy);  
		}
		else
		{
			get_xy(tx,ty,t.body[i],lastx,lasty);
		}
			lastx=tx;
			lasty=ty;
		if (mp[tx][ty]==2)
			mp[tx][ty]=0;	 
	}
}
int get_body_num(node &y)  //计算该状态特征值
{
	int i;
	int ret=0;
	for (i=1;i<=len;i++)
	{
	//	ret+=(y.body[i]-1)*  int(pow(4.0,i-1)); 下一行与本行等价
	 	ret+=(y.body[i]-1)*  1<<(2*(i-1)) ;  //位运算加速
	}
	return ret;
}

int main()
{
	int cnt=1;
	while(scanf("%d%d%d",&n,&m,&l)!=EOF)
	{
		if (!n&&!m&&!l) break;
		while(!qq.empty())
			qq.pop();
		node tm;
		memset(mp,0,sizeof(mp)); 
		 memset(state,0,sizeof(state));
		int i; 
		scanf("%d%d",&tm.hx,&tm.hy);
		for (i=0;i<=m+1;i++)   //lock the board 边界设置为障碍
			mp[0][i]=1;
		for (i=0;i<=m+1;i++)
			mp[n+1][i]=1;
		for (i=0;i<=n+1;i++)
			mp[i][0]=1;
		for (i=0;i<=n+1;i++)
			mp[i][m+1]=1;  
			int xx,yy;
		 
		int lastx,lasty;
		int delx,dely;
		for (i=1;i<=l-1;i++)
		{
			scanf("%d%d",&xx,&yy);
			if (i==1) 
			{
				delx=xx-tm.hx;
				dely=yy-tm.hy;
			}
			else
			{
				delx=xx-lastx;
				dely=yy-lasty;
			}
			tm.body[i]=trans(delx,dely);		//记录每个部分是在前一部分的什么方向
			lastx=xx;
			lasty=yy;
		}
		len=l-1;
		int num;	
	
		scanf("%d",&num);
		for (i=1;i<=num;i++) 
		{
			scanf("%d%d",&xx,&yy);
			mp[xx][yy]=1;			//屏蔽障碍
		}


	 
		tm.step=0;
		int ret=get_body_num(tm);
		state[tm.hx][tm.hy][ret]=true;	
		qq.push(tm); 
		int ans=-1;
		while(!qq.empty())
		{
			node t=qq.front();
			qq.pop();
			if (t.hx==1&&t.hy==1)		//结束条件
			{
				ans=t.step;
				break;
			}
			deal_map(t);		//设蛇身为障碍
			for (i=0;i<4;i++)
			{
				int x=t.hx+dirx[i];
				int y=t.hy+diry[i];
				if (mp[x][y]) continue;			//如果不可行 
				node res=move_to(t,x,y);		//移动后得到的新状态res
				int ret=get_body_num(res);		//计算是否出现过当前状态
				if (state[res.hx][res.hy][ret]==true)		
					continue;
				state[res.hx][res.hy][ret]=true;
				qq.push(res);				

			}
			un_deal_map(t);		//恢复地图
		}
		
		printf("Case %d: %d\n",cnt++,ans);
	}
	return 0;
	
}