NOIP[2013] 华容道

      一道能够让人深刻理解“状态“含义的图论题目。首先一个很简单的做法就是暴力BFS，以空格子为起点进行BFS，需要记录的有2个量：空格子的位置和目标棋子的位置，我们考虑这样的复杂度，由于是BFS，每个节点只会被访问一次，不同的节点数有n*m*n*m个，再算上q组询问，复杂度为q*n*m*n*m。可以拿到60分。

      其实我们有很多状态都是没有必要的，题目中的q组询问中的地图是不会改变的，于是我们考虑先预处理再搞，我们可以知道只有当一个棋子四周存在空格子时，这个棋子才能移动，那么我们可以考虑一个状态F[i][j][k][h]，表示从（i，j）这个位置，空格子在它的k方向，它要向h方向移动所需要的最少步数，这个我们可以提前用BFS求出，我们可以把格子（i，j）的K方向上是空格子当成一个状态，记作（i,j,k），那么对于一个F[i][j][k][h]来说，它就是连接（i,j,k）和（i+dx[h],j+dy[h],h的相反方向）的边的边权，这样我们先提前用BFS求出F[i][j][k][h]，复杂度是O(n*m*4*4)，然后对于每组询问，我们求出那个空格子到达起点格子四周的步数，然后连边，然后再跑一下spfa，复杂度大概是O(q*n*m*4*k)。

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;#define maxn 1000005int pre[maxn],last[maxn],other[maxn],len[maxn],dist[maxn];int qx[maxn],qy[maxn],n,m,q,tot,cnt,que[maxn],l,S,T;int f[31][31][5][5],g[31][31],vis[31][31],dis[32][32],id[31][31][5];bool flag[maxn];const int dx[5]={0,1,0,0,-1};const int dy[5]={0,0,1,-1,0};bool pd(int x,int y) {if (x<1||x>n||y<1||y>m) return 1;if (g[x][y]==0) return 1;return 0;}void connect(int x,int y,int z){l++;pre[l]=last[x];last[x]=l;other[l]=y;len[l]=z;}int bfs(int sx,int sy,int tx,int ty) {if (sx==tx&&sy==ty) return 0;int h=1,t=1;qx[1]=sx;qy[1]=sy;dis[sx][sy]=0;vis[sx][sy]=tot;while (h<=t) {int x=qx[h],y=qy[h];h++;for (int i=1;i<=4;i++) {int xx=x+dx[i];int yy=y+dy[i];if (pd(xx,yy)) continue;if (vis[xx][yy]!=tot) {dis[xx][yy]=dis[x][y]+1;vis[xx][yy]=tot;t++;qx[t]=xx;qy[t]=yy;if (xx==tx&&yy==ty) return dis[xx][yy];}}}return 1e7;}void spfa(int s){memset(dist,53,sizeof dist);memset(flag,0,sizeof flag);dist[s]=0;que[1]=s;int h=1,t=1;while (h<=t) {int u=que[h];h++;flag[u]=0;for (int p=last[u];p;p=pre[p]) {int v=other[p];if (dist[v]>dist[u]+len[p]) {dist[v]=dist[u]+len[p];if (!flag[v]){que[++t]=v;flag[v]=1;}}}}}int main(){scanf("%d%d%d",&n,&m,&q);for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) scanf("%d",&g[i][j]);for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) for (int k=1;k<=4;k++) id[i][j][k]=++cnt;memset(f,53,sizeof f);for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) if (g[i][j]) {for (int k=1;k<=4;k++) {if (pd(i+dx[k],j+dy[k])) continue;for (int h=1;h<=4;h++) {if (pd(i+dx[h],j+dy[h])) continue;vis[i][j]=++tot;f[i][j][k][h]=bfs(i+dx[k],j+dy[k],i+dx[h],j+dy[h])+1;if (f[i][j][k][h]>1e7-1) continue;//printf("%d %d %d %d %d\n",i,j,k,h,f[i][j][k][h]);connect(id[i][j][k],id[i+dx[h]][j+dy[h]][5-h],f[i][j][k][h]);}}}while (q--) {int ex,ey,tx,ty,sx,sy;scanf("%d%d%d%d%d%d",&ex,&ey,&sx,&sy,&tx,&ty);if (sx==tx&&sy==ty) {printf("0\n");continue;} S=++cnt;T=++cnt;for (int i=1;i<=4;i++) {vis[sx][sy]=++tot;if (pd(sx+dx[i],sy+dy[i])) continue;int temp=bfs(ex,ey,sx+dx[i],sy+dy[i]);if (temp>1e7-1) continue;connect(S,id[sx+dx[i]][sy+dy[i]][5-i],temp);}for (int i=1;i<=4;i++) {if (pd(tx+dx[i],ty+dy[i])) continue;connect(id[tx][ty][i],T,0); }spfa(S);if (dist[T]>1e7-1) printf("-1\n");else printf("%d\n",dist[T]+1);}return 0;}