【计算三维空间中两条直线间的最短距离】 问题描述: 已知三维

#include <cmath>
#include <cstdio>
#include <vector>
using namespace std;

struct Point {
    float x, y, z;
};

Point point2vector(const Point &point1, const Point &point2) {
    Point vec{point1.x - point2.x, point1.y - point2.y, point1.z - point2.z};
    return vec;
}

Point cross(const Point &point1, const Point &point2) {
    float i = point1.y * point2.z - point1.z * point2.y;
    float j = point1.z * point2.x - point1.x * point2.z;
    float k = point1.y * point2.x - point1.x * point2.y;
    return {i, j, k};
}

float dot(const Point &point1, const Point &point2) {
    return point1.x * point2.x + point1.y * point2.y + point1.z * point2.z;
}

float norm(const Point &point) {
    return sqrt(point.x * point.x + point.y * point.y + point.z * point.z);
}

int main() {
    vector<Point> points;
    float x, y, z;
    while (scanf("%f %f %f", &x, &y, &z) != EOF) {
        points.push_back({x, y, z});
    }
    Point A = points[0], B = points[1], C = points[2], D = points[3];
    Point vecAB = point2vector(B, A), vecCD = point2vector(D, C), vecAC = point2vector(C, A);
    Point vecN = cross(vecAB, vecCD);
    float dis = abs(dot(vecAC, vecN)) / norm(vecN);
    printf("%f", dis);
    return 0;
}