NOIP2013提高组D2T3 华容道

n<=30 * m<=30 的地图上，0表示墙壁，1表示可以放箱子的空地。q<=500次询问，每次问：当空地上唯一没有放箱子的空格子在(ex,ey)时，把位于(sx,sy)的箱子移动到(tx,ty)的最小步数。

第一档：n<=10,m<=10,不加剪枝地乱搞？？

第二档：n<=30,m<=30,q<=10，直接搜索，我觉得这60分够多了。。

所有档：

完成移动需要两个过程：空白格移动到指定箱子的周围（不能经过指定箱子），指定箱子在空白格的配合下逐渐逼近目标格子。

前者可以直接搜索得到，而后者比较麻烦。为了配合指定箱子，空白格必须时刻在指定箱子周围运动，也就是说，他会从指定箱子的上/下/左/右用最小代价移动到箱子的其他方位，然后让箱子朝那个方向前进一格。

这样的话，实际上是三元状态(x,y,z)，(x,y)是当前指定箱子的坐标，z是空格子相对他的方位，这样一些状态之间走来走去，最终走道(tx,ty,上下左右)，问最小代价。

可以发现状态之间走来走去的规则非常简单：一个是(x,y,上)转移成(x-1,y,下)这样的不同格子的相反空格方位的移动，另一种是(x,y,上)转移到(x,y,左)这样的同个格子、空格不同方位的转移。

同个格子不同方位间转移的代价可以bfs预处理，而，最后查最小代价只需要一次最短路。

 #include<stdio.h>
 #include<string.h>
 #include<math.h>
 #include<stdlib.h>
 #include<algorithm>
 #include<queue>
 //#include<iostream>
 using namespace std;
 
 //预处理1：play[i][j][0-3][0-3]表示格子(i,j)的上0下1左2右3到上下左右的不经过格(i,j)的最小移动步数
 //最短路
 
 bool isdigit(char c) {return c>='' && c<='';}
 int qread()
 {
     char c;int s=;while (!isdigit(c=getchar()));
     do s=s*+c-''; while (isdigit(c=getchar()));return s;
 }
 
 const int inf=0x3f3f3f3f;
 int n,m,Q;
 int mp[][],play[][][][];
 const int dx[]={-,,,},dy[]={,,-,};
 struct qnode1
 {
     int x,y,d;
 }que1[*];int head,tail;bool vis1[][];
 void pre(int x,int y)
 {
     mp[x][y]=;
     for (int i=;i<;i++)
         for (int j=;j<;j++)
         {
             if (i==j) {play[x][y][i][j]=;continue;}
             if (x+dx[i]< || x+dx[i]>n || y+dy[j]< || y+dy[j]>m ||
             !mp[x+dx[i]][y+dy[i]] || !mp[x+dx[j]][y+dy[j]])
             {play[x][y][i][j]=inf;continue;}
             play[x][y][i][j]=inf;
             head=tail=;
             que1[tail].d=,que1[tail].x=x+dx[i],que1[tail++].y=y+dy[i];
             memset(vis1,,sizeof(vis1));vis1[x+dx[i]][y+dy[i]]=;
             while (head!=tail)
             {
                 const int nx=que1[head].x,ny=que1[head].y,nd=que1[head].d;head++;
                 if (nx==x+dx[j] && ny==y+dy[j])
                 {
                     play[x][y][i][j]=nd;
                     break;
                 }
                 for (int k=;k<;k++)
                 {
                     int xx=nx+dx[k],yy=ny+dy[k];
                     if (xx< || xx>n || yy< || yy>m || !mp[xx][yy] || vis1[xx][yy]) continue;
                     vis1[xx][yy]=;
                     que1[tail].x=xx,que1[tail].y=yy;que1[tail++].d=nd+;
                 }
             }
         }
     mp[x][y]=;
 }
 int road(int sx,int sy,int tx,int ty)
 {
     if (tx==sx && ty==sy) return ;
     head=tail=;
     que1[tail].d=,que1[tail].x=sx,que1[tail++].y=sy;
     memset(vis1,,sizeof(vis1));vis1[sx][sy]=;
     while (head!=tail)
     {
         const int nx=que1[head].x,ny=que1[head].y,nd=que1[head].d;head++;
         if (nx==tx && ny==ty) return nd;
         for (int i=;i<;i++)
         {
             int xx=nx+dx[i],yy=ny+dy[i];
             if (xx< || xx>n || yy< || yy>m || !mp[xx][yy] || vis1[xx][yy]) continue;
             vis1[xx][yy]=;
             que1[tail].d=nd+;que1[tail].x=xx;que1[tail++].y=yy;
         }
     }
     return inf;
 }
 int dis[][][];bool vis[][][];
 struct qnode
 {
     int x,y,z,d;
     bool operator > (const qnode &b) const {return d>b.d;}
 };
 priority_queue<qnode,vector<qnode>,greater<qnode> > q;
 void work()
 {
     int ex=qread(),ey=qread(),sx=qread(),sy=qread(),tx=qread(),ty=qread();
     if (sx==tx && sy==ty)
     {
         puts("");
         return;
     }
     for (int i=;i<=n;i++)
         for (int j=;j<=m;j++)
             for (int k=;k<;k++)
                 dis[i][j][k]=inf;
     for (int i=;i<;i++)
     {
         int xx=sx+dx[i],yy=sy+dy[i];
         if (xx< || xx>n || yy< || yy>m || !mp[xx][yy]) continue;
         mp[sx][sy]=;
         dis[sx][sy][i]=road(ex,ey,sx+dx[i],sy+dy[i]);
         mp[sx][sy]=;
         if (dis[sx][sy][i]<inf) q.push((qnode){sx,sy,i,dis[sx][sy][i]});
     }
     memset(vis,,sizeof(vis));
     while (!q.empty())
     {
         const int x=q.top().x,y=q.top().y,z=q.top().z,d=q.top().d;q.pop();
         if (vis[x][y][z]) continue;
         vis[x][y][z]=;
         if (z== && x> && mp[x-][y])
             if (dis[x-][y][]>d+)
             {
                 dis[x-][y][]=d+;
                 q.push((qnode){x-,y,,d+});
             }
         if (z== && x<n && mp[x+][y])
             if (dis[x+][y][]>d+)
             {
                 dis[x+][y][]=d+;
                 q.push((qnode){x+,y,,d+});
             }
         if (z== && y> && mp[x][y-])
             if (dis[x][y-][]>d+)
             {
                 dis[x][y-][]=d+;
                 q.push((qnode){x,y-,,d+});
             }
         if (z== && y<m && mp[x][y+])
             if (dis[x][y+][]>d+)
             {
                 dis[x][y+][]=d+;
                 q.push((qnode){x,y+,,d+});
             }
         for (int i=;i<;i++) if (i!=z)
             if (dis[x][y][i]>d+play[x][y][z][i])
             {
                 dis[x][y][i]=d+play[x][y][z][i];
                 q.push((qnode){x,y,i,d+play[x][y][z][i]});
             }
     }
     int ans=min(dis[tx][ty][],min(dis[tx][ty][],min(dis[tx][ty][],dis[tx][ty][])));
     printf("%d\n",(ans)>=inf?-:ans);
 }
 int main()
 {
     n=qread(),m=qread(),Q=qread();
     for (int i=;i<=n;i++)
         for (int j=;j<=m;j++)
             mp[i][j]=qread();
     for (int i=;i<=n;i++)
         for (int j=;j<=m;j++)
             if (mp[i][j])
                 pre(i,j);
     while (Q--) work();
     return ;
 }

然而注意判断答案为0.。。。。。。没判被卡成65分 不如暴力