关闭

51nod 1264 线段相交(计算几何)

77人阅读 评论(0) 收藏 举报

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

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


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

#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;
		else
			cout<<"No"<<endl;
	}
}


0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:535次
    • 积分:46
    • 等级:
    • 排名:千里之外
    • 原创:5篇
    • 转载:1篇
    • 译文:0篇
    • 评论:0条
    文章存档