3

当我使用 gmp.h 头文件时。我需要一个函数,它接受 mpz_t 类型的输入并返回 mpz_t 类型。我是使用 gmp.h 的初学者所以,这里是我接近的代码的快照...

mpz_t sum_upto(mpz_t max)  
{    
    mpz_t sum;  
    mpz_init(sum);    
    mpz_init(result);  
    for(int i=0;i<=max-1;i++)    
        mpz_add_ui(sum,sum,pow(2,i));   
    return sum;   
}

但它会显示错误:

  1. pow 未在此范围内使用。”,尽管我在文件的最开头添加了 math.h。
  2. sum_upto 声明为返回数组的函数...
4

2 回答 2

8

使用 GMP 的函数约定可以在手册中找到。本质上,您必须遵循与 GMP 本身相同的约定 - 函数必须具有 void 返回类型,并且您必须提供将结果作为参数放入的值。

这是给出的示例:

 void foo (mpz_t result, const mpz_t param, unsigned long n)
 {
   unsigned long  i;
   mpz_mul_ui (result, param, n);
   for (i = 1; i < n; i++)
     mpz_add_ui (result, result, i*7);
 }

 int main (void)
 {
   mpz_t  r, n;
   mpz_init (r);
   mpz_init_set_str (n, "123456", 0);
   foo (r, n, 20L);
   gmp_printf ("%Zd\n", r);
   return 0;
 }
于 2012-11-15T10:54:19.380 回答
0

尝试以下操作:

mpz_t sum_upto(mpz_t max)
{
    mpz_t sum;
    mpz_init(sum);
    mpz_init(result);
    int val = 1;
    for(int i=0;i<=max-1;i++) {
        mpz_add_ui(sum,sum,val);
        val *= 2; //compiler should make a shift operation out of it
    }
    return sum;
}

此外,您可能会删除math.h标题。

于 2010-06-06T19:38:13.370 回答