HDU 1086 You can Solve a Geometry Problem too(规范相交模版)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1086 

代码：

 1 #include <iostream>
 2 #include <cstdio>
 3 using namespace std;
 4 const double eps=1e-10;
 5 struct point
 6 {
 7     double x;
 8     double y;
 9 };
10 struct line
11 {
12     point a;
13     point b;
14 }p[101];
15 double cross(point a,point b,point c)
16 {
17     //向量ac在ab的方向，顺时针为正
18     return (c.x-a.x)*(b.y-a.y)-(b.x-a.x)*(c.y-a.y);
19 }
20 bool cr(line u,line v)
21 {
22     double k1=cross(u.a,u.b,v.a);
23     double k2=cross(u.a,u.b,v.b);
24     double k3=cross(v.a,v.b,u.a);
25     double k4=cross(v.a,v.b,u.b);
26     if(k1*k2<eps&&k3*k4<eps)
27         return true;
28     return false;
29 }
30 int main()
31 {
32     int n,i,j,sum;
33     while(~scanf("%d",&n))
34     {
35         sum=0;
36         if(!n)
37             break;
38         for(i=0;i<n;i++)
39         {
40             scanf("%lf%lf%lf%lf",&p[i].a.x,&p[i].a.y,&p[i].b.x,&p[i].b.y);
41         }
42         for(i=0;i<n-1;i++)
43             for(j=i+1;j<n;j++)
44                 if(cr(p[i],p[j]))
45                     sum++;
46         printf("%d\n",sum);
47     }
48     return 0;
49 }