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])
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 var 2 sp:array[0..100,0..100]of longint; 3 a:array[0..20,0..100]of boolean; 4 first,next,last,w,q:array[0..10000]of longint; 5 dis:array[0..20]of longint; 6 flag,can:array[0..20]of boolean; 7 n,m,k,e,tot,l,r,head,tail:longint; 8 9 function min(x,y:longint):longint; 10 begin 11 if x<y then exit(x); 12 exit(y); 13 end; 14 15 procedure spfa; 16 var 17 i,j:longint; 18 begin 19 head:=1; 20 tail:=1; 21 q[1]:=1; 22 for i:=2 to m do 23 dis[i]:=21470000; 24 dis[1]:=0; 25 for i:=2 to m do 26 flag[m]:=false; 27 flag[1]:=true; 28 for i:=1 to m do 29 begin 30 can[i]:=true; 31 for j:=l to r do 32 if a[i,j] then can[i]:=false; 33 end; 34 while head<=tail do 35 begin 36 i:=first[q[head]]; 37 while i<>0 do 38 begin 39 if can[last[i]] then 40 if dis[last[i]]>dis[q[head]]+w[i] then 41 begin 42 dis[last[i]]:=dis[q[head]]+w[i]; 43 if flag[last[i]]=false then 44 begin 45 inc(tail); 46 flag[last[i]]:=true; 47 q[tail]:=last[i]; 48 end; 49 end; 50 i:=next[i]; 51 end; 52 flag[q[head]]:=false; 53 inc(head); 54 end; 55 sp[l,r]:=dis[m]; 56 end; 57 58 procedure insert(x,y,z:longint); 59 begin 60 inc(tot); 61 last[tot]:=y; 62 next[tot]:=first[x]; 63 first[x]:=tot; 64 w[tot]:=z; 65 end; 66 67 procedure init; 68 var 69 i,j,x,y,z:longint; 70 begin 71 read(n,m,k,e); 72 for i:=1 to e do 73 begin 74 read(x,y,z); 75 insert(x,y,z); 76 insert(y,x,z); 77 end; 78 read(e); 79 for i:=1 to e do 80 begin 81 read(x,y,z); 82 for j:=y to z do 83 a[x,j]:=true; 84 end; 85 for l:=1 to n do 86 for r:=l to n do 87 spfa; 88 end; 89 90 var 91 f:array[0..100,0..100]of longint; 92 93 procedure work; 94 var 95 i,j,l:longint; 96 begin 97 for i:=1 to n do 98 for j:=i to n do 99 f[i,j]:=sp[i,j]*(j-i+1); 100 for i:=1 to n-1 do 101 for j:=1 to n-i do 102 for l:=j to i+j-1 do 103 f[j,i+j]:=min(f[j,i+j],f[j,l]+k+f[l+1,i+j]); 104 write(f[1,n]); 105 end; 106 107 begin 108 init; 109 work; 110 end.