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