【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，注意考虑实数误差带来的问题。

题解:

         求出半平面交之后答案肯定是半平面交形成的函数折线与轮廓折线之间的最小距离.

         然后这个距离是一个分段函数.分段函数的最值在折点处取到.

         所以枚举上方的每个点和下方的每个点分别计算距离即可.

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 310
#define inf 210000000
using namespace std;
struct point{double x,y;}p[N],a[N];
struct line{
  point a,b;
  double ag;
}li[N],q[N];
int n,cnt,num,tot;
double ans;
int read(){
  int x(0);char ch=getchar();
  while (ch<'0'||ch>'9') ch=getchar();
  while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
  return x;
}
point operator-(point a,point b){
  point ans;
  ans.x=a.x-b.x;
  ans.y=a.y-b.y;
  return ans;
}
double operator*(point a,point b){
  double ans;
  ans=a.x*b.y-a.y*b.x;
  return ans;
}
bool cmp(line a,line b){
  if (a.ag==b.ag) return (a.b-a.a)*(b.b-a.a)>0;
  else return a.ag<b.ag;
}
void getline(){
  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+1;i++)
    li[++num].a=p[i-1],li[num].b=p[i];
  for (int i=1;i<=num;i++)
    li[i].ag=atan2(li[i].b.y-li[i].a.y,li[i].b.x-li[i].a.x);
  sort(li+1,li+num+1,cmp);
  for (int i=1;i<=num;i++){
    if (li[i].ag!=li[i-1].ag) cnt++;
    li[cnt]=li[i];
  }
}
point inter(line a,line b){
  point ans; 
  double p1=(b.b-a.a)*(a.b-a.a);
  double p2=(a.b-a.a)*(b.a-a.a);
  double p=p1/(p1+p2);
  ans.x=b.b.x+(b.a.x-b.b.x)*p;
  ans.y=b.b.y+(b.a.y-b.b.y)*p;
  return ans;
}
bool check(line a,line b,line c){
  point t=inter(a,b);
  return (c.b-c.a)*(t-c.a)<0;
}
void hpi(){
  int l=1,r=0;
  q[++r]=li[1];q[++r]=li[2];
  for (int i=3;i<=cnt;i++){
    while(l<r&&check(q[r-1],q[r],li[i])) r--;
    while(l<r&&check(q[l+1],q[l],li[i])) l++;
    q[++r]=li[i]; 
  }
  while(l<r&&check(q[r-1],q[r],q[l])) r--;
  while(l<r&&check(q[l+1],q[l],q[r])) l++; 
  for (int i=l;i<r;i++)
    a[++tot]=inter(q[i],q[i+1]);
}
void getans(){
  ans=1e60;
  for (int i=1;i<=tot;i++)
    for (int j=1;j<n;j++){
      point t;t.x=a[i].x;t.y=-1;
      if (p[j].x<=t.x&&p[j+1].x>=t.x)
        ans=min(ans,a[i].y-inter(line{p[j],p[j+1]},line{t,a[i]}).y);
    }
  for (int i=1;i<=n;i++)
    for (int j=1;j<tot;j++){
      point t;t.x=p[i].x;t.y=-1;
      if (a[j].x<=t.x&&a[j+1].x>=t.x)
        ans=min(ans,inter(line{a[j],a[j+1]},line{t,p[i]}).y-p[i].y);
    }  
}
int main(){
  //freopen("a.in","r",stdin);
  n=read();
  for (int i=1;i<=n;i++) p[i].x=read();
  for (int i=1;i<=n;i++) p[i].y=read();
  getline();    
  hpi();
  getans();
  printf("%.3lf\n",ans);
}