题意:给出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;
}