高斯消元模板

 

　 推荐一个讲得不错的博客：http://www.cnblogs.com/ECJTUACM-873284962/p/6880199.html　　

  洛谷模板题：https://www.luogu.org/problemnew/show/P3389

 

模板题代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define LL long long
#define RI register int
#define eps 1e-7 
using namespace std;
const int INF = 0x7ffffff ;
const int N = 100 + 10 ;

inline int read() {  // 快读 
    int k = 0 , f = 1 ; char c = getchar() ;
    for( ; !isdigit(c) ; c = getchar())
      if(c == '-') f = -1 ;
    for( ; isdigit(c) ; c = getchar())
      k = k*10 + c-'0' ;
    return k*f ;
}
int n ; double hh[N][N] ;

inline bool guass() {
    for(int i=0;i<n;i++) {  // 转换成上三角矩阵过程 
        int r = i ; double maxx = hh[i][i] ;
        for(int j=i+1;j<n;j++) {
            if(hh[j][i] > hh[r][i]) r = j ;
        }
        if(fabs(hh[r][i]) < eps) return 0 ;  // 因为有精度误差，所以除了x1有无数解时是严格等于0的，其他时候都只是接近0而已 
        if(r != i) for(int j=i;j<=n;j++) swap(hh[i][j],hh[r][j]) ;
        for(int j=i+1;j<n;j++) {
            double f = hh[j][i]/hh[i][i] ;
            for(int k=i;k<=n;k++) {
                hh[j][k] -= hh[i][k]*f ;
            }
        }
    }
    for(int i=n-1;i>=0;i--) {  // 回带过程 
        for(int j=i+1;j<n;j++) hh[i][n] -= hh[i][j]*hh[j][n] ;
        hh[i][n] /= hh[i][i] ;
    }
    return 1 ;
}

int main() {
    n = read() ;
    for(int i=0;i<n;i++) {
        for(int j=0;j<=n;j++)
          hh[i][j] = read() ;
    }
    if(!guass()) {
        printf("No Solution") ;
    } else {
        for(int i=0;i<n;i++) printf("%.2lf
",hh[i][n]) ;
    }
    return 0 ;
}

  大家看着代码照着样例手动模拟一下过程就能理解了。 （实践出真知啊qwq）