[NOIP2013]华容道

[NOIP2013]华容道

题目描述

小 B 最近迷上了华容道，可是他总是要花很长的时间才能完成一次。于是，他想到用编程来完成华容道：给定一种局面， 华容道是否根本就无法完成，如果能完成， 最少需要多少时间。

小 B 玩的华容道与经典的华容道游戏略有不同，游戏规则是这样的：

1. 在一个 n*m 棋盘上有 n*m 个格子，其中有且只有一个格子是空白的，其余 n*m-1个格子上每个格子上有一个棋子，每个棋子的大小都是 1*1 的；
2. 有些棋子是固定的，有些棋子则是可以移动的；
3. 任何与空白的格子相邻（有公共的边）的格子上的棋子都可以移动到空白格子上。

游戏的目的是把某个指定位置可以活动的棋子移动到目标位置。

给定一个棋盘，游戏可以玩 q 次，当然，每次棋盘上固定的格子是不会变的， 但是棋盘上空白的格子的初始位置、 指定的可移动的棋子的初始位置和目标位置却可能不同。第 i 次

玩的时候， 空白的格子在第 EXi 行第 EYi 列，指定的可移动棋子的初始位置为第 SXi 行第 SYi列，目标位置为第 TXi 行第 TYi 列。

假设小 B 每秒钟能进行一次移动棋子的操作，而其他操作的时间都可以忽略不计。请你告诉小 B 每一次游戏所需要的最少时间，或者告诉他不可能完成游戏。

输入输出格式

输入格式：

输入文件为 puzzle.in。

第一行有 3 个整数，每两个整数之间用一个空格隔开，依次表示 n、m 和 q；

接下来的 n 行描述一个 n*m 的棋盘，每行有 m 个整数，每两个整数之间用一个空格隔开，每个整数描述棋盘上一个格子的状态，0 表示该格子上的棋子是固定的，1 表示该格子上的棋子可以移动或者该格子是空白的。接下来的 q 行，每行包含 6 个整数依次是 EXi、EYi、SXi、SYi、TXi、TYi，每两个整数之间用一个空格隔开，表示每次游戏空白格子的位置，指定棋子的初始位置和目标位置。

输出格式：

输出文件名为 puzzle.out。

输出有 q 行，每行包含 1 个整数，表示每次游戏所需要的最少时间，如果某次游戏无法完成目标则输出−1。

输入输出样例

输入样例#1：

3 4 20 1 1 10 1 1 00 1 0 03 2 1 2 2 21 2 2 2 3 2

输出样例#1：

2-1

说明

【输入输出样例说明】

棋盘上划叉的格子是固定的，红色格子是目标位置，圆圈表示棋子，其中绿色圆圈表示目标棋子。

1. 第一次游戏，空白格子的初始位置是 (3, 2)（图中空白所示），游戏的目标是将初始位置在(1, 2)上的棋子（图中绿色圆圈所代表的棋子）移动到目标位置(2, 2)（图中红色的格子）上。

移动过程如下：

1. 第二次游戏，空白格子的初始位置是（1, 2）（图中空白所示），游戏的目标是将初始位置在（2, 2）上的棋子（图中绿色圆圈所示）移动到目标位置 (3, 2)上。

要将指定块移入目标位置，必须先将空白块移入目标位置，空白块要移动到目标位置，必然是从位置（2， 2）上与当前图中目标位置上的棋子交换位置，之后能与空白块交换位置的只有当前图中目标位置上的那个棋子，因此目标棋子永远无法走到它的目标位置， 游戏无

法完成。

【数据范围】

对于 30%的数据，1 ≤ n, m ≤ 10，q = 1；

对于 60%的数据，1 ≤ n, m ≤ 30，q ≤ 10；

对于 100%的数据，1 ≤ n, m ≤ 30，q ≤ 500。

题解：

这题难度给满分，终于不是一眼就出的模板题了，主要考察的是状态的表示方法和简化。

首先一看到这题肯定会想到bfs，尝试简化状态，可以发现由于每个状态可以划分为4个值来表示（白格子横纵作标，指定棋子的横纵坐标）可以证出时间复杂度为q（nm）^2，要超时的。

继续简化状态，发现只有空白格子与指定棋子相邻是棋子才可以移动，所以其他状态下指定棋子都无法移动，所以优化状态f[i][j][k]表示指定棋子坐标在i,j时白块在棋子的方向（上下左右）。

这样一来就无法用bfs了，因为4个状态无法一步就互相到达，这样我们可以用bfs先预处理出两个状态之间的距离，在用spfa跑最短路，时间复杂度4qnm。

代码：

#include<bits/stdc++.h>
#include<queue>
using namespace std;

const int max_n = 35;
const int inf = 1e6+7;

struct node
{
	int x,y;
	int t;
};

