#ifndef _GLIBCXX_NO_ASSERT
#include <cassert>
#endif
#include <cctype>
#include <cerrno>
#include <cfloat>
#include <ciso646>
#include <climits>
#include <clocale>
#include <cmath>
#include <csetjmp>
#include <csignal>
#include <cstdarg>
#include <cstddef>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>

//#ifdef __GXX_EXPERIMENTAL_CXX0X__
#include <ccomplex>
#include <cfenv>
#include <cinttypes>
#include <cstdbool>
#include <cstdint>
#include <ctgmath>
#include <cwchar>
#include <cwctype>
//#endif

#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <iterator>
#include <limits>
#include <list>
#include <locale>
#include <map>
#include <memory>
#include <new>
#include <numeric>
#include <ostream>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <typeinfo>
#include <utility>
#include <valarray>
#include <vector>

#include <array>
//#include <atomic>
#include <chrono>
//#include <condition_variable>
#include <forward_list>
//#include <future>
#include <initializer_list>
//#include <mutex>
#include <random>
#include <ratio>
#include <regex>
#include <system_error>
//#include <thread>
#include <tuple>
#include <typeindex>
#include <type_traits>
#include <unordered_map>
#include <unordered_set>


using namespace std;

typedef struct D3
{
	double x;
	double y;
	double z;
};

int main ()
{
	D3 P;
	int N;
	cin >> N;
	cin >> P.x >> P.y >> P.z;
	vector<D3>Q ( N );
	for( size_t i = 0; i < N; i++ )
	{
		cin >> Q[i].x >> Q[i].y >> Q[i].z;
	}
	double ans = 0.0;
	for( size_t i = 0; i < N; i++ )
	{
		for( size_t j = i + 1; j < N; j++ )
		{
			for( size_t k = j + 1; k < N; k++ )
			{
				D3 a;
				a.x = Q[j].x - Q[i].x;
				a.y = Q[j].y - Q[i].y;
				a.z = Q[j].z - Q[i].z;
				D3 b;
				b.x = Q[j].x - Q[k].x;
				b.y = Q[j].y - Q[k].y;
				b.z = Q[j].z - Q[k].z;

				D3 G;
				G.x = a.y*b.z - a.z*b.y;
				G.y = a.z*b.x - a.x*b.z;
				G.z = a.x*b.y - a.y*b.x;

				double S;
				S = sqrt ( G.x*G.x + G.y*G.y + G.z*G.z ) / 2;

				D3 c;
				c.x = Q[j].x - P.x;
				c.y = Q[j].y - P.y;
				c.z = Q[j].z - P.z;

				double V;
				V = ( G.x*c.x + G.y*c.y + G.z*c.z ) / 6;

				ans += abs ( V *3.0 / S );
			}
		}
	}
	cout << fixed << setprecision ( 12 ) << ans << endl;
	return 0;
}