本文主要是介绍HDU 1086 You can Solve a Geometry Problem too(判断线段相交），希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
题目地址：HDU 1086
就这么一道仅仅判断线段相交的题目写了2k多B的代码。。是不是有点浪费。。。但是我觉得似乎哪里也优化不了了。。。。
判断线段相交就是利用的叉积。假如现在两条线段分别是L1和L2，先求L1和L2两个端点与L1的某个端点的向量的叉积，如果这两个的叉积的乘积小于0的话，说明L1在是在L2两个端点之间的，但此时并不保证一定相交。此时需要用同样的方法去判断L2是否在L1的两个端点之间，如果L2也在L1的两个端点之间的话，那就足以说明L1与L2相交。但是这题还需要判断是否端点也相交，当时没想到这点，导致白白调了一段时间。。至于端点的判断，我也没想到什么好的方法。。就直接暴力判断4个端点是否是同一点的情况。。
搓代码如下：
#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <ctype.h>
#include <queue>
#include <map>
#include <set>
#include <algorithm>using namespace std;
#define eqs 1e-10
struct node
{double x, y;
} point[1000];
node xiang(node a, node b)
{node f1;f1.x=a.x-b.x;f1.y=a.y-b.y;return f1;
}
int cross(node a, node b)
{double c;c= a.x*b.y-a.y*b.x;if(c>0)return 1;else if(c==0)return 0;elsereturn -1;
}
int dcmp(double x, double y)
{if(fabs(x-y)<=eqs)return 1;return 0;
}
int main()
{int n, i, j;int c1, c2, c3, c4, ans;while(scanf("%d",&n)!=EOF&&n){ans=0;for(i=0; i<n; i++){scanf("%lf%lf%lf%lf",&point[2*i].x,&point[2*i].y,&point[2*i+1].x,&point[2*i+1].y);}for(i=0; i<n; i++){for(j=0; j<i; j++){c1=cross(xiang(point[2*i],point[2*i+1]),xiang(point[2*j],point[2*i+1]));c2=cross(xiang(point[2*i],point[2*i+1]),xiang(point[2*j+1],point[2*i+1]));c3=cross(xiang(point[2*j],point[2*j+1]),xiang(point[2*i],point[2*j+1]));c4=cross(xiang(point[2*j],point[2*j+1]),xiang(point[2*i+1],point[2*j+1]));if(c1*c2<0&&c3*c4<0)ans++;else if(dcmp(point[2*i].x,point[2*j].x)&&dcmp(point[2*i].y,point[2*j].y)){ans++;}else if(dcmp(point[2*i].x,point[2*j+1].x)&&dcmp(point[2*i].y,point[2*j+1].y)){ans++;}else if(dcmp(point[2*i+1].x,point[2*j+1].x)&&dcmp(point[2*i+1].y,point[2*j+1].y)){ans++;}else if(dcmp(point[2*i+1].x,point[2*j].x)&&dcmp(point[2*i+1].y,point[2*j].y)){ans++;}}}printf("%d\n",ans);}return 0;
}
这篇关于HDU 1086 You can Solve a Geometry Problem too(判断线段相交）的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
