luogu3389 【模板】高斯消元法

这里写链接内容
题目背景

Gauss消元
题目描述

给定一个线性方程组，对其求解
输入输出格式

输入格式：

第一行，一个正整数 nnn

第二至 n+1n+1n+1行，每行 n+1n+1n+1 个整数，为a1,a2⋯an a_1, a_2 \cdots a_na1​,a2​⋯an​ 和 bbb，代表一组方程。

输出格式：

共n行，每行一个数，第 iii行为 xix_ixi​ （保留2位小数）

如果不存在唯一解，在第一行输出”No Solution”.
输入输出样例
输入样例#1： 复制

3
1 3 4 5
1 4 7 3
9 3 2 2

输出样例#1： 复制

-0.97
5.18
-2.39

说明

1≤n≤100,∣ai∣≤104,∣b∣≤1041 \leq n \leq 100, \left | a_i \right| \leq {10}^4 , \left |b \right| \leq {10}^4 1≤n≤100,∣ai​∣≤104,∣b∣≤104

再次熟悉下高斯消元 辣鸡蒟蒻我的一些理解http://www.elijahqi.win/2017/12/23/bzoj1013-jsoi2008%e7%90%83%e5%bd%a2%e7%a9%ba%e9%97%b4%e4%ba%a7%e7%94%9f%e5%99%a8sphere/

#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
double matrix[110][110];
int n;
bool gauss(){
    int now=1,to;double t;
    for(int i=1;i<=n;++i){
        for (to=now;to<=n;++to) if (fabs(matrix[to][i])) break;
        if (to>n) continue;if (to!=now) for (int j=1;j<=n+1;++j) swap(matrix[to][j],matrix[now][j]);
        t=matrix[now][i];for (int j=1;j<=n+1;++j) matrix[now][j]/=t;
        for (int j=1;j<=n;++j){
            if(j==now) continue;t=matrix[j][i];
            for (int k=1;k<=n+1;++k) matrix[j][k]-=t*matrix[now][k];
        }
        ++now; 
    }if (now<=n) return 0;else return 1;
}
int main(){
    freopen("3389.in","r",stdin);
    scanf("%d",&n);
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n+1;++j) scanf("%lf",&matrix[i][j]);
    if(!gauss()) {printf("No Solution");return 0;}
    for (int i=1;i<=n;++i) printf("%.2f\n",matrix[i][n+1]);
    return 0;
}