1821: [JSOI2010]Group 部落划分 Group

2023-01-18 18:28:53

1821: [JSOI2010]Group 部落划分 Group

Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 1308 Solved: 627
[Submit][Status]

Description

聪聪研究发现，荒岛野人总是过着群居的生活，但是，并不是整个荒岛上的所有野人都属于同一个部落，野人们总是拉帮结派形成属于自己的部落，不同的部落之间则经常发生争斗。只是，这一切都成为谜团了——聪聪根本就不知道部落究竟是如何分布的。不过好消息是，聪聪得到了一份荒岛的地图。地图上标注了N个野人居住的地点（可以看作是平面上的坐标）。我们知道，同一个部落的野人总是生活在附近。我们把两个部落的距离，定义为部落中距离最近的那两个居住点的距离。聪聪还获得了一个有意义的信息——这些野人总共被分为了K个部落！这真是个好消息。聪聪希望从这些信息里挖掘出所有部落的详细信息。他正在尝试这样一种算法：对于任意一种部落划分的方法，都能够求出两个部落之间的距离，聪聪希望求出一种部落划分的方法，使靠得最近的两个部落尽可能远离。例如，下面的左图表示了一个好的划分，而右图则不是。请你编程帮助聪聪解决这个难题。

Input

第一行包含两个整数N和K(1<=N<=1000,1<k<=n)，分别代表了野人居住点的数量和部落的数量。接下来n行，每行包含两个正整数x,y，描述了一个居住点的坐标(0<="x," y<="10000)。" <="" div="">

Output

输出一行，为最优划分时，最近的两个部落的距离，精确到小数点后两位。

Sample Input

4 2
0 0
0 1
1 1
1 0

Sample Output

1.00

HINT

Source

JSOI2010第二轮Contest1

题解：萌萌哒最小生成树不解释（呵呵呵JSOI居然也出现过这样的题目）

 /**************************************************************

     Problem:

     User: HansBug

     Language: Pascal

     Result: Accepted

     Time: ms

     Memory: kb

 ****************************************************************/

 var

    i,j,k,l,m,n,tt:longint;

    a,d:array[..,..] of longint;

    b:array[..] of extended;

    c:array[..] of longint;

 function getfat(x:longint):longint;inline;

          begin

               if c[x]<>x then c[x]:=getfat(c[x]);

               getfat:=c[x];

          end;

 procedure merge(x,y:longint);inline;

           var a1,a2:longint;

           begin

                a1:=getfat(x);a2:=getfat(y);

                if a1=a2 then exit;

                dec(tt);c[a1]:=a2;

           end;

 function tog(x,y:longint):boolean;inline;

          begin

               exit(getfat(x)=getfat(y));

          end;

 procedure swap(var x,y:longint);inline;

           var z:longint;

           begin

                z:=x;x:=y;y:=z;

           end;

 procedure sort(l,r:longint);

           var i,j:longint;x,y:extended;

           begin

                i:=l;j:=r;x:=b[(l+r) div ];

                repeat

                      while b[i]<x do inc(i);

                      while b[j]>x do dec(j);

                      if i<=j then

                         begin

                              y:=b[i];

                              b[i]:=b[j];

                              b[j]:=y;

                              swap(a[i,],a[j,]);

                              swap(a[i,],a[j,]);

                              inc(i);dec(j);

                         end;

                until i>j;

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

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

           end;

 begin

      readln(n,m);

      for i:= to n do

          readln(d[i,],d[i,]);

      k:=;

      for i:= to n- do

          for j:=i+ to n do

              begin

                   inc(k);

                   a[k,]:=i;

                   a[k,]:=j;

                   b[k]:=sqrt(sqr(d[i,]-d[j,])+sqr(d[i,]-d[j,]));

              end;

      tt:=n;

      for i:= to n do c[i]:=i;

      sort(,k);

      i:=;

      while tt>m do

            begin

                 inc(i);

                 while tog(a[i,],a[i,]) do inc(i);

                 merge(a[i,],a[i,]);

            end;

      inc(i);

      while tog(a[i,],a[i,]) do inc(i);

      writeln(b[i]::);

 end.

码农公寓