BZOJ1601: [Usaco2008 Oct]灌水

1601: [Usaco2008 Oct]灌水

Time Limit: 5 Sec  Memory Limit: 162 MB
Submit: 1280  Solved: 839
[Submit][Status]

Description

Farmer John已经决定把水灌到他的n(1<=n<=300)块农田，农田被数字1到n标记。把一块土地进行灌水有两种方法，从其他农田饮水，或者这块土地建造水库。 建造一个水库需要花费wi(1<=wi<=100000),连接两块土地需要花费Pij(1<=pij<=100000,pij=pji,pii=0). 计算Farmer John所需的最少代价。

Input

*第一行：一个数n

*第二行到第n+1行：第i+1行含有一个数wi

*第n+2行到第2n+1行：第n+1+i行有n个被空格分开的数，第j个数代表pij。

Output

*第一行：一个单独的数代表最小代价.

Sample Input

4
5
4
4
3
0 2 2 2
2 0 3 3
2 3 0 4
2 3 4 0

Sample Output

9


输出详解：

Farmer John在第四块土地上建立水库，然后把其他的都连向那一个，这样就要花费3+2+2+2=9

HINT

 

Source

资格赛

题解：

刚开始看着觉得是网络流。。。YY了十几分钟一直过不了样例，后来一百度就惊呆了。。。

建立虚拟节点，向每个点连边，权值为在该点建立水库的费用，对这n+1个点求最小生成树。。。

虚拟节点的连通保证了至少有一个点建造了水库

漂亮！

代码：

var a,b,c,fa:array[0..100000] of longint;
    ans,i,j,n,tot:longint;
    procedure swap(var x,y:longint);
     var t:longint;
         begin
         t:=x;x:=y;y:=t;
         end;
function find(x:longint):longint;
 begin
 if fa[x]<>x then fa[x]:=find(fa[x]);
 exit(fa[x]);
 end;

procedure sort(l,r:longint);
 var i,j,x,y:longint;
 begin
 i:=l;j:=r;x:=c[(i+j)>>1];
 repeat
  while c[i]<x do inc(i);
  while c[j]>x do dec(j);
  if i<=j then
   begin
   swap(a[i],a[j]);swap(b[i],b[j]);swap(c[i],c[j]);
   inc(i);dec(j);
   end;
 until i>j;
 if i<r then sort(i,r);
 if j>l then sort(l,j);
 end;
procedure init;
 begin
   readln(n);
   for i:=1 to n do
    begin
      readln(c[i]);
      a[i]:=i;b[i]:=n+1;
    end;
   tot:=n;
   for i:=1 to n do
    begin
    for j:=1 to i-1 do
     begin
       inc(tot);
       read(c[tot]);a[tot]:=i;b[tot]:=j;
     end;
    readln;
    end;
 end;
procedure main;
 begin
 sort(1,tot);
 ans:=0;j:=1;
 for i:=1 to n+1 do fa[i]:=i;
 for i:=1 to n do
   begin
   while find(a[j])=find(b[j]) do inc(j);
   fa[find(a[j])]:=find(b[j]);inc(ans,c[j]);
   end;
 writeln(ans);
 end;

begin
 assign(input,'input.txt');assign(output,'output.txt');
 reset(input);rewrite(output);
 init;
 main;
 close(input);close(output);
end.