判断是否两条线段相交

/*    计算线段交
1、如何判断两条直线相交
跨立：是指某线段的两个端点分别处于另一线段所在的直线的两旁，
或者其中一个端点处于另一线段所在的直线上。
判断线段p1p2是否线段p3p4相交，只要判断下述两个条件是否同时成立
1、p1p2跨立线段p3p4所在的直线
2、p3p4跨立线段p1p2所在的直线
以p3p4跨立线段p1p2所在的直线实验为例，其设计思想是，从出发向另
一线段的两个端点p3和p4引出两条辅助线段p1p3,p1p4,然后计算两个叉积。
由(p3-p1)^(p2-p1)确定p1p3在p1p2的哪个方向；由(p4-p1)^(p2-p1)确定p1p4在p1p2的哪个方向
1、若是两个叉积正负号相同，都在p1p2的同边，说明p3p4不跨立p1p2所在的直线；
2、若是两个叉积正负号相同，分别在p1p2的异边，说明p3p4跨立p1p2所在的直线；
3、若是两个叉积有一个为0，都在p1p2的同边，说明p3p4跨立p1p2所在的直线。
*/
/* poj 2653
有n根不同长度的棍子，一次一根的随机扔在地上，找出没有其他棍子压在这样的棍子上面。
例子：
sample Input
5
1 1 4 2
2 3 3 1
1 -2.0 8 4
1 4 8 2
3 3 6 -2.0
0
Sample Output
Top sticks: 2, 4, 5.
*/
#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
using namespace std;
const double epsi=1e-10;
const int MAXN=100010;
struct Point
{
double x,y;
Point(double _x=0,double _y=0):x(_x),y(_y){}
Point operator -(const Point &op2) const
{
return Point(x-op2.x,y-op2.y);
}
double operator^(const Point &op2) const
{
return x*op2.y-y*op2.x;
}
};
inline int sign(const double &x)
{
if(x>epsi) return 1;
if(x<-epsi) return -1;
return 0;
}
inline double sqr( const double &x)
{
return x*x;
}
inline double dis(const Point &p1,const Point &p2)
{
return sqrt(sqr(p1.x-p2.x)+sqr(p1.y-p2.y));
}
inline double mul(const Point &p0,const Point &p1,const Point &p2)
{
return (p1-p0)^(p2-p1);
}
inline int cross(const Point &p1,const Point &p2,const Point &p3,const Point &p4)
{//判断p1p2是否跨立与p3p4
double a1=mul(p1,p2,p3),a2=mul(p1,p2,p4);
if(sign(a1)==0&&sign(a2)==0) return 2;
if(sign(a1)==sign(a2)) return 0;
return 1;
}
int n;
Point p1[MAXN],p2[MAXN];
int main()
{
int test=0;
while(scanf("%d",&n)==1)
{
if(n==0) break;
printf("Top sticks:");
bool flag1=false;
int i,j;
for(i=1;i<=n;i++)
{
scanf("%lf%lf%lf%lf",&p1[i].x,&p1[i].y,&p2[i].x,&p2[i].y);
}
for(i=1;i<=n;i++)
{
bool flag2=false;
for(j=i+1;j<=n;j++)
{
if(cross(p1[i],p2[i],p1[j],p2[j])==1&&cross(p1[j],p2[j],p1[i],p2[i])==1)
{
flag2=true;break;
}
}
if(flag2==false&&flag1==true) printf(",");
if(flag2==false) printf(" %d",i),flag1=true;
}
printf(".\n");

}
return 0;
}