2022春每日一题：Day 12

题目：[SDOI2006]线性方程组

显然，高斯消元模板题

代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 105
#define eps 1e-6
using namespace std;
int n;
double a[N][N];
int guass()
{
	int r=0;
	for(int c=0;c<n;c++)
	{
		int t=r;
		for(int i=r;i<n;i++)
		    if(fabs(a[i][c])>fabs(a[t][c]))
		        t=i;
		if(fabs(a[t][c])<eps)
		    continue;
		for(int i=c;i<=n;i++)
		    swap(a[t][i],a[r][i]);
		for(int i=n;i>=c;i--)
		    a[r][i]/=a[r][c];
		for(int i=r+1;i<n;i++)
		    if(fabs(a[i][c])>eps)
		        for(int j=n;j>=c;j--)
		            a[i][j]-=a[r][j]*a[i][c];
		++r;
	}
	if(r<n)
	{
		for(int i=r;i<n;i++)
		    if(fabs(a[i][n])>eps)
		        return 2;
		return 1;
	}
	for(int i=n-1;i>=0;i--)
	    for(int j=i+1;j<n;j++)
	        a[i][n]-=a[i][j]*a[j][n];
	return 0;
}
int main()
{
	scanf("%d",&n);
	for(int i=0;i<n;i++)
	    for(int j=0;j<=n;j++)
	        scanf("%lf",&a[i][j]);
	int t=guass();
	if(!t)
	{
		for(int i=0;i<n;i++)
		{
			if(fabs(a[i][n])<eps)
			    a[i][n]=fabs(a[i][n]);
			printf("x%d=%.2lf\n",i+1,a[i][n]);
		}
	}
	else if(t==1)
	    puts("0");
	else
	    puts("-1");
	return 0;
}