首页 > 代码库 > USACO6.5-Closed Fences:计算几何

USACO6.5-Closed Fences:计算几何

Closed Fences

A closed fence in the plane is a set of non-crossing, connected linesegments with N corners (3 < N < 200). The corners or verticesare each distinct and are listed in counter-clockwise order in an array{xi, yi}, i in (1..N).

Every pair of adjacent vertices defines a side of the fence. Thus{xi yi xi+1 yi+1} is a sideof the fence for all i in (1..N). For our purposes, N+1 = 1, so thatthe first and last vertices making the fence closed.

Here is a typical closed fence and a point x,y:

                         * x3,y3                 x5,y5  /     x,y *          *   /                     / \ /                      /   *                  x6,y6*   x4,y4                     |                              |                          x1,y1*----------------* x2,y2

Write a program which will do the following:

  • Test an ordered list of vertices {xi,yi}, i in (1..N) to see if the array is a valid fence.
  • Find the set of fence sides that a person (with no height) who is standing in the plane at position (x,y) can "see" when looking at the fence. The location x,y may fall anywhere not on the fence.

A fence side can be seen if there exists a ray that connects (x,y) and any point on the side, and the ray does not intersect any other side of the fence. A side that is parallel to the line of sight is not considered visible. In the figure, above the segments x3,y3-x4,y4; x5,y5-x6,y6; and x6-y6-x1,y1 are visible or partially visible from x,y.

PROGRAM NAME: fence4

INPUT FORMAT

Line 1:N, the number of corners in the fence
Line 2:Two space-separated integers, x and y, that are the location of the observer. Both integers will fit into 16 bits.
Line 3-N+2:A pair of space-separated integers denoting the X,Y location of the corner. The pairs are given in counterclockwise order. Both integers are no larger than 1000 in magnitude.

NOTE: I have added anNew test case #12 for this task. Let me know if you think it‘s wrong. Rob Be sure to include USACO in your mail subject!

SAMPLE INPUT (file fence4.in)

135 50 07 05 27 55 73 54 91 82 50 9-2 70 3-3 1 

OUTPUT FORMAT

If the sequence is not a valid fence, the output is a single line containing the word "NOFENCE".

Otherwise, the output is a listing of visible fence segments, one per line, shown as four space-separated integers that represent the two corners. Express the points in the segment by showing first the point that is earlier in the input, then the point that is later. Sort the segments for output by examining the last point and showing first those points that are earlier in the input. Use the same rule on the first of the two points in case of ties.

SAMPLE OUTPUT (file fence4.out)

70 0 7 05 2 7 57 5 5 75 7 3 5-2 7 0 30 0 -3 10 3 -3 1


一道计算几何,需要细心做,特判几种情况,一种是点在某段篱笆的平行线上,这种情况是不能算看到的,这也是我第一次WA的原因。。。。T.T
一种是将端点沿着篱笆移动一段很小的距离,这样就避免了篱笆在端点处相交的情况。
还有就是按题目要求来排序。
这题要是没数据的话,我肯定找不到我错在哪的 o(>﹏<)o


Executing...
   Test 1: TEST OK [0.003 secs, 3508 KB]
   Test 2: TEST OK [0.003 secs, 3508 KB]
   Test 3: TEST OK [0.003 secs, 3508 KB]
   Test 4: TEST OK [0.005 secs, 3508 KB]
   Test 5: TEST OK [0.003 secs, 3508 KB]
   Test 6: TEST OK [0.014 secs, 3508 KB]
   Test 7: TEST OK [0.008 secs, 3508 KB]
   Test 8: TEST OK [0.008 secs, 3508 KB]
   Test 9: TEST OK [0.011 secs, 3508 KB]
   Test 10: TEST OK [0.008 secs, 3508 KB]
   Test 11: TEST OK [0.003 secs, 3508 KB]
   Test 12: TEST OK [0.003 secs, 3508 KB]

