poj2074 Line of Sight【计算几何】

题目大意：

题意：给出一个房子（看成线段）的端点坐标，和一条观光区的两端坐标，给出一些障碍物（看成线段）的两端坐标（都平行于x轴）。问在路上能看到完整房子的最大连续长度是多长。
【注意】
1. 输入的线段(x1,x2,y)保证x1 < <script type="math/tex" id="MathJax-Element-8"><</script>x2
2. 输入的home一定在观光区上面
3. 输入的障碍物不一定在home和观光区之间
4. 如果输入的障碍物与home同y轴，不管重合与部分重合还是不重合，不算遮住。
5. 如果输入的障碍物与观光区重合，同理。

解题思路：

求出每个障碍在观光区上所阻挡的视野区间，然后去掉被覆盖的小区间求相邻两个区间之间的间隙最大值即可。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
#include<ctime>
#include<vector>
#include<queue>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
    if(c=='-')c=getchar(),f=-1;
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=1005;
const double eps=1e-8;
int n;
double ans;
struct node
{
    double x1,x2,y;
    inline friend bool operator <(const node &a,const node &b)
    {
        if(a.x1==b.x1)return a.x2>b.x2;
        return a.x1<b.x1;
    }
}a[N],h,p;

double calc(double x1,double y,double x2)
{
    if(fabs(x1-x2)<eps)return x1;
    double k=(h.y-y)*1.0/(x2-x1);
    double m=y-k*x1;
    double res=(p.y-m)*1.0/k;
    return res;
}

int main()
{
    //freopen("lx.in","r",stdin);
    double x1,x2,y;
    while(scanf("%lf%lf%lf",&h.x1,&h.x2,&h.y))
    {
        if(!h.x1&&!h.x2&&!h.y)break;
        ans=0;
        scanf("%lf%lf%lf",&p.x1,&p.x2,&p.y);
        n=getint();
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf",&x1,&x2,&y);
            if(y>=h.y||y<=p.y)
            {
                i--,n--;
                continue;
            }
            a[i].x1=calc(x1,y,h.x2);
            a[i].x2=calc(x2,y,h.x1);
        }
        ++n,a[n].x1=a[n].x2=p.x1;
        ++n,a[n].x1=a[n].x2=p.x2;
        sort(a+1,a+n+1);
        int top=1;
        for(int i=2;i<=n;i++)
            if(a[i].x2>a[top].x2)a[++top]=a[i];
        n=top;
        int i=1;
        while(a[i].x2<p.x1)i++;
        for(i;i<n;i++)ans=max(ans,a[i+1].x1-a[i].x2);
        if(!ans)puts("No View");
        else printf("%0.2f\n",ans);
    }
    return 0;
}