poj 3714 最近点对

题意:

给出n个a类点,n个b类点,求a类点到b类点的最近距离。

限制:

1 <= n <= 1e5

0 <= x,y <= 1e9

思路:

点分治

/*poj 3714
  题意:
  给出n个a类点,n个b类点,求a类点到b类点的最近距离。
  限制:
  1 <= n <= 1e5
  0 <= x,y <= 1e9
  思路:
  点分治
 */
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const double INF=1e30;
struct Point{
	double x,y;
	int s;
	Point(){}
	Point(double _x,double _y){ x=_x; y=_y; }
	Point operator - (Point p){ return Point(x-p.x,y-p.y); }
	double dot(Point p){ return x*p.x+y*p.y; }
};
const int N=1e5+5;
Point p[2*N];
int t[2*N];
bool cmpxy(Point a,Point b){
	if(a.x==b.x) return a.y<b.y;
	return a.x<b.x;
}
bool cmpy(int a,int b){
	return p[a].y<p[b].y;
}
double dist(Point a,Point b){
	return sqrt((a-b).dot(a-b));
}
double gao(int l,int r){
	double ret=INF;
	if(l==r) return ret;
	if(l+1==r){
		if(p[l].s!=p[r].s) return dist(p[l],p[r]);
		else return ret;
	}
	int mid=(l+r)>>1;
	double d=min(gao(l,mid),gao(mid+1,r));
	int cnt=0;
	for(int i=l;i<r;++i){
		if(fabs(p[mid].x-p[i].x)<d)
			t[cnt++]=i;
	}
	sort(t,t+cnt,cmpy);
	for(int i=0;i<cnt;++i)
		for(int j=i+1;j<cnt && p[t[j]].y-p[t[i]].y<d;++j)
			if(p[t[i]].s!=p[t[j]].s) d=min(d,dist(p[t[i]],p[t[j]]));
	return d;
}
int main(){
	int T;
	scanf("%d",&T);
	int n;
	while(T--){
		scanf("%d",&n);
		for(int i=0;i<n;++i){
			scanf("%lf%lf",&p[i].x,&p[i].y);
			p[i].s=0;
		}
		for(int i=n;i<2*n;++i){
			scanf("%lf%lf",&p[i].x,&p[i].y);
			p[i].s=1;
		}
		sort(p,p+2*n,cmpxy);
		printf("%.3f\n",gao(0,2*n-1));
	}
	return 0;
}

时间: 2024-10-10 12:44:03

poj 3714 最近点对的相关文章

poj 3714 Raid 分治法求平面最近点对

题意: 给平面上的n个点,求两点间的最短距离. 分析: 分治法,保存点用vector会tle... 代码: //poj 3714 //sep9 #include <iostream> #include <algorithm> #include <cmath> using namespace std; const double INF=1e50; struct P { double x,y; int type; }p[240000],b[240000]; bool cmp

POJ 3714 Raid 最近对点题解

本题是一般最近对点求解,稍微增加点限定:有两个集合点,要求不同集合中的点的最近对. 那么就增加一个判断,如果是同一个集合中的点,那么就返回最大值,其他和一般的最近对点解法一样. 注意:本题数据有重合点,那么就要防止分类的时候溢出. Geeks上的最近对的程序是无法处理有重合点的情况的. #include <stdio.h> #include <stdlib.h> #include <float.h> #include <math.h> #include &l

POJ 3714 Raid(平面最近点对)

解题思路: 分治法求平面最近点对,点分成两部分,加个标记就好了. #include <iostream> #include <cstring> #include <cstdlib> #include <cstdio> #include <cmath> #include <vector> #include <queue> #include <algorithm> #include <iomanip>

平面最近点对 ZOJ 2107 POJ 3714

分治法求最近点对,模板题 第二题稍微判断一下即可 总结一下分治法基本写法: 第一部分:边界判断 第二部分:递归函数 第三部分:区间合并 第一题: #include<stdio.h> #include<math.h> #include<algorithm> using namespace std; const double eps=1e-8; struct Point { double x,y; }p[100005],tmp[100005]; int n; double

POJ 3714 Raid 平面最近点对

题目大意:给出两个集合的点,问这两个集合之间最近的点对的距离是多少. 思路:先要知道平面最近点对的分治算法,剩下的就简单了,只需要在更新答案的时候判断一下两个点是否属于两个集合就可以了. 分治算法总是十分神奇的. 对于平面最近点对,首先按照x坐标排序,然后递归进行分治,每次分治时,先获得分治得到的结果,然后按照这个结果来计算本区间.由于我们只需要计算答案小于这个结果的点对就行了,其中(l,mid)和(mid + 1,r)我们已经得到答案,只需要统计一个点在(l,mid),另一个点在(mid +

POJ 3714 分治/求平面最近点对

第一次见这种问题直接懵圈...没想到分治法这么强大,借鉴了lyd的代码: 代码如下 #include<cstdio> #include<algorithm> #include<cstring> #include<cmath> #include<iostream> using namespace std; const int maxn=200010; struct point{ int x,y,z; }; point a[maxn],b[maxn]

【POJ 3714】 Raid

[题目链接] http://poj.org/problem?id=3714 [算法] 分治求平面最近点对 [代码] #include <algorithm> #include <bitset> #include <cctype> #include <cerrno> #include <clocale> #include <cmath> #include <complex> #include <cstdio> #

计算几何学习8

由于poj炸了 而题单上有很多poj的题 就先开始第二部分了 学习了两个固定算法 最小圆覆盖和平面上最近点对 平面上最近点对采用的是分治的思想 把一个x有序的序列分成A,B左右两部分 当得到A内最近点对距离,B类最近点对距离后 先更新大序列的答案ans A,B间最近点对的产生 显然在x坐标距离mid点不超过ans的区间内产生 我们把这段区间拿出来,对y排序 再逐对更新 值得注意的地方 1)为了确保速度 可以在第二次对y排序的时候采取对标号排序的方式 2)在枚举点对更新的时候注意在y差值 >= a

计算几何题目分类

转载 一.基础题目 1.1 有固定算法的题目 A, 最近点对问题最近点对问题的算法基于扫描线算法.ZOJ 2107    Quoit Design    典型最近点对问题POJ    3714    Raid    变种最近点对问题 B,最小包围圆最小包围圆的算法是一种增量算法,期望是O(n).ZOJ    1450    Minimal Circle  HDU    3007    Buried memory C,旋转卡壳POJ 3608    Bridge Across Islands