Poj-3304 Segments(计算几何—— 判断线段相交）

this -> (m≧w≦m)

yada 终于写完了

思路 ↓

投影直线上有公共点的话就是说这条直线存在一个垂线经过所有线段 就是与所有线段均相交
只要枚举这些线段端点所确定的直线 是否存在一条与所有线段均相交即可
判断线段相交用的是叉乘 首先用两点p1,p2确定了一条直线 在用p1,p2分别与计算线段两个端点计算叉乘即可
叉乘之积>0就说明线段两端点在直线的同侧 也就是直线不经过此线段
注意 ：
本题是double类型
判断是否有重点 重点定义为坐标非常接近的两点 <1e-8即可
判断叉乘之积时 只要<1e-8 即可

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
#define maxn 105
using namespace std;
const double eps = 1e-9;//趋向于0的数
int n;
/*
int sgn(double x)
{
    if(fabs(x) < eps)return 0;
    if(x < 0) return -1;
    return 1;
}
int dcmp(double x,double y)
{
	if (fabs(x - y) < eps) return 0;
	else if (x > y) return 1;
	else return -1;
}
*/
struct Point{
	double x;
	double y;
};
Point p1[maxn],p2[maxn]; //两端点 

double Cross(Point v0,Point v1,Point v2)//叉积 
{
	return (v1.x-v0.x)*(v2.y-v0.y)-(v1.y-v0.y)*(v2.x-v0.x);
}

double Dis(Point p1,Point p2)//两点距离 
{
	return sqrt((p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));
}
bool judge(Point t1,Point t2)//t1 t2两端点连成的直线 如果与其他所有线段都相交，返回1 
{
	if(Dis(t1,t2) < eps)
	return 0;
	for(int i = 0;i < n;i++)
	{
		if(Cross(t1,t2,p1[i])*Cross(t1,t2,p2[i]) > eps)
			return 0;
	}
	return 1;
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&n);
		for(int i = 0;i < n;i++)
		{
			scanf("%lf%lf%lf%lf",&p1[i].x,&p1[i].y,&p2[i].x,&p2[i].y);
			
		}
		if(n == 1)
		{
			printf("Yes!\n");
			continue;	
		}
		int ans = 0;
		for(int i = 0;i < n;i++)
		{
			if(ans) break;
			for(int j = i+1;j < n;j++)
			{
				if(judge(p1[i],p1[j])||judge(p1[i],p2[j])||judge(p2[i],p1[j])||judge(p2[i],p2[j]))
				{
					ans = 1;
				}
			}
		}
		if(ans)
		printf("Yes!\n");
		else
		printf("No!\n");

	}
	return 0;
}