凡人升天传5——【模板】高斯消元法

题目来源：【模板】高斯消元法 - 洛谷

题目背景

Gauss 消元

本题数据很水，且在题目的数据范围下可以构造数据使高斯消元产生极大的精度误差，因此本题仅供参考。

题目描述

给定一个线性方程组，对其求解

输入格式

第一行，一个正整数 nn

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

输出格式

共 nn 行，每行一个数，第 ii 行为 x_ixi​ （保留 2 位小数）

如果不存在唯一解，在第一行输出 No Solution.

输入输出样例

输入 #

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

输出 #

-0.97
5.18
-2.39

题解

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll read(){
	ll ans = 0,f = 1;
	char ch = getchar();
	while(!isdigit(ch)){
		f *= (ch == '-') ? -1 : 1;
		ch = getchar();
	}
	while(isdigit(ch)){
		ans = (ans << 3) + (ans << 1) + (ch ^ 48);
		ch = getchar();
	}
	return ans * f;
}
const int MAXN = 1e4+5;
const double eps = 1e-7;
double g[MAXN][MAXN];
int n;
int main(){
    n = read();
    for (register int i = 0;i < n; i ++){
        for (register int j = 0;j <= n; j ++)
            g[i][j] = read();
    }
    for (register int i = 0;i < n; i ++){
        int p = i;
        for (register int j = i;j < n; j ++)
            if (fabs(g[j][i]-g[p][i]) <= eps) 	
				p = j;
        for (register int j = 0;j <= n; j ++){
            double t = g[i][j];
            g[i][j] = g[p][j];
            g[p][j] = t;
        }
        if (fabs(g[i][i]) <= eps){
            printf("No Solution");
            return 0;
        }
        for (register int j = i + 1;j <= n; j ++) 	
			g[i][j] /= g[i][i];
        for (register int j = 0;j < n; j ++)
            if (i != j)
                for (register int k = i + 1;k <= n; k ++) 
					g[j][k] -= g[j][i] * g[i][k];
    }
    for (register int i = 0;i < n; i ++) 
		printf("%.2lf\n",g[i][n]);
    return 0;
}