社区讨论

90pts WA on #3求助

P2455[SDOI2006] 线性方程组参与者 2已保存回复 2

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
2 条
当前快照
1 份
快照标识符
@lo7y8y7p
此快照首次捕获于
2023/10/27 09:43
2 年前
此快照最后确认于
2023/10/27 09:43
2 年前
查看原帖
CPP
#include<iostream>
#include<cstdio>
#include<cmath>
#define int long long
using namespace std;
int n;
double a[101][102];
const double eps=1e-7;
signed main()
{
	scanf("%lld",&n);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n+1;j++) scanf("%lf",&a[i][j]);
	}
	for(int i=1;i<=n;i++)
	{
		int maxn=i;
		for(int j=i+1;j<=n;j++) 
		{
			if(fabs(a[j][i])>fabs(a[maxn][i])) maxn=j;
		}
		for(int j=1;j<=n+1;j++) swap(a[i][j],a[maxn][j]);
		if(fabs(a[i][i])<eps)
		{
			continue;
//			if(fabs(a[i][n+1])<eps) printf("0");
//			else printf("-1");
//			return 0;
		}
		for(int j=1;j<=n;j++) 
		{
			if(i!=j) 
			{
				double K=a[j][i]/a[i][i];
				for(int k=i+1;k<=n+1;k++) a[j][k]-=a[i][k]*K;
			}
		}
	}
	for(int i=1;i<=n;i++) 
	{
		if(fabs(a[i][i])<eps)
		{
			if(fabs(a[i][n+1])>eps&&fabs(a[i][i])<eps) printf("-1");
			else if(fabs(a[i][n+1])<eps&&fabs(a[i][i])<eps) printf("0");
			return 0;
		}
	}
	for(int i=1;i<=n;i++)
	{
		if(fabs(a[i][n+1]/a[i][i])<eps) printf("x%lld=0.00\n",i);
		else printf("x%lld=%.2lf\n",i,a[i][n+1]/a[i][i]);
	}
	return 0;
} 

回复

2 条回复,欢迎继续交流。

正在加载回复...