POJ 3304 Segments

标签：计算几何

//按题目意思，所有线段按某方向进行投影，有一个公共点，即为找一条直线过所有线段
//我们可以先找到一条过所有线段的直线，再将该直线进行偏移直到该直线经过任两个端点，则问题转化为对任两个端点形成的直线是否过所有线段 
//像这种先找到我们需要的答案，再将答案进行偏移的思想是计算几何中一个常见的思想。 
#include<stdio.h>
#include<math.h>
const double eps=1e-8;
struct Point {
	double x,y;
	Point(){
	}
	Point(double xx,double yy){
		x=xx;y=yy;
	}
	
}p[205];
int T,n;
int cross(Point a,Point b,Point c,Point d){
	double tmp1=(c.x-a.x)*(b.y-a.y)-(c.y-a.y)*(b.x-a.x);
	double tmp2=(d.x-a.x)*(b.y-a.y)-(d.y-a.y)*(b.x-a.x);
	if(fabs(tmp1*tmp2)<eps) return 0;
	else if(tmp1*tmp2>eps) return -1;
	else return 1;
}
int judge(Point a,Point b){
	if(fabs(a.x-b.x)<eps&& fabs(a.y-b.y)<eps) return 0;//这个一定要加 
	for(int i=0;i<n;i++){
		if(cross(a,b,p[i*2],p[i*2+1])==-1) return 0;
	}
	return 1;
}
int main() {
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
#endif
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);
		for(int i=0;i<n;i++){
			scanf("%lf%lf%lf%lf",&p[i*2].x,&p[i*2].y,&p[i*2+1].x,&p[i*2+1].y);
		}
		int flag=0;
		for(int i=0;i<2*n;i++){
			for(int j=i+1;j<2*n;j++){
				if(judge(p[i],p[j])) flag=1;
			}
		}
		if(flag) printf("Yes!\n");
		else printf("No!\n");
	}
}

POJ 3304 Segments

标签：计算几何

原文地址：http://blog.csdn.net/lj94093/article/details/45309309