poj 2074 Line of Sight 计算几何

/**
大意：给定一个建筑--水平放置，给定n个障碍物， 给定一条街道，从街道上能看到整个建筑的最长的连续的区域
思路： 分别确定每一个障碍物所确立的盲区，即----建筑物的终点与障碍物的起点的连线，建筑物的起点与障碍物的终点的连线。。这段区域即为盲区，，，有多个盲区，需要去重。。
学习之处： 1、障碍物在输入时去重，非常巧妙
                   2、 盲区的去重。与重组，。。巧妙。。
                    3、 寻找最大连续区域，，
**/
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
const double eps = 1e-8;
double  max(double  a,double b){
    return a>b?a:b;
}

struct point{
    double x,y;
};

struct segment{
    point s,e;
};

struct line{
    double a,b,c;
};

struct data{
    double l,r;
}a[100];

line getline(point p1,point p2){
    line tmp;
    tmp.a = p1.y-p2.y;
    tmp.b = p2.x-p1.x;
    tmp.c = p1.x*p2.y-p2.x*p1.y;
    return tmp;
}

double getInterx(line l1,line l2){
    return (l1.b*l2.c-l2.b*l1.c)/(l1.a*l2.b-l2.a*l1.b);
}

bool cmp(data a,data b){
    return a.l<b.l;
}

int main()
{
    int n,i,j;
    double x1,x2,y;
    segment hou,pro,obs[100];
    while(cin>>x1>>x2>>y){
        if(x1==0&&x2==0&&y==0)
            break;
        hou.s.x = x1,hou.s.y = y;
        hou.e.x = x2,hou.e.y = y;
        cin>>x1>>x2>>y;
        pro.s.x = x1,pro.s.y = y;
        pro.e.x = x2,pro.e.y = y;
        cin>>n;
        for(i=0;i<n;i++){
            cin>>x1>>x2>>y;
            obs[i].s.x = x1,obs[i].s.y = y;
            obs[i].e.x = x2,obs[i].e.y =y;
            if(y>hou.s.y||y<pro.s.y){
                n--;
                i--;
            }
        }
        line l_pro = getline(pro.s,pro.e);
        for(i=j=0;i<n;i++,j++){
            line l1,l2;
            l1 = getline(hou.s,obs[i].e);
            l2 = getline(hou.e,obs[i].s);
            a[j].r = getInterx(l_pro,l1);
            a[j].l = getInterx(l_pro,l2);
            if(a[j].r<pro.s.x||a[j].l>pro.e.x){
                j--;
                continue;
            }
            if(a[j].l<pro.s.x)     //对于盲区超出街道的区域。
                a[j].l = pro.s.x;
            if(a[j].r>pro.e.x)
                a[j].r = pro.e.x;
        }
        n = j;
        if(n==0){     //若没有障碍物在建筑物与街道之间，，则输出整个街道的长度
            printf("%.2lf
",pro.e.x-pro.s.x);
            continue;
        }
        sort(a,a+n,cmp);
        double l[100],r[100];
        l[0] = a[0].l,r[0] = a[0].r;
        for(j=0,i=1;i<n;i++){     //盲区的重组与去重
            if(r[j]<a[i].l){
                j++;
                l[j] = a[i].l;
                r[j] = a[i].r;
            }else
            r[j] = max(r[j],a[i].r);
        }
        double ans = max(l[0]-pro.s.x,pro.e.x-r[j]);          //寻找最大值。
        for(i=0;i<j;i++){
            ans = max(ans,l[i+1]-r[i]);
        }
        if(ans<eps)
            printf("No View
");
        else
            printf("%.2lf
",ans);
    }
    return 0;
}