题意:给出许多线段,问能否找到一条直线,使得所有线段在它上面的投影存在交点。
题解:转自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;
}