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>='0' && c<='9';}
int qread()
{
    char c;int s=0;while (!isdigit(c=getchar()));
    do s=s*10+c-'0'; while (isdigit(c=getchar()));return s;
}

const int inf=0x3f3f3f3f;
int n,m,Q;
int mp[33][33],play[33][33][4][4];
const int dx[]={-1,1,0,0},dy[]={0,0,-1,1};
struct qnode1
{
    int x,y,d;
}que1[33*33];int head,tail;bool vis1[33][33];
void pre(int x,int y)
{
    mp[x][y]=0;
    for (int i=0;i<4;i++)
        for (int j=0;j<4;j++)
        {
            if (i==j) {play[x][y][i][j]=0;continue;}
            if (x+dx[i]<1 || x+dx[i]>n || y+dy[j]<1 || 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=0;
            que1[tail].d=0,que1[tail].x=x+dx[i],que1[tail++].y=y+dy[i];
            memset(vis1,0,sizeof(vis1));vis1[x+dx[i]][y+dy[i]]=1;
            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=0;k<4;k++)
                {
                    int xx=nx+dx[k],yy=ny+dy[k];
                    if (xx<1 || xx>n || yy<1 || yy>m || !mp[xx][yy] || vis1[xx][yy]) continue;
                    vis1[xx][yy]=1;
                    que1[tail].x=xx,que1[tail].y=yy;que1[tail++].d=nd+1;
                }
            }
        }
    mp[x][y]=1;
}
int road(int sx,int sy,int tx,int ty)
{
    if (tx==sx && ty==sy) return 0;
    head=tail=0;
    que1[tail].d=0,que1[tail].x=sx,que1[tail++].y=sy;
    memset(vis1,0,sizeof(vis1));vis1[sx][sy]=1;
    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=0;i<4;i++)
        {
            int xx=nx+dx[i],yy=ny+dy[i];
            if (xx<1 || xx>n || yy<1 || yy>m || !mp[xx][yy] || vis1[xx][yy]) continue;
            vis1[xx][yy]=1;
            que1[tail].d=nd+1;que1[tail].x=xx;que1[tail++].y=yy;
        }
    }
    return inf;
}
int dis[33][33][4];bool vis[33][33][4];
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("0");
        return;
    }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            for (int k=0;k<4;k++)
                dis[i][j][k]=inf;
    for (int i=0;i<4;i++)
    {
        int xx=sx+dx[i],yy=sy+dy[i];
        if (xx<1 || xx>n || yy<1 || yy>m || !mp[xx][yy]) continue;
        mp[sx][sy]=0;
        dis[sx][sy][i]=road(ex,ey,sx+dx[i],sy+dy[i]);
        mp[sx][sy]=1;
        if (dis[sx][sy][i]<inf) q.push((qnode){sx,sy,i,dis[sx][sy][i]});
    }
    memset(vis,0,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]=1;
        if (z==0 && x>1 && mp[x-1][y])
            if (dis[x-1][y][1]>d+1)
            {
                dis[x-1][y][1]=d+1;
                q.push((qnode){x-1,y,1,d+1});
            }
        if (z==1 && x<n && mp[x+1][y])
            if (dis[x+1][y][0]>d+1)
            {
                dis[x+1][y][0]=d+1;
                q.push((qnode){x+1,y,0,d+1});
            }
        if (z==2 && y>1 && mp[x][y-1])
            if (dis[x][y-1][3]>d+1)
            {
                dis[x][y-1][3]=d+1;
                q.push((qnode){x,y-1,3,d+1});
            }
        if (z==3 && y<m && mp[x][y+1])
            if (dis[x][y+1][2]>d+1)
            {
                dis[x][y+1][2]=d+1;
                q.push((qnode){x,y+1,2,d+1});
            }
        for (int i=0;i<4;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][0],min(dis[tx][ty][1],min(dis[tx][ty][2],dis[tx][ty][3])));
    printf("%d
",(ans)>=inf?-1:ans);
}
int main()
{
    n=qread(),m=qread(),Q=qread();
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            mp[i][j]=qread();
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            if (mp[i][j])
                pre(i,j);
    while (Q--) work();
    return 0;
}

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