All tests OK.

  1 /*  2 TASK:fence4  3 LANG:C++  4 */  5   6 #include <iostream>  7 #include <vector>  8 #include <stdlib.h>  9 #include <stdio.h> 10 #include <cmath> 11 using namespace std; 12  13 //定义点 14 struct Point 15 { 16     double x; 17     double y; 18 }; 19 typedef Point Vector; 20 bool operator == (const Point& p1, const Point& p2) 21 { 22     return fabs(p1.x - p2.x)<1e-12 && fabs(p1.y - p2.y)<1e-12; 23 } 24 typedef struct Point point; 25 Vector operator - (const Point& A, const Point& B) 26 { 27     return Vector{A.x-B.x, A.y-B.y}; 28 } 29 double Cross(const Vector& A, const Vector& B) 30 { 31     return A.x*B.y - A.y*B.x; 32 } 33 //叉积 34 double multi(point p0, point p1, point p2) 35 { 36     return (p1.x - p0.x )*( p2.y - p0.y )-( p2.x -p0.x )*( p1.y -p0.y ); 37 } 38 //点积 39 double Dot(Vector A,Vector B) 40 { 41     return A.x*B.x+A.y*B.y; 42 } 43  44 //相交返回true,否则为false, 接口为两线段的端点 45 bool isIntersected(point s1,point e1, point s2,point e2) 46 { 47     if(s1==s2 || s1==e2 || e1==s2 || e1==e2) 48     { 49         return false; 50     } 51     return (max(s1.x,e1.x) >=min(s2.x,e2.x)) && 52            (max(s2.x,e2.x)>=min(s1.x,e1.x)) && 53            (max(s1.y,e1.y) >=min(s2.y,e2.y)) && 54            (max(s2.y,e2.y) >=min(s1.y,e1.y)) && 55            (multi(s1,s2,e1)*multi(s1,e1,e2)>0) && 56            (multi(s2,s1,e2)*multi(s2,e2,e1)>0); 57 } 58  59 point ps[300]; 60 double x,y; 61 int n ; 62 bool check(int idx) 63 { 64     point pos=point {x,y}; 65     // 特判pos与线段idx平行的情况 66     if(fabs(Cross(ps[idx]-pos,ps[idx+1]-pos))<1e-13) 67     { 68         return false; 69     } 70  71     bool flag=false; 72     for(int i=0; i<n; i++) 73     { 74         // 将ps[idx]移动一小段距离 75         if(i==idx) 76             continue; 77         double dx=(ps[idx+1].x-ps[idx].x)*1e-8; 78         double dy=(ps[idx+1].y-ps[idx].y)*1e-8; 79  80         if(isIntersected(pos,point {ps[idx].x+dx,ps[idx].y+dy},ps[i],ps[i+1])) 81             flag=true; 82     } 83     if(!flag) 84         return true; 85     for(int i=0; i<n; i++) 86     { 87         if(i==idx) 88             continue; 89  90  91         double dx=(ps[idx].x-ps[idx+1].x)*1e-8; 92         double dy=(ps[idx].y-ps[idx+1].y)*1e-8; 93  94         if(isIntersected(pos,point {ps[idx+1].x+dx,ps[idx+1].y+dy},ps[i],ps[i+1])) 95             return false; 96     } 97     return true; 98 } 99 100 101 int main()102 {103     freopen("fence4.in","r",stdin);104     freopen("fence4.out","w",stdout);105 106     cin>>n;107 108     scanf("%lf%lf",&x,&y);109     for(int i=0; i<n; i++)110     {111         scanf("%lf%lf",&ps[i].x,&ps[i].y);112     }113     ps[n]=ps[0];114 115     // 判断是否符合116     for(int i=1; i<n; i++)117     {118         for(int j=0; j<i-1; j++)119         {120             if(isIntersected(ps[j],ps[j+1],ps[i],ps[i+1]))121             {122                 puts("NOFENCE");123                 exit(0);124             }125         }126     }127 128     vector<int> ans;129     for(int i=0; i<n; i++)130     {131         if(check(i))132             ans.push_back(i);133     }134 135     int sz=ans.size();136     if(sz>=2 && ans[sz-1]==n-1 && ans[sz-2]==n-2)137         swap(ans[sz-1],ans[sz-2]);138     printf("%d\n",sz);139     for(int i=0; i<sz; i++)140     {141         if(ans[i]==n-1)142             printf("%.0f %.0f %.0f %.0f\n",ps[ans[i]+1].x,ps[ans[i]+1].y,ps[ans[i]].x,ps[ans[i]].y);143         else144             printf("%.0f %.0f %.0f %.0f\n",ps[ans[i]].x,ps[ans[i]].y,ps[ans[i]+1].x,ps[ans[i]+1].y);145     }146 147     return 0;148 }

 





USACO6.5-Closed Fences:计算几何