malloc: corrupted top size错误排查：方阵求逆C程序问题

错误根源 #

malloc: corrupted top size 是堆内存被越界写入破坏导致的，核心问题有两处：

1. matrix_for_det 函数的越界写入 #

生成子矩阵时，每处理完一行后未重置列索引jn，导致下一行的写入位置直接延续上一行的jn值，超出子矩阵的边界，直接破坏堆结构的元数据，后续malloc调用会直接崩溃。

比如处理3x3矩阵的2x2子矩阵时，第一行处理完jn变为2，第二行开始时jn未归零，会直接写入res[1][2]，而子矩阵仅2列，触发越界写入。

2. allocate_matrix 的内存对齐隐患 #

一次性分配指针数组和数据区时，直接将matrix + n强转为double*可能违反double类型的内存对齐要求（通常需8字节对齐）。若系统中double**的大小与double的对齐要求不匹配，会导致数据区地址未对齐，写入数据时触发未定义行为，甚至破坏堆。

修复方案 #

修复 matrix_for_det 函数 #

将jn的声明移至i循环内部，每次处理新行时重置为0：

void matrix_for_det(double **matrix, int n, int m, double **res, int y, int x) {
    int in = 0;
    for (int i = 0; i < n; i++) {
        if (i == y) continue;
        int jn = 0; // 每次处理新行时重置列索引
        for (int j = 0; j < m; j++) {
            if (j == x) continue;
            res[in][jn] = matrix[i][j];
            jn++;
        }
        in++;
    }
}

修复 allocate_matrix 函数（更安全的实现） #

改为分开分配指针数组和数据区，彻底避免对齐问题，同时用calloc直接初始化0值：

double **allocate_matrix(int n, int m) {
    double **matrix = malloc(n * sizeof(double *));
    if (!matrix) return NULL;
    double *data = calloc(n * m, sizeof(double));
    if (!data) {
        free(matrix);
        return NULL;
    }
    for (int i = 0; i < n; i++) {
        matrix[i] = data + m * i;
    }
    return matrix;
}

注意：这种方式释放内存时需要先释放数据区，再释放指针数组，需修改所有free(matrix)的位置：

• main函数：

  free(matrix[0]);
  free(matrix);
  

• transpose函数：

  free((*matrix)[0]);
  free(*matrix);
  

• det函数：

  free(next[0]);
  free(next);
  

• invert函数：

  free(next[0]);
  free(next);
  // ...
  free(new[0]);
  free(new);
  

• input函数：

  free((*matrix)[0]);
  free(*matrix);
  

额外说明 #

• invert函数中调用transpose后交换n和m，但由于程序仅处理方阵，n与m始终相等，不会影响功能；若后续扩展到非方阵需调整逻辑。
• det函数的递归实现对于大矩阵效率极低，属于算法优化问题，并非本次内存错误的根源。

内容的提问来源于stack exchange，提问作者andrey-dru-mel