int b[4][4]={{1,0},{0,1},{-1,0},{0,-1}};
int point[max_n*max_n*4],nxt[max_n*max_n*4],v[max_n*max_n*4],val[max_n*max_n*4];
int dis[max_n][max_n];
int a[max_n][max_n],num[max_n][max_n][4],dist[max_n*max_n*4];
bool vis[max_n*max_n*4];
int n,m,q,x2,y2,x3,y3,x4,y4,tot;

inline int read() 
{ 
	int ans=0;char c;bool flag=false;
	while( (c=getchar())==' ' || c=='\n' || c=='\r');
	if(c=='-') flag=true;else ans=c-'0';
	while((c=getchar()) >='0' && c<='9') ans=ans*10+c-'0';
	return ans*(flag?-1:1);
} 

inline void addedge(int x,int y,int z)
{
	++tot; nxt[tot]=point[x]; 
	point[x]=tot; v[tot]=y;
	val[tot]=z;
}

inline int get_u(int k)
{
	if(k==0) return 2;
	if(k==1) return 3;
	if(k==2) return 0;
	if(k==3) return 1;
}

inline int bfs(int sx,int sy,int ex,int ey,int k)
{
	memset(dis,-1,sizeof(dis));
	queue<node> q;
	node tmp;
	tmp.x=sx; tmp.y=sy; tmp.t=0;
	q.push(tmp);
	while(!q.empty())
	{
		node now=q.front(); q.pop();
		for(int i=0; i<4; ++i)
		{
			int xx=now.x+b[i][0],yy=now.y+b[i][1];
			if(!a[xx][yy] || dis[xx][yy]!=-1) continue;
			dis[xx][yy]=now.t+1;
			if(xx==ex && yy==ey) continue; 
			node tmp;
			tmp.x=xx; tmp.y=yy; tmp.t=dis[xx][yy];
			q.push(tmp);
		}
	}
	if(k==4) return dis[ex][ey];
	for(int i=0; i<4; ++i)
	{
		int xx=ex+b[i][0],yy=ey+b[i][1];
		if(xx==sx && yy==sy) continue;
		if(dis[xx][yy]!=-1) 
		  addedge(num[ex][ey][k],num[ex][ey][i],dis[xx][yy]);
	}
	addedge(num[ex][ey][k],num[sx][sy][get_u(k)],1);
}

inline int spfa(int x2,int y2,int x3,int y3,int x4,int y4)
{
	queue<int> q;
	memset(dist,0x7,sizeof(dist));
	memset(vis,0,sizeof(vis));
	for(int i=0; i<4; ++i)
	{
		int xx=x3+b[i][0],yy=y3+b[i][1];
		a[x3][y3]=0;
		int tmp=bfs(x2,y2,xx,yy,4);
		if(xx==x2 && yy==y2) tmp=0;
		a[x3][y3]=1;
		if(tmp!=-1)
		{
			vis[num[x3][y3][i]]=true;
			q.push(num[x3][y3][i]);
			dist[num[x3][y3][i]]=tmp;
		}
	}
	while(!q.empty())
	{
		int now=q.front(); q.pop();
		vis[now]=false;
		for(int i=point[now]; i!=-1; i=nxt[i])
		  if(dist[v[i]]>dist[now]+val[i])
		  {
		  	dist[v[i]]=dist[now]+val[i];
		  	if(!vis[v[i]]) q.push(v[i]);
		  	vis[v[i]]=true;
		  }
	}
	int res=inf;
	for(int i=0; i<4; ++i)
	  res=min(res,dist[num[x4][y4][i]]);
	if(res>=inf) return -1;
	else return res;
}

inline void init()
{
	int tot=0;
	for(int i=1; i<=n; ++i)
	  for(int j=1; j<=m; ++j)
	    for(int k=0; k<4; ++k)
	      num[i][j][k]=++tot;
	memset(point,-1,sizeof(point));
	memset(nxt,-1,sizeof(nxt));
	tot=-1;
}

int main() 
{ 
	n=read(); m=read(); q=read();
	for(int i=1; i<=n; ++i)
	  for(int j=1; j<=m; ++j)
	    a[i][j]=read();
	init();
	for(int i=1; i<=n; ++i)
	  for(int j=1; j<=m; ++j)
	  {
	  	if(!a[i][j]) continue; 
	  	if(a[i+1][j]) bfs(i+1,j,i,j,0);
	  	if(a[i][j+1]) bfs(i,j+1,i,j,1);
	  	if(a[i-1][j]) bfs(i-1,j,i,j,2);
	  	if(a[i][j-1]) bfs(i,j-1,i,j,3);
	  }
	for(int i=1; i<=q; ++i)
	{
		x2=read(); y2=read(); x3=read();
		y3=read(); x4=read(); y4=read();
		int tmp1,tmp2;
		tmp1=spfa(x2,y2,x3,y3,x4,y4);
		if(x3==x4 && y3==y4) printf("0\n");
		else printf("%d\n",tmp1);
	}
}