判断直线与线段相交 POJ 3304 Segments

题意:在二维平面中,给定一些线段,然后判断在某直线上的投影是否有公共点。

转化,既然是投影,那么就是求是否存在一条直线L和所有的线段都相交。

证明:

下面给出具体的分析:先考虑一个特殊的情况,即n=1的时候,如下图,线段AB在直线L上的投影为线段A'B',则过任意介于A'B'之间的点C'做直线L的垂线必交线段AB与一点C;反之,过线段AB之间任意一点C做直线L的垂线,垂足必定落在A'B'之间。

判断直线与线段相交 POJ 3304 Segments

不难将此结论推广到n条线段的情况,假设存在一满足题意的直线L,则设点A为各个线段在L上投影的公共点,那么过A做一条直线L的垂线L',则L'必定与n条线段都相交;反之,过所有线段做一个直线L1使其与n条线段均相交,做直线L1的垂线L2,容易发现垂足即为所有n条线段投影的公共点。

                                ---------摘自http://blog.csdn.net/Once_HNU/article/details/6327906

  判断直线和线段相交,这个和判断线段和线段相交差不多,从直线上任取两个点,和线段的端点比较,如果两个在同一侧,那么就不相交,反之,就相交。关键问题就是怎么找这个直线, 假设存在一条直线跟所有线段都相交,我们可以让这个直线旋转,条件是旋转之后仍然相交,所有一定有个临界条件,而且这个临界条件一定发生在线段的端点处, 所以,可以通过枚举所有 线段端点的方式来求这个直线。

还有需要注意的是:

1.  n = 1, n = 2时要特判,因为2条线段一定可以找第三条直线与他们相交

2. 如果枚举线段端点的时候,两个点的距离小于10^-8,那么就跳过此直线。由于计算机在浮点类型的计算中1e-8会当成三点共线来算,所以不能用它来当成要求的那条直线。

代码如下:

/*************************************************************************
> File Name: poj_3304.cpp
> Author:
> Mail:
> Created Time: 2015年04月01日 星期三 19时01分05秒
************************************************************************/ #include<iostream>
#include <cstdio>
#include <math.h>
#define EPS 1e-8
using namespace std;
struct point{
double x, y;
};
const int N = ;
point p[N];
int T, n;
//直线的叉积,来判断点在直线ab的哪一侧
double get_direction(point a, point b, point c)
{
point t1, t2;
t1.x = c.x - a.x; t1.y = c.y - a.y;
t2.x = b.x - a.x; t2.y = b.y - a.y;
return (t1.x * t2.y - t1.y * t2.x);
}
//算两点之间的距离
double get_distance(int i, int j)
{
return sqrt((p[j].x - p[i].x) * (p[j].x - p[i].x) + (p[j].y - p[i].y) * (p[j].y - p[i].y));
}
int main()
{
scanf("%d", &T);
while (T--)
{
double d1, d2;
bool flag = false;
scanf("%d", &n);
for (int i = ; i < * n; i++)
scanf("%lf %lf", &p[i].x, &p[i].y);
if (n == || n == )//特判
{
puts("Yes!");
continue;
}
for (int i = ; i < * n; i++)
{
for (int j = i + ; j < * n; j++)//这两层for是枚举所有的线段端点所在直线
{
if (get_distance(i, j) < EPS)//精度
continue;
bool tmp_flag = true;
for (int k = ; k < * n; k++)
{
d1 = get_direction(p[i], p[j], p[k]);
d2 = get_direction(p[i], p[j], p[++k]);
if (d1 * d2 > )
{
tmp_flag = false;
break;
}
}
if (tmp_flag)//如果找到直接跳出循环
{
flag = true;
break;
}
}
if (flag)//优化,找到之后直接跳出循环
break;
}
if (flag)
puts("Yes!");
else
puts("No!");
} return ;
}
上一篇:POJ 3304 Segments (直线和线段相交判断)


下一篇:程序员求爱的创意程序^_^