bzoj3538[Usaco2014 Open]Dueling GPS*
bzoj3538[Usaco2014 Open]Dueling GPS
题意:
给你一个N个点的有向图,设定初始位置为1,结束位置为n。有两个GPS定位系统,分别认为经过边i的时间为Pi,和Qi.每走一条边的时候,如果一个系统认为走的这条边不是它认为的最短路,就会受到警告一次。如果走的这条边都不在两个系统认为的最短路范围内,就会受到2次警告。求最少需要受到多少次警告。n≤10000,边数≤50000
题解:
分别按两个GPS的边权求最短路,然后枚举每条边,把该边的边权变为其警告数,然后再求一次最短路。判断该边是不是最短路的条件是是否dis[es[i].f]+es[i].w==dis[es[i].t]。
代码:
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 10010 7 #define INF 0x3fffffff 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 f,t,w,n;}es[maxn*5]; int g[maxn],ess; 17 void pe(int f,int t,int w){es[++ess]=(e){f,t,w,g[f]}; g[f]=ess;} 18 int n,m,a[2][maxn*5],b[maxn*5],d[2][maxn]; queue<int>q; bool inq[maxn]; 19 void spfa(int s,int o){ 20 while(!q.empty())q.pop(); memset(inq,0,sizeof(inq)); inc(i,1,n)d[o][i]=INF; 21 q.push(s); inq[s]=1; d[o][s]=0; 22 while(!q.empty()){ 23 int x=q.front(); q.pop(); inq[x]=0; 24 for(int i=g[x];i;i=es[i].n)if(d[o][x]+es[i].w<d[o][es[i].t]){ 25 d[o][es[i].t]=d[o][x]+es[i].w; if(!inq[es[i].t])q.push(es[i].t),inq[es[i].t]=1; 26 } 27 } 28 } 29 int main(){ 30 n=read(); m=read(); inc(i,1,m){int x=read(),y=read(),z=read(); pe(y,x,z); a[0][i]=z; a[1][i]=read();} 31 spfa(n,0); inc(i,1,m)es[i].w=a[1][i]; spfa(n,1); 32 inc(i,1,m){inc(j,0,1)if(d[j][es[i].f]+a[j][i]>d[j][es[i].t])b[i]++; es[i].w=b[i];} 33 spfa(n,0); printf("%d",d[0][1]); return 0; 34 }
20160909