HDU1411 欧拉四面体

用向量解决：

三角形面积：S=1/2*|x1*y2-x2*y1|;      (粗体表示向量)

 

三棱锥体积：V=1/6*(OA*OB)*OC                                                               

    

不知道哪里去找的代码，毕竟很线性代数矩阵什么的很头疼，晚上再手动比划下。

#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
int main()
{
    double l,n,a,m,b,c;
    while(scanf("%lf%lf%lf%lf%lf%lf",&l,&n,&a,&m,&b,&c)!=EOF){
        double x,y;
        x=4*a*a*b*b*c*c-a*a*(b*b+c*c-m*m)*(b*b+c*c-m*m)-b*b*(c*c+a*a-n*n)*(c*c+a*a-n*n);
        y=c*c*(a*a+b*b-l*l)*(a*a+b*b-l*l)-(a*a+b*b-l*l)*(b*b+c*c-m*m)*(c*c+a*a-n*n);
        double num=(sqrt(x-y)/12);
        printf("%.4lf\n",num);
    }
    return 0;
}