1003: [ZJOI2006]物流运输trans

Description

很久以前，在一个遥远的星系，一个黑暗的帝国靠着它的超级武器统治者整个星系。某一天，凭着一个偶然的机遇，一支反抗军摧毁了帝国的超级武器，并攻下了星系中几乎所有的星球。这些星球通过特殊的以太隧道互相直接或间接地连接。 但好景不长，很快帝国又重新造出了他的超级武器。凭借这超级武器的力量，帝国开始有计划地摧毁反抗军占领的星球。由于星球的不断被摧毁，两个星球之间的通讯通道也开始不可靠起来。现在，反抗军首领交给你一个任务：给出原来两个星球之间的以太隧道连通情况以及帝国打击的星球顺序，以尽量快的速度求出每一次打击之后反抗军占据的星球的连通快的个数。（如果两个星球可以通过现存的以太通道直接或间接地连通，则这两个星球在同一个连通块中）。
Input

输入文件第一行包含两个整数，N (1 <= N <= 2M) 和M (1 <= M <= 200,000)，分别表示星球的数目和以太隧道的数目。星球用0~N-1的整数编号。接下来的M行，每行包括两个整数X, Y，其中（0<=X<>Y
Output

输出文件的第一行是开始时星球的连通块个数。接下来的N行，每行一个整数，表示经过该次打击后现存星球的连通块个数。
Sample Input
8 13
0 1
1 6
6 5
5 0
0 6
1 2
2 3
3 4
4 5
7 1
7 2
7 6
3 6
5
1
6
3
5
7
Sample Output
1
1
1
2
3
3

额，不会写，看题解

先预处理出sp[l,r]表示l到r天都能用的最短路

然后dp，f[i,j]表示i到j天的最小花费

f[i,j]=min(f[i,k]+cost+f[k+1,j])

var
    sp:array[0..100,0..100]of longint;
    a:array[0..20,0..100]of boolean;
    first,next,last,w,q:array[0..10000]of longint;
    dis:array[0..20]of longint;
    flag,can:array[0..20]of boolean;
    n,m,k,e,tot,l,r,head,tail:longint;

function min(x,y:longint):longint;
begin
    if x<y then exit(x);
    exit(y);
end;

procedure spfa;
var
    i,j:longint;
begin
    head:=1;
    tail:=1;
    q[1]:=1;
    for i:=2 to m do
      dis[i]:=21470000;
    dis[1]:=0;
    for i:=2 to m do
      flag[m]:=false;
    flag[1]:=true;
    for i:=1 to m do
      begin
        can[i]:=true;
        for j:=l to r do
          if a[i,j] then can[i]:=false;
      end;
    while head<=tail do
      begin
        i:=first[q[head]];
        while i<>0 do
          begin
            if can[last[i]] then
            if dis[last[i]]>dis[q[head]]+w[i] then
            begin
              dis[last[i]]:=dis[q[head]]+w[i];
              if flag[last[i]]=false then
              begin
                inc(tail);
                flag[last[i]]:=true;
                q[tail]:=last[i];
              end;
            end;
            i:=next[i];
          end;
        flag[q[head]]:=false;
        inc(head);
      end;
    sp[l,r]:=dis[m];
end;

procedure insert(x,y,z:longint);
begin
    inc(tot);
    last[tot]:=y;
    next[tot]:=first[x];
    first[x]:=tot;
    w[tot]:=z;
end;

procedure init;
var
    i,j,x,y,z:longint;
begin
    read(n,m,k,e);
    for i:=1 to e do
      begin
        read(x,y,z);
        insert(x,y,z);
        insert(y,x,z);
      end;
    read(e);
    for i:=1 to e do
      begin
        read(x,y,z);
        for j:=y to z do
          a[x,j]:=true;
      end;
    for l:=1 to n do
      for r:=l to n do
        spfa;
end;

var
    f:array[0..100,0..100]of longint;

procedure work;
var
    i,j,l:longint;
begin
    for i:=1 to n do
      for j:=i to n do
        f[i,j]:=sp[i,j]*(j-i+1);
    for i:=1 to n-1 do
      for j:=1 to n-i do
        for l:=j to i+j-1 do
          f[j,i+j]:=min(f[j,i+j],f[j,l]+k+f[l+1,i+j]);
    write(f[1,n]);
end;

begin
    init;
    work;
end.