538: Gaussian elimination

解法：

#include <bits/stdc++.h>
using namespace std;

const int N=1e3+3;
double a[N][N];

int main(int argc, char** argv) {
	int n,m;
	cin>>n>>m;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
			cin>>a[i][j];
	
	for (int k=1;k<=n;k++){
		for (int i=1;i<=n;i++){
			if (i==k) continue;
			double p=a[i][k]/a[k][k];
			for (int j=1;j<=m;j++){
				a[i][j]-=a[k][j]*p;
			}
		}
	}
	for (int i=1;i<=n;i++){
		printf("%.2lf ",a[i][m]/a[i][i]);
	}
	return 0;
}