[bzoj2574] [Poi1999]Store-Keeper

　　坑啊。。

　　膜了半天byvoid大爷的题解。https://www.byvoid.com/blog/poi-1999-mag/?replytocom=1335/

　　一开始从人的位置bfs一波，看看能走到初始包裹哪些方向上。要注意不能穿过初始包裹...byvoid的标程在处理穿过包裹什么的地方有问题...但数据略弱。

　　之后如题解所述。

　　发现自己不会求点双连通分量QAQ。。那部分就抄标程了...

　　大概就是记录一下当前走过的边，遇到割点就一直出栈...和求强连通分量差不多，就是一个点可能在多个点双里面..这个比较烦

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 const int maxn=;
 const int xx[]={-,,,},yy[]={,,-,};
 struct zs{int too,pre;}e[maxn<<];int tot,last[maxn];
 int too[],pre[],la[maxn],tt,st1[maxn],top,cnt,st2[maxn],uuu[][][][];bool uu[maxn];
 struct poi{int x,y;}S,T,M;
 struct zs1{int x,y,dir;}dl[maxn<<];
 int dfn[maxn],low[maxn],tim;
 bool gd[maxn],u[maxn];
 int id[][],dis[][][],q[maxn];
 char mp[][];
 int i,j,k,n,m,man,st,ed;

 inline void ins(int a,int b){
     int i=la[a];
     while(i&&too[i]!=b)i=pre[i];
     if(i)return;
     too[++tt]=b,pre[tt]=la[a],la[a]=tt;
 }
 inline bool same(int a,int b){
     int i;bool flag;
     for(i=la[a];i;i=pre[i])uu[too[i]]=;
     for(i=la[b];i&&!uu[too[i]];i=pre[i]);
     flag=i!=;
     for(i=la[a];i;i=pre[i])uu[too[i]]=;
     return flag;
 }
 inline bool canget(int x,int y,int dir1,int dir2){
     if(!gd[id[x][y]])return ;
     if(uuu[x][y][dir1][dir2]<)
         return (uuu[x][y][dir1][dir2]=same( id[x+xx[dir1]][y+yy[dir1]] , id[x+xx[dir2]][y+yy[dir2]] ));
     else return uuu[x][y][dir1][dir2];
 }
 void tarjan(int x,int fa){
     dfn[x]=low[x]=++tim;int too,u,v;
     for(int i=last[x];i;i=e[i].pre)if(dfn[e[i].too]<dfn[x]){
         too=e[i].too;
         if(!dfn[too]){
             st1[++top]=x,st2[top]=too;
             tarjan(too,x),low[x]=min(low[x],low[too]);
             if(low[too]>=dfn[x]){
                 gd[x]=;
                 cnt++;
                 do{
                     u=st1[top],v=st2[top],
                     ins(u,cnt),ins(v,cnt),top--;
                 }while((u!=x||v!=too)&&(u!=too||v!=x));
             }
         }
         else low[x]=min(low[x],dfn[e[i].too]);
     }
 }

 inline void insert(int a,int b){
     e[++tot].too=b,e[tot].pre=last[a],last[a]=tot,
     e[++tot].too=a,e[tot].pre=last[b],last[b]=tot;
 }

 inline bool check(int x,int y){return x>&&y>&&x<=n&&y<=m&&mp[x][y]!='S';}

 inline void bfs(int s){
     int l=,r=,i,now;q[]=s,u[s]=;
     while(l<r)
         for(i=last[now=q[++l]];i;i=e[i].pre)if(!u[e[i].too]&&e[i].too!=st)
             u[e[i].too]=,q[++r]=e[i].too;
 }
 inline int run(){
     int l=,r=,i,nx,ny,ndir,ndis,x,y;
     for(i=;i<;i++){
         x=S.x+xx[i],y=S.y+yy[i];
         if(check(x,y)&&u[id[x][y]])dl[++r]=(zs1){S.x,S.y,i},dis[S.x][S.y][i]=;//,printf("st:%d,%d\n",x,y);
     }
     while(l<r){
         l++,nx=dl[l].x,ny=dl[l].y,ndir=dl[l].dir,ndis=dis[nx][ny][ndir];
         x=nx+xx[ndir^],y=ny+yy[ndir^];//printf("case:%d,%d  dis:%d  dir:%d\n",nx,ny,ndis,ndir);
         if(!check(x,y))continue;
         if(x==T.x&&y==T.y)return ndis;
         for(i=;i<;i++)
             if(check(x+xx[i],y+yy[i])&&!dis[x][y][i]&&canget( x,y,ndir,i ))
                 dis[x][y][i]=ndis+,dl[++r]=(zs1){x,y,i};
     }
     return -;
 }
 int main(){
     memset(uuu,,sizeof(uuu));
     scanf("%d%d",&n,&m);int tmp=;tot=;
     for(i=;i<=n;i++){
         scanf("%s",mp[i]+);
         for(j=;j<=m;j++)if(mp[i][j]!='S'){
             id[i][j]=++tmp;
             for(k=;k<;k+=){
                 int x=i+xx[k],y=j+yy[k];
                 if(check(x,y))
                     insert(id[i][j],id[x][y]);//,printf("(%d,%d)-->(%d,%d)\n",x,y,i,j);
             }
             if(mp[i][j]=='M')man=tmp,M=(poi){i,j};
             if(mp[i][j]=='P')st=tmp,S=(poi){i,j};
             if(mp[i][j]=='K')ed=tmp,T=(poi){i,j};
         }
     }//return 233;
     for(i=;i<=tmp;i++)if(!dfn[i])tarjan(i,);
 //  for(i=1;i<=n;i++)for(j=1;j<=m;j++)if(id[i][j]){
 //      if(mp[i][j]!='S')printf("(%d,%d)  %d   fa:%d\n",i,j,id[i][j],getfa(id[i][j]));
 //      if(gd[id[i][j]])printf("        gd:  (%d,%d)\n",i,j);
 //  }
     bfs(man);
     int ans=run();
     if(ans<)puts("NO");else printf("%d\n",ans);
 }