首页 > 代码库 > bzoj3575 [Hnoi2014]道路堵塞

bzoj3575 [Hnoi2014]道路堵塞

A国有N座城市,依次标为1到N。同时,在这N座城市间有M条单向道路,每条道路的长度是一个正整数。现在,A国交通部指定了一条从城市1到城市N的路径,并且保证这条路径的长度是所有从城市1到城市N的路径中最短的。不幸的是,因为从城市1到城市N旅行的人越来越多,这条由交通部指定的路径经常发生堵塞。现在A国想知道,这条路径中的任意一条道路无法通行时,由城市1到N的最短路径长度是多少。

Input

输入文件第一行是三个用空格分开的正整数N、M和L,分别表示城市数目、单向道路数目和交通部指定的最短路径包含多少条道路。
按下来M行,每行三个用空格分开的整数a、b和c,表示存在一条由城市a到城市b的长度为c的单向道路。这M行的行号也是对应道路的编号,即其中第1行对应的道路编号为1,第2行对应的道路编号为2,…,第M行对应的道路编号为M。最后一行为L个用空格分开的整数sp(1)…,,sp(L),依次表示从城市1到城市N的由交通部指定的最短路径上的道路的编号。

Output

输出文件包含L行,每行为一个整数,第i行(i=1,2…,,L)的整数表示删去编号为sp(i)的道路后从城市1到城市N的最短路径长度。如果去掉后没有从城市1到城市N的路径,则输出一1。

Sample Input

4 5 2
1 2 2
1 3 2
3 4 4
3 2 1
2 4 3
1 5

Sample Output

6
6

HINT

100%的数据满足2<N<100000,1<M<200000。所用道路长度大于0小于10000。

数据已加强By Vfleaking

 

 

这题我是看别人题解的,感觉复杂度很玄学。

传送门:

http://www.cnblogs.com/mmlz/p/4297696.html

http://www.cnblogs.com/MashiroSky/p/6440116.html

 1 program road(input,output);
 2 const
 3   inf=1000000000;
 4 type
 5   etype=record
 6      t,w,next:longint;
 7   end;
 8 var
 9   a,b,c,d,g,q,dis:array[0..100010]of longint;
10   h,p:array[0..1000000]of longint;
11   inq,ing:array[0..100010]of boolean;
12   e:array[0..200020]of etype;
13   n,m,l,i,j,k,head,tail,t,x,y,z,cnt,tot:longint;
14 procedure add(x,y,w:longint);
15 begin
16   e[i].t:=y;e[i].w:=w;e[i].next:=a[x];a[x]:=i;
17 end;
18 procedure swap(var a,b:longint);
19 begin
20    t:=a;a:=b;b:=t;
21 end;
22 procedure ins(x,y:longint);
23 begin
24    inc(tot);p[tot]:=x;h[tot]:=y;
25    k:=tot;
26    while (h[k]<h[k>>1]) and (k>1) do begin swap(h[k],h[k>>1]);swap(p[k],p[k>>1]);k:=k>>1; end;
27 end;
28 procedure del;
29 begin
30    p[1]:=p[tot];h[1]:=h[tot];dec(tot);
31    k:=1;
32    while true do
33       begin
34          if k+k>tot then break;
35          if k+k=tot then
36             begin
37                if h[k+k]<h[k] then begin swap(h[k],h[k+k]);swap(p[k],p[k+k]); end;
38                break;
39             end;
40          if (h[k]<=h[k+k]) and (h[k]<=h[k+k+1]) then break;
41          if h[k+k]<h[k+k+1] then begin swap(h[k],h[k+k]);swap(p[k],p[k+k]);k:=k+k; end
42          else begin swap(h[k],h[k+k+1]);swap(p[k],p[k+k+1]);k:=k+k+1; end;
43       end;
44 end;
45 begin
46    assign(input,road.in);assign(output,road.out);reset(input);rewrite(output);
47    readln(n,m,l);
48    fillchar(a,sizeof(a),0);
49    for i:=1 to m do begin readln(x,y,z);add(x,y,z); end;
50    fillchar(c,sizeof(c),0);c[1]:=1;
51    for i:=1 to l do begin read(b[i]);c[e[b[i]].t]:=i+1; end;
52    d[l+1]:=0;
53    for i:=l downto 1 do d[i]:=d[i+1]+e[b[i]].w;
54    dis[1]:=0;for i:=2 to n do dis[i]:=inf;
55    fillchar(ing,sizeof(ing),false);
56    fillchar(inq,sizeof(inq),false);
57    b[0]:=0;e[0].t:=1;
58    tot:=0;
59    for i:=1 to l do
60       begin
61          head:=0;tail:=1;cnt:=0;q[1]:=e[b[i-1]].t;inq[q[1]]:=true;
62          while head<>tail do
63             begin
64                inc(head);if head>100000 then head:=1;
65                j:=a[q[head]];
66                while j<>0 do
67                   begin
68                      if j<>b[i] then
69                         begin
70                            if dis[q[head]]+e[j].w<dis[e[j].t] then
71                               begin
72                                  dis[e[j].t]:=dis[q[head]]+e[j].w;
73                                  if c[e[j].t]>=c[e[b[i]].t] then
74                                     begin
75                                        if not ing[e[j].t] then begin inc(cnt);g[cnt]:=e[j].t;ing[e[j].t]:=true; end;
76                                     end
77                                  else if not inq[e[j].t] then
78                                     begin
79                                        inc(tail);if tail>100000 then tail:=1;
80                                        q[tail]:=e[j].t;inq[e[j].t]:=true;
81                                     end;
82                               end;
83                         end;
84                      j:=e[j].next;
85                   end;
86                inq[q[head]]:=false;
87             end;
88          for j:=1 to cnt do begin ins(c[g[j]],dis[g[j]]+d[c[g[j]]]);ing[g[j]]:=false; end;
89          while (p[1]<c[e[b[i]].t]) and (tot>0) do del;
90          if tot=0 then writeln(-1) else writeln(h[1]);
91          dis[e[b[i]].t]:=dis[e[b[i-1]].t]+e[b[i]].w;
92       end;
93    close(input);close(output);
94 end.

 

bzoj3575 [Hnoi2014]道路堵塞