JZOJ 3808 【NOIP2014模拟8.25】道路值守

道路值守

Description

有一个有$N$个点，$M$条边的图，问对于该图内的每一对点对($x$,$y$)（$x$≠$y$），问$x$,$y$之间的所有的最短路径覆盖了多少条边。

Data Constraint

$N$<=$500$，任意两个点之间最多只有一条边连接。

Solution

先用$Floyd$最短路算法求出任意两点之间的最短距离。
（下面的$F_{X,Y}$表示$X$,$Y$之间的最短距离，$B_{X,Y}$表示$X$,$Y$之间的连的边的长度）
然后枚举一个点$X$，再枚举$Y$,$Z$，若$F_{X,Y}$+$B_{Y,Z}$=$F_{X,Z}$则说明$X$到$Z$之间的某条最短路径之中覆盖住了($Y$,$Z$)这一条边，$S_Z$++。（$S_i$表示以$i$结尾的上述的那种边的数量）
$S$数组处理好以后，再枚举$K_1$,$K_2$，若$F_{X,K_1}$+$F_{K_1,K_2}$=$F_{X,K_2}$则说明$X$到$K_2$的所有最短路径覆盖了$S_{K_1}$条边（易证），$Ans_{(X,K_2)}$+=$S_{K_1}$，如此这般即可求出答案。

Code(Pacal)

var
    n,m,j,k,l,i,o,x,y,c,kk,ans:longint;
    be,ne,bh,bq:array[0..10000000] of longint;
    bz:array[0..1000] of boolean;
    dd,u,p,qy:qword;
    way:array[0..1000,0..5000] of qword;
    dy:array[0..270000,1..2] of qword;
    ok:array[0..250000] of qword;
    f,bc:array[0..500,0..500] of longint;
    dl,mi,sd:array[0..1000] of longint;
function min(a,b:longint):longint;
    begin
        if a<b then exit(a) else exit(b);
    end;
begin
    readln(n,m);
    for i:=1 to n do
    begin
        for l:=1 to n do f[i,l]:=maxlongint div 2;
        f[i,i]:=0;
    end;
    for i:=1 to n do
    begin
        for l:=1 to n do bc[i,l]:=maxlongint div 2+1;
        bc[i,i]:=0;
    end;
    for i:=1 to m do
    begin
        readln(x,y,c);
        inc(o); ne[o]:=be[x];
        be[x]:=o; bh[o]:=y; bq[o]:=c;
        inc(o); ne[o]:=be[y];
        be[y]:=o; bh[o]:=x; bq[o]:=c;
        f[x,y]:=min(f[x,y],c);  f[y,x]:=f[x,y];
        bc[x,y]:=c;  bc[y,x]:=c;
    end;
    for k:=1 to n do
    for i:=1 to n do
    for l:=1 to n do
    if f[i,k]+f[k,l]<f[i,l] then f[i,l]:=f[i,k]+f[k,l];
    for kk:=1 to n-1 do
    begin
        u:=1; p:=1;
        for i:=1 to n do sd[i]:=0;
        for k:=1 to n do
        for i:=1 to n do
        if i<>k then
        if f[kk,k]+bc[k,i]=f[kk,i] then inc(sd[i]);
        for i:=kk+1 to n do
        begin
            ans:=0;
            for l:=1 to n do
            if f[kk,l]+f[l,i]=f[kk,i] then ans:=ans+sd[l];
            write(ans,' ');
        end;
    end;
end.