用 mpz_t 制作帕斯卡三角形

Question

嘿，我正在尝试将我编写的用于生成表示 Pascal 三角形的 long 数组的函数转换为返回 mpz_t 数组的函数。但是，使用以下代码：

mpz_t* make_triangle(int rows, int* count) {
//compute triangle size using 1 + 2 + 3 + ... n = n(n + 1) / 2
*count = (rows * (rows + 1)) / 2;
mpz_t* triangle = malloc((*count) * sizeof(mpz_t));

//fill in first two rows
mpz_t one;
mpz_init(one);
mpz_set_si(one, 1);
triangle[0] = one; triangle[1] = one; triangle[2] = one;

int nums_to_fill = 1;
int position = 3;
int last_row_pos;
int r, i;
for(r = 3; r <= rows; r++) {
    //left most side
    triangle[position] = one;
    position++;

    //inner numbers
    mpz_t new_num;
    mpz_init(new_num);
    last_row_pos = ((r - 1) * (r - 2)) / 2;
    for(i = 0; i < nums_to_fill; i++) {
        mpz_add(new_num, triangle[last_row_pos + i], triangle[last_row_pos + i + 1]);
        triangle[position] = new_num;
        mpz_clear(new_num);
        position++;
    }
    nums_to_fill++;

    //right most side
    triangle[position] = one;
    position++;
}

return triangle;
}

我收到错误消息：在设置三角形位置的所有行的赋值中类型不兼容（即：triangle[position] = one;）。

有谁知道我可能做错了什么？

Answer 1

mpz_t 被定义为 struct __mpz_struct 长度为 1 的数组，这会阻止赋值。这样做是因为正常的 C 赋值是浅拷贝，并且各种 gmp 数字类型存储指向 数组的指针需要深度复制的“肢体”。您需要使用 mpz_set 或 < a href="http://gmplib.org/manual/Simultaneous-Integer-Init-_0026-Assign.html#Simultaneous-Integer-Init-_0026-Assign" rel="nofollow noreferrer">mpz_init_set （甚至 mpz_init_set_si）来分配 MP 整数，确保在使用前者之前初始化目标。

另外，您应该为每个 mpz_init 最多调用一次 mpz_clear （在这方面它们类似于 malloc 和 free，并且出于相同的原因）。通过在外部循环中调用 mpz_init(new_nom) 在内循环中的 mpz_clear(new_num) ，您将引入一个错误，当您检查 mpz_clear(new_num) 的结果时，该错误将会很明显。代码>make_triangle。但是，您甚至不需要 new_num；初始化triangle的下一个元素并将其用作mpz_add的目的地。

    mpz_init(triangle[position]);
    mpz_add(triangle[position++], triangle[last_row_pos + i], triangle[last_row_pos + i + 1]);

小数值优化：您可以使用加法和减法来更新 last_row_pos，而不是使用两个减法（乘法和除法）。看看你是否能弄清楚如何做。