挺好的一道计算几何题目

挺好的一道计算几何题目
题目大意是：给一条线段代表房子，给一条线段代表路，给一些障碍物，求在路上能完全看到房子的最长连续长度
题目中所有线段都是和x轴平行的
两个难点
1、利用相似三角形，求出再房子和障碍物连线在路上的交点

2、判断线段的相交方式，并切割线段

这些问题都解决了之后，就是写代码了

  1 /* 
  2  * SOUR:pku 2074
  3  * ALGO:computional geometry
  4  * DATE: Thu, 15 Oct 2009 23:22:48 +0800
  5  * COMM:3
  6  * */
  7 #include<iostream>
  8 #include<cstdio>
  9 #include<cstdlib>
 10 #include<cstring>
 11 #include<algorithm>
 12 #include<vector>
 13 #include<cassert>
 14 #include<cmath>
 15 using namespace std;
 16 typedef long long LL;
 17 const int maxint = 0x7fffffff;
 18 const long long max64 = 0x7fffffffffffffffll;
 19 template<class T>void show(T a, int n){for(int i=0;i<n;++i)cout<<a[i]<<' ';cout<<endl;}
 20 template<class T>void show(T a,int r,int l){for(int i=0;i<r;++i)show(a[i],l);cout<<endl;}
 21 #define pr(x) fprintf(stderr, x)
 22 /* #define pr(x) for(;0;) */
 23 const int N = 512;
 24 const double eps = 1e-7;
 25 struct NODE {
 26     double x1,x2,y;
 27     NODE(){}
 28     NODE(double a,double b){ x1 = a,x2 = b; }
 29 }h,ob,root;
 30 vector<NODE> g[N];
 31 int n,top,pt;
 32 
 33 
 34 struct point_t{
 35     double x,y;
 36     point_t(){}
 37     point_t(double a,double b){
 38         x = a,y = b;
 39     }
 40 };
 41 point_t operator +(point_t a,point_t b) { return point_t(a.x + b.x,a.y + b.y); }
 42 point_t operator -(point_t a,point_t b) { return point_t(a.x - b.x,a.y - b.y); }
 43 point_t operator *(point_t a,double b) { return point_t(a.x * b,a.y * b); }
 44 point_t operator /(point_t a,double b) { return point_t(a.x / b,a.y / b); }
 45 double sqr(double x) {return x * x;}
 46 double dist(point_t a) { return sqrt(sqr(a.x) + sqr(a.y)); }
 47 double dist(point_t a,point_t b) { return dist(a-b); }
 48 
 49 void cut(int idx,double v1,double v2)
 50 {
 51     int i,j,k;
 52     for(i = 0;i < g[idx-1].size();i++) {
 53         if(v1 <= g[idx-1][i].x1 && v2 >= g[idx-1][i].x2) {
 54         }else if(v1 <= g[idx-1][i].x1 && v2 > g[idx-1][i].x1 && v2 <= g[idx-1][i].x2) {
 55             g[idx].push_back(NODE(v2,g[idx-1][i].x2));
 56         }else if(v1 >= g[idx-1][i].x1 && v1 < g[idx-1][i].x2 && v2 >= g[idx-1][i].x2) {
 57             g[idx].push_back(NODE(g[idx-1][i].x1,v1));
 58         }else if(v1 >= g[idx-1][i].x1 && v2 <= g[idx-1][i].x2) {
 59             g[idx].push_back(NODE(g[idx-1][i].x1,v1));
 60             g[idx].push_back(NODE(v2,g[idx-1][i].x2));
 61         }else {
 62             g[idx].push_back(g[idx-1][i]);
 63         }
 64     }
 65 }
 66 
 67 
 68 point_t cacu(point_t a,point_t b)
 69     //b为ob,a为h
 70 {
 71     assert(h.y - root.y >= 0);
 72     assert(h.y - ob.y >= 0);
 73     return a + (b - a) * (h.y - root.y) / (h.y - ob.y);
 74 }
 75 
 76 int main()
 77 {
 78     int i,j,k;
 79     while(scanf("%lf%lf%lf",&h.x1,&h.x2,&h.y) && (h.x1 || h.x2 || h.y)) {
 80         scanf("%lf%lf%lf",&root.x1,&root.x2,&root.y);
 81         scanf("%d",&n);
 82         for(i = 0;i < N;i++) {
 83             g[i].clear();
 84         }
 85         g[0].push_back(root);
 86         for(i = 1,j = 0;i <= n;i++) {
 87             scanf("%lf%lf%lf",&ob.x1,&ob.x2,&ob.y);
 88             if(ob.y < h.y && ob.y > root.y) {
 89                 point_t v1 = cacu(point_t(h.x2,h.y) , point_t(ob.x1,ob.y));
 90                 point_t v2 = cacu(point_t(h.x1,h.y) , point_t(ob.x2,ob.y));
 91                 //printf("<%f,%f> <%f,%f>\n",v1.x,v1.y,v2.x,v2.y);
 92                 assert(v1.x <= v2.x && fabs(v1.y - v2.y) < eps);
 93                 cut(++j,v1.x,v2.x);
 94             }
 95         }
 96         double res = 0;
 97         for(i = 0;i < g[j].size();i++) {
 98             res = max(res,g[j][i].x2 - g[j][i].x1);
 99         }
100         if(res < eps) {
101             printf("No View\n");
102         }else {
103             printf("%.2f\n",res);
104         }
105     }
106     return 0;
107 }
108 
109