HDU 1086 判断线段相交

/*
  判断线段相交利用叉乘，
  因为 两向量叉乘==两向量构成的平行四边形(以两向量为邻边)的面积 , 所以上面的公式也不难理解.
而且由于向量是有方向的, 所以面积也是有方向的, 通常我们以逆时针为正, 顺时针为负数.


改良算法关键点就是:
如果"线段ab和点c构成的三角形面积"与"线段ab和点d构成的三角形面积" 构成的三角形面积的正负符号相异,
那么点c和点d位于线段ab两侧.


注意面积等于0也是交点
*/




#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
using namespace std;
const int maxn=105;


struct node
{
    double x1,y1,x2,y2;
}p[maxn];




int slove(int i,int j)
{
    int flag1=0;
    int flag2=0;
    double acd=(p[j].x1-p[i].x1)*(p[j].y2-p[i].y1)-(p[j].x2-p[i].x1)*(p[j].y1-p[i].y1);//这里要自己去理解，利用叉乘求面积
    double bcd=(p[j].x1-p[i].x2)*(p[j].y2-p[i].y2)-(p[j].x2-p[i].x2)*(p[j].y1-p[i].y2);
    if(acd*bcd<=0)
        flag1=1;
    double dab=(p[i].x1-p[j].x2)*(p[i].y2-p[j].y2)-(p[i].x2-p[j].x2)*(p[i].y1-p[j].y2);
    double cad=(p[i].x1-p[j].x1)*(p[i].y2-p[j].y1)-(p[i].x2-p[j].x1)*(p[i].y1-p[j].y1);
    if(dab*cad<=0)
        flag2=1;
    if((flag1+flag2)==2)//这里要满足两个限制条件，就是既要你在他两侧，也要他在你两侧
        return 1;
    else
        return 0;
}
int main()
{
    int n;
    while(~scanf("%d",&n),n)
    {
        for(int i=1; i<=n; i++)
            scanf("%lf%lf%lf%lf",&p[i].x1,&p[i].y1,&p[i].x2,&p[i].y2);
        int sum=0;
        for(int i=1; i<n; i++)
        {
            for(int j=i+1; j<=n; j++)
            {
                if(slove(i,j))
                    sum++;
            }
        }
        printf("%d\n",sum);
    }
    return 0;
}