【Vijos P1052】贾老二算算术高斯消元傻题

最新推荐文章于 2019-03-09 18:21:00 发布

空灰冰魂

最新推荐文章于 2019-03-09 18:21:00 发布

阅读量1.6k

点赞数 1

分类专栏：高斯消元文章标签： Vijos P1052 贾老二算算术高斯消元傻题

本文链接：https://blog.youkuaiyun.com/Vmurder/article/details/43445323

版权

高斯消元专栏收录该内容

10 篇文章

订阅专栏

本文详细解析了使用C++实现高斯消元法解决线性方程组的过程，包括算法原理、代码实现及测试验证，旨在为学习算法和C++编程的读者提供实用指导。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <stdio.h>
int main()
{
	puts("转载请注明出处谢谢");
	puts("http://blog.youkuaiyun.com/vmurder/article/details/43445323");
}

题解：

虽然是说有唯一解，但是依然需要在消到某一层时判0/非0，然后适当交换位置。

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 105
#define eps 1e-3
using namespace std;
int n;
double a[N][N],ans[N];
void Gauss()
{
	int i,j,k;
	for(i=1;i<=n;i++)
	{
		for(j=i;j<=n&&!a[j][i];j++);
		if(i!=j)for(k=i;k<=n+1;k++)swap(a[i][k],a[j][k]);
		for(j=i+1;j<=n;j++)
		{
			for(k=n+1;k>=i;k--)
			{
				a[j][k]-=a[i][k]/a[i][i]*a[j][i];
			}
		}
	}
	for(i=n;i;i--)
	{
		ans[i]=a[i][n+1];
		for(j=i+1;j<=n;j++)ans[i]-=a[i][j]*ans[j];
		ans[i]/=a[i][i];
	}
}
int main()
{
	freopen("test.in","r",stdin);
	int i,j,k;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
		for(j=1;j<=n+1;j++)
			scanf("%lf",&a[i][j]);
	Gauss();
	for(i=1;i<=n;i++)
	{
		printf("%d",(int)(ans[i]+eps));
		if(i!=n)printf(" ");
	}
	return 0;
}