luogu P3389 【模板】高斯消元法

题目传送门:https://www.luogu.org/problemnew/show/P3389



题意:

给你一个n元一次方程,求着n个未知数的解。



思路:

gauss消元。



代码:

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
	int n;
	double a[110][110];
void gauss()
{
	for(int i=1;i<=n;i++)
	{
		int p=i;
		for(int j=i+1;j<=n;j++)
			if(fabs(a[p][i])<fabs(a[j][i])) p=j;
		if(fabs(a[p][i])<1e-8)
		{
			printf("No Solution");
			return;
		}
		for(int j=1;j<=n+1;j++)
			swap(a[i][j],a[p][j]);
		for(int j=i+1;j<=n;j++)
		{
			double d=a[j][i]/a[i][i];
			for(int k=i;k<=n+1;k++)
				a[j][k]-=d*a[i][k];
		}		
	}
	for(int i=n;i>=1;i--)
	{
		for(int j=i+1;j<=n;j++)
			a[i][n+1]-=a[j][n+1]*a[i][j];
		a[i][n+1]/=a[i][i];
	}
	for(int i=1;i<=n;i++)
		printf("%.2lf\n",a[i][n+1]);
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n+1;j++)
			scanf("%lf",&a[i][j]);
	gauss();
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值