【NOIP2013】华容道
看别人的代码然后被坑了一下午+一晚上，睡一觉第二天醒悟过来打表过了
果然打表才是正确的调试方法，跟踪什么的去屎(╯‵□′)╯︵┻━┻
原题：
小 B 最近迷上了华容道，可是他总是要花很长的时间才能完成一次。于是，他想到用编程来完成华容道：给定一种局面，华容道是否根本就无法完成，如果能完成，最少需要多少时间。
小 B 玩的华容道与经典的华容道游戏略有不同，游戏规则是这样的：
1. 在一个 n*m 棋盘上有 n*m 个格子，其中有且只有一个格子是空白的，其余 n*m-1 个格子上每个格子上有一个棋子，每个棋子的大小都是 1*1 的；
2. 有些棋子是固定的，有些棋子则是可以移动的；
3. 任何与空白的格子相邻（有公共的边）的格子上的棋子都可以移动到空白格子上。
游戏的目的是把某个指定位置可以活动的棋子移动到目标位置。
给定一个棋盘，游戏可以玩 q 次，当然，每次棋盘上固定的格子是不会变的，但是棋盘上空白的格子的初始位置、指定的可移动的棋子的初始位置和目标位置却可能不同。第 i 次玩的时候，空白的格子在第 EX[i] 行第 EY[i] 列，指定的可移动棋子的初始位置为第 SX[i] 行第 SY[i] 列，目标位置为第 TX[i] 行第 TY[i] 列。
假设小 B 每秒钟能进行一次移动棋子的操作，而其他操作的时间都可以忽略不计。请
你告诉小 B 每一次游戏所需要的最少时间，或者告诉他不可能完成游戏。
1 ≤ n, m ≤ 30，q ≤ 500。
这题暴力分很实惠，很容易写出一个bfs能拿到60-70
想要A掉怎么办呐
正解就要魔性p♂l♀a♂y了
用一个状态b[i][j][k][h]表示第i行第j列的格子，空格在k方向，这个格子要移动到h方向需要的最小步数
这个步数怎么求呐
可以用一个bfs，求出空格从k方向到h方向需要的步数，然后需要移动的格子和h方向的空格交换即可，可以证明这样做是最优的（逃
然后一个三维的点c[i][j][k]表示第i行第j列的格子，空格在k方向，上面用bfs求出步数↑后就可以把c[i][j][k]与c[i+k方向上x的变化量][j+k方向上y的变化量][h的反方向]连上权值为步数的边
因为求出步数的方法是先把空格移到h方向，在把格子和空格交换↑，所以边的汇点的方向是h的反方向
然后每次询问的时候，可以直接套用上面的bfs↑求出把空白格移动到要移动的方格附近需要的步数，四个方向都搞一下，建一个源点s，把s与c[i][j][k]连上权值为步数的边
然后再建一个汇点t，把t和c[目标节点x][目标节点y][四个方向]连上权值为0的边
然后就可以开心地搞最短路辣
需要注意的的方：
最好判断一下要移动的格子和目标格子是一个东西以及要移动的格子或目标格子不合法的情况
连边的时候需要注意什么时候坐标要加上方向，什么时候坐标是原格子
小技巧：
可以用一个计数，在输入数据的时候给每个c[i][j][k]都标上号，连边的时候就很方便
把x和y的变化量分别设成{1,-1,0,0}，{0,0,1,-1}，这样直接用k^1就可以快速取到反向边
代码：
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #include<cstring> 5 #include<cmath> 6 using namespace std; 7 const int oo=168430090; 8 int fx[4]={1,-1,0,0},fy[4]={0,0,1,-1}; 9 struct ddd{int next,y,value;}e[510000];int LINK[110000],ltop=0; 10 inline void insert(int x,int y,int z){e[++ltop].next=LINK[x];LINK[x]=ltop;e[ltop].y=y;e[ltop].value=z;} 11 int m,n,q,a[50][50]; 12 int hash[50][50][4],hash_cnt=0; 13 int dui[510000],tou=0; 14 int bf[50][50]; 15 int s,t; 16 int sf[11000]; 17 bool visited[11000]; 18 int bfs(int sx,int sy,int tx,int ty){ 19 memset(bf,10,sizeof(bf)); 20 dui[tou=1]=(sx-1)*n+sy-1; bf[sx][sy]=0; 21 for(int k=1;k<=tou;k++){ 22 int x=dui[k]/n+1,y=dui[k]%n+1; 23 for(int i=0;i<4;i++)if(a[x+fx[i]][y+fy[i]] && bf[x+fx[i]][y+fy[i]]>bf[x][y]+1){ 24 bf[x+fx[i]][y+fy[i]]=bf[x][y]+1; 25 //if(x+fx[i]==tx && y+fy[i]==ty) return bf[tx][ty]; 26 dui[++tou]=(x+fx[i]-1)*n+y+fy[i]-1; 27 } 28 } 29 return bf[tx][ty]; 30 } 31 int spfa(){ 32 memset(sf,10,sizeof(sf)); 33 memset(visited,0,sizeof(visited)); 34 dui[tou=1]=s; sf[s]=0; 35 for(int k=1;k<=tou;k++){ 36 visited[dui[k]]=false; 37 for(int i=LINK[dui[k]];i;i=e[i].next) 38 if(sf[e[i].y]>sf[dui[k]]+e[i].value){ 39 sf[e[i].y]=sf[dui[k]]+e[i].value; 40 if(!visited[e[i].y]) dui[++tou]=e[i].y,visited[e[i].y]=true; 41 } 42 } 43 return (sf[t]<oo) ? sf[t] : -1; 44 } 45 int main(){//freopen("ddd.in","r",stdin); 46 cin>>m>>n>>q; 47 for(int i=1;i<=m;i++) for(int j=1;j<=n;j++){ 48 scanf("%d",&a[i][j]); 49 for(int k=0;k<4;k++) hash[i][j][k]=++hash_cnt; 50 } 51 for(int i=1;i<=m;i++) 52 for(int j=1;j<=n;j++)if(a[i][j]) 53 for(int k=0;k<4;k++)if(a[i+fx[k]][j+fy[k]]){ 54 insert(hash[i][j][k],hash[i+fx[k]][j+fy[k]][k^1],1); 55 //cout<<i<<" "<<j<<" "<<k<<" "<<i+fx[k]<<" "<<j+fy[k]<<" "<<(k^1)<<" "<<1<<endl; 56 for(int q=0;q<4;q++)if(q!=k && a[i+fx[q]][j+fy[q]]){ 57 a[i][j]=0; 58 int bu=bfs(i+fx[k],j+fy[k],i+fx[q],j+fy[q])+1; 59 if(bu<oo){ insert(hash[i][j][k],hash[i+fx[q]][j+fy[q]][q^1],bu); /*cout<<i<<" "<<j<<" "<<k<<" "<<i+fx[q]<<" "<<j+fy[q]<<" "<<(q^1)<<" "<<bu<<endl;*/} 60 a[i][j]=1; 61 } 62 } 63 int sx,sy,tx,ty,kx,ky; 64 while(q --> 0){//趋向于 65 scanf("%d%d%d%d%d%d",&kx,&ky,&sx,&sy,&tx,&ty); 66 if(sx==tx && sy==ty){ printf("0\n"); continue;} 67 if(!a[sx][sy] || !a[tx][ty]){ printf("-1\n"); continue;} 68 s=++hash_cnt;t=++hash_cnt; 69 a[sx][sy]=0; 70 for(int k=0;k<4;k++)if(a[sx+fx[k]][sy+fy[k]]){ 71 int bu=bfs(kx,ky,sx+fx[k],sy+fy[k]); 72 if(bu<oo) insert(s,hash[sx][sy][k],bu); 73 } 74 a[sx][sy]=1; 75 for(int k=0;k<4;k++)if(a[tx+fx[k]][ty+fy[k]]) 76 insert(hash[tx][ty][k],t,0); 77 printf("%d ",spfa()); 78 } 79 return 0; 80 }
