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。

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

1 个答案: