线性方程组[高斯消元]

​​传送门​​

高斯消元模板

#include<bits/stdc++.h>
#define N 150
using namespace std;
double a[N][N]; int n;
void gauss(){
  for(int i=1;i<=n;i++){
    for(int j=i;j<=n;j++)
      if(fabs(a[j][i])>1e-8)
        for(int k=1;k<=n+1;k++) swap(a[i][k],a[j][k]);
    if(fabs(a[i][i])<1e-8 && fabs(a[i][n+1])>1e-8){printf("-1"); return;}
    if(fabs(a[i][i])<1e-8 && fabs(a[i][n+1])<1e-8){printf("0"); return;}
    for(int j=1;j<=n;j++){
      if(i==j) continue;
      double now = a[j][i] / a[i][i];
      for(int k=1;k<=n+1;k++) a[j][k] -= a[i][k] * now;
    }
  }
  for(int i=1;i<=n;i++) printf("x%d=%0.2lf\n",i,a[i][n+1]/a[i][i]);
}
int main(){
  cin>>n;
  for(int i=1;i<=n;i++)
    for(int j=1;j<=n+1;j++)
      scanf("%lf",&a[i][j]);
  gauss(); return 0;
}