题目：[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;

}

巴特西

2022春每日一题：Day 12

题目：[SDOI2006]线性方程组

最新文章

热门文章