usaco 3.4 Closed Fences 计算几何

题目给定一个多边形,要求

1)判断是否简单多边形

方法如下

 

对于临边：假设p1,p2和p2,p3为临边,如果p1,p2,p3重合并且p1,p3在p2的同侧,在存在交叉边,可以通过叉积cross(p1,p2,p3)==0 和点击dot(p1,p2,p3) < 0来判断.

 

对于不相邻的边：如果规范相交，则存在交叉边.

 

2)判断观察者是否能看到，题目的意思是只要能看到一点就满足条件，但是如果线段和观察者共线，则不能看到

 

方法:对于每个顶点p,将其分别做左右微小移动得到点pl,pr, 假设观察者为o,作射线(o,pl)和(o,pr),对于两射线分别求出与其相交并且举例最近的线段,就是能找到观察到的线段

 

 

3)输出要求先按第二点排序,再第一点,实际上,只需要把最后一条线段的两点交换顺序,再把最后两条线段交换顺序,输出即可

 

1. #include <iostream>
2. #include <cmath>
3. #include <algorithm>
4. using namespace std;
6. /*
7. PROG: fence4
8. LANG: C++
9. ID: heben991
10. */
11. const int N = 500;
12. typedef double T;
14. const T inf = 1e15, eps = 1e-5, pi = acos(-1.0), C = 1e8;
17. struct point
18. {
19.     T x, y;
20.     point(T a=0, T b=0)
21.     {
22.         x=a,y=b;
23.     }
24. }p[N], o;
26. int n;
27. bool see[N];
30. T dot(point a, point b, point c)
31. {
32.     return (b.x-a.x)*(c.x-a.x) + (b.y-a.y)*(c.y-a.y);
33. }
35. T cross(point a, point b)
36. {
37.     return a.x*b.y - a.y*b.x;
38. }
40. T cross(point a, point b, point c)
41. {
42.     return (b.x-a.x)*(c.y-a.y) - (c.x-a.x)*(b.y-a.y);
43. }
45. int sign(T x)
46. {
47.     if( x == 0 ) return 0;
48.     if( x > 0 )return 1;
49.     return -1;
50. }
52. bool checkssint(point a, point b, point c, point d)
53. {
54.     return max(a.x,b.x) > min(c.x,d.x)
55.     && max(a.y,b.y) > min(c.y,d.y)
56.     && max(c.x,d.x) > min(a.x,b.x)
57.     && max(c.y,d.y) > min(a.y,b.y)
58.     && sign(cross(a,b,c))*sign(cross(a,b,d)) < 0
59.     && sign(cross(c,d,a))*sign(cross(c,d,b)) < 0;
60. }
62. bool segsegint(point p1, point p2, point p3, point p4, point &p)
63. {
64.     if( !checkssint(p1,p2,p3,p4) ) return 0;
66.     double d, d1, d2;
67.     d = (p1.x-p2.x)*(p4.y-p3.y) - (p1.y-p2.y)*(p4.x-p3.x);
68.     d1 = cross(p3,p4)*(p1.x-p2.x) - cross(p1,p2)*(p3.x-p4.x);
69.     d2 = cross(p1,p2)*(p4.y-p3.y) - cross(p3,p4)*(p2.y-p1.y);
70.     p.x = d1/d;
71.     p.y = d2/d;
72.     return 1;
73. }
75. double dist2(point a, point b)
76. {
77.     return (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y);
78. }
80. void nearest(T xx, T yy)
81. {
82.     point pt(xx,yy), inter;
83.     double dt, dm = inf;
84.     int near = -1, i, j, k;
86.     pt.x += (pt.x-o.x)*C;
87.     pt.y += (pt.y-o.y)*C;
89.     for(i = 0; i < n; ++i)
90.     {
91.         if(!segsegint(o,pt,p[i],p[i+1],inter))continue;
92.         dt = dist2(inter, o);
93.         if(dt < dm)
94.         {
95.             dm = dt;
96.             near = i;
97.         }
98.     }
100.     if(near != -1) see[near] = 1;
101. }
104. int main()
105. {
106.     int i, j, k;
108.     freopen("fence4.in", "r", stdin);
109.     freopen("fence4.out","w",stdout);
111.     scanf("%d%lf%lf", &n, &o.x, &o.y);
112.     for(i = 0; i < n; ++i) scanf("%lf%lf", &p[i].x, &p[i].y);
113.     p[n] = p[0];
114.     p[n+1] = p[1];
116.     for(i = 0; i < n; ++i)
117.     {
118.         if( fabs( cross(p[i],p[i+1],p[i+2]) ) <= eps
119.             && dot(p[i],p[i+1],p[i+2]) < 0 )
120.         {
121.             printf("---%d/n", i);
122.             break;
123.         }
125.         for(j = i+2; j < n; ++j)
126.         if((j+1)%n != i && checkssint(p[i],p[i+1],p[j],p[j+1]))
127.         {
128.             printf("------%d %d/n", i, j);
129.             break;
130.         }
132.         if(j < n)break;
133.     }
135.     if(i < n)
136.     {
137.         puts("NOFENCE");
138.         return 0;
139.     }
141.     for(i = 0; i < n; ++i)
142.     {
143.         nearest(p[i].x-10.0*eps, p[i].y-pi*eps);
144.         nearest(p[i].x+10.0*eps, p[i].y+10*eps);
145.     }
147.     int cnt = 0;
148.     for(i = 0; i < n; ++i) if(see[i]) ++cnt;
150.     printf("%d/n", cnt);
152.     for(i = 0; i < n-2; ++i)
153.     if(see[i])
154.     {
155.         printf("%.0lf %.0lf %.0lf %.0lf/n",p[i].x,p[i].y,p[i+1].x,p[i+1].y);
156.     }
158.     if(see[i=n-1])
159.     {
160.         printf("%.0lf %.0lf %.0lf %.0lf/n",p[i+1].x,p[i+1].y,p[i].x,p[i].y);
161.     }
163.     if(see[i=n-2])
164.     {
165.         printf("%.0lf %.0lf %.0lf %.0lf/n",p[i].x,p[i].y,p[i+1].x,p[i+1].y);
166.     }
167.     return 0;
168. }