C - 递归算法计算二维矩阵的行列式

Question

我想根据辅因子的定义及其在计算行列式中的用途递归计算二维矩阵的行列式。我的代码是：

double **supr_mat(int size, double mat[size][size], int nb_colonne)
{
    double **temp;
    *temp = malloc((size - 1) * sizeof(double*));
    if(*temp == NULL)
        return 0;
    for(int i = 0; i < size - 1; i++)                  
    {
        temp[i] = malloc((size - 1)*sizeof(double));
        if(temp[i] == NULL)
        {
            for(int k = 0; k < i; k++)
                free(temp[i]);
            free(*temp);
            return 0;
        }
        for(int j = 0; j < size -1; j++)
            temp[i][j] = 0;
    }
    for(int i = 0; i < size - 1; i++)               //suppresses line, column.
    {
        for(int j = 0; j < size - 1; j++)
        {
            if(j<nb_colonne)
                temp[i][j] = mat[i+1][j];
            else
                temp[i][j] = mat[i+1][j+1];
        }
    }
    return temp;
}

void destroy_supr_mat(double** temp, int size)
{
    for(int i = 0; i < size; i++)
        free(temp[i]);
    free(*temp);
}

double determinant(int size, double mat[size][size])
{
    double det = 0;
    double temp[size - 1][size - 1];
    if(size == 1)
        return mat[0][0];
    for(int i = 0; i < size; i++)
    {
        temp = supr_mat(size, mat, i);
        if(i%2==0)
            det += (-1)*mat[0][i]*determinant(size - 1, temp);
        else
            det += 1*mat[0][i]*determinant(size - 1, temp);
    }
    free(temp, size - 1);
    return det;
}

我收到以下错误：assignment to expression with array type: temp = supr_mat(size, mat, i); 我不明白为什么会出现这样的错误，因为我的函数返回的是地址而不是数组。编译也给了我一些关于指针类型的警告，所以我想我误解了一些东西。任何解释将不胜感激。我也假设递归算法是正确的，但我不确定。

Answer 1

double**与double[size - 1][size - 1]不同。它的实现不是指针上的指针。它实际上“等效于”...double*（实现为直接指向 2D 数据的指针，而不是指向数据指针的指针）。所以你必须改变你的代码来适应这个问题。

您也可以就地使用 VLA，您不必分配任何东西（不确定 size-1 是否有效，但您明白了）：

void supr_mat(int size, double mat[size][size], int nb_colonne, double temp[size-1][size-1])
{
for(int i = 0; i < size - 1; i++)               //suppresses line, column.
{
    for(int j = 0; j < size - 1; j++)
    {
        if(j<nb_colonne)
            temp[i][j] = mat[i+1][j];
        else
            temp[i][j] = mat[i+1][j+1];
    }
}
}