把矩阵化成行最简型（伪）

三种初等行变换

因为最后显示出来矩阵行还是乱的（雾
所以是伪
但是俺觉得已经不影响使用了呢（笑

输入：

行数 列数
矩阵

#include<iostream>
#include<algorithm>
using namespace std;
double mt[100][100];
int n,m;
void num_multiply_line(int a,double x)
{
	for(int j=0;j<m;j++)
		mt[a][j]*=x;
	return ;
}
void line_plus_line(int a,double x,int b)
{
	for(int j=0;j<m;j++)
		mt[a][j]+=x*mt[b][j];
	return ;
}
int ffnz(int a)
{
	for(int j=0;j<m;j++)
		if(mt[a][j]!=0)
			return j;
	return m;
}
void solve()
{
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<m;j++)
			printf("%.2f  ",mt[i][j]+1e-6);
		cout<<endl;
	}
	cout<<endl;
	for(int i=0;i<n;i++)  //遍历每一行 
	{
		int j=ffnz(i);
		if(j==m)
			continue;
		for(int k=i+1;k<n;k++)   //让首非零元这一列，他下面的全变成零，首非零元为mt[i][j],他下面每一行都进行 line_plus_line
			line_plus_line(k,-mt[k][j]/mt[i][j],i);
	}
	for(int i=0;i<n;i++)
	{
		int j=ffnz(i);
		if(j==m)
			continue;
		if(mt[i][j]!=1)
			num_multiply_line(i,1/mt[i][j]);
		for(int k=0;k<n;k++)
			if(mt[k][j]!=0&&k!=i)
				line_plus_line(k,-mt[k][j]/mt[i][j],i);
	}
	return ;
}
int main()
{
	cin>>n>>m;
	for(int i=0;i<n;i++)
		for(int j=0;j<m;j++)
			scanf("%lf",&mt[i][j]);
	solve();
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<m;j++)
			printf("%.2f  ",mt[i][j]+1e-6);
		cout<<endl;
	}
}
/*
3 4
1 0 2 -1
2 0 3 1
3 0 4 3

4 5
2 1 8 3 7
2 -3 0 7 -5
3 -2 5 8 0
1 0 3 2 0
*/