============
*** TODO ***
============
Technical:
 - Use bigints (maybe GMP or another similar library)
 - Parallelize computations (MPI and/or OpenMP)