欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

[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");
    }
}
相关标签: 计算几何 叉积