Segments POJ 3304_叉乘+直线与线段相交+思维

 题目链接：Segments

题意：给出n条线段。问是否存在这样一条直线，使得各线段到这段直线的投影至少有一个交点。

题解：假设存在题目所说的这条直线，作这条直线的垂线，
由于交点在这条直线上，因为投影是线段根据垂线映射下来的一部分，
要使各线段的投影至少有一个交点，那么作经过这个交点的垂线，必定与所有线段相交。
1. 所以，只要找出一条直线(即上述的垂线)能与所有线段都相交即可
2. 该直线(垂线)至少经过所有线段中的两个端点
    理由：假设所有线段的投影重合部分(这个部分以下简称C)只有一个点，那么只有两种情况：

        (1)至少存在一条线段,与垂线平行(此时这条线段就有2个端点了)
        (2)至少存在一条线段,只有端点P1的投影在C处，现在假设除了这一条，其他线段投影重合部分不止一个点。但同时也因为C处只有一个点，那么可以反证至少还存在第二条线段的端点P2的投影刚好在了重合处(否则第二条线段的投影就会越过C处，使得第一条线段的投影不止一个点，那么C也从而不止一个点，就违背了前面的假设)
        综上，以上两种情况，我们都能找到经过垂线的两个端点(P1和P2)。

所以，枚举所有点，从中取两个点作为我们所要找的垂线，利用叉乘判断这条垂线是否能与所有线段都相交即可(相交判断：以直线为基准，线段中的一个点在直线的左边，另一个点在直线的右边)

3. 注意点（坑点）
    （1）考虑重合的点：如果两个点小于题目精度，则是重合点
    （2）叉乘判断那里，注意要以你利用2个端点来确定的直线为基准，不能以线段为基准，否则叉乘会有错，通过简单画图就可以印证

#include <algorithm>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <set>
#include<map>
#include<cmath>
#include <vector>
#define EPS 1e-8//精度
using namespace std;
struct Edge {
    double x1, y1, x2, y2;
    Edge() : x1(0), y1(0), x2(0), y2(0) {}
    Edge(double _x1, double _y1, double _x2, double _y2)
        : x1(_x1), y1(_y1), x2(_x2), y2(_y2) {}
};
struct Point {
    double x, y;
};
bool judge(Edge a,Edge b){//判断直线与线段是否相交
    double x1,x2,y1,y2;
    x1=(b.x2-b.x1),y1=(b.y2-b.y1);
    x2=(a.x1-b.x1),y2=(a.y1-b.y1);
    double res1=x1*y2-x2*y1;
    x2 = (a.x2 - b.x1), y2 = (a.y2 - b.y1);
    double res2=x1*y2-x2*y1;
    if(res1*res2<=EPS)
        return 1;//相交/重合
    return 0;//不相交

}
bool judgePoint(Point a,Point b){//判断点是否重合，小于精度则重合
    double res=(a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
    res=sqrt(res);
    if(res<=EPS)
        return 1;//同一个点
    return 0;
}
Point point[220];
Edge edge[110];
int main()
{
    ios::sync_with_stdio(false);
    int t,n;
    cin>>t;
    while(t--){
        memset(point,0,sizeof(point));
        memset(edge,0,sizeof(edge));
        cin>>n;
        for(int i=0;i<2*n;i++){
            double a,b,c,d;
            cin>>a>>b>>c>>d;
            point[i].x=a;point[i].y=b;point[++i].x=c;point[i].y=d;
            int t=(i-1)/2;
            edge[t].x1=a,edge[t].y1=b,edge[t].x2=c,edge[t].y2=d;
        }
        bool flag = 1;
        /* 找垂线。枚举两个点组成直线 */
        for(int i=0;i<2*n;i++){
            for(int j=i+1;j<2*n;j++){
                flag=1;
                if(judgePoint(point[i],point[j]))  continue;//坑点：判断重合点
                Edge temp(point[i].x,point[i].y,point[j].x,point[j].y);
                for(int z=0;z<n;z++){//取一条线段
                    bool res=judge(edge[z],temp);
                    if(!res){
                        flag=0;
                        break;
                    }
                }
                if(flag){
                    cout << "Yes!"<<endl;
                    i=n*n;
                    break;
                }
            }
        }
        if(!flag){
            cout << "No!"<<endl;
        }
    }
    return 0;
}