51nod 1264 线段相交（计算几何）

来源：互联网发布：英国经济学博士知乎编辑：程序博客网时间：2024/05/23 11:47

今天开始写第一个博客，一点一点累计吧

其实这道题算是基础题，也只是简单地用到了计算几何的知识，我就写了一个比较笨的方法，就是算出两条线的方程，例如算出第一条线的方程，再判断另外两个点的是否在这条线的两边，如果不是的话就是不相交，如果是的活就判断另外一条线，如果都分布在两边的话就是相交的了，当然，判断过程中，端点就在另一条线上必相交。

下面是源程序，以后再更新一下伪代码吧

#include<bits/stdc++.h>using namespace std;const double mn=1e-6;int main(){int m;cin>>m;while(m--){int x[4],y[4],i;for(i=0;i<4;i++)cin>>x[i]>>y[i];bool f1=0,f2=0,ff=0;if(x[0]==x[1]){if((x[2]==x[0]&&y[2]<=max(y[0],y[1])&&y[2]>=min(y[0],y[1]))||(x[3]==x[0]&&y[3]<=max(y[0],y[1])&&y[3]>=min(y[0],y[1])))f1=1,ff=1;else if(x[2]<x[0]^x[3]<x[0])f1=1;}else{double k,b;k=(double)(y[0]-y[1])/(x[0]-x[1]);b=y[0]-k*x[0];//cout<<k<<' '<<b<<endl;if((fabs(y[2]-k*x[2]-b)<mn&&x[2]<=max(x[0],x[1])&&x[2]>=min(x[0],x[1]))||(fabs(y[3]-k*x[3]-b)<mn&&x[3]<=max(x[0],x[1])&&x[3]>=min(x[0],x[1])))f1=1,ff=1;else if((y[2]-k*x[2]-b<0)^(y[3]-k*x[3]-b<0))f1=1;}if(x[2]==x[3]){if((x[0]==x[2]&&y[0]<=max(y[2],y[3])&&y[0]>=min(y[2],y[3]))||(x[1]==x[2]&&y[1]<=max(y[2],y[3])&&y[1]>=min(y[2],y[3])))f2=1,ff=1;else if(x[1]<x[2]^x[0]<x[2])f2=1;}else{double k,b;k=(double)(y[2]-y[3])/(x[2]-x[3]);b=y[2]-k*x[2];//cout<<k<<' '<<b<<endl;if((fabs(y[0]-k*x[0]-b)<mn&&x[0]<=max(x[2],x[3])&&x[0]>=min(x[2],x[3]))||(fabs(y[1]-k*x[1]-b)<mn&&x[1]<=max(x[2],x[3])&&x[1]>=min(x[2],x[3])))f2=1,ff=1;else if((y[0]-k*x[0]-b<0)^(y[1]-k*x[1]-b<0))f2=1;}if(ff==1)cout<<"Yes"<<endl;else if(f1&&f2)cout<<"Yes"<<endl;elsecout<<"No"<<endl;}}

0 0