最短网络

题目描述

农民约翰被选为他们镇的镇长！他其中一个竞选承诺就是在镇上建立起互联网，并连接到所有的农场。当然，他需要你的帮助。

约翰已经给他的农场安排了一条高速的网络线路，他想把这条线路共享给其他农场。为了用最小的消费，他想铺设最短的光纤去连接所有的农场。

你将得到一份各农场之间连接费用的列表，你必须找出能连接所有农场并所用光纤最短的方案。

每两个农场间的距离不会超过100000

 

INPUT FORMAT

第一行：农场的个数，N（3<=N<=100）。

第二行..结尾: 后来的行包含了一个N*N的矩阵,表示每个农场之间的距离。理论上，他们是N行，每行由N个用空格分隔的数组成，实际上，他们限制在80个字符，因此，某些行会紧接着另一些行。当然，对角线将会是0，因为不会有线路从第i个农场到它本身。

 

SAMPLE INPUT

4

0 4 9 21

4 0 8 17

9 8 0 16

21 17 16 0

 

OUTPUT FORMAT

只有一个输出，其中包含连接到每个农场的光纤的最小长度。

 

 

SAMPLE OUTPUT

28

思路：Kruskal(克鲁斯卡尔)算法+并查集

代码：

var

        a:array[1..10001,1..3]of longint;

        n,m,i,j,k,x,y,z,s,ls,h:longint;

        b:array[0..10001]of longint;

        f:array[0..1001,0..1001]of longint;

procedure qs(l,r:longint);  //快排

var

        i,j,m:longint;

        t:array[1..3]of longint;

begin

        i:=l;

        j:=r;

        m:=a[(l+r) div 2,1];

        repeat

                while a[i,1]<m do inc(i);

                while a[j,1]>m do dec(j);

                        if i<=j then

                        begin

                                t:=a[i];

                                a[i]:=a[j];

                                a[j]:=t;

                                inc(i);

                                dec(j);

                        end;

        until i>j;

        if l<j then qs(l,j);

        if i<r then qs(i,r);

end;

function shang(k:longint):longint;  //查父亲节点（并查集优化）

var

        ls:longint;

begin

        if b[k]=k then exit(k);

        ls:=shang(b[k]);

        b[k]:=ls;

        exit(ls);

end;

begin

        readln(m);

        for i:=1 to m do

        for j:=1 to m do

        read(f[i,j]);

        k:=0;

        for i:=1 to m do

        for j:=1 to i-1 do

        if f[i,j]<>0 then

        begin

                inc(K);

                a[k,1]:=f[i,j];

                a[k,2]:=i;

                a[k,3]:=j;

        end;

        qs(1,k);

        h:=k;

        k:=0;

        for i:=0 to m do

        b[i]:=i;

        for i:=1 to h do

        begin

                if shang(a[i,2])<>shang(a[i,3]) then

                begin

                        b[b[a[i,2]]]:=b[a[i,3]];

                        s:=s+a[i,1];

                        inc(k);

                        if k=m-1 then

                        break;

                end;

        end;

        writeln(s);

end.