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;}