Segments POJ - 3304

http://poj.org/problem?id=3304

利用叉积可以判断一个点在直线的哪一侧 对于两线段确定的直线(u1 v1) (u2 v2)

如果u1与v1在(u2 v2)不同侧或u2与v2在(u1 v1)不同侧 则两直线相交

对于这道题来说 是否存在一条直线 使所有线段在该直线上的投影都相交于一点 这个问题等价于是否存在一条直线与所有线段相交 因为如果存在这样一条直线的话 该直线的垂直线便满足所有投影相交 且此两直线交点即为所有投影相交点

还有要注意的是 判断线段与线段或线段与直线相交时 只求一次叉积即可（如果是线段与直线 必须以直线为基准） 但如果是直线与直线 则必须求两次

 

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const double eps=1e-8;
const int maxn=1e2+10;

struct node
{
    double x1,y1,x2,y2;
};

node seg[maxn];
double x[2*maxn],y[2*maxn];
int n;

double getval(double x1,double y1,double x2,double y2)
{
    return x1*y2-x2*y1;
}

bool judge(double x1,double y1,double x2,double y2)
{
    double res1,res2;
    int i;
    if(fabs(x1-x2)<eps&&fabs(y1-y2)<eps) return 0;
    for(i=1;i<=n;i++){
        res1=getval(seg[i].x1-x1,seg[i].y1-y1,x2-x1,y2-y1);
        res2=getval(seg[i].x2-x1,seg[i].y2-y1,x2-x1,y2-y1);
        if(res1*res2>eps) return 0;
    }
    return 1;
}

int main()
{
    int t,i,j,flag;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        for(i=1;i<=n;i++){
            scanf("%lf%lf%lf%lf",&seg[i].x1,&seg[i].y1,&seg[i].x2,&seg[i].y2);
            x[2*i-1]=seg[i].x1,y[2*i-1]=seg[i].y1;
            x[2*i]=seg[i].x2,y[2*i]=seg[i].y2;
        }
        if(n<=2) printf("Yes!\n");
        else{
            flag=0;
            for(i=1;i<=2*n;i++){
                for(j=i+1;j<=2*n;j++){
                    if(judge(x[i],y[i],x[j],y[j])) flag=1;
                }
            }
            if(flag) printf("Yes!\n");
            else printf("No!\n");
        }
    }
    return 0;
}