bzoj1038: [ZJOI2008]瞭望塔

bzoj1038: [ZJOI2008]瞭望塔

Description

　　致力于建设全国示范和谐小村庄的H村村长dadzhi，决定在村中建立一个瞭望塔，以此加强村中的治安。我们
将H村抽象为一维的轮廓。如下图所示 我们可以用一条山的上方轮廓折线(x1, y1), (x2, y2), …. (xn, yn)来描
述H村的形状，这里x1 < x2 < …< xn。瞭望塔可以建造在[x1, xn]间的任意位置, 但必须满足从瞭望塔的顶端可
以看到H村的任意位置。可见在不同的位置建造瞭望塔，所需要建造的高度是不同的。为了节省开支，dadzhi村长
希望建造的塔高度尽可能小。请你写一个程序，帮助dadzhi村长计算塔的最小高度。

Input

　　第一行包含一个整数n，表示轮廓折线的节点数目。接下来第一行n个整数, 为x1 ~ xn. 第三行n个整数，为y1
~ yn。

Output

　　仅包含一个实数，为塔的最小高度，精确到小数点后三位。

Sample Input

【输入样例一】
6
1 2 4 5 6 7
1 2 2 4 2 1
【输入样例二】
4
10 20 49 59
0 10 10 0

Sample Output

【输出样例一】
1.000
【输出样例二】
14.500

HINT

N ≤ 300，输入坐标绝对值不超过106，注意考虑实数误差带来的问题。

分析

主要是理解题意
瞭望塔是建在那个轮廓上的，并且要看到每个山顶。
其实就是在轮廓直线半平面交上一点到轮廓最短距离
距离是一堆线性的分段函数。
所以极值一定在分界点上取到。
所以把半平面交上的点投影到山轮廓，再把山轮廓上的点投影到半平面交上即可。
哦还有这是一道半半平面交的题
具体做法参照：bzoj1007: [HNOI2008]水平可见直线 关于一类半半平面交题目的解法研究

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
#define inf 1000000000
#define ll long long
using namespace std;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
double ans=1e60;
int n,cnt,top,tot;
struct P{double x,y;}p[1005],a[1005];
struct L{P a,b;double slop;}l[1005],q[1005];
P operator-(P a,P b)
{
    P t;t.x=a.x-b.x;t.y=a.y-b.y;return t;
}
double operator*(P a,P b)
{
    return a.x*b.y-a.y*b.x;
}
inline bool operator<(L a,L b)
{
    if(a.slop!=b.slop)return a.slop<b.slop;
    else return (a.b-a.a)*(b.b-a.a)>0;
}
P inter(L a,L b)
{
    double k1,k2,t;
    k1=(b.b-a.a)*(a.b-a.a);
    k2=(a.b-a.a)*(b.a-a.a);
    t=k1/(k1+k2);
    P ans;
    ans.x=b.b.x+(b.a.x-b.b.x)*t;
    ans.y=b.b.y+(b.a.y-b.b.y)*t;
    return ans;
}
bool jud(L a,L b,L t)
{
    P p=inter(a,b);
    return (t.b-t.a)*(p-t.a)<0;
}
void hpi()
{
    int L=1,R=0;tot=0;
    for(int i=1;i<=cnt;i++)
    {
        if(l[i].slop!=l[i-1].slop)tot++;
        l[tot]=l[i];
    }
    cnt=tot;
    q[++R]=l[1];q[++R]=l[2];
    for(int i=3;i<=cnt;i++)
    {
        while(L<R&&jud(q[R-1],q[R],l[i]))R--;
        while(L<R&&jud(q[L+1],q[L],l[i]))L++;
        q[++R]=l[i];
    }
    while(L<R&&jud(q[R-1],q[R],q[L]))R--;
    while(L<R&&jud(q[L+1],q[L],q[R]))L++;
    tot=0;
    for(int i=L;i<R;i++)
        a[++tot]=inter(q[i],q[i+1]);
}
void pre()
{
    p[0].x=p[1].x;p[0].y=100001;
    p[n+1].x=p[n].x;p[n+1].y=100001;
    for(int i=1;i<=n;i++)
    {
        l[++cnt].a=p[i-1];l[cnt].b=p[i];
        l[++cnt].a=p[i];l[cnt].b=p[i+1];
    }
    for(int i=1;i<=cnt;i++)
    {
        l[i].slop=atan2(l[i].b.y-l[i].a.y,l[i].b.x-l[i].a.x);
    }
    sort(l+1,l+cnt+1);
}
void getans()
{
    for(int k=1;k<=tot;k++)
        for(int i=1;i<n;i++)
        {
            P t;t.x=a[k].x;t.y=-1;
            if(a[k].x>=p[i].x&&a[k].x<=p[i+1].x)
                ans=min(ans,a[k].y-inter((L){p[i],p[i+1]},(L){t,a[k]}).y);
        }
    for(int k=1;k<=n;k++)
        for(int i=1;i<tot;i++)
        {
            P t;t.x=p[k].x;t.y=-1;
            if(p[k].x>=a[i].x&&p[k].x<=a[i+1].x)
                ans=min(ans,inter((L){a[i],a[i+1]},(L){t,p[k]}).y-p[k].y);
        }
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
        p[i].x=read();
    for(int i=1;i<=n;i++)
        p[i].y=read();
    pre();
    hpi();
    getans();
    printf("%.3lf",ans);
    return 0;
}