挺好的一道计算几何题目

挺好的一道计算几何题目
题目大意是：给一条线段代表房子，给一条线段代表路，给一些障碍物，求在路上能完全看到房子的最长连续长度
题目中所有线段都是和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
;
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