1

私は現在libgmpを学んでおり、そのために素因数を見つける小さなプログラムを書いています。私のプログラムは、さまざまな量の mpz_t 整数 (特定の数値の素因数) で配列を埋める関数を呼び出します。これを返す必要があります。最後の要素を NULL に設定することを計画しているので、関数が検出した mpz_t 整数の数を知っています。

私の問題は、mpz_t 整数へのポインターの配列でダブル フリー エラーが発生することです。私の問題を説明するサンプルコードをいくつか書きました：

#include <stdlib.h>
#include <stdio.h>
#include <gmp.h>

int main(void)
{
    mpz_t *p = malloc(5*sizeof(mpz_t*));
    mpz_init_set_ui(p[0], 2UL);
    mpz_init_set_ui(p[1], 5UL);
    gmp_printf("%Zd %Zd\n", p[0], p[1]);
    mpz_clear(p[0]);
    mpz_clear(p[1]);
    free(p);
    return 0;
}

2 と 5 は stdout に出力されるので、割り当ては問題ないようです。しかし、私は以下のダブルフリーエラーを取得しています:

2 5
*** glibc detected *** ./lol: double free or corruption (out): 0x08e20020 ***
======= Backtrace: =========
/lib/libc.so.6(+0x6b6c1)[0xb77126c1]
/lib/libc.so.6(+0x6cf18)[0xb7713f18]
/lib/libc.so.6(cfree+0x6d)[0xb7716f8d]
/usr/lib/libgmp.so.3(__gmp_default_free+0x1d)[0xb77f53fd]
/usr/lib/libgmp.so.3(__gmpz_clear+0x2c)[0xb77ff08c]
./lol[0x80485e3]
/lib/libc.so.6(__libc_start_main+0xe6)[0xb76bdb86]
./lol[0x80484e1]

私はまだポインターに完全に慣れており、gccはエラーを出しませんが、これは間違っていると確信しており、次のようなことをする必要があります

mpz_init_set_ui(*p[0], 2UL);

それ以外の：

mpz_init_set_ui(p[0], 2UL);

しかし、それは私にコンパイラエラーを与えます

test.c:8: error: incompatible type for argument 1 of ‘__gmpz_init_set_ui’
/usr/include/gmp.h:925: note: expected ‘mpz_ptr’ but argument is of type ‘__mpz_struct’

とにかく、私の質問は次のとおりです。

  1. mpz_init_set_ui() 呼び出しでポインターを逆参照する必要があると確信していますが、なぜそれが間違っているのですか?
  2. これを行うより良い方法はありますか？連結リストを使用する必要がありますか?(連結リストはまだ学習していません。これには配列が最適だと思いますが、本当に難しい場合は教えてください) 3. 構造体を作成した方がよいでしょうか配列へのポインターと、配列内の要素の量を含む別の変数を使用して、代わりにそれへのポインターを返しますか?

プラットフォームは、関連する場合に備えて Linux 32 ビットです。

これが私が今持っているコードです。これを変更したいのですが、スタックで mpz_t の配列を宣言します。しかし、 main() を関数にしたい:

#include <stdio.h>
#include <stdlib.h>
#include "prime.h"

#define MAXFACTORS 100

int main(void)
{
    mpz_t numToFactor, factor;
    mpz_t result;/* used to pass return values from getPrimeFactor() */
    mpz_t primeFactors[MAXFACTORS];

    mpz_init_set_str(numToFactor, "18 446 744 073 709 551 615 436 457 568", 10);
    mpz_init(factor);
    mpz_init(result);

    int pFLen = 0;
    mpz_init(primeFactors[pFLen]);

    getPrimeFactor(numToFactor, result);
    mpz_set(factor, result);
    while(mpz_cmp_ui(factor, 0UL))
    {
        mpz_set(primeFactors[pFLen], factor);
        pFLen++;
        if(pFLen == MAXFACTORS)
        {
            puts("Ran out of space to store prime factors, quitting...");
        }
        mpz_init(primeFactors[pFLen]);

        mpz_divexact(factor, numToFactor, factor);
        mpz_set(numToFactor, factor);

        getPrimeFactor(factor, result);
        mpz_set(factor, result);
    }
    mpz_set(primeFactors[pFLen], numToFactor);
    pFLen++;

    int i;
    for(i = 0; i < pFLen; i++)
    {
        gmp_printf("%Zd ", primeFactors[i]);
    }
    puts("");

    mpz_clear(numToFactor);
    mpz_clear(factor);
    return 0;
}

よろしくお願いします。

4

2 に答える 2

3

この行

 mpz_t *p = malloc(5*sizeof(mpz_t*));

おそらくあなたの悩みの原因です。5 s ではなく、sへの 5 つのポインターに十分なスペースを割り当てました。mpz_t のサイズによっては、配列の末尾を超えて書き込んでいる可能性があります。mpz_tmpz_t

あなたは言いたくなるでしょう

 mpz_t *p = malloc(5*sizeof(mpz_t));

5 つの mpz_t の配列を割り当てます。

于 2009-12-12T03:00:40.873 に答える
2

あなたの質問のほんの一部

mpz_t *p = ...;

pですpointer to mpz_t。p[0]( と同じ*p) はmpz_t、そのままp[1]( と同じ*(p + 1)) です ...

mpz_init_set_ui(*p[1], 5UL); /* error */

p[1]ですmpz_t。逆参照することはできません（と思います）。
必要に応じて、次の構文を使用できます

mpz_init_set_ui(*(p + 1), 5UL);
于 2009-12-12T03:00:55.903 に答える