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