判断是否两条线段相交

/*    计算线段交
  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;
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

u014068781

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值