Как вычислить факториал большого числа?

#include <mpirxx.h>
int main ()
{
mpz_class a;
mpz_fac_ui(a,1000);
cout << "result is " << a << "\n";
return 0;
}

#include <iostream>
#include <vector>
#include <sstream>
 
void addDigit(std::string &number, int digit)
{
    number.insert(0, std::string(1, '0' + digit));
}
 
void multiply(std::string &multiplicand, std::string multiplier)
{
    std::vector<std::string > addendums;
 
    std::string::reverse_iterator i, j;
    int memorized = 0;
    int shift = 0;
 
    for(i = multiplicand.rbegin(); i != multiplicand.rend(); ++i)
    {
        std::string addend = "";
        for(j = multiplier.rbegin(); j != multiplier.rend(); ++j)
        {
            int result = (*i - '0') * (*j - '0');
            result += memorized;
            addDigit(addend, result % 10);
            memorized = result / 10;
        }
 
        if(memorized > 0)
            addDigit(addend, memorized);
 
        memorized = 0;
        addend.append(std::string(shift, '0'));
        ++shift;
        addendums.push_back(addend);
    }
 
    int maxaddend = addendums.back().size();
 
    for(std::string &x: addendums)
    {
        shift = maxaddend - x.size();
        x = std::string(shift, '0') + x;
    }
 
    std::string answer = "";
    memorized = 0;
 
    for(unsigned int digit_index = maxaddend; digit_index > 0; --digit_index)
    {
        int result = 0;
        for(unsigned int addend_index = 0; addend_index < addendums.size(); ++addend_index)
        {
            std::string addend = addendums[addend_index];
            result += (addend[digit_index - 1] - '0');
        }
 
        result += memorized;
        addDigit(answer, result % 10);
        memorized = result / 10;
    }
 
    if(memorized > 0)
        addDigit(answer, memorized);
 
    multiplicand = answer;
    std::cout << answer;
}
 
int main()
{
    std::string multiplicand = "1";
 
    for(int i = 1; i <= 1000; ++i)
    {
        std::ostringstream oss;
        oss << i;
        std::string multiplier = oss.str();
        std::cout << i << "!: ";
        multiply(multiplicand, multiplier);
        std::cout << std::endl;
    }
 
    return 0;
}