高斯-约旦消元法、列主元消去法及其差异

最新推荐文章于 2024-05-03 21:03:50 发布

小鱼yn

最新推荐文章于 2024-05-03 21:03:50 发布

阅读量9.5k

点赞数 8

分类专栏：小东西

本文链接：https://blog.csdn.net/qq_43803508/article/details/115601829

版权

文章目录

高斯消元法求解线性方程组：

对增广矩阵做相应的初等行变换，将右边系数矩阵化为单位矩阵后，最右边一列就是原线性方程组的解。

例题：P3389 【模板】高斯消元法

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