1264 线段相交解题报告（51nod）（有关叉乘 ）

推荐博客：https://www.cnblogs.com/fantastic123/p/9038411.html

给出平面上两条线段的两个端点，判断这两条线段是否相交（有一个公共点或有部分重合认为相交）。 如果相交，输出"Yes"，否则输出"No"。

Input

第1行：一个数T，表示输入的测试数量(1 <= T <= 1000)
第2 - T + 1行：每行8个数，x1,y1,x2,y2,x3,y3,x4,y4。(-10^8 <= xi, yi <= 10^8)
(直线1的两个端点为x1,y1 | x2, y2,直线2的两个端点为x3,y3 | x4, y4)

Output

输出共T行，如果相交输出"Yes"，否则输出"No"。

Input示例

2
1 2 2 1 0 0 2 2
-1 1 1 1 0 0 1 -1

Output示例

Yes
No

---

我开始时这样想的 

在图中，若两线段相交，这四个向量（a1,a2,a3,a4）会有一定的关系 ,即若相交，则（a1+a2）*(a3+a4)<0-------------结果这种方法不行

其实发现只要四个端点形成的四边形是凸四边形（这个该怎么判断呢，），则两线段相交。

别人是用叉积算的：

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
struct Point
{
    double x;
    double y;
};
Point a,b,c,d;
//叉积
double mult(Point a, Point b, Point c)
{
	return a.x*b.y+a.y*c.x+b.x*c.y-b.y*c.x-a.y*b.x-a.x*c.y;//三角形的有向面积
	//等价(a.x-c.x)*(b.y-c.y)-(b.x-c.x)*(a.y-c.y);---》这个是两向量的叉乘，外的积减内的积,大小为  两边的乘积乘以他们的夹角
}


bool intersect(Point aa, Point bb, Point cc, Point dd)
{

	if(max(aa.x,bb.x)<min(cc.x,dd.x)||max(aa.y,bb.y)<min(cc.y,dd.y)||max(cc.x,dd.x)<min(aa.x,bb.x)||max(cc.y,dd.y)<min(aa.y,bb.y))
        return false;
    if(mult(aa,bb,cc)*mult(aa,bb,dd)>0||mult(cc,dd,aa)*mult(cc,dd,bb)>0)
        return false;
    return true;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y,&d.x,&d.y);
        if(intersect(a,b,c,d))
            printf("Yes\n");
        else
            printf("No\n");
    }
    return 0;
}

与https://blog.csdn.net/weixin_42373330/article/details/81782965有些关联

---

上面的叉乘是计算其数值大小的，其正负号表示方向，但不知到具体的向量坐标。有时我们需要知道叉乘后的向量，可以这样算：百度了下是这样的（看了看发现和线代里的代数余子式是一样的，可能整个行列式为零

x       y       z

X1   Y1     Z1

X2   Y2     Z2

 

）

设a=(X1,Y1,Z1),b=(X2,Y2,Z2),
a×b=（Y1Z2-Y2Z1,Z1X2-Z2X1,X1Y2-X2Y1）

下面的cross_product就是（下面的代码拷自51nod 关鑫）

#include <stdio.h>

struct  p3d_t {
        double x, y, z;
        p3d_t  operator - ( p3d_t & o ) {
                p3d_t  ans;
                ans.x = this->x - o.x;
                ans.y = this->y - o.y;
                ans.z = this->z - o.z;
                return ans;
        }
};

double  dot_product( const p3d_t & a, const p3d_t & b ) {
        return a.x * b.x + a.y * b.y + a.z * b.z;
}

p3d_t   cross_product( const p3d_t & a, const p3d_t & b ) {
        p3d_t  ans;
        ans.x = a.y * b.z - a.z * b.y;
        ans.y = a.z * b.x - a.x * b.z;
        ans.z = a.x * b.y - a.y * b.x;
        return ans;
}

int main() {
        p3d_t  p[4];
        int T;
        for( scanf( "%d", &T ); T--; ) {
                for( int i = 0; i < 4; ++i ) scanf( "%lf%lf%lf", &p[i].x, &p[i].y, &p[i].z );
                puts( dot_product( p[3] - p[0], cross_product( p[2] - p[0], p[1] - p[0] ) ) == 0.0 ? "Yes" : "No" );
        }
        return 0;
}