最近点对,zoj 2107,hdu1007

本文分享了《计算机算法设计与分析》中最近点对算法的修正版本,并提供了一段通过ZOJ平台验证的C++代码实现。修正了书中33页的一个坐标轴错误,并详细展示了算法的具体实现细节。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

今天看《计算机算法设计与分析》第二版里的最近点对,


按着写了,可以过hud的,不能过zoj的,郁闷死了。


查了半天,发现书上印错了,33页第一行第一个Y改成X


下面是可以过zoj的代码:


#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <math.h>
#include <string.h>
using namespace std;
#define N 120005
struct point
{
	double x,y;
	int id;
}aa[N],bb[N],cc[N];
bool cmp_x(point a,point b)
{
	return a.x<b.x;
}
bool cmp_y(point a,point b)
{
	return a.y<b.y;
}
double dis(point a,point b)
{
	double dx=a.x-b.x,dy=a.y-b.y;
	return sqrt(dx*dx+dy*dy);
}

void merge(point c[],point d[],int l,int m,int r)
{
	int i=l,j=m+1,k=l;
	while(i<=m&&j<=r)
		if(c[i].y<=c[j].y)d[k++]=c[i++];
		else d[k++]=c[j++];
	if(i>m)
		while(j<=r)
			d[k++]=c[j++];
	else 
		while(i<=m)
			d[k++]=c[i++];
}
double closest(point X[],point Y[],point Z[],int l,int r)
{
	if(r-l==1)return dis(X[l],X[r]);
	if(r-l==2)
	{
		double d1=dis(X[l],X[r]);
		double d2=dis(X[l+1],X[r]);
		double d3=dis(X[l],X[l+1]);
		if(d1<d2&&d1<d3)return d1;
		if(d2<d3)return d2;
		return d3;
	}
	int m=(l+r)/2;
	int f=l,g=m+1;
	int i,j,k;
	double d1,d2,dm;
	for(i=l;i<=r;i++)
		if(Y[i].id>m)Z[g++]=Y[i];
		else Z[f++]=Y[i];
	d1=closest(X,Z,Y,l,m);
	d2=closest(X,Z,Y,m+1,r);
	dm=min(d1,d2);
	merge(Z,Y,l,m,r);
	k=l;
	for(i=l;i<=r;i++)
		if(fabs(X[m].x-Y[i].x)<dm)Z[k++]=Y[i];
	for(i=l;i<k;i++)
	{
		for(j=i+1;j<k&&Z[j].y-Z[i].y<dm;j++)
		{
			double tmp=dis(Z[i],Z[j]);
			if(tmp<dm)dm=tmp;
		}
	}
	return dm;
}
int main()
{
	int n,i;
	while(scanf("%d",&n)&&n)
	{
		for(i=0;i<n;i++)
			cin>>aa[i].x>>aa[i].y;
		sort(aa,aa+n,cmp_x);
		for(i=0;i<n;i++)
		{
			aa[i].id=i;
			bb[i]=aa[i];
		}
		//memcpy(b,a,n*sizeof(a[0]));
		sort(bb,bb+n,cmp_y);
		printf("%.2lf\n",closest(aa,bb,cc,0,n-1)/2);
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值