POJ 1039——计算几何初步（判断线段与直线相交)
此题太恶心了。。。。。（开始的时候思路错了，没有考虑到交点可能是拐点，这是其一。其二，就是精度。。。。。。。尼玛啊，这个害我WA了N次）
唉～～～～～精度很重要啊。
还有一点就是先判断枚举的直线是否合理，即是否在管道内。
然后对于每个拐点映射到直线上的Y值是否在管道里，不是则求出X值。
接下来枚举就行。
View Code
1 #include<iostream>
2 #include<cstdio>
3 #define EPS 1e-8
4 using namespace std;
5
6 struct point {
7 double x, y;
8 point(double _x, double _y):x(_x),y(_y){};
9 point(){};
10 point(point &temp)
11 {
12 x = temp.x;
13 y = temp.y;
14 }
15 }pipe[20][2];
16
17 double cross(point &p, point &q)
18 {
19 return q.y*p.x - q.x*p.y;
20 }
21
22 bool line_law(point &st, point &ed, int n)
23 {
24 double line_k, line_b;
25 line_k = (ed.y - st.y)/(ed.x - st.x);
26 line_b = ed.y - line_k*ed.x;
27 double x = pipe[0][0].x;
28 double y = line_k*x + line_b;
29 if (y < pipe[0][0].y + EPS && pipe[0][1].y - EPS <= y)return true;
30 return false;
31 }
32
33 bool check(point &st, point &ed, double &MAX, int n)
34 {
35 double line_k, line_b;
36 line_k = (ed.y - st.y)/(ed.x - st.x);
37 line_b = ed.y - line_k*ed.x;
38 int flag;
39 if (!line_law(st,ed,n))return false;
40 double MIN = 1000000;
41 for (int i(0); i<n; ++i) {
42 double x = pipe[i][0].x;
43 double y = line_k*x + line_b;
44 if (y < pipe[i][0].y + EPS && y >= pipe[i][1].y - EPS)continue;
45 double seg_k, seg_b;
46 seg_k = (pipe[i][0].y - pipe[i-1][0].y)/(pipe[i][0].x - pipe[i-1][0].x);
47 seg_b = pipe[i][0].y - pipe[i][0].x*seg_k;
48 if (y < pipe[i][1].y - EPS) {
49 seg_k = (pipe[i][1].y - pipe[i-1][1].y)/(pipe[i][1].x - pipe[i-1][1].x);
50 seg_b = pipe[i][1].y - pipe[i][1].x*seg_k;
51 }
52 MIN = (line_b - seg_b)/(seg_k - line_k);
53 if (MIN != 1000000)break;
54 }
55 if (MIN != 1000000) {
56 MAX = MAX > MIN ? MAX:MIN;
57 return false;
58 } else {
59 return true;
60 }
61 }
62
63 int main()
64 {
65 int t;
66 while (scanf("%d",&t) && t) {
67 double len = -1000000;
68 for (int i(0); i<t; ++i) {
69 scanf("%lf%lf",&pipe[i][0].x,&pipe[i][0].y);
70 pipe[i][1].x = pipe[i][0].x;
71 pipe[i][1].y = pipe[i][0].y - 1;
72 }
73 bool flag = false;
74 for (int i(0); i<t; ++i) {
75 for (int k(0); k<2; ++k) {
76 for (int j(0); j<t; ++j) {
77 for (int l(0); l<2; ++l) {
78 if (i != j) {
79 if (check(pipe[i][k],pipe[j][l],len,t)) {
80 flag = true;
81 }
82 }
83 if (flag)break;
84 }
85 if (flag)break;
86 }
87 if (flag)break;
88 }
89 if (flag)break;
90 }
91 if (flag) {
92 cout<<"Through all the pipe."<<endl;
93 } else {
94 printf("%.2lf\n",len);
95 }
96 }
97 return 0;
98 }
