c - GMP GNU 代码中的错误？

Question

当我运行此代码时，出现错误“分段错误（核心转储）”。但是不使用 mpz_get_str 就可以了。

#include<stdio.h>
#include<gmp.h>
int main()
{
char *A;
mpz_t P,Q;
mpz_init(P);
mpz_init(Q);
mpz_set_str(P,"1201858877187548528922917",10);
mpz_set_str(Q,"1248833599132922783100713",10);

mpz_t PHI,E,K,d;
mpz_init(PHI);
mpz_init(E);
mpz_init(K);
mpz_init(d);

mpz_sub_ui(PHI,P,1);
mpz_sub_ui(K,Q,1);
mpz_mul(PHI,PHI,K);
mpz_set_str(E,"37",10);
mpz_gcd(d,E,PHI);
gmp_printf("d=%Zd\n",d);
mpz_invert(d,E,PHI);
gmp_printf("d=%Zd\n",d);
A=mpz_get_str(A,2,d);

mpz_mul(K,d,E);
mpz_sub_ui(K,K,1);
mpz_divexact(K,K,PHI);



}

Answer 1

您应该将您的初始化char *A为 NULL 或 0，否则（如果它不是 NULL 或 0 偶然），GMP 将尝试写入随机内存地址，这将导致段错误。

来自 GMP 文档： http: //gmplib.org/manual/Converting-Integers.html

如果 str 为 NULL，则使用当前分配函数分配结果字符串（请参阅自定义分配）。该块将是 strlen(str)+1 个字节，这对于字符串和空终止符来说已经足够了。

如果 str 不为 NULL，它应该指向一个足够大的存储块，即 mpz_sizeinbase (op, base) + 2。两个额外的字节用于可能的减号和空终止符。