【BOI2012】Mobile jzoj 3636 计算几何

题目

给出一个端点分别为(0.0)和(L,0)的线段，按照x坐标不减给出n个整点。问线段上的点中，距离最近整点距离的最大值。
n<=1000000,1<=L<=10^9，坐标绝对值不超过10^9。

分析

出题人丧心病狂的用100个数据和0.5s的时限来卡二分……

以下摘抄自题解:

假设所有n个点x坐标均不相同。如果有两个点x坐标相同，那么显然只需要保留y坐标绝对值较小的一个即可。

如图，接下来对于相邻顺序的两个点，作两点连线段的中垂线。显然，中垂线左侧的点距离x坐标较小的点较近，中垂线右侧的点距离x坐标较大的点较近。
因此按x坐标的顺序从小到大扫描n个整点，用一个栈保存可能成为最近点的整点。同时维护栈中两点连线段中垂线与给出线段交点x坐标单调递增。每当加入一个点前，判断是否破坏了栈的单调性，如果不满足单调性则弹出栈顶元素直到单调性被满足，再将新点加入栈顶。最后答案即为线段的端点以及所有中垂线交点与他们最近点的距离的最大值。最后总的时间复杂度为O(n)。

code

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<string>
#include<algorithm>
#define maxn 1200005  
#define maxm 300005 
#define INF 0x7fffffff;

using namespace std;

double x[maxn],y[maxn];
int n,L;

int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}

double coordinate(int a,int b)
{
    double x1=(x[a]+x[b])*0.5,y1=(y[a]+y[b])*0.5;
    if (y[a]==y[b]) return x1;
    double k=(double)(-1)/(double)((double)(y[a]-y[b])/(double)(x[a]-x[b]));
    double B=(double)y1-(double)x1*k;
    return -(B/k);
}

double dis(double a,int b)
{
    return sqrt((a-x[b])*(a-x[b])+y[b]*y[b]);
}

int stack[maxn];
int top=0;
int start=0;

int main()
{
    freopen("mobile.in","r",stdin);
    freopen("mobile.out","w",stdout);
    scanf("%d%d",&n,&L);
    for (int i=1;i<=n;i++)
    {
        int xx,yy;
        read(xx),read(yy);
        x[i]=(double)xx; y[i]=(double)yy;
    }
    stack[++top]=1;
    for (int i=2;i<=n;i++)
    {
        if (x[i]==x[stack[top]]) 
            if (abs(y[i])>abs(y[stack[top]]))
                continue;
            else top--;
        while ((top>1)&&(coordinate(stack[top],stack[top-1])>coordinate(stack[top],i))) top--;
        stack[++top]=i;
    }
    start=1;
    for (int i=1;i<=top;i++)
    {
        if (coordinate(stack[i],stack[i+1])<0) start++;
        else break;
    }
    for (int i=top;i>=start;i--)
    {
        if (coordinate(stack[i],stack[i-1])>L) top--;
        else break;
    }
    double ans=dis(0,stack[start]);
    for (int i=start+1;i<=top;i++)
    {
        double k=coordinate(stack[i],stack[i-1]);
        double distance=dis(k,stack[i]);
        if (distance>ans) ans=distance;
    }
    ans=max(ans,dis(L,stack[top]));
    printf("%.6f",ans);
}