hdu 1086(判断线段相交)

传送门:You can Solve a Geometry Problem too

题意:给n条线段,判断相交的点数。

分析:判断线段相交模板题,快速排斥实验原理就是每条线段代表的向量和该线段的一个端点与 另一条线段的两个端点构成的两个向量求叉积,如果线段相交那么另一条线段两个端点必定在该线段的两边,则该线段代表的向量必定会顺时针转一遍逆时针转一遍,叉积必定会小于等于0,同样对另一条线段这样判断一次即可。

#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
#define N 1010
#define PI acos(-1.0)
using namespace std;
const double eps=1e-;
int sgn(double x)
{
if(fabs(x)<eps)return ;
if(x<)return -;
return ;
}
struct Point
{
double x,y;
Point(){}
Point(double _x,double _y):x(_x),y(_y){}
Point operator -(const Point &a)const
{
return Point(x-a.x,y-a.y);
}
double operator *(const Point &a)const
{
return a.x*x+a.y*y;
}
double operator ^(const Point &a)const
{
return x*a.y-y*a.x;
}
};
struct Line
{
Point s,e;
Line(){}
Line(Point _s,Point _e):s(_s),e(_e){}
};
bool segcrossseg(Line l1,Line l2)
{
return
max(l1.s.x,l1.e.x)>=min(l2.s.x,l2.e.x)&&
max(l2.s.x,l2.e.x)>=min(l1.s.x,l1.e.x)&&
max(l1.s.y,l1.e.y)>=min(l2.s.y,l2.e.y)&&
max(l2.s.y,l2.e.y)>=min(l1.s.y,l1.e.y)&&
sgn((l1.s-l1.e)^(l2.s-l1.e))*sgn((l1.s-l1.e)^(l2.e-l1.e))<=&&
sgn((l2.s-l2.e)^(l1.s-l2.e))*sgn((l2.s-l2.e)^(l1.e-l2.e))<=;
}
Line L[];
int main()
{
int n;
while(scanf("%d",&n),n)
{
for(int i=;i<=n;i++)
{
double a,b,c,d;
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
L[i]=Line(Point(a,b),Point(c,d));
}
int ans=;
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
if(segcrossseg(L[i],L[j]))ans++;
printf("%d\n",ans);
}
}
上一篇:Windows下通过VNC远程访问Linux服务器,并实现可视化


下一篇:css3背景总结与解析