POJ 3304 Segments (直线与线段相交）

最新推荐文章于 2024-03-22 03:00:00 发布

原创最新推荐文章于 2024-03-22 03:00:00 发布 · 656 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

计算几何

本文探讨了一种方法来确定是否存在一条直线，使得所有给定的线段在其投影上有交点。通过枚举线段端点并判断线段与直线的交点，最终找到了满足条件的直线。

题意：给出许多线段，问能否找到一条直线，使得所有线段在它上面的投影存在交点。

题解：转自http://blog.sina.com.cn/s/blog_6635898a0100n2lv.html

1：把问题转化为是否存在一条直线与每条线段都有交点。证明：若存在一条直线l和所有线段相交，作一条直线m和l垂直，则m就是题中要求的直线，所有线段投影的一个公共点即为垂足。

2：枚举两两线段的各一个端点，连一条直线，再判断剩下的线段是否都和这条直线有交点。证明：若有l和所有线段相交，则可保持l和所有线段相交，左右平移l到和某一线段交于端点停止（“移不动了”）。然后绕这个交点旋转。也是转到“转不动了”（和另一线段交于其一个端点）为止。这样就找到了一个新的l满足题意，而且经过其中两线段的端点。

#include <cmath>
#include <iostream>
using namespace std;

#define MAX 300
#define eps 1e-8

struct Segement
{
	double x1, y1, x2, y2;
} seg[MAX];

struct Point
{
	double x, y;
} p[MAX];

bool samePoint ( Point & a, Point & b )
{
	if ( (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y) < eps * eps )
		return true;
	return false;
}

bool solve ( Point & a, Point & b, int n )
{
	double t1, t2;
	for ( int i = 1; i <= n; i++ )
	{
		t1 = (a.x-b.x)*(seg[i].y1-b.y) - (seg[i].x1-b.x)*(a.y-b.y);
		t2 = (a.x-b.x)*(seg[i].y2-b.y) - (seg[i].x2-b.x)*(a.y-b.y);
		if ( t1 * t2 > eps ) return false;
	}
	return true;
}

int main()
{
	int t, n;
	scanf("%d",&t); 
	while ( t-- )
	{
		scanf("%d",&n);
		int i, j = 0;
		for ( i = 1; i <= n; i++ )
		{
			scanf("%lf %lf %lf %lf",&seg[i].x1,&seg[i].y1,&seg[i].x2,&seg[i].y2);
			p[++j].x = seg[i].x1;
			p[j].y = seg[i].y1;
			p[++j].x = seg[i].x2;
			p[j].y = seg[i].y2;
		}

		bool flag = false;
		for ( i = 1; i < 2 * n; i++ )
		{
			for ( j = i + 1; j <= 2 * n; j++ )
			{
				if ( ! samePoint ( p[i], p[j] ) )
				    flag = solve ( p[i], p[j], n );
				if ( flag ) break;
			}
			if ( flag ) break;
		}
		if ( flag )
			printf("Yes!\n");
		else
			printf("No!\n");
	}
	return 0;
}

标签

#struct #n2