高斯消元与行列式求值 part1

时间:2023-03-08 18:16:41

两道模板题,思路与算法却是相当经典。

先说最开始做的行列式求值,题目大致为给一个10*10的行列式,求其值

具体思路(一开始看到题我的思路):

1.暴算,把每种可能组合试一遍,求逆序数,做相应加减运算,一看就知道不是正解。

2.暴算2.0 用递归和代数余子式计算,但同样需计算逆序数。

3.一看就是知道是正解高斯消元,将行列式消为上三角,将对角线相乘。

看看代码:

#include<bits/stdc++.h>      //喜闻乐见的万用头
using namespace std;
int n;
double a[][]; //因需处理小数,开double
int t=;
bool no;
double tim=; //储存约去系数int main(){
cin>>n;
for(int i=;i<=n;i++){
for(int j=;j<=n;j++) scanf("%lf",&a[i][j]); //由于数据较少(10*10),不用快读
}
for(int i=;i<=n;i++){
if(!a[i][i]){ //判断当前处理数据是否为0,point
t=i;
while(!a[t][i]&&t<=n) t++;
if(t==n+){no=; continue;} //筛选首项不为0的行,进行交换,别忘记*=-1,最少wa4个
for(int j=;j<=n;j++) swap(a[i][j],a[t][j]);
tim*=-;
}
double x=a[i][i]; //将该行非0首项(之前的运算保证其不为0且为第一的非0数)
for(int j=i;j<=n;j++) a[i][j]/=x; //将每个a[i][i]除为1
tim*=x; //将消去系数累乘保存
for(int j=i+;j<=n;j++){
x=a[j][i];
for(int k=;k<=n;k++) a[j][k]-=x*a[i][k]; //将每行首项消为0,其余数做相同运算
}
/*for(int j=1;j<=n;j++){ //输出看看行列式消得对不对
for(int k=1;k<=n;k++) cout<<a[j][k]<<" ";
cout<<endl;
}*/
}
printf("%0.0lf",tim*a[n][n]); //输出系数与最后一项的乘积,取整
return ;
}
p.s.:“no”好像在这里没有用,不知为啥就写上了。。。
//以下是例题没有的特殊数据
/*8
5 5 10 9 5 9 10 4
3 3 6 1 4 6 7 10
1 1 2 10 8 9 8 7
6 2 3 4 8 3 6 9
1 2 6 3 2 7 8 9
9 5 4 5 1 7 3 10
2 4 6 10 10 5 7 8
4 5 6 10 4 7 5 2*/

point:(以后把需要较大量文字叙述的重点用文下注释解释)

如果a[i][i]为0,则在将要进行的化简运算中会出现  n/0 情况,

对于此类数字输出为“nan”(同学说这是暗示这道题难。。。)“not a number”

而对于此处出现的0,根据高斯消元的相关理论应该往下换,毕竟要组成下三角需把0“沉下去”,所以与下面换行。

p.s.:这道题给了我一点启示:

当初同学最快做出来得了11分,我初次尝试得了64分,这给了我做出来的动力,非常大的动力

我便梦想总有一天我要在全56级初学者之前a掉它,自此,我在家打开的窗口不再是虐杀原形,而是c++,为此,我改变了3次算法,代码重构好多次,提交。。。10几次是有了吧,乘积由64提到95,再提到96,最后查出致命错误终是a了,并达成了梦想。。。不总结了,只是希望以后万念俱灰的自己看到这个能。。。cheer下吧