HDU 1086 判定直线相交

来源：互联网发布：js中字符串比较编辑：程序博客网时间：2024/06/03 19:37

主要用叉积来求，感觉计算几何用叉积的地方很多

计算几何这块，黑书上讲得比较好，可以看下，直线相交，有规范相交和非规范相交，求交点，求交点个数等类型

#include<iostream>                                                   #include<stdio.h>                                                    #include<math.h>                                                     struct point                                                         {                                                                    double x,y;                                                          }a[110];                                                                                                                                  double det(double x1,double y1,double x2,double y2)                  {                                                                         return x1*y2-x2*y1;                                             }                                                                                                                                         double cross(point a,point b,point c)                                {                                                                         return det(b.x-a.x,b.y-a.y,c.x-a.x,c.y-a.y);                    }                                                                                                                                         int dblcmp(double d)                                                 {                                                                        if(fabs(d)<1e-6)                                                     return 0;                                                            return d>0?1:-1;                                                 }                                                                                                                                         int seg(point a,point b,point c,point d)                             {                                                                       if(cross(a,c,d)*cross(b,c,d)>0)                                      return 0;                                                            if(cross(c,a,b)*cross(d,a,b)>0)                                      return 0;                                                            return 1;                                                         }                                                                                                                                         int main()                                                           {                                                                      int n;                                                               while(scanf("%d",&n)==1&&n)                                          {                                                                       int i,j;                                                             for(i=0;i<2*n;i+=2)                                                  {                                                                         scanf("%lf%lf%lf%lf",&a[i].x,&a[i].y,&a[i+1].x,&a[i+1].y);         }                                                                    int count=0;                                                         for(j=0;j<2*n;j+=2)                                                  {                                                                    for(i=j+2;i<2*n;i+=2)                                                {                                                                      if(seg(a[j],a[j+1],a[i],a[i+1]))                                     count++;                                                           }                                                                    }                                                                    printf("%d\n",count);                                             }                                                                    return 0;                                                          }