高斯消元法求解线性方程组:
对增广矩阵做相应的初等行变换,将右边系数矩阵化为单位矩阵后,最右边一列就是原线性方程组的解。
1、高斯-约旦(G-J)消元法
直接消成对角矩阵
(1)从第一行开始,先顺序寻找到主元(不能为0),并将主元通过行变换移动到主对角线上
(2)将主元所在的行内的所有元素除以主元,使得主元化为1(也可以放在第四部之后)
(3)其他行减去 主元所在行 乘以一定的倍数,使得主元所在的列内除主元外的其他元素为0
(4)在剩下的矩阵范围内寻找下一个主元,重复以上步骤
#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
using namespace std;
//#define int long long
#define ll long long
#define uql unsigned long long
#define pii pair<int,int>
#define mid ((l + r)>>1)
#define chl (root<<1)
#define chr (root<<1|1)
#define lowbit(x) ( x&(-x) )
const int manx = 100+10;
const int manx2 = 4e7 + 10;
const ll INF = 1e18;
const ll mod = 1000000007;
const double eps=1e-10;
//运算过程中精度丢失导致本来应当为0的值的机器表示不为0,所以用abs(x)<eps 来判0
int n;
double f[manx][manx];
int Gauss()
{
for(int i=1;i<=n;i++){
for(int j=i;j<=n;j++){
/********对应步骤(1)********/
if(fabs(f[j][i])>eps){
swap(f[i],f[j]);
break;
}
}
if