poj 3130 How I Mathematician Wonder What You Are! 半平面交判断核存在问题

http://poj.org/problem?id=3130

和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 }