我应该如何使用 mpfr::mpfr_fac_ui 函数？

Question

我尝试在互联网上找到一个如何使用mpfr::mpfr_fac_ui 的示例，但我无法做到，所以我决定在这里提问。

我有自己的迭代阶乘

boost::multiprecision::mpfr_float factorial(int start, int end)
{
    boost::multiprecision::mpfr_float fact = 1;

    for (; start <= end; ++start)
        fact *= start;

    return fact;
}

但我想试试内置的factorial。

我不知道我做错了什么，因为当我像这样测试它时

mpfr_t test;
mpfr_init2(test, 1000);

std::cout << mpfr_fac_ui(test, 5, MPFR_RNDN) << std::endl;
std::cout << factorial(1, 5) << std::endl;

mpfr_fac_ui 不返回任何错误（返回 0），test 是 0 而应该是 120。

是我做错了什么还是我遗漏了什么？

Answer 1

在 C 语言中，我得到 120 的预期值：

#include <stdio.h>
#include <mpfr.h>

int main (void)
{
  mpfr_t test;
  mpfr_init2 (test, 1000);
  mpfr_fac_ui (test, 5, MPFR_RNDN);
  mpfr_printf ("%Rg\n", test);
  mpfr_clear (test);
  return 0;
}

在您的程序中，您没有显示如何打印test 的值。你所做的就是打印mpfr_fac_ui的返回值，即0。