bzoj1038 [ZJOI2008]瞭望塔（半平面交）

bzoj1038 [ZJOI2008]瞭望塔

原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=1038

题意：
村中要建立一个瞭望塔，我们将H村抽象为一维的轮廓。如下图所示 我们可以用一条山的上方轮廓折线(x1, y1), (x2, y2), …. (xn, yn)来描述H村的形状，这里x1 < x2 < …< xn。瞭望塔可以建造在[x1, xn]间的任意位置, 但必须满足从瞭望塔的顶端可以看到H村的任意位置。
可见在不同的位置建造瞭望塔，所需要建造的高度是不同的。请你写一个程序，帮助村长计算塔的最小高度。

数据范围
N ≤ 300，输入坐标绝对值不超过1e6

题解：
首先，瞭望塔塔顶的的位置必须在半平面交交出的这个范围内。

但是考虑到地面是不平的，不能直接找最低点。

分段一次函数的极值只能在边界或分段点取到

于是枚举这些“折点”算距离。

注意村子的一些点可能在HPI交出的平面的第一个点左，最后一个点右，因此要首位先加入直向下和直向上，两个半平面。

代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#define Poi Vec
using namespace std;
const double EPS=1e-10;
const int N=1600;
struct Vec
{
    double x,y;
    Vec(){}
    Vec(double x,double y):x(x),y(y){}
    double ang(){return atan2(y,x);}
}P[N],pp[N];
Vec operator+(const Vec &A,const Vec &B){return Vec(A.x+B.x,A.y+B.y);}
Vec operator-(const Vec &A,const Vec &B){return Vec(A.x-B.x,A.y-B.y);}
Vec operator*(const Vec &A,double s){return Vec(A.x*s,A.y*s);}
Vec operator/(const Vec &A,double s){return Vec(A.x/s,A.y/s);}
double cross(const Vec &A,const Vec &B){return A.x*B.y-A.y*B.x;}
double dot(const Vec &A,const Vec &B){return A.x*B.x+A.y*B.y;}
struct Line
{
    Poi p;
    Vec u;
    double ang;
    Line(){}
    Line(Poi p,Vec u):p(p),u(u),ang(u.ang()){}
}L[N],Q[N];
bool cmp(const Line &A,const Line &B){return A.ang<B.ang;}
int n,m,top=0,h,t;
double ans;
int sign(double x)
{if(fabs(x)<=EPS) return 0; else return x<0?-1:1;}
bool onleft(Poi p,Line l) {return sign(cross(l.u,p-l.p))>=0;}
Poi intersig(Poi P,Vec u,Poi Q,Vec v)
{
    double t=cross(Q-P,v)/cross(u,v);
    return P+u*t;
}
void HPI()
{
    sort(L,L+n,cmp);
    h=0,t=0; Q[0]=L[0];
    for(int i=1;i<n;i++)
    {
        while(h<t&&!onleft(P[t-1],L[i])) t--;
        while(h<t&&!onleft(P[h],L[i])) h++;
        Q[++t]=L[i];
        if(sign(cross(Q[t].u,Q[t-1].u))==0)
        {t--; if(onleft(L[i].p,Q[t])) Q[t]=L[i];}
        if(h<t) P[t-1]=intersig(Q[t-1].p,Q[t-1].u,Q[t].p,Q[t].u);
    }
    while(h<t&&!onleft(P[t-1],Q[h])) t--; 
    top=0;
    for(int i=h;i<t;i++) 
    P[++top]=intersig(Q[i].p,Q[i].u,Q[i+1].p,Q[i+1].u);
}
void getans()
{
    ans=1e60; Vec vv=Vec(0,1);
    for(int i=1;i<=m;i++)
    {
        Poi v=Poi(pp[i].x,-1); 
        for(int j=1;j<top;j++)
        {
            if(P[j].x<=pp[i].x&&P[j+1].x>=pp[i].x)
            ans=min(ans,intersig(P[j],P[j+1]-P[j],v,vv).y-pp[i].y);
        }
    }
    for(int i=1;i<=top;i++)
    { 
        Poi v=Poi(P[i].x,-1); 
        for(int j=1;j<m;j++)
        {
            if(pp[j].x<=P[i].x&&P[i].x<=pp[j+1].x)
            ans=min(ans,P[i].y-intersig(v,P[i]-v,pp[j],pp[j+1]-pp[j]).y);
        }
    }
}
void init()
{
    pp[0].x=pp[1].x; pp[0].y=1000005;
    pp[n+1].x=pp[n].x; pp[n+1].y=1000005; 
    for(int i=0;i<=n;i++) L[i]=Line(pp[i],pp[i+1]-pp[i]); 
    n++;
}
int main()
{
    scanf("%d",&n); m=n;
    for(int i=1;i<=n;i++) scanf("%lf",&pp[i].x);
    for(int i=1;i<=n;i++) scanf("%lf",&pp[i].y);
    init();
    HPI();
    getans();
    printf("%0.3lf\n",ans);
    return 0;
}