高斯约旦消元法

求解线性方程组。

首先，我们要知道 \(\frac{x}{0.01}\) 和 \(\frac{x}{100}\) 的精度是 \(\frac{x}{100}\) 要大的，后面的数要小，所以精度大。

然后就是简单的模拟消元的过程了，高斯约旦消元法最后的解是一个对角线矩阵。

我构造的是一个对角线矩阵，所以循环从 \(i + 1\) 开始，然后我们可以找到最大值，找最大值是精度的保障。

AC code:

/*
	Work by: TLE_Automation
*/
#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
#define int long long
using namespace std;

const int N = 1e6 + 10;
const int MAXN = 2e5 + 10;

inline char readchar() {
	static char buf[100000], *p1 = buf, *p2 = buf;
	return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2) ? EOF : *p1++;
}

inline int read() {
#define readchar getchar
	int res = 0, f = 0;char ch = readchar();
	for(; !isdigit(ch); ch = readchar()) if(ch == '-') f = 1;
	for(; isdigit(ch); ch = readchar()) res = (res << 1) + (res << 3) + (ch ^ '0');
	return f ? -res : res;
}

inline void print(int x) {
	if (x < 0 ) putchar('-'), x = -x;
	if (x > 9 ) print(x / 10);
	putchar(x % 10 + '0');
}

double a[102][102];
int n;

void Gauss() {	
	for(int i = 1; i <= n; i++) {
		int Max = i;
		for(int j = i + 1; j <= n; j++) if(fabs(a[j][i]) > fabs(a[Max][i])) Max = j;
		for(int j = 1; j <= n + 1; j++) swap(a[i][j], a[Max][j]);
		if(a[i][i] == 0) puts("No Solution"), exit(0);
		for(int j = 1; j <= n; j++) {
			if(j == i) continue;
			double tmp = a[j][i] / a[i][i];
			for(int k = 1; k <= n + 1; k++) a[j][k] -= a[i][k] * tmp; 
		}
	}
}
signed main() {
	n = read();
	for(int i = 1; i <= n; i++) {
		for(int j = 1; j <= n + 1; j++) {
			scanf("%lf", &a[i][j]);
		}
	}
	Gauss();
	for(int i = 1; i <= n; i++) {
		printf("%.2lf\n",a[i][n + 1] / a[i][i]);
	}
}