【计算几何】You can Solve a Geometry Problem too—hdoj 1086（判断线段相交）

You can Solve a Geometry Problem too—hdoj 1086（线段相交）

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

让你求相交点的个数，其实题目是简单的，就是考察了一个知识点，判断线段相交与否。判断线段相交的有很多种方法，最容易想到的就是求出点坐标，但代码实现起来有点麻烦。我学到了一个更简便的方法，在此分享。
网上搜到的其他方法，链接：http://blog.csdn.net/rickliuxiao/article/details/6259322
在此就介绍其上链接的第2种思路：
(1)首先判断x，y范围来确定其有没可能相交
快速排斥试验
设以线段 P1P2 为对角线的矩形为R， 设以线段 Q1Q2 为对角线的矩形为T，如果R和T不相交，显然两线段不会相交。
(2)其次如果两直线相交，则俩直线互相跨立对方
跨立试验
如果两线段相交，则两线段必然相互跨立对方。

如上图：如果向量P2Q1和P2Q2在向量P2P1的两侧即可证明Q1Q2 跨立 P1P2，
设向量P2Q1为向量a，向量P2Q2为向量b，向量P2P1为向量c，
如果 （a×c）*（c×b）>= 0，那么可以证明Q1Q2 跨立 P1P2。
但是不能证明P1P2 跨立 Q1Q2，因为P2P1是以向量计算，可以伸长和缩短就不符合跨立，
如果我们知道P1P2是直线，那么我们只需知道Q1Q2 跨立 P1P2就可以证明直线与线段相交，
但P1P2是线段的话，我们需要两线段分别跨立对方才可以判断其相交。

ac代码

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const int N = 1010;
const double esp = 1e-3;

struct point   //点坐标
{
    double x;
    double y;
};
struct V   //线段两点式
{
    point start;
    point end;
}v[N];
V change(V a)   //线段两点式起点化为 0，0
{
    V b;
    b.start.x = 0;
    b.start.y = 0;
    b.end.x = a.end.x - a.start.x;
    b.end.y = a.end.y - a.start.y;
    return b;
}
double cross_mul(V a,V b)
{
    double result;
    V aa = change(a);
    V bb = change(b);
    result = aa.end.x*bb.end.y - aa.end.y*bb.end.x;
    return result;
}
bool Across(V a,V b)
{
    V aa,bb,cc;
    aa.start = a.start;
    aa.end = b.start;
    bb.start = a.start;
    bb.end = b.end;
    cc = a;
    if(cross_mul(aa,cc) * cross_mul(cc,bb) < 0)
        return false;
    a = b;
    b = cc;
    aa.start = a.start;
    aa.end = b.start;
    bb.start = a.start;
    bb.end = b.end;
    cc = a;
    if(cross_mul(aa,cc) * cross_mul(cc,bb) < 0)
        return false;
    return true;
}

int main()
{
    int n;
    while(~scanf("%d",&n)&&n)
    {
        int sum = 0;
        for(int i=0;i<n;i++)
            scanf("%lf%lf%lf%lf",&v[i].start.x,&v[i].start.y,&v[i].end.x,&v[i].end.y);
        for(int i=0;i<n;i++)
        {
            for(int j=i+1;j<n;j++)
            {
                if(Across(v[i],v[j]))
                    sum++;
            }
        }
        printf("%d\n",sum);
    }
    return 0;
}