poj 1474 Video Surveillance 半平面交
http://poj.org/problem?id=1474
和前面做过的两道一样，又是一题判断多边形是否存在核问题，照样粘贴模板。
点是顺时针给出。
1 #include<stdio.h> 2 #include<stdlib.h> 3 #include<cmath> 4 #include<algorithm> 5 using namespace std; 6 const double MAX =100000000; 7 const double pi =acos(-1.0); 8 const double eps=1e-8; 9 int m,s; 10 struct node 11 { 12 double x,y; //注意类型 13 }tr[110],p[110],q[110]; 14 int sig(double k) 15 { 16 return (k<-eps)?-1:(k>eps); 17 } 18 void interect(node x,node y,double a,double b,double c) 19 { 20 double u=fabs(a*x.x+b*x.y+c); 21 double v=fabs(a*y.x+b*y.y+c); 22 q[++s].x=(x.x*v+y.x*u)/(u+v); 23 q[s].y=(x.y*v+y.y*u)/(u+v); 24 } 25 //利用半平面切割 26 void cut(double a,double b,double c) 27 { 28 s=0; 29 int i; 30 for(i=1;i<=m;i++) //遍历所有顶点是否能观察到该边 31 { 32 if(sig(a*p[i].x+b*p[i].y+c)>=0)//因为线段是顺时针给出的，如果是逆时针就是<=0 33 { 34 q[++s]=p[i]; //若是则存储 35 } 36 else 37 { 38 if(sig(a*p[i-1].x+b*p[i-1].y+c)>0)//逆时针就是<0 39 interect(p[i-1],p[i],a,b,c); 40 if(sig(a*p[i+1].x+b*p[i+1].y+c)>0)//逆时针就是<0 41 interect(p[i+1],p[i],a,b,c); 42 } 43 } 44 //最后的p数组存放半平面的点集合 45 for(i=1;i<=s;i++) 46 p[i]=q[i]; 47 p[s+1]=p[1],p[0]=p[s]; 48 m=s; 49 } 50 51 52 int main() 53 { 54 int n,i,j,t=1; 55 while(scanf("%d",&n),n) 56 { 57 for(i=0;i<n;i++) 58 { 59 scanf("%lf%lf",&tr[i].x,&tr[i].y); 60 p[i+1]=tr[i]; //初始化边界 61 } 62 tr[n]=tr[0]; 63 p[n+1]=p[1];p[0]=p[n]; 64 m=n; 65 double a,b,c; 66 for(i=0;i<n;i++) 67 { 68 a=tr[i+1].y-tr[i].y; //计算出相邻两点所在直线ax+by+c=0 69 b=tr[i].x-tr[i+1].x; 70 c=tr[i+1].x*tr[i].y-tr[i].x*tr[i+1].y; 71 cut(a,b,c); 72 } 73 printf("Floor #%d\n",t++); 74 if(!m) puts("Surveillance is impossible.");//这里如果有一个点，或者一条线段都可以，所以判断m是不是等于0就行了，不用判断面积 75 else 76 puts("Surveillance is possible."); 77 printf("\n"); 78 } 79 return 0; 80 }
