三维几何（three-dimensional geometry）.docVIP

三维几何（three-dimensional geometry） # includes stdio.h # includes string.h # includes math.h # includes iostream # algorithm includes # vector includes using namespace std; # define mp make _ pair # define pb push _ back # define eps 1 - 8 # define pi acos (- 1.0) int sig (double) { return to - eps? - 1: - eps. } struct tpoint { double x, y, z; tpoint () {} tpoint (double x0, y0 is double, double 20): x (x0), y (y0), z (20) {} void in () { scanf (% lf% lf% lf, x y z); } tpoint operator * (double h) { return tpoint (x * h, y * h, z * h); } tpoint operators / (double h) { return tpoint (x / h, y / h, z / h); } double operator ^ (tpoint h) {/ / 点乘 return x + y * * h.x h.y + z * h.z; } tpoint operator + (tpoint h) { return tpoint (x + y + h.x, h.y, z + h.z); } tpoint operator - (tpoint h) { return tpoint (x - h.x, y - h.y, z - h.z); } tpoint operator * (tpoint b) {/ / 叉乘 return tpoint (y * b.z - b.y * z, z * b.x - b.z * x, x * b.y - b.x * y); } double len2 () const { return x * x + y * y + z * z; } double len () const { return sqrt (len2 ()); } tpoint turnlen (double l) const {/ / 改变长度 double r = l / len (); return tpoint (x r, y * r, z * r). } }; struct tline { tpoint, b; tline () {} tline (tpoint tpoint a0, b0), (a0), b (b0) {} }; struct tplane { tpoint a, b, c, or; tplane () {} tplane (tpoint tpoint a0, b0, c0 tpoint) { a = a0. b = b0. c = c0. or = (b - a) * (c - a); } tplane (double a0, b0 c0 double, double, double d0) / / ax + by + cz + d = 0, 构成点法式 { or = tpoint (a0, b0, c0). if (sig (a0). = 0) a = tpoint ((b0 - c0 - d0) / a0,1,1); else if (sig (b0). = 0) a = (1, tpoint (a0 - c0 - d0) / b0,1); else if (sig (c0). = 0) a = tpoint (1.1, (a0 - b0 - d0) / c0). } }; double area (tpoint, tpoint b) / / 两点之间距离 { return (a - b).len (); } double area (tpoint, tpoint b, tpoint c) / / 三角形面积 * 2 { return ((b - a) * (c -).len (); } double volume (tpoint, tpoint b, tpoint c, tpoint d) / / 四面体有向体积 * 6 { return (b - a) * (c - a) ^ (d -). } bool is _ conan (tpoint, tpoin