hdu1086 You can Solve a Geometry Problem too【排斥实验+跨立实验】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1086
题意：给你n条线段，让你输出n条线段的交点个数，相同的交点也要计算
解析：由于不用求具体交点，所以可以用跨立实验和排斥实验判断是否有交点
排斥实验：
两条线段各自有一个矩形，若两个矩形相交，则这两条线有可能相交
跨立实验：
通过叉乘判断两条线段是否相交，也就是一条线段的两个端点对另一条线段的叉乘，如果异号则有可能相交

#include <cmath>#include <algorithm>#include <iostream>#include <cstdio>#include <vector>#include <cstring>using namespace std;const int maxn = 1e5;const double eps = 1e-5;struct point{    double x;    double y;    point() {}    point(double _x,double _y)    {        x = _x;        y = _y;    }};struct line{    point a;    point b;}a[maxn];double x_mul(point p0,point p1,point p2){    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);}bool judge(line p1,line p2){    //排斥实验    bool flag1 = max(p1.a.x,p1.b.x)>=min(p2.a.x,p2.b.x)&&                 max(p2.a.x,p2.b.x)>=min(p1.a.x,p1.b.x)&&                 max(p1.a.y,p1.b.y)>=min(p2.a.y,p2.b.y)&&                 max(p2.a.y,p2.b.y)>=min(p1.a.y,p1.b.y);    //跨立实验    bool flag2 = (x_mul(p1.a,p2.a,p2.b)*x_mul(p1.b,p2.a,p2.b)<=0)&&                 (x_mul(p2.a,p1.a,p1.b)*x_mul(p2.b,p1.a,p1.b)<=0);    return flag1&&flag2;}int main(void){    int n;    while(~scanf("%d",&n)&&n)    {        for(int i=0;i<n;i++)        {            double x1,y1,x2,y2;            scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);            a[i].a = point(x1,y1);            a[i].b = point(x2,y2);        }        int ans = 0;        for(int i=0;i<n;i++)        {            for(int j=i+1;j<n;j++)            {                if(judge(a[i],a[j]))                    ans++;            }        }        printf("%d\n",ans);    }    return 0;}