部落划分 (二分+并查集/kruskal)

问题:

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

解:

再次明确几点重要的东西 别搞忘了

1:处理集合问题一定要用并查集

2: 最近的最远 最远的最近 不是二分就是生成树

考试的时候想到正解的 突然卡住 难受QWQ

$1$

二分一个最小距离

暴力枚举节点之间的距离,小于枚举的答案就合并成一个集合 

                    两个部落的距离只会比这个更小所以必须合并

大于的话是另一个集合(考试的时候我就是在这里卡住了 我还以为他会被别的点更新 但实际上是不会的 因为暴力枚举的是所有点的距离)

并查集套用即可

//
#include<bits/stdc++.h>
using namespace std;
#define maxnn 100000
#define ll long long 
ll n,k;
ll f[maxnn];
ll x[maxnn],y[maxnn];
double dist(ll x1,ll x2,ll y1,ll y2)
{
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
int gf(int v)
{
    if(f[v]==v) return v;
    else return f[v]=gf(f[v]);
}
void merge(int v,int u)
{
    int f1=gf(v);
    int f2=gf(u);
    if(f1!=f2) f[f1]=f2;
}
bool isok(double ttt)
{
    for(int i=1;i<=n;i++)
    {
        f[i]=i;
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(i!=j)
            {
                double dis=dist(x[i],x[j],y[i],y[j]);
                if(dis<ttt) merge(i,j);
            }
        }
        
    }
    int sec=0;
        for(int i=1;i<=n;i++)
    {
        if(f[i]==i) sec++;
    }
    return sec>=k;
}
int main()
{
    
    cin>>n>>k;
    double  l=0,r=12000,eps=1e-1;
    for(int i=1;i<=n;i++)
    {
        scanf("%lld%lld",&x[i],&y[i]);
    }
    int cnt=0;
    while(cnt<=60) 
    {
        cnt++;
        double mid=(l+r)/2;
        if(isok(mid)) l=mid;
        else r=mid;
    }
    printf("%.2lf",l);
}

$2$ kruskal 

注意到最大距离最小  不只可以二分 还可以跑生成树 

生成树也是处理集合问题的利器

按照边从小到大排序 

优先合并小的 直到联通块(集合)的值刚好 ==k;

输出下一条边

来自hyhjulao 的code:

#include<bits/stdc++.h>
using namespace std;
template<typename T>
void rin(T &t)
{t=0;int k=1;char c=getchar();while(!isdigit(c)){if(c=='-')k=-1;c=getchar();}while(isdigit(c)){t=t*10+c-'0';c=getchar();}t*=k;}
const int maxn=1e3+5;
struct node
{
    int a,b;
    double c;
}e[maxn*maxn];
int n,k,fa[maxn],tot;
bool cmp(node x,node y)
{return x.c<y.c;}
typedef pair<int,int> P;
P NODE[maxn];
inline double dis(int i,int j)
{return sqrt((NODE[i].first-NODE[j].first)*(NODE[i].first-NODE[j].first)+(NODE[i].second-NODE[j].second)*(NODE[i].second-NODE[j].second));}
int getfa(int v)
{return fa[v]==v?v:fa[v]=getfa(fa[v]);}
int main()
{
    rin(n);rin(k);
    for(int i=1;i<=n;i++)
    {
        rin(NODE[i].first);rin(NODE[i].second);
        fa[i]=i;
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<i;j++)
    {
        e[++tot].a=i;e[tot].b=j;
        e[tot].c=dis(i,j);
    }
    sort(e+1,e+1+tot,cmp);
    int cnt=n,fx,fy;
    double ans;
    for(int i=1;i<=tot;i++)
    {
        fx=getfa(e[i].a);
        fy=getfa(e[i].b);
        if(fx!=fy)
        {
            if(cnt==k){ans=e[i].c;break;}
            else cnt--;
            fa[fx]=fy;
        }
    }
    printf("%.2lf",ans);

}

猜你喜欢

转载自www.cnblogs.com/OIEREDSION/p/11370731.html