[JSOI2010]Group部落划分

15 篇文章 0 订阅
2 篇文章 0 订阅

这道题可以贪心做

每次把距离最近的点对合并到一起

直到不能合并……

我2b的想到什么最近点对之类的……

一看n=1000,其实就是暴力求每两个点之间的距离,排个序并查集就行……


//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
#include<set>
#include<map>
using namespace std;
//Macro
#define	rep(i,a,b)	for(int i=a,tt=b;i<=tt;++i)
#define	drep(i,a,b)	for(int i=a,tt=b;i>=tt;--i)
#define	erep(i,e,x)	for(int i=x;i;i=e[i].next)
#define	irep(i,x)	for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define	read()	(strtol(ipos,&ipos,10))
#define	sqr(x)	((x)*(x))
#define	pb	push_back
#define	PS	system("pause");
typedef	long long	ll;
typedef	pair<int,int>	pii;
const int oo=~0U>>1;
const double inf=1e100;
const double eps=1e-6;
string name="", in=".in", out=".out";
//Var
struct D
{
	double x,y;
}dot[1008];
struct E
{
	double dis;
	int i,j;
	bool operator<(const E &o)const{return dis<o.dis;}
}edge[1000008];
struct U
{
	int f[1008];
	void set(){rep(i,1,1000)f[i]=i;}
	int find(int x){return f[x]==x?x:f[x]=find(f[x]);}
	void Union(int x,int y){f[x]=y;}
}U;
int n,k,tot;
double Dist(D i,D j){return sqrt(sqr(i.x-j.x)+sqr(i.y-j.y));}
void Init()
{
	scanf("%d%d",&n,&k);
	rep(i,1,n)scanf("%lf%lf",&dot[i].x,&dot[i].y);
	rep(i,1,n)rep(j,i+1,n)edge[++tot].dis=Dist(dot[i],dot[j]),edge[tot].i=i,edge[tot].j=j;
	sort(edge+1,edge+1+tot);
}
void Work()
{
	U.set();
	int fx,fy;
	rep(i,1,tot)
	{
		fx=edge[i].i,fy=edge[i].j;
		fx=U.find(fx),fy=U.find(fy);
		if(fx!=fy)
		{
			if(n>k)n--,U.Union(fx,fy);
			else{printf("%.2lf\n",edge[i].dis);return;}
		}
	}
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
	Init();
	Work();
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值