poj 2074 Line of Sight 计算几何

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

 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[];
     while(cin>>x1>>x2>>y){
         if(x1==&&x2==&&y==)
             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=;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=;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==){     //若没有障碍物在建筑物与街道之间，，则输出整个街道的长度
             printf("%.2lf\n",pro.e.x-pro.s.x);
             continue;
         }
         sort(a,a+n,cmp);
         double l[],r[];
         l[] = a[].l,r[] = a[].r;
         for(j=,i=;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[]-pro.s.x,pro.e.x-r[j]);          //寻找最大值。
         for(i=;i<j;i++){
             ans = max(ans,l[i+]-r[i]);
         }
         if(ans<eps)
             printf("No View\n");
         else
             printf("%.2lf\n",ans);
     }
     return ;
 }