poj 3340 Segments（计算几何，直线跟线段的交点）

原题友情链接快速通道~ poj 3304

大概说下题意：给你一些线段的信息，要求是判断是否存在一条直线使得这些线段投影到直线上至少存在一个公共交点。

思路延伸：倘若存在一条直线满上上面要求，过公共点做该直线的垂线，这垂线必定与所有的线段都相交！由此转化为是否存在一条直线与所有线段都相交。

两点确定一条直线，则该直线至少相较于线段所有端点中的两个。

到这里，思路就很明了了。 依次枚举出所有端点的两两组合可能，判断由这两个端点确定的直线是否与所有线段都相交。

要注意的是判断重点。

具体实现看代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#define eps 1e-8
#define zero(x)(((x) > 0 ?(x):-(x)) < eps)/*浮点数判零*/
using namespace std;

const int Maxsize = 110;

typedef struct
{
	double x,y;
}Point;

typedef struct
{
	Point a,b;
}Line;

Line line[Maxsize];
int nsegment;

double CrossProduct(Point p1, Point p2, Point p0)
{
    return ( p1.x - p0.x )*( p2.y - p0.y )-( p2.x - p0.x )*( p1.y - p0.y );
}

double dis(Point a,Point b)/*两点间距离，判断是否存在重点*/
{
	return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y));
}

bool judge(Point A,Point B)/*判断A,B两个点确定的直线是否全部穿过其他线段*/
{
	if(dis(A,B)<eps)
	{
		return false;
	}
	for(int i = 0 ; i < nsegment ; i++)
	{
		if(CrossProduct(A,B,line[i].a) * CrossProduct(A,B,line[i].b) > eps)return false;
	}
	return true;
}

int main()
{
	int ncase;
	scanf("%d",&ncase);
	while(ncase--)
	{
		scanf("%d",&nsegment);
		for(int i = 0 ; i < nsegment ; i++)
		{
			scanf("%lf %lf %lf %lf",&line[i].a.x,&line[i].a.y,&line[i].b.x,&line[i].b.y);
		}
		bool flag = false;
		for(int i = 0 ; i < nsegment ; i++)
		{
			for(int j = 0 ; j < nsegment && !flag; j++)/*枚举所有可能的点的直线组合*/
			{
				if(judge(line[i].a,line[j].a)||judge(line[i].b,line[j].b)||judge(line[i].b,line[j].a)||judge(line[i].b,line[j].b))
				{
					flag = true;
				}
			}
		}
		if(flag || nsegment == 1)/*一条线段的特殊情况*/
		{
			printf("Yes!\n");
		}
		else
		{
			printf("No!\n");
		}
	}
	return 0;
}