【模板】高斯消元

---

CODE

#include<bits/stdc++.h>
using namespace std;
const double eps=1e-10;
double uu,a[52][52],b[52];
int n,l[52];bool pd;
inline void zzd(int &maxx,int i,int cnt){
	for(int j=cnt+1;j<=n;++j){//找系数最大值 
		if(fabs(a[j][i])>fabs(a[maxx][i]))
			maxx=j;
	}
	return ;
}
signed main(void){
	scanf("%d",&n);
	for(int i=1;i<=n;++i){
		for(int j=1;j<=n;++j){
			scanf("%lf",&a[i][j]);
		}
		scanf("%lf",b+i);
	}int cnt=1;
	for(int i=1;i<=n;++i){//枚举列 
		int maxx=cnt;
		zzd(maxx,i,cnt);
		if(fabs(a[maxx][i])<eps) continue;
		for(int j=1;j<=n;++j) swap(a[cnt][j],a[maxx][j]);//交换 
		swap(b[cnt],b[maxx]);
		uu=a[cnt][i];
		for (int j = i; j <= n; j++)
			a[cnt][j]/=uu;
		b[cnt]/=uu;
		for(int j=1;j<=n;++j){
			if(j!=cnt){
				uu=a[j][i];
				for(int k=i;k<=n;++k)
					a[j][k]-=a[cnt][k]*uu;
				b[j]-=b[cnt]*uu;
			}
		}
		l[cnt]=i;++cnt;
	}
	for(int i=1;i<=n;++i){
		if(fabs(a[i][l[i]])<eps){
			if(fabs(b[i])>=eps) return puts("-1"),0;
			pd=1;continue;
		}
	}
	if(pd) return puts("0")&0;
	double ans;
	for(int i=1;i<=n;++i)
	{
		ans=b[i]/a[i][i];
		if(ans==0) printf("x%d=0\n",i);
		else printf("x%d=%.2lf\n",i,ans);
	}
	return 0;
} 

__EOF__

本文作者：GOD_HJ
本文链接：https://www.cnblogs.com/GOD-HJ/p/17321232.html
关于博主：评论和私信会在第一时间回复。或者直接私信我。
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！