COGS1442 [NOIP2013] 华容道

传送门

这道题我做了三天，非常chishi（其实还是自己太弱）
第一眼看到这题：hash！（于是写了双hash只得了50，改成单hash得40。。。）
too naive
noip day2T3怎么会有这么水
正解：强行建立与白块位置有关的奇怪的图，预处理后对于每个询问进行SPFA
考场上肯定是写不出来的，只能拿60分部分分（我也不知道为什么我的暴力常数如此之大导致我少得了10分）
分析题目：白块肯定是先要移动到目标棋子旁边，然后在目标棋子周围转悠，并且带着目标棋子一步步来到目标位置。
我们可以看到，对于每一个测试点，其中的询问都在同一个图上，所以我们可以bfs预处理出对于每一个可移动的点，白块在它的上下左右（如果可以的话）之间移动所需要的最少步数。
对于每个询问，我们可以将它拆成两个问题：将白块移动到目标棋子周围，和从目标棋子到目标位置。将白块移动到目标棋子周围，我们可以再进行一次bfs，计算图中每一个点到目标棋子四周的最少移动次数。其中白块所有的移动不能经过目标棋子所在位置。将目标棋子移动到目标位置，我们可以做一遍SPFA，利用我们预处理出的数据连边即可。
为了便于理解，发一下自己调试的时候画的图

代码非常难看，大改了好几次，，

CODE：

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int N=1e4+10;
const int INF=1e9;
const int step[5][2]={0,0,1,0,0,1,0,-1,-1,0};
struct edge
{
    int nxt,to,dis;
}a[N<<1];
int s[35][35];
int head[N],dis[N],d[35][35][5][5],dist[35][35];
bool b[N],vis[35][35];
int n,m,T,e1,e2,s1,s2,t1,t2,ans,num;
queue<int>q,q1,q2;
inline void add(int x,int y,int dis)
{
    a[++num].nxt=head[x],a[num].to=y,a[num].dis=dis,head[x]=num;
    a[++num].nxt=head[y],a[num].to=x,a[num].dis=dis,head[y]=num;
}
inline int pos(int i,int j,int k){return ((i-1)*m+j-1)*4+k;}
inline bool check(int x,int y){return x>=1&&x<=n&&y>=1&&y<=m&&s[x][y];}
inline int calc(int x,int y)
{
    int ans=0;
    for(int i=1;i<=4;i++)
      if(check(x+step[i][0],y+step[i][1])) ans++;
    return ans;
}
inline int judge(int x,int y,int X,int Y)
{
    for(int i=1;i<=4;i++)
      if(x+step[i][0]==X&&y+step[i][1]==Y) return i;
    return 0;
}
inline void bfs()
{
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++)
      {
        if(!s[i][j]) continue;
        int tot=calc(i,j);
        for(int t=1;t<=4;t++)
        {
            int x=i+step[t][0],y=j+step[t][1];
            if(!check(x,y)) continue;
            memset(vis,0,sizeof(vis));
            memset(dist,0x3f,sizeof(dist));
            while(!q1.empty()) q1.pop();
            while(!q2.empty()) q2.pop();
            q1.push(x),q2.push(y);
            vis[x][y]=1;dist[x][y]=0;
            vis[i][j]=1;
            int num=0;
            while(!q1.empty())
            {
                int x1=q1.front(),y1=q2.front();
                q1.pop(),q2.pop();
                for(int k=1;k<=4;k++)
                {
                    int X=x1+step[k][0],Y=y1+step[k][1],tmp;
                    if(!check(X,Y)||vis[X][Y]) continue;
                    vis[X][Y]=1;dist[X][Y]=dist[x1][y1]+1;
                    q1.push(X),q2.push(Y);
                    if(tmp=(judge(i,j,X,Y)))
                    {
                        num++;
                        d[i][j][t][tmp]=dist[X][Y];
                        if(num==tot-1) break;
                    }
                }
                if(num==tot-1) break;
            }
        }
      }
}
inline void solve()
{
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++)
        if(s[i][j]) for(int k=1;k<=4;k++)
        {
            int x=i+step[k][0],y=j+step[k][1];
            if(!check(x,y)) continue;
            for(int l=k+1;l<=4;l++)
              if(d[i][j][k][l]<INF)
                add(pos(i,j,k),pos(i,j,l),d[i][j][k][l]);
            add(pos(i,j,k),pos(x,y,5-k),1);
        }
}
inline void bfs2(int x,int y)
{
    memset(dist,0x3f,sizeof(dist));
    memset(vis,0,sizeof(vis));
    while(!q1.empty()) q1.pop();
    while(!q2.empty()) q2.pop();
    q1.push(x),q2.push(y);
    dist[x][y]=0;vis[x][y]=1;vis[s1][s2]=1;
    while(!q1.empty())
    {
        int x1=q1.front(),y1=q2.front();
        q1.pop(),q2.pop();
        for(int i=1;i<=4;i++)
        {
            int X=x1+step[i][0],Y=y1+step[i][1];
            if(!check(X,Y)||vis[X][Y]) continue;
            q1.push(X),q2.push(Y);
            dist[X][Y]=dist[x1][y1]+1;vis[X][Y]=1;
        }
    }
}
inline int SPFA()
{
    if(s1==t1&&s2==t2) return 0;
    int ans=INF;
    memset(dis,0x3f,sizeof(dis));
    for(int i=1;i<=4;i++)
    {
        int x=s1+step[i][0],y=s2+step[i][1];
        if(check(x,y)) q.push(pos(s1,s2,i)),dis[pos(s1,s2,i)]=dist[x][y];
    }
    while(!q.empty())
    {
        int tmp=q.front();q.pop();
        b[tmp]=0;
        for(int i=head[tmp];i;i=a[i].nxt)
          if(dis[a[i].to]>dis[tmp]+a[i].dis)
          {
            dis[a[i].to]=dis[tmp]+a[i].dis;
            if(!b[a[i].to]) b[a[i].to]=1,q.push(a[i].to);
          }
    }
    for(int i=1;i<=4;i++)
    {
        int x=t1+step[i][0],y=t2+step[i][1];
        if(check(x,y)) ans=min(ans,dis[pos(t1,t2,i)]);
    }
    return ans<INF?ans:-1;
}
int main()
{
    scanf("%d%d%d",&n,&m,&T);
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++)
        scanf("%d",&s[i][j]);
    memset(d,0x3f,sizeof(d));
    bfs();solve();
    while(T--)
    {
        scanf("%d%d%d%d%d%d",&e1,&e2,&s1,&s2,&t1,&t2);
        bfs2(e1,e2);
        printf("%d\n",SPFA());
    }
    return 0;
}