洛谷 P3389 高斯消元模板题

最新推荐文章于 2020-02-07 15:21:14 发布

shiyongyang

最新推荐文章于 2020-02-07 15:21:14 发布

阅读量256

点赞数

分类专栏：高斯消元

本文链接：https://blog.csdn.net/shiyongyang/article/details/78113744

版权

高斯消元专栏收录该内容

2 篇文章 0 订阅

订阅专栏

高斯消元模板（求多元一次线性方程的解）

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define N 105
using namespace std;
int n;
double date[N][N],ans[N];
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
	    for(int j=1;j<=n+1;j++)
	    {
	    	scanf("%lf",&date[i][j]);
	    	date[i][j]*=1.00;//转化为浮点数 
		}	
		double tmp=date[i][1];
		for(int j=1;j<=n+1;j++)
		    date[i][j]=date[i][j]*1.0/tmp;//将第一项化为1 
	}
	for(int i=1;i<=n;i++)
	{
		if(date[i][i]==0)
		{
			printf("No Solution\n");return 0;
		}
		for(int j=i+1;j<=n+1;j++) date[i][j]/=date[i][i];//将对角线化为1 
		date[i][i]=1.0;
		for(int j=i+1;j<=n;j++)//将下三角化为0 
		{
			for(int k=i+1;k<=n+1;k++)
			     date[j][k]-=date[i][k]*date[j][i];
			date[j][i]=0;
		}  
	}
	for(int i=n;i;i--)//求ans； 
	{
		double opt=date[i][n+1];
		for(int j=i+1;j<=n;j++)
		     opt-=ans[j]*date[i][j];
		ans[i]=opt;
	}
	for(int i=1;i<=n;i++) printf("%.2lf\n",ans[i]);
	return 0;
 }