POJ 1039 Pipe

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

题意：是否存在一条光线通过弯曲的管道，不行的话输出能到达最远的横坐标

思路：枚举每两个点，然后判断

 

#include<cstdio>

#include<cstring>

#include<cmath>

#include<queue>

#include<algorithm>

const int maxn = 1e2 + 10;

const double eps = 1e-8;

const double INF = (double)1e9;

using namespace std;

//加法误差

double add(double a, double b) {

    if(fabs(a + b) < eps) return 0;

    return a + b;

}

//向量运算

struct P {

    double x, y;

    P() {}

    P(double x, double y) : x(x), y(y) {}

    void input() {

        scanf("%lf %lf", &x, &y);

    }

    P operator + (P p) {

        return P(add(x, p.x), add(y, p.y));

    }

    P operator - (P p) {

        return P(add(x, -p.x), add(y, -p.y));

    }

    P operator * (double d) {

        return P(x * d, y * d);

    }

    //点乘

    bool operator != (P p) {

        return add(x, -p.x) || add(y, -p.y);

    }

    double dot(P p) {

        return add(x * p.x, y * p.y);

    }

    //叉乘

    double det(P p) {

        return add(x * p.y, -y * p.x);

    }

    //是否平行

    bool ispar(P q) {

        P p(x, y);

        return p.det(q) == 0;

    }

    //点是否在线段上

    bool onseg(P p1, P p2) {

        P p(x, y);

        return (p - p1).det(p - p2) == 0 && (p - p1).dot(p - p2) <= 0;

    }

};

//直线p1-p2 和 q1-q2 交点(有时需要特判线段共线)

P in(P p1, P p2, P q1, P q2) {

    return p1 + (p2 - p1) * ((q2 - q1).det(q1 - p1) / (q2 - q1).det(p2 - p1));

}

 

struct Seg {

    P p, q;

    Seg() {}

    Seg(P p, P q) : p(p), q(q) {}

};

Seg line[4], seg[maxn];

P f[5];

 

int main() {

    int T, n;

    while(scanf("%d", &n) && n) {

        for(int i = 0; i < n; i++) {

            seg[i].p.input();

            seg[i].q.x = seg[i].p.x; seg[i].q.y = seg[i].p.y - 1;

        }

        int sign = 1;

        double maxx = -INF;

        for(int i = 0; i < n - 1; i++) {

            for(int j = i + 1; j < n; j++) {

                //四个方向 + 四条直线

                line[0] = Seg(seg[j].p, seg[i].p);

                line[1] = Seg(seg[j].p, seg[i].q);

                line[2] = Seg(seg[j].q, seg[i].p);

                line[3] = Seg(seg[j].q, seg[i].q);

                for(int k = 0; k < 4; k++) {

                    sign = 1;

                    for(int ind = 0; ind < n; ind++) {

                        if(ind == i || ind == j) continue;

                        P ja = in(seg[ind].p, seg[ind].q, line[k].p, line[k].q);

                        double lf = (seg[ind].p - ja).dot(seg[ind].q - ja);

                        if(lf > 0) {

                            if(ind) {

                                P p1 = in(seg[ind].p, seg[ind - 1].p, line[k].p, line[k].q);

                                P p2 = in(seg[ind].q, seg[ind - 1].q, line[k].p, line[k].q);

                                if(p1.onseg(seg[ind].p, seg[ind - 1].p)) maxx = max(maxx, p1.x);

                                if(p2.onseg(seg[ind].q, seg[ind - 1].q)) maxx = max(maxx, p2.x);

                            } else maxx = max(maxx, seg[ind].p.x);

                            sign = 0; break;

                        }

                    }

                    if(sign) break;

                }

                if(sign) break;

            }

            if(sign) break;

        }

        if(sign) printf("Through all the pipe.\n");

        else printf("%.2f\n", maxx);

    }

    return 0;

}