bzoj1681[Usaco2005 Mar]Checking an Alibi 不在场的证明*

bzoj1681[Usaco2005 Mar]Checking an Alibi 不在场的证明

题意：

给个点集，求无向有权图中点集里的哪些点到点1的距离小于等于M。点集内点数≤100，图中点数≤500，边数≤1000。

题解：

spfa。

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#define inc(i,j,k) for(int i=j;i<=k;i++)
#define maxn 510
#define INF 0x3ffffff
using namespace std;

inline int read(){
    char ch=getchar(); int f=1,x=0;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return f*x;
}
struct e{int t,w,n;}; e es[maxn*4]; int g[maxn],ess;
void pe(int f,int t,int w){
    es[++ess]=(e){t,w,g[f]}; g[f]=ess; es[++ess]=(e){f,w,g[t]}; g[t]=ess;
}
int d[maxn],f,p,c,m,cow[maxn],tot; bool inq[maxn];
queue <int> q;
void spfa(){
    while(!q.empty())q.pop(); memset(inq,0,sizeof(inq)); inc(i,1,f)d[i]=INF;
    q.push(1); inq[1]=1; d[1]=0;
    while(!q.empty()){
        int x=q.front(); q.pop(); inq[x]=0;
        for(int i=g[x];i;i=es[i].n)if(d[x]+es[i].w<d[es[i].t]){
            d[es[i].t]=d[x]+es[i].w;
            if(!inq[es[i].t])inq[es[i].t]=1,q.push(es[i].t);
        }
    }
}
int main(){
    f=read(); p=read(); c=read(); m=read(); inc(i,1,p){int x=read(),y=read(),z=read(); pe(x,y,z);}
    inc(i,1,c)cow[i]=read(); spfa(); inc(i,1,c)if(d[cow[i]]<=m)tot++; printf("%d
",tot);
    inc(i,1,c)if(d[cow[i]]<=m)printf("%d
",i); return 0;
}

20160803