Вывести последовательность из цифр огромного числа

#include <iostream>
#include <math.h>
 
using namespace std;
double factorial(int to);
double exponential(int to, int basis);
 
int main()
{
    int i=1;
    double mantiss, whole;
 
    double toOutput = exponential(40, 2);
    
    cout << toOutput << '\n';
 
    while(true)
        {
            toOutput = toOutput / 10;
            mantiss = modf(toOutput, &whole);
            cout << (mantiss * 10) << '\n';
            toOutput = toOutput - mantiss;
            if(toOutput<1)
                break;
            i++;
        }
    cout << '\n' << i << '\n';
    return 0;
}
double factorial(int to)
{
    double factor = 1;
    for (int i = 1; i <= to; i++)
    {
        factor = factor * i;
    }
    return factor;
}
 
double exponential(int to, int basis)
{
    double expon = 1;
    for (int i = 1; i <= to; i++)
    {
        expon *= basis;
    }
    return expon;
}

#include <iostream>
#include <boost/multiprecision/cpp_int.hpp>
 
namespace mp = boost::multiprecision;
 
mp::cpp_int fact(unsigned n)
{
   mp::cpp_int res = 1;
   for (unsigned i = 2; i <= n; ++i)
      res *= i;
   return res;
}
 
int main()
{
   std::cout << "5!   = " << fact(5) << std::endl;
   std::cout << "100! = " << fact(100) << std::endl;
   std::cout << "1000! = " << fact(1000) << std::endl;
 
   std::cout<<'\n';
   return 0;
}