2458: [BeiJing2011]最小三角形 分治

Description
Xaviera现在遇到了一个有趣的问题。
平面上有N个点，Xaviera想找出周长最小的三角形。
由于点非常多，分布也非常乱，所以Xaviera想请你来解决这个问题。
为了减小问题的难度，这里的三角形也包括共线的三点。

题解：

分治，利用分治得到的答案缩小选点的范围，时间复杂度未知。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=200010;
struct Point{double x,y;}p[Maxn],tmp[Maxn];
bool cmpx(Point a,Point b){return a.x<b.x;}
bool cmpy(Point a,Point b){return a.y<b.y;}
int n;double ans=2147483647.0;
double dis(int a,int b){return sqrt((p[a].x-p[b].x)*(p[a].x-p[b].x)+(p[a].y-p[b].y)*(p[a].y-p[b].y));}
double _dis(int a,int b){return sqrt((tmp[a].x-tmp[b].x)*(tmp[a].x-tmp[b].x)+(tmp[a].y-tmp[b].y)*(tmp[a].y-tmp[b].y));}
void solve(int l,int r)
{
    if(l+2>r)return;
    if(l+2==r){ans=min(ans,dis(l,l+1)+dis(l+1,l+2)+dis(l,l+2));return;}
    int mid=l+r>>1,cnt=0,j=1;
    solve(l,mid),solve(mid+1,r);double lim=ans/2.0;
    for(int i=l;i<=r;i++)if(abs(p[i].x-p[mid].x)<=lim)tmp[++cnt]=p[i];
    sort(tmp+1,tmp+1+cnt,cmpy);
    for(int i=1;i<=cnt;i++)
    {
        while(abs(tmp[i].y-tmp[j].y)<=lim&&j<cnt)j++;
        for(int k=i+1;k<=j;k++)
        for(int l=i+1;l<k;l++)
        ans=min(ans,_dis(i,k)+_dis(i,l)+_dis(k,l));
    }
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%lf%lf",&p[i].x,&p[i].y);
    sort(p+1,p+1+n,cmpx);
    solve(1,n);
    printf("%.6lf",ans);
}