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

program road(input,output);
const
  inf=1000000000;
type
  etype=record
     t,w,next:longint;
  end;
var
  a,b,c,d,g,q,dis:array[0..100010]of longint;
  h,p:array[0..1000000]of longint;
  inq,ing:array[0..100010]of boolean;
  e:array[0..200020]of etype;
  n,m,l,i,j,k,head,tail,t,x,y,z,cnt,tot:longint;
procedure add(x,y,w:longint);
begin
  e[i].t:=y;e[i].w:=w;e[i].next:=a[x];a[x]:=i;
end;
procedure swap(var a,b:longint);
begin
   t:=a;a:=b;b:=t;
end;
procedure ins(x,y:longint);
begin
   inc(tot);p[tot]:=x;h[tot]:=y;
   k:=tot;
   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;
end;
procedure del;
begin
   p[1]:=p[tot];h[1]:=h[tot];dec(tot);
   k:=1;
   while true do
      begin
         if k+k>tot then break;
         if k+k=tot then
            begin
               if h[k+k]<h[k] then begin swap(h[k],h[k+k]);swap(p[k],p[k+k]); end;
               break;
            end;
         if (h[k]<=h[k+k]) and (h[k]<=h[k+k+1]) then break;
         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
         else begin swap(h[k],h[k+k+1]);swap(p[k],p[k+k+1]);k:=k+k+1; end;
      end;
end;
begin
   assign(input,'road.in');assign(output,'road.out');reset(input);rewrite(output);
   readln(n,m,l);
   fillchar(a,sizeof(a),0);
   for i:=1 to m do begin readln(x,y,z);add(x,y,z); end;
   fillchar(c,sizeof(c),0);c[1]:=1;
   for i:=1 to l do begin read(b[i]);c[e[b[i]].t]:=i+1; end;
   d[l+1]:=0;
   for i:=l downto 1 do d[i]:=d[i+1]+e[b[i]].w;
   dis[1]:=0;for i:=2 to n do dis[i]:=inf;
   fillchar(ing,sizeof(ing),false);
   fillchar(inq,sizeof(inq),false);
   b[0]:=0;e[0].t:=1;
   tot:=0;
   for i:=1 to l do
      begin
         head:=0;tail:=1;cnt:=0;q[1]:=e[b[i-1]].t;inq[q[1]]:=true;
         while head<>tail do
            begin
               inc(head);if head>100000 then head:=1;
               j:=a[q[head]];
               while j<>0 do
                  begin
                     if j<>b[i] then
                        begin
                           if dis[q[head]]+e[j].w<dis[e[j].t] then
                              begin
                                 dis[e[j].t]:=dis[q[head]]+e[j].w;
                                 if c[e[j].t]>=c[e[b[i]].t] then
                                    begin
                                       if not ing[e[j].t] then begin inc(cnt);g[cnt]:=e[j].t;ing[e[j].t]:=true; end;
                                    end
                                 else if not inq[e[j].t] then
                                    begin
                                       inc(tail);if tail>100000 then tail:=1;
                                       q[tail]:=e[j].t;inq[e[j].t]:=true;
                                    end;
                              end;
                        end;
                     j:=e[j].next;
                  end;
               inq[q[head]]:=false;
            end;
         for j:=1 to cnt do begin ins(c[g[j]],dis[g[j]]+d[c[g[j]]]);ing[g[j]]:=false; end;
         while (p[1]<c[e[b[i]].t]) and (tot>0) do del;
         if tot=0 then writeln(-1) else writeln(h[1]);
         dis[e[b[i]].t]:=dis[e[b[i-1]].t]+e[b[i]].w;
      end;
   close(input);close(output);
end.