gmp中mpz_t是指针?结构体?或者是其它什么东西?

没——问题

如果是个结构体,那怎么这样改变变量的值呢:

1. mpz_t a,b,c;
2. mpz_add(a,b,c);

复制代码

此时a就等于b+c了 可如果mpz_t声明的是指针,那它的空间是什么时候分配的呢?

mathe

指针。 mpz_init初始化

没——问题

可是mpz_init(a);的参数也是mpz_t型的,这样的话mpz_init()得到的就是一个指针的副本,它如何做到使原指针指向自己开辟的空间呢? 比如下面的代码是会出错的

1. #include <stdio.h>
2. #include <stdlib.h>
4. void f(int *x)
5. {
6. x=(int *)malloc(sizeof(int));
7. }
8. int main()
9. {
10. int *a;
11. f(a);
12. *a=3;
13. printf("%d",*a);
14. return 0;
15. }

复制代码

mathe

有几种不同的做法，可以通过宏来做到

没——问题

感谢mathe的解答

wayne

5# 没——问题 请教一个问题，自己分配空间自己再free 有什么好处吗。 ================================================== GMP里面有一大堆的 初始化，赋值的函数可以用的，如 mpz_init，mpz_inits， 甚至，初始化和赋值合二为一，一个函数就能搞定的，比如： mpz_init_set ，mpz_init_set_ui，mpz_init_set_si，mpz_init_set_d，mpz_init_set_str 我感觉GMP提供的这些函数都已经很充分了。

没——问题

回wayne兄... 没想自己弄,只是用着用着对mpz_t的性质起了疑心,翻代码没翻出头绪就来请教了

G-Spider

7# 没——问题 从gmp.h的代码片段

1. #ifdef __GMP_SHORT_LIMB
2. typedef unsigned int mp_limb_t;
3. typedef int mp_limb_signed_t;
4. #else
5. #ifdef _LONG_LONG_LIMB
6. typedef unsigned long long int mp_limb_t;
7. typedef long long int mp_limb_signed_t;
8. #else
9. typedef unsigned long int mp_limb_t;
10. typedef long int mp_limb_signed_t;
11. #endif
12. #endif
13. typedef unsigned long int mp_bitcnt_t;
15. /* For reference, note that the name __mpz_struct gets into C++ mangled
16. function names, which means although the "__" suggests an internal, we
17. must leave this name for binary compatibility. */
18. typedef struct
19. {
20. int _mp_alloc; /* Number of *limbs* allocated and pointed
21. to by the _mp_d field. */
22. int _mp_size; /* abs(_mp_size) is the number of limbs the
23. last field points to. If _mp_size is
24. negative this is a negative number. */
25. mp_limb_t *_mp_d; /* Pointer to the limbs. */
26. } __mpz_struct;
28. #endif /* __GNU_MP__ */
31. typedef __mpz_struct MP_INT; /* gmp 1 source compatibility */
32. typedef __mpz_struct mpz_t[1];

复制代码

可以看出定义: mpz_t a,b,c; 就是定义结构体类型的指针变量。 内存管理通过三个函数malloc(), realloc(),free() __mpz_struct结构体中，_mp_d为数据块指针，_mp_alloc为分配的数据块大小，_mp_size为实际有效数据大小和正负空信息。 init(); 根据_mp_alloc大小（初始化为1）通过malloc()分配空间，并赋块指针给_mp_d，并置_mp_size=0 mpn为通用模块（给mpz,mpf,mpq调用），mpz为整形大数模块，mpf为浮点大数模块，mpq为有理大数模块。