高斯消元

高斯消元求解线性方程组

算法步骤：
枚举每一列c，找出每一列当中中绝对值最大的那一行，然后将这一行换到当前没有固定的行当中最上面的那一行（即当前所在的行，第一列最上面的那一行就是第一行，第二列就是第二行），再把这一行的第一个数变成1，然后把这一行下面的所有行的第c列变成0，即从这一行的最后一个数开始枚举，一直到当前所在的列，将每一个数都除以第一个数即可，（第一个行从第一列开始枚举，第二行从第二列开始枚举，以此类推…，只需要在每次枚举列的时候，r++即可）
​

输入：

3

1.00 2.00 -1.00 -6.00

2.00 1.00 -3.00 -9.00

-1.00 -1.00 2.00 7.00

输出：1.00

-2.00

3.00

代码实例：

#include<iostream>
#include<algorithm>
#include<iostream>
#include<algorithm>
#include<vector>
#include<set>
#include<stack>
#include<queue>
#include<string>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<ctype.h>
#include<iomanip>
#include<fstream>
using namespace std;
#define endl '\n'
#define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
typedef unsigned long long ULL;
typedef pair<int, int> PII;
typedef long long ll;
const int N = 1e3 + 10, null = 0x3f3f3f3f;
const double eps=1e-6;
double a[N][N];
int n;
int gaoss()
{
	int r,c;//r表示枚举的行，c表示枚举的列
	for(c=0,r=0;c<n;c++)//枚举每一列，c从第一列开始枚举，r从第一行开始枚举，每次循环r++。总之可以理解为第一列从第一行开始枚举，第二列从第二行开始枚举，即每次循环一次，r++
	{
		int t=r;//t表示每一列当中绝对值最大的那一行，初始化为r
		for(int i=r+1;i<n;i++)//枚举当前没有确定的行当中的所有行
		{
			if(fabs(a[i][c])>fabs(a[t][c])) t=i;//找到绝对值最大的那一行，更新t的值
		}
		if(fabs(a[t][c])<eps) continue;//如果绝对值最大的那一行的第c列值为0，则不用进行处理，跳出该循环
		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];//将这一行的第一个数变为1，方法为从右到左，将这一行的每一个数均除以第一个数
		for(int i=r+1;i<n;i++)//将下面所有行的第c列置为0
		{
			if(fabs(a[i][c])>eps)//如果该行的第c列不为0的话，则从这一行的最后一个数开始枚举，将其置为0
			{
				for(int j=n;j>=c;j--)
				{
					a[i][j]-=a[r][j]*a[i][c];
				}
			}
		}
		r++;//每次循环r++，保证第一列是从第一行开始寻找绝对值最大的行，第二列从第二行开始寻找绝对值最大的行。
	}
	if(r<n) //如果r小于n的话，说明要么方程组存在无穷多解，要么方程组无解，分为两种情况
	{
		for(int i=r;i<n;i++)//从第r行开始枚举，到最后一行
		{
			if(fabs(a[i][n])>eps) return 2;//如果该行的最后一列值不为0的话，即方程左边为0，右边不为0，说明方程组无解
		}
		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;
}
int main()
{
	ios;
	//输入初始化
	cin>>n;
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<=n;j++)
		{
			cin>>a[i][j];
		}
	}
	int t=gaoss();
	if(t==0)//唯一解的情况
	{
		for(int i=0;i<n;i++)
		{
			cout<<fixed<<setprecision(2)<<a[i][n]<<endl;//输出方程组的解，同时保留两位小数
		}
	}
	else if(t==1) cout<<"方程存在无穷多解"<<endl;
	else cout<<"方程组无解"<<endl;
	return 0;
}

​