BZOJ 1821 JSOI2010 部落划分 Group prim

Description

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

Input

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

Output

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

Sample Input

Sample Output

1.00

题解

最小生成树求第K大边

代码

 #include <cstdio>

 #include <cmath>

 #include <cstring>

 #include <ctime>

 #include <iostream>

 #include <algorithm>

 #include <set>

 #include <vector>

 #include <queue>

 #include <typeinfo>

 #include <map>

 #include <stack>

 typedef long long ll;

 #define inf 0x7fffffff

 using namespace std;

 inline ll read()

 {

     ll x=,f=;

     char ch=getchar();

     while(ch<''||ch>'')

     {

         if(ch=='-')f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

     {

         x=x*+ch-'';

         ch=getchar();

     }

     return x*f;

 }

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

 double mp[][];

 double lowcost[];

 int used[];

 int n;

  int k;

 int kk;

 double  ans[];

 void prim()

 {

    int x=;

    double mn;

    memset(used,,sizeof(used));

    memset(lowcost,,sizeof(lowcost));

    for(int i=;i<=n;i++)

     lowcost[i]=mp[x][i];

    used[x]=;

    for(int i=;i<=n;i++)

    {

        mn=inf*1.0;

        for(int j=;j<=n;j++)

        {

            if(!used[j]&&lowcost[j]<mn)

            {

                x=j;

                mn=lowcost[j];

            }

        }

        ans[++kk]=mn;

        used[x]=;

        for(int j=;j<=n;j++)

        {

            if(!used[j]&&mp[x][j]<lowcost[j])

             lowcost[j]=mp[x][j];

        }

    }

 }

 int main()

 {

  double x[];

  double y[];

     scanf("%d%d",&n,&k);

     for(int i=;i<=n;i++)

     {

         scanf("%lf%lf",&x[i],&y[i]);

     }

     for(int i=;i<=n;i++)

     {

         for(int j=;j<=n;j++)

         {

             if(i==j)mp[i][j]=;

             else

             {

                 mp[i][j]=sqrt((y[i]-y[j])*(y[i]-y[j])+(x[i]-x[j])*(x[i]-x[j]));

             }

         }

     }

     kk=;

     prim();

     sort(ans+,ans++kk);

     printf("%.2f\n",ans[kk-k+]);

     return ;

 }

巴特西

BZOJ 1821 JSOI2010 部落划分 Group prim

Description

最新文章

热门文章