Вычисление числа пи свыше 4 миллионов знаков.

#include <iostream>
#include <fstream>
#include <string>
#include <limits>
#include <boost/multiprecision/mpfr.hpp>
 
namespace mp = boost::multiprecision;
constexpr unsigned prec = std::numeric_limits<unsigned>::max();
using big_float = mp::number<mp::mpfr_float_backend<prec>>;
 
// вывод big_float в виде форматированной строки разделенной по sizeOfLine символов
void print(const std::string& number, unsigned sizeOfLine, std::ostream& os = std::cout) {
    unsigned size = number.size() - sizeOfLine;
    unsigned i = 0;
    for (; i < size; i += sizeOfLine) {
        for (unsigned j = 0; j < sizeOfLine; ++j) {
            os << number[i+j];
        }
        os << "\n";
    }
    os << number.substr(i);
}
 
int main()
{
    size_t count = 100;
 
    big_float a = 1, two = 2, an, bn;
    big_float b = a / mp::sqrt(two);
    big_float t = 0.25;
    big_float p = 1;
 
    while (count--)
    {
        an = (a + b) / 2;
        bn = mp::sqrt(a * b);
        t -= p * mp::pow((a - an), 2);
        p *= 2;
 
        a = an; b = bn;
    }
    big_float pi = mp::pow((a + b), 2) / (4 * t);
 
    std::ofstream ofs("PI.txt");
    print(pi.str(), 100, ofs);
 
    std::cout << "\nDone.\n";
    return 0;
}