• BZOJ1821 [JSOI2010] Group 部落划分 Group


    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1821

    Description

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

    Input

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

    Output

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

    贪心,把边权小的边划分到部落内部,第n - k + 1条边的边权就是答案。

     1 #include <iostream>
     2 #include <cstdio>
     3 #include <algorithm>
     4 #include <cstring>
     5 #include <cmath>
     6 #define rep(i,l,r) for(int i=l; i<=r; i++)
     7 #define clr(x,y) memset(x,y,sizeof(x))
     8 #define travel(x) for(Edge *p=last[x]; p; p=p->pre)
     9 using namespace std;
    10 const int INF = 0x3f3f3f3f;
    11 const int maxn = 1010;
    12 inline int read(){
    13     int ans = 0, f = 1;
    14     char c = getchar();
    15     for(; !isdigit(c); c = getchar())
    16     if (c == '-') f = -1;
    17     for(; isdigit(c); c = getchar())
    18     ans = ans * 10 + c - '0';
    19     return ans * f;
    20 }
    21 struct Edge{
    22     int from,to;
    23     double cost;
    24     inline bool operator < (const Edge &_Tp) const{
    25         return cost < _Tp.cost;
    26     }
    27 }edge[500010];
    28 struct Point{
    29     int x,y;
    30 }t[maxn];
    31 int n,k,cnt=0,tot=0,fa[maxn];
    32 inline double getdis(Point a,Point b){
    33     return sqrt((a.x-b.x) * (a.x-b.x) + (a.y-b.y) * (a.y-b.y));
    34 }
    35 inline void addedge(int x,int y,double dis){
    36     edge[++cnt].from = x; edge[cnt].to = y; edge[cnt].cost = dis;
    37 }
    38 int getfa(int x){
    39     return fa[x] == x ? x : fa[x] = getfa(fa[x]);
    40 }
    41 int main(){
    42     n = read(); k = read();
    43     rep(i,1,n) t[i].x = read(), t[i].y = read(), fa[i] = i;
    44     rep(i,1,n-1) rep(j,i+1,n) addedge(i,j,getdis(t[i],t[j]));
    45     sort(edge+1,edge+cnt+1);
    46     k = n - k;
    47     rep(i,1,cnt){
    48         int a = getfa(edge[i].from), b = getfa(edge[i].to);
    49         if (a == b) continue;
    50         fa[a] = b;
    51         if (++tot > k){
    52             printf("%.2lf
    ",edge[i].cost);
    53             return 0;
    54         }
    55     }
    56     return 0;
    57 }
    View Code
  • 相关阅读:
    com.alibaba.fastjson.JSONObject之对象与JSON转换方法
    Spring最新下载教程(转)
    Docker安装出现的问题
    使用 rsync 的一些记录
    manjaro 相关软件配置
    nginx开启gzip
    Navicat Premium 15 永久破解激活工具及安装教程(亲测可用)
    基于Spring Boot+Shiro+Vue+Element实现的后台权限管理系统
    在线激活win11、win10、win8/8.1、win7和office2021/2019/2016/2013/2010神器
    SQLyog的下载、安装、破解、配置(MySQL可视化工具安装)
  • 原文地址:https://www.cnblogs.com/jimzeng/p/bzoj1821.html
Copyright © 2020-2023  润新知