Question

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

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];
    }
}
}

C-计算2D矩阵行列式的递归算法

1 个答案: