高斯消元

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int eps=1e-6;// 因为double储存有误差所以用eps;
int n;double a[110][110];
int gauss();
int main()
{
	scanf("%d",&n);
	for(int i=0;i<n;i++)
		for(int j=0;j<n+1;j++)
			scanf("%lf",&a[i][j]);
			
	int t=gauss();
	if(t==1){
		printf("No Solution");
	}else{
		for(int i=0;i<n;i++){
			printf("%0.2lf\n",a[i][n]);
		}
	}
	///
	return 0;
}
/
int gauss(){
	int c,r;//C表示每一列   R表示每一行
	for(c=0,r=0;c<n;c++){
		int t=r;
		for(int i=r;i<n;i++)
			if(fabs(a[i][c])>fabs(a[t][c]))
				t=i;// 每次都去找到最大的
		if(fabs(a[t][c])<eps)continue;//  此处则说明该元系数为0；就可能存在吴姐或无数姐情况
		/*
			如果说三个方程；其中一个未知数的系数都为0；
			以下两种可能：
				1.舍去那个未知数，如果说对于剩下的两个未知数方程有唯一解则出现了0x=0；
				  就是无数姐；
				2.舍去之后因为那啥吴姐，所以吴姐；
		*/
		for(int i=c;i<=n;i++)swap(a[t][i],a[r][i]);//  交换
		for(int i=n;i>=c;i--)a[r][i]/=a[r][c];// 系数归一
		for(int i=r+1;i<n;i++)
			if(fabs(a[i][c])>eps)// 如果不为0
				for(int j=n;j>=c;j--)
					a[i][j]-=a[r][j] * a[i][c];//  那就变成0
		r++;//下一行
	}
	
	if(r<n)return 1;// 吴姐或无穷解
	for(int i=n-1;i>=0;i--)
		for(int j=i+1;j<n;j++)
			a[i][n]-=a[i][j]*a[j][n];   //结果是一个阶梯，根据梯形来处理；
	return 0;
}

这个模板是徐倞舟老师教的；感谢徐倞老师；