#include <bits/stdc++.h>
using namespace std;
mt19937 rnd(1);
#define N 1005
using db=long double;
const db PI=acos(-1);
const db eps=1e-8;
const db inf=1e15;
int sgn(db x)
{
	if(x<-eps) return -1;
	else if(x>eps) return 1;
	else return 0;
}
struct p3
{
	db x,y,z;
	p3(db a=0,db b=0,db c=0) {x=a,y=b,z=c;}
	db norm() const {return x*x+y*y+z*z;}
	db abs() const {return sqrt(norm());}
};
p3 operator + (p3 x,p3 y) {return p3(x.x+y.x,x.y+y.y,x.z+y.z);}
p3 operator - (p3 x,p3 y) {return p3(x.x-y.x,x.y-y.y,x.z-y.z);}
p3 operator * (p3 x,db y) {return p3(x.x*y,x.y*y,x.z*y);}
p3 operator / (p3 x,db y) {return p3(x.x/y,x.y/y,x.z/y);}
p3 operator * (p3 x,p3 y) {return p3(x.y*y.z-x.z*y.y,x.z*y.x-x.x*y.z,x.x*y.y-x.y*y.x);}
db operator % (p3 x,p3 y) {return x.x*y.x+x.y*y.y+x.z*y.z;}
struct plane
{
	p3 n; db d; // n dot x = d
	plane() {}
	plane(p3 a,p3 b,p3 c) {n=(c-a)*(b-a),d=n%a;}
	db side(p3 x) const {return n%x-d;}
	db dist(p3 x) const {return side(x)/n.abs();}
	p3 proj(p3 w) const {return w-n*(side(w)/n.abs()/n.abs());}
};
struct face
{
	int a[3];
	plane p;
};
int vis[N][N];
vector<face> f;
void convex(vector<p3> a)
{
	auto getface=[&](int i,int j,int k)->face{
		return {{i,j,k},plane(a[i],a[j],a[k])};
	};
	f={getface(0,1,2),getface(0,2,1)};
	vector<face> tmp[2];
	for(int i=3;i<(int)a.size();i++)
	{
		for(auto x:f)
		{
			if(x.p.dist(a[i])<-eps) tmp[0].push_back(x);
			else
			{
				tmp[1].push_back(x);
				for(int t:{0,1,2}) vis[x.a[t]][x.a[(t+1)%3]]=i;
			}
		}
		for(auto x:tmp[1])
		{
			for(int t:{0,1,2})
			{
				if(vis[x.a[t]][x.a[(t+1)%3]]==i&&vis[x.a[(t+1)%3]][x.a[t]]!=i) tmp[0].push_back(getface(x.a[t],x.a[(t+1)%3],i));
			}
		}
		f=tmp[0],tmp[0].clear(),tmp[1].clear();
	}
}
p3 bary(p3 A,p3 B,p3 C,db a,db b,db c)
{
	return (A*a+B*b+C*c)/(a+b+c);
}
struct circle:p3{db r;};
p3 circumcenter(p3 A,p3 B,p3 C)
{
	db a=(B-C).norm(),b=(C-A).norm(),c=(A-B).norm();
	return bary(A,B,C,a*(b+c-a),b*(c+a-b),c*(a+b-c));
}
circle getc(p3 A,p3 B,p3 C)
{
	p3 o=circumcenter(A,B,C);
	// printf("* %Lf %Lf %Lf\n",(A-o).abs(),(B-o).abs(),(C-o).abs());
	return {o,(A-o).abs()};
}
circle cir(p3 A,p3 B)
{
	return {(A+B)/2,(A-B).abs()/2};
}
int in(circle x,p3 y) {return ((p3)x-y).abs()<=x.r+eps;}
db mec(vector<p3> a)
{
	shuffle(a.begin(),a.end(),rnd);
	circle o={a[0],0};
	for(int i=1;i<(int)a.size();i++)
	{
		if(in(o,a[i])) continue;
		o=cir(a[i],a[0]);
		for(int j=1;j<i;j++)
		{
			if(in(o,a[j])) continue;
			o=cir(a[j],a[i]);
			for(int k=0;k<j;k++)
			{
				if(in(o,a[k])) continue;
				o=getc(a[i],a[j],a[k]);
			}
		}
	}
	return o.r;
}
int n;
p3 a[N];
db solve(face f)
{
	// printf("solving: %d %d %d\n",f.a[0],f.a[1],f.a[2]);
	db h=0;
	for(int i=0;i<n;i++) h=max(h,-f.p.dist(a[i]));
	vector<p3> vec;
	for(int i=0;i<n;i++)
	{
		p3 c=f.p.proj(a[i]);
		// printf("%d : %Lf %Lf %Lf\n",i,c.x,c.y,c.z);
		vec.push_back(c);
	}
	db r=mec(vec);
	// printf("%Lf %Lf\n",r,h);
	return PI*r*r*h;
}
void qaq()
{
	p3 A(0,0,10),B(2,4,10),C(4,1,10);
	p3 o=circumcenter(A,B,C);
	printf("%Lf %Lf %Lf\n",o.x,o.y,o.z);
}
signed main() {
#ifdef xay5421
	freopen("b.in","r",stdin);
#endif
	cin>>n;
	for(int i=0;i<n;i++) cin>>a[i].x>>a[i].y>>a[i].z;
	convex(vector<p3>(a,a+n));
	db ans=inf;
	/*
	for(auto x:f)
	{
		printf("%d %d %d : (%Lf,%Lf,%Lf) %Lf\n",x.a[0],x.a[1],x.a[2],x.p.n.x,x.p.n.y,x.p.n.z,x.p.d);
		for(int j=0;j<n;j++) printf("%d : %Lf\n",j,x.p.dist(a[j]));
	}
	*/
	for(auto x:f) ans=min(ans,solve(x));
	printf("%.10Lf\n",ans);
	return 0;
}