POJ 3304 Segments

题目链接: http://poj.org/problem?id=3304

------------------------------------------------------------------------------------------------------

由于几乎没写过计算几何 只要$WA$了第一反应是调$eps$ 然而这题坑点并不是调$eps$

对于一条直线是否穿过一条线段 如果线段比较短的话 只要看看线段两端点是否都贴在直线旁就好

然而如果直线的方向向量由两个非常接近的点确定 那么与其他的向量做叉积的时候几乎都是$0$

所以我们枚举确定直线方向的两端点的时候要把这一种情况跳过

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const double eps = 1e-8;
const int N = 110;
struct node
{
double x, y;
}a[N << 1];
int t, n, n2;
double cross(node &aa, node &bb, node &cc)
{
return (bb.x - aa.x) * (cc.y - aa.y) - (cc.x - aa.x) * (bb.y - aa.y);
}
bool check(node &aa, node &bb, node &cc, node &dd)
{
double t1, t2;
t1 = cross(aa, bb, cc);
t2 = cross(aa, bb, dd);
return (t1 <= eps && t2 >= -eps) || (t1 >= -eps && t2 <= eps);
}
int main()
{
scanf("%d", &t);
while(t--)
{
scanf("%d", &n);
for(int i = 1; i <= n; ++i)
scanf("%lf%lf%lf%lf", &a[i * 2 - 1].x, &a[i * 2 -1].y,
&a[i * 2].x, &a[i * 2].y);
n2 = n << 1;
bool gain = 0;
for(int i = 1; i < n2 && !gain; ++i)
for(int j = i + 1; j <= n2 && !gain; ++j)
{
if(abs(a[i].x - a[j].x) + abs(a[i].y - a[j].y) < eps * 2)
continue;
for(int k = 1; k < n2; k += 2)
if(!check(a[i], a[j], a[k], a[k + 1]))
break;
else if(k == n2 - 1)
{
gain = 1 ;
break;
}
}
if(gain)
puts("Yes!");
else
puts("No!");
}
return 0;
}