[POJ3304]Segments(计算几何-叉积)
程序员文章站
2022-04-01 14:45:50
...
题目:
题意:
给出一些线段,问是否存在一条直线,使所有线段在直线上的射影至少有一个公共点。
题解:
如果所有线段在直线上的射影至少有一个公共点的话,那么过这个点做这条直线的垂线,垂线一定与所有线段都相交
问题可以转化为判断是否存在一条直线与所有线段都相交
这个点的数量很少啊,而且怎么想都是跟端点有关吧,那就。。枚举端点看看这条直线是不是过所有直线?
用到了直线与线段的交点啊
利用叉积的性质
在直线上任取两个点组成一个向量,再以其中的某一个点为起点到线段的两个端点分别组成两个向量,判断这两个向量和直线上向量的叉积是否同号,同号则不相交,异号则相交
注意如果两个端点是一个点的话就不能形成一条直线,这个时候要直接退出
代码:
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int n;
const double eps=1e-9;
int dcmp(double x)
{
if (x<=eps && x>=-eps) return 0;
return (x>0)?1:-1;
}
struct po{
double x,y;
po(double X=0,double Y=0){x=X;y=Y;}
}dian[300];
struct line{
po x,y;
line(po X=po(0,0),po Y=po(0,0)){x=X;y=Y;}
}d[105];
po operator -(po x,po y){return po(x.x-y.x,x.y-y.y);}
bool operator ==(po x,po y){return x.x==y.x&&x.y==y.y;}
double cj(po x,po y){return x.x*y.y-x.y*y.x;}
bool check(po x,po y)
{
if (x==y) return 0;///
for (int i=1;i<=n;i++)
if (dcmp(cj(y-x,x-dian[i*2-1]))==dcmp(cj(y-x,x-dian[i*2])) && cj(y-x,x-dian[i*2-1])) return 0;
return 1;
}
int main()
{
int T;scanf("%d",&T);
while (T--)
{
scanf("%d",&n);
int s=0;
for (int i=1;i<=n;i++)
{
double x1,y1,x2,y2;
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
d[i]=line(po(x1,y1),po(x2,y2));
dian[++s]=po(x1,y1); dian[++s]=po(x2,y2);
}
bool fff=0;
for (int i=1;i<=s&&!fff;i++)
for (int j=i+1;j<=s&&!fff;j++)
if (check(dian[i],dian[j])){fff=1; break;}
if (fff) printf("Yes!\n");
else printf("No!\n");
}
}
推荐阅读
-
2020 年百度之星·程序设计大赛 - 测试赛 度度熊保护村庄 计算几何(叉积)+floyd(最小环)
-
牛客多校3 - Operation Love(几何+叉积确定三点顺逆)
-
2020 年百度之星·程序设计大赛 - 测试赛 度度熊保护村庄 计算几何(叉积)+floyd(最小环)
-
计算几何 - 点积与叉积
-
计算几何之用叉乘求多边形面积
-
1016 - 计算几何之直线与线段的交 - Segments(POJ 3304)
-
[POJ2653]Pick-up sticks(计算几何-叉积)
-
[POJ3304]Segments(计算几何-叉积)
-
POJ 1696 Space Ant 【几何叉积】
-
[POJ1269]Intersecting Lines(计算几何-叉积)