空间四面体体积公式

美少女zss

已于 2022-07-09 21:10:39 修改

阅读量1.3k

点赞数 2

分类专栏：重要的小操作几何文章标签： c++

于 2021-10-23 13:56:23 首次发布

本文链接：https://blog.csdn.net/m0_52398496/article/details/120920148

版权

重要的小操作同时被 2 个专栏收录

17 篇文章

订阅专栏

几何

5 篇文章

订阅专栏

int T,n,m;
int va[N];

db volume(db x1,db y1,db z1,db x2,db y2,db z2,db x3,db y3,db z3,db x4,db y4,db z4);

signed main()
{
	IOS;
	double x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4;
    cin>>x1>>y1>>z1;
    cin>>x2>>y2>>z2;
    cin>>x3>>y3>>z3;
    cin>>x4>>y4>>z4;
    double s = volume(x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4);
    printf("%.10lf\n",s);
}

db volume(db x1,db y1,db z1,db x2,db y2,db z2,db x3,db y3,db z3,db x4,db y4,db z4)
{
    db A1, A2, A3, A4, value;
    A1 = x2 * (y3 * z4 - y4 * z3) - x3 * (y2 * z4 - y4 * z2) + x4 * (y2 * z3 - y3 * z2);
    A2 = x1 * (y3 * z4 - y4 * z3) - x3 * (y1 * z4 - y4 * z1) + x4 * (y1 * z3 - y3 * z1);
    A3 = x1 * (y2 * z4 - y4 * z2) - x2 * (y1 * z4 - y4 * z1) + x4 * (y1 * z2 - y2 * z1);
    A4 = x1 * (y2 * z3 - y3 * z2) - x2 * (y1 * z3 - y3 * z1) + x3 * (y1 * z2 - y2 * z1);
    value = (A1 - A2 + A3 - A4) / 6;
    return fabs(value);
}