poj1475 Pushing Boxes[双重BFS(毒瘤搜索题)]
地址。
很重要的搜索题。★★★
吐槽：算是写过的一道码量比较大的搜索题了，细节多，还比较毒瘤。虽然是一遍AC的，其实我提前偷了大数据，但是思路还是想了好长时间，照理说想了半小时出不来，我就会翻题解，但是这次总觉得自己快想出来了结果磕了两个小时才想到正解上。嘛主要还是我太弱了QwQ。
题目内容很简单。给一张图20*20，推箱子，求推动次数最少的路径，有推动次数相同方案时输出走路次数最少的路径。
乍一看这题觉得很棘手。箱子在动，人还要来回跑来跑去，虽然可以感觉出来是搜索题，但怎么搜？DFS效率太低明显不行的，先排除。那就剩BFS。但是，两个对象，还要维护步数和推数两个玩意儿，还是很棘手。模拟了一下人玩推箱子的过程，发现有一个规律：除了第一次人是奔着箱子去的，其他时候总是人推了箱子后要么继续推要么到箱子另一个方向继续推的。所以时时刻刻人和箱子总是黏在一起的，或者走一点路到另一个方向再“黏”上去。这样就可以把人和箱子作为整体看了。然后就有了记录推数的数组。$pus[x][y][dir]$表示人在坐标$(x,y)$处准备推他上/下/左/右边$(dir=0～3)$的箱子前已推的次数,$walk[x][y][dir]$维护当下已行步数。队列维护一个推数单调增的队列以保证最先求得答案。每次取队头先推一下，记下当前已推次数，拓展的时候就枚举方向看可不可以从箱子新位置四周推动。可以就将这种状态压入队尾，准备下次开推。以上是粗略描述。下面记一下诸多细节。
- 队列存的是三元组x,y,dir. 每次用最新的推动情况来拓展后面的状况。有两种情形：1、继续向前推：没什么好说的；2、到箱子其他方向推：也需判断能不能推动，可不可以到，可以的话就把现在的步数和推数赋值给他，进队。
- 当考虑到箱子其他方向上推时，因为需要了解此处到彼处的最短路程，每次，箱子的位置不能通行，所以每个点在每次拓展四个方向前要先跑一次图。因为状态最多$20*20*4=1600$个（每个点每个朝向），所以总共跑图不会昊太多时间。
- 既然提到队列中状态数量不会超过1600，就要避免重复或无效推动。即我推过来又推回去，或者现在推的状态以前已经拓展过了等等。处理方法：拓展时看四个方向的推法有没有被算过，有的话是不是比当前推动的次数少，少的话就不用了。次数相同表明是我现在同一层拓展出来的产物，那就要看谁到这个状态的步行次数少了。如果没算过，直接赋值入队。又因为该状态是在下一层的，这一层所有状态都有可能去更新到他，所以该状态在取出前一定是最优解。这样同时保证了不会推来推去和冗余重复，后面的效仿前面的。
- 注意队列只是推数单调，步数并不单调。所以我找到目标后要等这一层全拓展完看哪个步数少。
- 有可能我推了之后的状态以前有过次数更少的，直接不算了（见line48,如转圈）
- 寻路过程还具体看code。其实框架还基本差不多。
- 关于输出路径：行走过程不太好记，所以索性不记了，只记箱子推动的路径（code里面的las），las[x][y][dir]表示将要向dir推的当下坐标上一次是在哪儿向哪边推的。这个想想就明白了。中间的行走路径就再跑一边BFS即可，所以是一个推动和走路穿插的串。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<cmath> 5 #include<algorithm> 6 #include<queue> 7 #define mst(a,x) memset(a,x,sizeof a) 8 #define dbg(x) cerr<<#x<<" = "<<x<<endl 9 #define ddbg(x,y) cerr<<#x<<" = "<<x<<" "<<#y<<" = "<<y<<endl 10 using namespace std; 11 typedef long long ll; 12 typedef pair<int,int> pii; 13 template<typename T>inline char MIN(T&A,T B){return A>B?A=B,1:0;} 14 template<typename T>inline char MAX(T&A,T B){return A<B?A=B,1:0;} 15 template<typename T>inline T _min(T A,T B){return A<B?A:B;} 16 template<typename T>inline T _max(T A,T B){return A>B?A:B;} 17 template<typename T>inline T read(T&x){ 18 x=0;int f=0;char c;while(!isdigit(c=getchar()))if(c=='-')f=1; 19 while(isdigit(c))x=x*10+(c&15),c=getchar();return f?x=-x:x; 20 } 21 const int N=20+2,INF=0x3f3f3f3f; 22 const int dx[4]={-1,1,0,0},dy[4]={0,0,-1,1}; 23 const char d[4]={'N','S','W','E'}; 24 int mp[N][N],dis[N][N],n,m,bx,by,sx,sy,tx,ty; 25 #define xx x+dx[i] 26 #define yy y+dy[i] 27 queue<pii> Q; 28 inline void find_way(int X,int Y,int bX,int bY){ 29 mp[bX][bY]=1;mst(dis,INF);dis[X][Y]=0,Q.push(make_pair(X,Y)); 30 while(!Q.empty()){ 31 int x=Q.front().first,y=Q.front().second;Q.pop(); 32 for(register int i=0;i<4;++i)if(!mp[xx][yy]&&dis[xx][yy]>dis[x][y]+1)dis[xx][yy]=dis[x][y]+1,Q.push(make_pair(xx,yy)); 33 } 34 mp[bX][bY]=0; 35 } 36 #undef xx 37 #undef yy 38 int pus[N][N][4],walk[N][N][4],vis[N][N][4],pre[N][N][4],flag,ansx,ansy,ansdir,answ; 39 pair<pii,int> las[N][N][4]; 40 inline void start_pushing(){ 41 find_way(sx,sy,bx,by);queue<pair<pii,int> > q; 42 mst(pus,INF),mst(walk,INF),mst(vis,0);answ=INF;ansx=ansy=ansdir=0;flag=-1; 43 for(register int i=0,xx=bx+dx[0],yy=by+dy[0];i<4;++i,xx=bx+dx[i],yy=by+dy[i])if(dis[xx][yy]<INF&&!mp[bx-dx[i]][by-dy[i]]) 44 q.push(make_pair(make_pair(xx,yy),i^1)),pus[xx][yy][i^1]=0,walk[xx][yy][i^1]=dis[xx][yy],las[xx][yy][i^1]=make_pair(make_pair(0,0),0); 45 while(!q.empty()){ 46 int x=q.front().first.first,y=q.front().first.second,dir=q.front().second;q.pop(); 47 int pushx=pus[x][y][dir]+1,walkx=walk[x][y][dir]; 48 if(~flag&&pushx>flag)return;if(pus[x+dx[dir]][y+dy[dir]][dir]<pushx)continue; 49 int kx=x+2*dx[dir],ky=y+2*dy[dir]; 50 if(kx==tx&&ky==ty){flag=pushx;if(MIN(answ,walkx))ansx=x,ansy=y,ansdir=dir;continue;} 51 find_way(x+dx[dir],y+dy[dir],kx,ky); 52 for(register int i=0,xx=kx+dx[0],yy=ky+dy[0];i<4;++i,xx=kx+dx[i],yy=ky+dy[i])if(dis[xx][yy]<INF&&!mp[kx-dx[i]][ky-dy[i]]) 53 if(pus[xx][yy][i^1]>=pushx&&walk[xx][yy][i^1]>=walkx+dis[xx][yy]){ 54 pus[xx][yy][i^1]=pushx,walk[xx][yy][i^1]=walkx+dis[xx][yy],las[xx][yy][i^1]=make_pair(make_pair(x,y),dir); 55 if(!vis[xx][yy][i^1])q.push(make_pair(make_pair(xx,yy),i^1)),vis[xx][yy][i^1]=1; 56 } 57 } 58 } 59 char ans[2400];int len,T,ent;pii Pre[N][N]; 60 inline void print_way(int fx,int fy,int lx,int ly,int bX,int bY){ 61 if(fx==lx&&fy==ly)return; 62 mp[bX][bY]=1;mst(dis,INF);dis[fx][fy]=0,Pre[fx][fy]=make_pair(0,0);queue<pii> q;q.push(make_pair(fx,fy)); 63 while(!q.empty()){ 64 int x=q.front().first,y=q.front().second;q.pop();if(x==lx&&y==ly)break; 65 for(register int i=0;i<4;++i)if(!mp[x+dx[i]][y+dy[i]]&&dis[x+dx[i]][y+dy[i]]>dis[x][y]+1) 66 dis[x+dx[i]][y+dy[i]]=dis[x][y]+1,q.push(make_pair(x+dx[i],y+dy[i])),Pre[x+dx[i]][y+dy[i]]=make_pair(x,y); 67 }mp[bX][bY]=0;int x=Pre[lx][ly].first,y=Pre[lx][ly].second; 68 while(x)ans[++len]=(x^lx?(x==lx-1?'s':'n'):(y==ly-1?'e':'w')),lx=x,ly=y,x=Pre[lx][ly].first,y=Pre[lx][ly].second; 69 } 70 inline void print(){ 71 len=0;while(ansx){ 72 ans[++len]=d[ansdir]; 73 int x=las[ansx][ansy][ansdir].first.first,y=las[ansx][ansy][ansdir].first.second,dir=las[ansx][ansy][ansdir].second; 74 if(x)print_way(x+dx[dir],y+dy[dir],ansx,ansy,x+(dx[dir]<<1),y+(dy[dir]<<1));else break; 75 ansx=x,ansy=y,ansdir=dir; 76 } 77 print_way(sx,sy,ansx,ansy,bx,by); 78 for(register int i=len;i;--i)printf("%c",ans[i]);puts("");puts(""); 79 } 80 81 int main(){//freopen("test.in","r",stdin);freopen("test.out","w",stdout); 82 while(read(n),read(m),n&&m){ 83 char s[24];//if(ent)puts("");else ent=1; 84 for(register int i=1;i<=n;++i){ 85 scanf("%s",s+1);mp[i][0]=mp[i][m+1]=1; 86 for(register int j=1;j<=m;++j)switch(s[j]){ 87 case '#':mp[i][j]=1;break; 88 case '.':mp[i][j]=0;break; 89 case 'S':mp[i][j]=0,sx=i,sy=j;break; 90 case 'T':mp[i][j]=0,tx=i,ty=j;break; 91 case 'B':mp[i][j]=0,bx=i,by=j;break; 92 } 93 } 94 for(register int j=0;j<=m+1;++j)mp[0][j]=mp[n+1][j]=1; 95 start_pushing();printf("Maze #%d\n",++T);if(~flag)print();else printf("Impossible.\n\n"); 96 } 97 return 0; 98 }
