#include <bits/stdc++.h> using namespace std; #include <boost/multiprecision/cpp_dec_float.hpp> namespace mp = boost::multiprecision; typedef mp::number<mp::cpp_dec_float<5>> d; int main(){ int N; cin>>N; vector<d> P(3); for(int i=0;i<3;i++){ cin>>P[i]; } vector<vector<d>> Q(N,vector<d>(3)); for(int i=0;i<N;i++){ for(int j=0;j<3;j++){ cin>>Q[i][j]; Q[i][j]-=P[j]; } } vector<d> a(3),b(3),law(3); d length; d product; d ans=0; for(int i=0;i<N;i++){ for(int j=i+1;j<N;j++){ for(int k=j+1;k<N;k++){ for(int l=0;l<3;l++){ a[l]=Q[i][l]-Q[k][l]; b[l]=Q[j][l]-Q[k][l]; } length=0; for(int l=0;l<3;l++){ law[l]=a[(l+1)%3]*b[(l+2)%3]-a[(l+2)%3]*b[(l+1)%3]; length+=law[l]*law[l]; } length=mp::sqrt(length); product=0; for(int l=0;l<3;l++){ product+=law[l]*Q[k][l]; } ans+=abs(product/length); } } } cout << std::fixed << std::setprecision(15) << ans << endl; }