RunningPhoton's Nightmare

1.orzzzzz，不知道BFS为啥不对！！！！非要改成SPFA。。。

2.经典套路，对时间装置重置点，起点和终点都跑一遍BFS，建一张图，然后在图上从起点到终点跑一遍SPFA，答案就出来了。

3.建图的时候两点之间的距离小于K才放进去。

4.QAQ，终于搞明白了。。BFS跑单位边权的图返回的才一定是最优解。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
typedef pair<int,int> P;

const int INF=100000000;

char map[605][605];
int N,M;
int sx,sy;
int gx,gy;

int T,k,tot,cnt;
int d[605][605],head[200],vis[200],dis[200],p[200][2];
int dx[4]={1,0,-1,0},dy[4]={0,1,0,-1};

struct edge{
    int to,next;
    int w;
}e[40000];

void e_init(){
    tot=0;
    memset(head,-1,sizeof(head));
}

void addedge(int u,int v,int w){
    e[tot].to=v;
    e[tot].w=w;
    e[tot].next=head[u];
    head[u]=tot++;
}

void init(int np,int m_sx,int m_sy){
    queue<P> q;
    for(int i=0;i<N;i++){
        for(int j=0;j<M;j++)
            d[i][j]=INF;
    }
    q.push(P(m_sx,m_sy));
    d[m_sx][m_sy]=0;
    
    while(q.size()){
        P p=q.front();
        q.pop();
        for(int i=0;i<4;i++){
            int nx=p.first+dx[i],ny=p.second+dy[i];
            if(0<=nx&&nx<N&&0<=ny&&ny<M&&map[nx][ny]!='W'&&d[nx][ny]==INF){
                 q.push(P(nx,ny));
                 d[nx][ny]=d[p.first][p.second]+1;
            }
        }
    }
    for(int i=0;i<=cnt;i++){
        int mx=p[i][0],my=p[i][1];
        if(i==np||d[mx][my]>=k) continue;
        addedge(np,i,d[mx][my]);
    }
}
void SPFA(int st){
    for(int i=1;i<=cnt;i++){
        dis[i]=INF;
        vis[i]=0;
    }
    dis[st]=0;
    vis[st]=1;
    queue<int> q;
    q.push(st);
    while(q.size()){
        int v=q.front();
        q.pop();
        for(int i=head[v];i!=-1;i=e[i].next){
            int t=e[i].to;
            if(dis[t]>dis[v]+e[i].w){
                dis[t]=dis[v]+e[i].w;
                if(!vis[t]){
                    q.push(t);
                    vis[t]=1;
                }
            }
        }
        vis[v]=0;
    }    
}
/*
int BFS(int st,int en){
    queue<int> q;
    memset(dis,0,sizeof(dis));
    memset(vis,0,sizeof(vis));
    
    q.push(st);
    vis[st]=1;
    while(q.size()){
        int v=q.front();
        q.pop();
        for(int i=head[v];i!=-1;i=e[i].next){
            int t=e[i].to;
            if(!vis[t]){
                vis[t]=1;
                dis[t]=dis[v]+e[i].w;
                if(t==en) return dis[en];
                q.push(t);
            }
        }
    }
    return -1;
}
*/
int main()
{   scanf("%d",&T);
    while(T--){
        cnt=0;
        e_init();
        scanf("%d%d%d",&N,&M,&k);
        memset(p,0,sizeof(p));
        for(int i=0;i<N;i++){
            for(int j=0;j<M;j++){
                cin>>map[i][j];
                if(map[i][j]=='S') sx=i,sy=j;
                if(map[i][j]=='E') gx=i,gy=j;
                if(map[i][j]=='R'){
                    cnt++;
                    p[cnt][0]=i;
                    p[cnt][1]=j;
                }
            }
            getchar(); 
        }
        cnt++;
        p[0][0]=sx,p[0][1]=sy,p[cnt][0]=gx,p[cnt][1]=gy;
        for(int i=0;i<=cnt;i++) 
            init(i,p[i][0],p[i][1]);
        SPFA(0);   
        if(dis[cnt]==INF) printf("Poor RunningPhoton!
");
        else printf("%d
",dis[cnt]);
    }
    return 0;
}