HDU 1007 Quoit Design 最近点对
2015-09-11 07:21
351 查看
题意:懒得去找链接了。
方法:递归最近点对。
解析:
题中要求最近点对的距离除以2。
如何求最近点对呢?
有这么一种递归的方法。
首先我们按x排序所有的点。
分别递归左半边的最近点对距离以及右半边最近点对距离。
直至点的个数小于三再回溯。
但是这不够,显然有一种情况是左右半边的两个点间距离最短。
如何求这个呢?
这也就是递归回溯的一个合并过程。
首先我们假设左边求出来的距离为d1,右边求出来的距离为d2。
那么我们把左右半边中x距离相距小于Min(d1,d2)的点都推入数组中。
然后按y坐标排序,枚举扫一下y距离小于ans的点对即可。
这复杂度会不会爆炸呢?
不会,算导上给出了我们最后的枚举的期望枚举点个数为7,所以只是一个常数,并不必在意。
所以这么写的复杂度是
T(n)=T(n2)+O(nlogn)=O(nlog2n)T(n)=T(\frac{n}{2})+O(nlogn)=O(nlog^2n)
可以过。
但是算导上给出了更优越的O(nlogn)算法。
其优越性体现在将如上的算法每一次递归排序y这一步省略。
它是优先处理出一个按照y排序的数组,每一次扫这个数组的话,就可以得到一个y已经排好序的数组。
所以复杂度掉到了O(nlogn)
代码:
#include <cmath> #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define N 100100 #define INF 0x7f7f7f7f using namespace std; int n; struct Point { double x,y; friend istream& operator >> (istream &_,Point &a) { scanf("%lf%lf",&a.x,&a.y); return _; } Point(){} Point(double _x,double _y):x(_x),y(_y){} Point operator - (const Point &a) {return Point(x-a.x,y-a.y);} double operator * (const Point &a) {return x*a.x+y*a.y;} friend bool operator < (Point a,Point b) { return a.x<b.x; } }pt ,a ; bool cmp(Point a,Point b) { return a.y<b.y; } double dis(Point &a,Point &b) { return sqrt((a-b)*(a-b)); } double find(int l,int r) { if(l==r)return INF; if(l+1==r)return dis(pt[l],pt[r]); int mid=(l+r)>>1; double ans=min(find(l,mid),find(mid+1,r)); int tot=0; for(int i=l;i<=r;i++) { if(pt[i].x>=pt[mid].x-ans&&pt[i].x<=pt[mid].x+ans) a[++tot]=pt[i]; } sort(a+1,a+tot+1,cmp); for(int i=1;i<=tot;i++) { for(int j=i+1;j<=tot;j++) { if(dis(a[i],a[j])>=ans)break; ans=min(ans,dis(a[i],a[j])); } } return ans; } int main() { while(~scanf("%d",&n)&&n!=0) { for(int i=1;i<=n;i++)cin>>pt[i]; sort(pt+1,pt+n+1); printf("%.2lf\n",find(1,n)/2); } }
相关文章推荐
- OSChina 周五乱弹 —— 买个 6s 送老师
- C++Primer第五版 2.3.3节练习
- Servlet3.0引入的新特性
- 位运算
- fastcgi配置
- 不翻墙如何调试Spark源码?
- Leetcode 82. Remove Duplicates from Sorted List II
- ref out
- PHP-P函数打印数据
- AndroidManifest.xml文件的作用和简单使用
- PHP-C函数读取配置信息
- LeetCode Binary Tree Maximum Path Sum
- LeetCode Binary Tree Maximum Path Sum
- Unity5内部渲染的优化3:移除固定功能
- Leetcode: Read N Characters Given Read4
- CentOS上解压ZIP乱码的解决办法
- 遍历二叉树的各种操作(非递归遍历) C++
- CentOS上解压ZIP乱码的解决办法
- 11g之前 修改主机名后如何启动使用ASM存储的数据库
- Leetcode NO.220 Contains Duplicate III