我是编程新手,我一直在寻找一种找到矩阵行列式的方法。我在网上找到了这段代码,但我很难理解这里的算法。我对递归的基础没有问题,但继续和主循环我很难理解。非常感谢任何可以向我解释该算法的人。
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;
}
}
该算法使用分而治之的方法来解决问题(找到 N*N 矩阵的行列式)。
该算法使用递归模式,这是分而治之的方法之一。您可以通过注意算法在第三个条件语句中调用自身来发现这一点。
每个递归算法都有一个退出条件,它是代码中的第一个 if 语句。它们还包含一个部分,它是最方便问题的解决方案,或者是主要大问题的原子问题,而这些问题一开始就很难解决。原子问题或最分割的问题可以很容易地解决,因为您可以看到代码的第二个 if 语句。在你的例子中,它实际上是在求解 2*2 矩阵的行列式。
代码中最重要的部分是进行除法的部分(这也是递归的!),这也有点具有挑战性。
这部分是征服其中任何一个的关键。通过做一些回溯和数值示例,您可以找到答案:
det = det + a[0][p] * pow(-1,p) * determ(temp,n-1);
对于最后的建议,尝试使用仅需要一次除法的 3*3 矩阵。
祝你好运。
这本书非常适合开始学习和理解算法 http://mitpress.mit.edu/books/introduction-algorithms
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)