使用GMP整数时出现分段错误(核心转储)

Question

我在C++中编写了一个简单的程序,使用gmp库,它打印出数字模数,其数字由1组成(如1,11,111,1111,1111,...)和2001; 问题是当程序达到23位数1时,我得到一个错误说Segmantation fault(core dumped).你能指出问题出在哪里吗？这是我的代码:

#include <iostream>
#include <stdio.h>
#include <string>
#include <gmpxx.h>

int main(int argc, char** args){

    mpz_t currNumber;
    mpz_init(currNumber);   
    mpz_set_str(currNumber, "1", 10);

    while(mpz_sizeinbase(currNumber, 10) < 24){
        char* digits =  mpz_get_str(nullptr, 10, currNumber);
        strcat(digits, "1");
        mpz_set_str(currNumber, digits, 10);
        digits = nullptr;


        mpz_t r;
        mpz_init(r);    
        mpz_set_str(r, "1", 20);

        mpz_t divisor;
        mpz_init(divisor);  
        mpz_set_str(divisor, "2001", 20);


        mpz_mmod(r, currNumber, divisor);

        std::cout << "====>" << currNumber << " mod(2001) = " << r << "\n\n\n";

        //Clean up
        mpz_clear(r);
        mpz_clear(divisor);     
    }

    std::cout << "Went until " << mpz_sizeinbase(currNumber, 10) << " digits !" << "\n";

    ///Clean up
    mpz_clear(currNumber);

    return 0;
}

Answer 1

第一个明显的错误是:

    char* digits =  mpz_get_str(nullptr, 10, currNumber);
    strcat(digits, "1");

由mpz_get_str分配的缓冲区不能可靠地为您将一个额外字符连接到其内容上.

我想你可以用:

    char* digits =  mpz_get_str(nullptr, 10, currNumber);
    std::string more_digits = std::string(digits) + "1";
    free(digits);
    mpz_set_str(currNumber, more_digits.c_str(), 10);

1)由于你使用的是C++,你应该使用std :: string进行大多数字符串操作,而不需要学习使用C字符串的麻烦.
2)如果我正确理解mpz_get_str(尽管从未使用过它)你需要释放它分配的缓冲区(就像我在这个建议的代码中所做的那样)以避免内存泄漏.