P3389 【模板】高斯消元法 常规写法+约旦消元

把方程化简成   简化梯形矩阵 然后回带即可

#include<bits/stdc++.h>
using namespace std;
const int M=100+7;
double mp[M][M];
double ans[M];
double eps=1e-7;
int main()
{
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n+1;j++)
	scanf("%lf",&mp[i][j]);
	for(int i=1;i<=n;i++)//处理第i个方程
	{
		int r=i;
		for(int j=i+1;j<=n;j++)
			if(fabs(mp[r][i])<fabs(mp[j][i]))r=j;
		if(fabs(mp[r][i])<eps)//方程存在零系数，无解或无穷多解，反之方程化简后不存在零系数，一定存在唯一解 
		{
		//	cout<<mp[r][i]<<"  "<<r<<"  "<<i<<endl;
			printf("No Solution\n");
			return 0;
		 } 
		 if(i!=r)swap(mp[i],mp[r]);//每次都从剩下的方程种选取当前未知数系数最大的方程，优化精度
		 double div=mp[i][i];
		 for(int j=i;j<=n+1;j++)mp[i][j]/=div;
		 for(int j=i+1;j<=n;j+