BZOJ2458: [BeiJing2011]最小三角形

这题和求平面最近点对有点像，这种题有个分治套路
先把所有点按x排序，我们用(l,r)代表(l,r)内周长最小的三角形
那么当我们处理到(l,r)时，先递归下去处理(l,mid),(mid+1,r)
令ans=min( (l,mid), (mid+1,r) )
那么考虑 (l,mid) 中的点和 (mid+1,r)中的点合并成三角形
对于1个点在左侧，2个点在右侧的情况
我们枚举左侧这个点，那么右侧可能成为答案的点一定在这个蓝色区域内

否则答案不可能比ans优
如果是在这个区域内，画一画就会发现这里面的点数很少，否则(mid+1,r)就会< ans
所以在(l,r)返回时按y排序，合并时就可以在右侧维护两个指针扫过去，这个区域内的点暴力枚举和左边的点组三角形更新ans

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

inline void read(int &x)
{
    char c; int f=1;
    while(!((c=getchar())>='0'&&c<='9')) if(c=='-') f=-1;
    x=c-'0';
    while((c=getchar())>='0'&&c<='9') (x*=10)+=c-'0';
    if(f==-1) x=-x;
}
const int maxn = 410000;
const ll sqr(ll x){return x*x;}

int n;
struct point{int x,y;}p[maxn];
inline bool cmpx(const point x,const point y){return x.x<y.x;}

double dis(point x,point y){return sqrt((double)sqr(x.x-y.x)+(double)sqr(x.y-y.y));}
double sum(point x,point y,point z)
{
    double re=dis(x,y)+dis(y,z)+dis(z,x);
    return re;
}

double ans;
int ld[maxn],lu[maxn],rd[maxn],ru[maxn];
point s1[maxn],s2[maxn];
int n1,n2;
void solve(const int l,const int r)
{
    if(r-l+1<=2)
    {
        if(l!=r&&p[l].y>p[r].y) swap(p[l],p[r]);
        return;
    }
    int mid=l+r>>1,mi=p[mid].x;
    solve(l,mid); solve(mid+1,r);

    n1=n2=0;
    double hf=ans/2.0;
    for(int i=l;i<=mid;i++)
        if((double)mi-p[i].x<=hf) s1[++n1]=p[i];
    for(int i=mid+1;i<=r;i++)
        if((double)p[i].x-mi<=hf) s2[++n2]=p[i];

    if(n1&&n2)
    {
        ld[0]=lu[0]=1;
        for(int i=1;i<=n1;i++)
        {
            ld[i]=ld[i-1],lu[i]=lu[i-1];
            while(ld[i]<=n2&&(double)s1[i].y-s2[ld[i]].y>=ans) ld[i]++;
            while(lu[i]<n2&&(double)s2[lu[i]+1].y-s1[i].y<ans) lu[i]++;
        }
        rd[0]=ru[0]=1;
        for(int i=1;i<=n2;i++)
        {
            rd[i]=rd[i-1],ru[i]=ru[i-1];
            while(rd[i]<=n1&&(double)s2[i].y-s1[rd[i]].y>=ans) rd[i]++;
            while(ru[i]<n1&&(double)s1[ru[i]+1].y-s2[i].y<ans) ru[i]++;
        }
        for(int i=1;i<=n1;i++)
        {
            for(int j=ld[i];j<lu[i];j++) for(int k=j+1;k<=lu[i];k++)
                ans=min(ans,sum(s1[i],s2[j],s2[k]));
        }
        for(int i=1;i<=n2;i++)
        {
            for(int j=rd[i];j<ru[i];j++) for(int k=j+1;k<=ru[i];k++)
                ans=min(ans,sum(s2[i],s1[j],s1[k]));
        }
    }
    n1=n2=0;
    for(int i=l;i<=mid;i++) s1[++n1]=p[i];
    for(int i=mid+1;i<=r;i++) s2[++n2]=p[i];
    int pos=r;
    while(n1|n2)
    {
        if(!n1) p[pos]=s2[n2--];
        else if(!n2) p[pos]=s1[n1--];
        else p[pos]=s1[n1].y>s2[n2].y?s1[n1--]:s2[n2--];
        pos--;
    }
}

int main()
{

    read(n);
    for(int i=1;i<=n;i++) read(p[i].x),read(p[i].y);
    ans=sum(p[1],p[2],p[3]);

    sort(p+1,p+n+1,cmpx);
    solve(1,n);
    printf("%.6lf\n",ans);

    return 0;
}