Problem Description
杭州电子科技大学即将迎来50周年的校庆,作为校庆委员会成员的我被上级要求设计一座神秘的建筑物来迎合校庆,因此我苦思冥想了一个月,终于设计出了一套方案,这座建筑物有点象古老埃及的金字塔,不过这个神秘建筑的根基是三角形的而不是矩形的,从数学的专业角度来讲,它是四面体。当我打算上交我的设计图纸的时候发现,我不知道怎么计算这个神秘建筑的体积(我知道这座建筑的各边的尺寸),于是我找来了聪明的你来帮助我解决这个难题。
Input
输入文件包含6个不超过1000的实数,每个数之间用空格隔开。每个数代表金字塔ABCD的一条棱边长度,棱边排序如下:AB,AC,AD,BC,BD,CD。
Output
输出数据应是一个实数,表示金字塔的体积,精确到4位小数。
Sample Input
2 2 2 2 2 2
Sample Output
0.9428
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1411
题解:已知任意四面体P-ABC,记PA=a,PB=b,PC=c,cos角APB=x,cos角BPC=y,cos角CPA=z,则:V=1/6*abc*sqrt(1+2xyz-x^2-y^2-z^2)
#include <iostream>
#include <algorithm>
#include <vector>
#include <cstdio>
#include <cstring>
#include <string>
#include <queue>
#define s(x) (x*x)
using namespace std;
int main()
{
double ab,ac,ad,bc,bd,cd;
while(~scanf("%lf %lf %lf %lf %lf %lf",&ab,&ac,&ad,&bc,&bd,&cd))
{
double x = (s(ab)+s(ac)-s(bc))/ab/ac/2.0;
double y = (s(ac)+s(ad)-s(cd))/ac/ad/2.0;
double z = (s(ad)+s(ab)-s(bd))/ad/ab/2.0;
double tmp = sqrt(1.0+2*x*y*z-s(x)-s(y)-s(z));
double v = ab*ac*ad*tmp/6;
printf("%.4f\n",v);
}
return 0;
}