最近点对问题模板题
在与联盟的战斗中屡战屡败后,帝国撤退到了最后一个据点。
依靠其强大的防御系统,帝国击退了联盟的六波猛烈进攻。
经过几天的苦思冥想,联盟将军亚瑟终于注意到帝国防御系统唯一的弱点就是能源供应。
该系统由N个核电站供应能源,其中任何一个被摧毁都会使防御系统失效。
将军派出了N个特工进入据点之中,打算对能源站展开一次突袭。
不幸的是,由于受到了帝国空军的袭击,他们未能降落在预期位置。
作为一名经验丰富的将军,亚瑟很快意识到他需要重新安排突袭计划。
他现在最想知道的事情就是哪个特工距离其中任意一个发电站的距离最短。
你能帮他算出来这最短的距离是多少吗?
输入格式
输入中包含多组测试用例。
第一行输入整数T,代表测试用例的数量。
对于每个测试用例,第一行输入整数N。
接下来N行,每行输入两个整数X和Y,代表每个核电站的位置的X,Y坐标。
在接下来N行,每行输入两个整数X和Y,代表每名特工的位置的X,Y坐标。
输出格式
每个测试用例,输出一个最短距离值,结果保留三位小数。
每个输出结果占一行。
数据范围
1≤N≤1000001≤N≤100000
0≤X,Y≤10000000000≤X,Y≤1000000000
样例
输入样例:
2
4
0 0
0 1
1 0
1 1
2 2
2 3
3 2
3 3
4
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
输出样例:
1.414
0.000
作者:秦淮岸灯火阑珊
链接:https://www.acwing.com/solution/acwing/content/826/
来源:AcWing
(我还不能写出这么厉害的题解,只能先搬一下,以后自己看看)
算法分析:这是一道经典的最近点对问题的模板,这里略微讲述一下这种问题的解法,首先呢,我们将这些点的xx坐标,为第一关键字,yy坐标为第二关键字,从小到大排序,接着我们取一个中点midmid点,将midmid点左边的点,统统归为平面d1d1,然后midmid点右边的点统统归为平面d2d2如下图所示,
至于midmid点可以随便放到一个平面即可,然后我们设ans1ans1为平面d1d1的最近距离,ans2ans2为平面d2d2的最近距离,那么最终答案就为ans=min(ans1,ans2)ans=min(ans1,ans2),看上取很合理,但是实际上还有一种情况,那就是比如说上图中的2和4,是不是距离更近呢,所以我们还要考虑一种非常恶心的方案,就是左边一个点,右边一个点,而且这两个点距离可能是最小的方案。
最难点就是这个方案,因为有很多高大上的优化,但是实际上玄学优化,就可以让我们很容易地理解这种方法,而且时间复杂度不是很高,满足一般用法。
我们知道,当前局部最优解为ans,那么最终全局的答案肯定是<=ans,这一点很明确,既然如此,那么最终答案的点对,就只可能在下图中的区域里面出现,因为如果不在这个区域,那么肯定是大于ans的。
)
既然如此,那么我们就可以只需要在这个区域里面寻找答案即可,至于寻找的过程,我们还要注意一点,那就是我们需要将这个区间里面的数,将yy坐标排序,从小到大,这样便于我们接下来的break剪枝,具体实现看代码。
#include
#include
include
#include
using namespace std;
const int N=100100;
struct node{
double x,y;
int f;
}a[N<<1],t[N<<1];
#define INF 1e8
int n,m;
int cmp1(node a,node b)
{
if(a.x==b.x)
return a.y>1;\\找到中间的数
if(r-l==0) \\只有一个数的时候返回无限大
return INF;
if(r-l==1)\\有两个数的时候就返回他俩的距离
return dist(a[l],a[r]);
double ans=min(slove(l,mid),slove(mid+1,r));\\通过递归,找到按照x坐标排序,两边距离里的最小值,
int cnt=0;
\\还有第三种情况就是在中间位置,两个点里的最近,先找出在中间
与上面最小值之间的点,添加到另一个数组里。
for(int i=l;i<=r;i++)
{
if(a[i].x>=a[mid].x-ans&&a[i].x<=a[mid].x+ans)
t[++cnt]=a[i];
}
sort(t+1,t+cnt+1,cmp2);
\\按照y坐标的大小排序
for(int i=1;i<=cnt;i++)
{
for(int j=1;j<=cnt;j++)
{ \\这个是吧每一个点都当做中心点,另一个点y坐标大于中心点y坐标与最小值的值的时候,就退出,剪枝。
if(t[j].y>=t[i].y+ans)
break;
ans=min(ans,dist(t[i],t[j]));
}
}
return ans;
}
int main()
{ cin>>m;
while(m--)
\\将这些坐标都输入到一个结构体数组里,按照x坐标的从小到大排序
{
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>a[i].x>>a[i].y;
a[i].f=1;
}
for(int i=n+1;i<=n<<1;i++)
{
cin>>a[i].x>>a[i].y;
a[i].f=0;
}
n<<=1;
sort(a+1,a+n+1,cmp1);
printf("%.3lf\n", slove(1,n));
}
return 0;
}