高斯消元法

题目传送门sx
高斯消元原理
x+2y+3z=11①
2x+y+4z=11②
3x+4y+z=19③
⑴
用①式分别把②③式中的x消去；
再用②式把③式中的y消去；
⑵
得到的③式可直接得到z的值；
然后再把z带到②式，得到y；
再把y，z带到①式；
⑴
x+2y+3z=11①
-3y-2z=-11②
-20/3z=-20/3（懒得弄分数了）
⑵
回带，解得
x=2,y=3,z=1;
放到代码中用同样的操作原理
这是上面那道题的代码

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int n;
double a[111][111],ans[111];
const double eps=1e-9;
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n+1;j++)
	    scanf("%lf",&a[i][j]);
	for(int i=1;i<=n;i++){
		int r=i;
		for(int j=i+1;j<=n;j++)
		  if(fabs(a[r][i])<fabs(a[j][i]))//减小误差，找到最大的
		    r=i;
		if(fabs(a[r][i])<eps) {cout<<"No Solution";return 0;}
		if(i!=r) swap(a[r],a[i]);//交换一下
		double div=a[i][i];
		for(int j=i;j<=n+1;j++)
		  a[i][j]/=div;//为了方便下面倍数关系，将i处未知数系数化为1
		for(int j=i+1;j<=n;j++){
			div=a[j][i];            //该式是i式的div倍；
			for(int k=i;k<=n+1;k++)//两式做减法运算
			  a[j][k]-=div*a[i][k];
		}
	}
	//回带
	ans[n]=a[n][n+1];
	for(int i=n-1;i;i--){
		ans[i]=a[i][n+1];
		for(int j=i+1;j<=n;j++)
		  ans[i]-=a[i][j]*ans[j];
	}
	for(int i=1;i<=n;i++)
	  printf("%.2lf\n",ans[i]);
	return 0;
}