NxN矩阵行列式递归问题

我目前正在尝试编写一个程序来找到NxN矩阵的行列式，但我有一个N大于2的递归问题。基本上我可以告诉它，它没有这样做，它只是运行一次使用函数我的调试选项显示该函数在列中运行，但顺序永远不会下降，然后无论如何，它都会为我的行列式提供零。 香港专业教育学院曾试图寻找任何想法，因为我做错了什么，但我似乎找不到任何答案，我甚至找到了与我做基本相同的事情的例子，使用它们给我零，不管怎么样，所以我非常困惑:(。如果有人可以快速查看我的代码并告诉我哪里是一个白痴，我会非常感激！（抱歉格式化，它在我的编辑器中看起来不错但我似乎无法得到它在这里）

码：

#include #include #include double det(double **mat, int order); int main (int argc, char* argv[]) { FILE* input; int row,column,N; double **matrix; N=3; matrix=(double**)malloc(N*sizeof(double)); input=fopen("matrix.dat", "r"); if(input !=(FILE*) NULL) { for(row=0; row<N; row++) { matrix[row]=(double*)malloc(N*sizeof(double)); } for(row=0; row<N; row++) { printf("| "); for(column=0; column<N; column++) { fscanf(input,"%lf ", &matrix[row][column]); printf("%g ", matrix[row][column]); } if(row != (N/2)) { printf("|\n"); } else { printf("|= %lf \n", det(matrix, N) ); } } return(EXIT_SUCCESS); } else { printf("*********************ERROR*********************\n"); printf("** Cannot open input file 'matrix.dat' make **\n"); printf("** sure file is present in working directory **\n"); printf("***********************************************\n"); return(EXIT_FAILURE); } } double det(double **mat, int order) { int debug; double cofact[order], determinant, **temp; determinant = 0; debug=0; if(order==1) { determinant=mat[0][0]; if(debug==1) { printf("order 1 if\n"); } } else if(order==2) { determinant= ((mat[0][0]*mat[1][1])-(mat[0][1]*mat[1][0])); if(debug==1) { printf("order 2 if\n"); } } else { int column, rowtemp, coltemp, colread; for (column=0; column<order; column++) { /* Now create an array of size N-1 to store temporary data used for calculating minors */ temp= malloc((order-1)*sizeof(*temp)); for(rowtemp=0; rowtemp<(order-1); rowtemp++) { /* Now asign each element in the array temp as an array of size N-1 itself */ temp[rowtemp]=malloc((order-1)*sizeof(double)); } for(rowtemp=1; rowtemp<order; rowtemp++) { /* We now have our empty array, and will now fill it by assinging row and collumn values from the original mat with the aprroriate elements excluded */ coltemp=0; for(colread=0; colread<order; colread++) { /* When the collumn of temp is equal to the collumn of the matrix, this indicates this row should be exlcuded and is skiped over */ if(colread==column) { continue; } temp[rowtemp-1][coltemp] = mat[rowtemp][colread]; coltemp++; } } if(debug==1) { printf("column =%d, order=%d\n", column, order); } determinant+=(mat[0][column]*(1 - 2*(column & 1))*det(temp, order-1)); } } return(determinant); } 

Answer 1

 temp= (double **)malloc((order-1)*sizeof(double)); 

只要sizeof(double*) <= sizeof(double)通常是32位或64位系统就是这种情况sizeof(double*) <= sizeof(double) ，这不会导致崩溃，但这在概念上是错误的。 你正在为double*数组分配空间，所以因子应该是sizeof(double*)或者更好，因为当类型改变时它是不变的， sizeof *temp ，

 temp = malloc((order-1) * sizeof *temp); 

（而且你不需要在C中malloc的结果，最好不要这样做，因为强制转换可以隐藏错误，比如忘记#include 。）

这同样适用于分配

 matrix=(double**)malloc(N*sizeof(double)); 

在main 。

在计算行列式时，

  for(coltemp=0; coltemp 

您在列中循环两次，对于coltemp == 0循环，然后在内循环中， coltemp按order-1递增order-1次，因此内循环在第一次运行时使用coltemp == order-1 。 然后coltemp在循环中再次递增多次，并且您正在写出已分配内存的范围。

应删除外环， coltemp = 0; 内循环是你需要的。

 pow(-1,column)) 

不是确定标志的好方法。

 (1 - 2*(column & 1)) 

比召唤pow更快。

最后，在main

 for(row=0; row 

您正在N/2行打印行列式，这看起来很不错，但此时您尚未扫描整个矩阵，因此N/2 + 1到N-1包含未初始化的数据，这不太可能全是零。 如果将if (row != N/2)更改为if (row != N-1) ，它将起作用，但是，处理它的正确方法是将矩阵的扫描与计算和打印分开。 这些都是独立的操作，应该在各自独立的function中处理。