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

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

题意:

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

题解:

spfa。

代码:

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 #include <queue>
 5 #define inc(i,j,k) for(int i=j;i<=k;i++)
 6 #define maxn 510
 7 #define INF 0x3ffffff
 8 using namespace std;
 9 
10 inline int read(){
11     char ch=getchar(); int f=1,x=0;
12     while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
13     while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
14     return f*x;
15 }
16 struct e{int t,w,n;}; e es[maxn*4]; int g[maxn],ess;
17 void pe(int f,int t,int w){
18     es[++ess]=(e){t,w,g[f]}; g[f]=ess; es[++ess]=(e){f,w,g[t]}; g[t]=ess;
19 }
20 int d[maxn],f,p,c,m,cow[maxn],tot; bool inq[maxn];
21 queue <int> q;
22 void spfa(){
23     while(!q.empty())q.pop(); memset(inq,0,sizeof(inq)); inc(i,1,f)d[i]=INF;
24     q.push(1); inq[1]=1; d[1]=0;
25     while(!q.empty()){
26         int x=q.front(); q.pop(); inq[x]=0;
27         for(int i=g[x];i;i=es[i].n)if(d[x]+es[i].w<d[es[i].t]){
28             d[es[i].t]=d[x]+es[i].w;
29             if(!inq[es[i].t])inq[es[i].t]=1,q.push(es[i].t);
30         }
31     }
32 }
33 int main(){
34     f=read(); p=read(); c=read(); m=read(); inc(i,1,p){int x=read(),y=read(),z=read(); pe(x,y,z);}
35     inc(i,1,c)cow[i]=read(); spfa(); inc(i,1,c)if(d[cow[i]]<=m)tot++; printf("%d\n",tot);
36     inc(i,1,c)if(d[cow[i]]<=m)printf("%d\n",i); return 0;
37 }

 

20160803

posted @ 2016-08-06 14:11  YuanZiming  阅读(224)  评论(0编辑  收藏  举报