7

我是编程新手，我一直在寻找一种方法来找到矩阵的行列式。我在网上找到了这段代码，但我无法理解这里的算法。我对 recursion 的基础没有任何问题，但是我无法理解 continue 和 main 循环。非常感谢任何可以向我解释算法的人。

int determ(int a[MAX][MAX],int n) {
  int det=0, p, h, k, i, j, temp[MAX][MAX];
  if(n==1) {
    return a[0][0];
  } else if(n==2) {
    det=(a[0][0]*a[1][1]-a[0][1]*a[1][0]);
    return det;
  } else {
    for(p=0;p<n;p++) {
      h = 0;
      k = 0;
      for(i=1;i<n;i++) {
        for( j=0;j<n;j++) {
          if(j==p) {
            continue;
          }
          temp[h][k] = a[i][j];
          k++;
          if(k==n-1) {
            h++;
            k = 0;
          }
        }
      }
      det=det+a[0][p]*pow(-1,p)*determ(temp,n-1);
    }
    return det;
  }
}
4

2 回答 2

7

该算法使用分治法来解决问题（找到 N*N 矩阵的行列式）。

该算法使用递归模式，这是一种分而治之的方法。您可以通过注意算法在第三个条件语句中调用自身来发现这一点。

每个递归算法都有一个退出条件，这是代码中的第一个 if 语句。它们还包含一个部分，该部分是最方便的问题的解决方案，或者是最初难以解决的主要大问题的原子问题。原子问题或最分裂的问题可以很容易地解决，因为您可以看到代码的第二个 if 语句。在您的情况下，它实际上是在解决 2*2 矩阵的行列式。

代码中最重要的部分要理解哪个也有一点挑战性，那就是你进行除法的部分（这也是递归的！）。这部分是征服的关键。通过做一些回溯和数值示例，您可以找到它：

det = det + a[0][p] * pow(-1,p) * determ(temp,n-1);

最后的建议是尝试一个 3*3 矩阵，它只需要一个除法。祝你好运。

这本书是开始学习和理解算法的好书

于 2014-01-19T18:34:22.330 回答
0 
#include <iostream>

using std::cin;
using std::cout;
using std::endl;

int **submatrix(int **matrix, unsigned int n, unsigned int x, unsigned int y) {
    int **submatrix = new int *[n - 1];
    int subi = 0;
    for (int i = 0; i < n; i++) {
        submatrix[subi] = new int[n - 1];
        int subj = 0;
        if (i == y) {
            continue;
        }
        for (int j = 0; j < n; j++) {
            if (j == x) {
                continue;
            }
            submatrix[subi][subj] = matrix[i][j];
            subj++;
        }
        subi++;
    }
    return submatrix;
}

int determinant(int **matrix, unsigned int n) {
    int det = 0;
    if (n == 2) {
        return matrix[0][0] * matrix[1][1] - matrix[1][0] * matrix[0][1];
    }
    for (int x = 0; x < n; ++x) {
        det += ((x % 2 == 0 ? 1 : -1) * matrix[0][x] * determinant(submatrix(matrix, n, x, 0), n - 1));
    }

    return det;
}

int main() {
    int n;
    cin >> n;
    int **matrix = new int *[n];
    for (int i = 0; i < n; ++i) {
        matrix[i] = new int[n];
        for (int j = 0; j < n; ++j) {
            cin >> matrix[i][j];
        }
    }

    cout << determinant(matrix, n);

    return 0;
}
于 2021-03-22T10:39:55.803 回答