求外接圆面积

【题目描述】
已知三角形的三个顶点的坐标，求三角形外接圆的面积。　
【输入格式】
共三行。每行三个实数。第一行是x1,y1,z1，第二行x2,y2,z2,第三行x3,y3,z3。
-10000.0<=x1,x2,x3,y1,y2,y3<=10000.0
【输出格式】
共一行，以(x1,y1,z1),(x2,y2,z2),(x3,y3,z3)这三点为顶点的三角形的外接圆的面积。精确到小数点后三位。
【样例输入】
0 0 0
3 0 0
0 4 0
【样例输出】
19.635
【分析】
设a,b,c分别为三角形的三边长，则可得出圆的半径为

再求出圆的面积即可。
注意，只乘3.14是不行的窝~

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<math.h>
using namespace std;
double sqr(double x){ return x*x; }
int main(){
    double x1,y1,z1,x2,y2,z2,x3,y3,z3;
    cin>>x1>>y1>>z1;
    cin>>x2>>y2>>z2;
    cin>>x3>>y3>>z3;
    double a=sqrt(sqr(x1-x2)+sqr(y1-y2)+sqr(z1-z2));
    double b=sqrt(sqr(x1-x3)+sqr(y1-y3)+sqr(z1-z3));
    double c=sqrt(sqr(x3-x2)+sqr(y3-y2)+sqr(z3-z2));
    double r=a*b*c/sqrt((a+b+c)*(a+b-c)*(a-b+c)*(-a+b+c));
    double s=3.1415926535*sqr(r);
    printf("%0.3f",s);
}