poj 3304 Segments

此题求解的是是否有一条直线可以穿过所有的线段，线段的断点可以刚好在这条直线上。

方法就是枚举不在同一线段上的两端点所在的直线，然后与所有的线段求叉乘。

#include <iostream>
#include <cstdio>
#include <cmath>
#define MAX 105
#define INF 99999999999999

struct node1
{
	double k;
	double b;
}segment[MAX];

struct node2
{
	double x;
	double y;
}endpoint[MAX*2];

int n;
const double eps=1e-8;
const double ep=1e-15;

double min(double a,double b)
{
	return a<b?a:b;
}
double max(double a,double b)
{
	return a>b?a:b;
}

bool search()
{
	int i,j,k;
	double a,b;
	bool flag=false;
	double x,y;
	for(i=1 ; i<=2*n&&!flag ; i++)
	{
		for(j=i+1;j<=2*n &&!flag ;j++)
		{
			flag=true;
			if(fabs(endpoint[j].x-endpoint[i].x)<eps)
			{
				a=INF;
				b=0;
			}
			else
			{
				a = ( endpoint[j].y - endpoint[i].y ) / ( endpoint[j].x - endpoint[i].x );
				b=endpoint[i].y-a*endpoint[i].x;
			}
			for(k=1;k<=n;k++)
			{
				if(fabs(a-segment[k].k)<ep)
				{
					if(a==INF)
					{
						if(fabs(endpoint[i].x-endpoint[2*k].x)>eps)
						{
							flag=false;
							break;
						}
					   else if( min(endpoint[i].y,endpoint[j].y)-max(endpoint[2*k].y,endpoint[2*k-1].y)>eps 
						      || min(endpoint[2*k].y,endpoint[2*k-1].y)-max(endpoint[i].y,endpoint[j].y)>eps )
					 	{
							flag=false;
							break;
						}
					}
					else if(a==0)
					{
						if(endpoint[i].y !=endpoint[2*k].y)
						{
							flag=false;
							break;
						}
						else if(min(endpoint[i].x,endpoint[j].x)-max(endpoint[2*k].x,endpoint[2*k-1].x)>eps 
							    ||min(endpoint[2*k].x,endpoint[2*k-1].x)-max(endpoint[i].x,endpoint[j].x)>eps)
						{
							flag=false;
							break;
						}
					}
					else if(fabs(b-segment[k].b)>eps)
					{
					
						flag=false;
						break;
					}
				}
				else
				{
					if(segment[k].k==INF)
					{
						y=a*(endpoint[2*k-1].x)+b;
						if(y-max(endpoint[2*k-1].y,endpoint[2*k].y)>eps || min(endpoint[2*k-1].y,endpoint[2*k].y)-y>eps)
						{
							flag=false;
							break;
						}
					}
					else if(a==INF)
					{
						if(endpoint[i].x-max(endpoint[2*k-1].x,endpoint[2*k].x)>eps || min(endpoint[2*k-1].x,endpoint[2*k].x)-endpoint[i].x>eps)
						{
							flag=false;
							break;
						}
					}
					else
					{
						x=(segment[k].b-b)/(a-segment[k].k);
						y=a*x+b;
						if(x-max(endpoint[2*k-1].x,endpoint[2*k].x)>eps || min(endpoint[2*k-1].x,endpoint[2*k].x)-x>eps )
						{
							flag=false;
								break;
						}
					}
				}
			}
			if(flag==true)
				return true;
		}
	}
	return false;
}

int main()
{
	freopen("test.txt","r",stdin);
	int T;
	bool flag;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&n);
	
		for(int i=1;i<=n*2;i+=2)
		{
			scanf("%lf%lf" ,&(endpoint[i].x),&(endpoint[i].y));
			scanf("%lf%lf" ,&endpoint[i+1].x,&endpoint[i+1].y);
		
			if(endpoint[i+1].x == endpoint[i].x)
			{
				segment[(i+1)/2].k=INF;
				segment[(i+1)/2].b=0;
			}
			else
			{
				segment[(i+1)/2].k = ( endpoint[i+1].y - endpoint[i].y ) / ( endpoint[i+1].x - endpoint[i].x ) ;
				segment[(i+1)/2].b = endpoint[i].y - segment[(i+1)/2].k * endpoint[i].x ;
			}
	
		}

		flag=search();
		if(flag)
			printf("Yes!\n");
		else
			printf("No!\n");

	}
	return 0;
}