数论--高斯消元

高斯消元

时间复杂度O(n^3)
此图来自Acwing灰之魔女

算法步骤：

枚举每一列c；
1.找到当前列绝对值最大的一行
2. 把这一行换到最上面（未确定阶梯型的行，并不是第一行）
3.将该行的第一个数变成 1 （其余所有的数字依次跟着变化）
4. 将下面所有行的当且列的值变成 0；

注释没有引起来，看着舒服点；
#include<iostream>
#include <cmath>
using namespace std;
const int N=1010;
const double age=1e-6;
int n;
double a[N][N];
int gauss()
{
    int c,r;
    for(c=0,r=0;c<n;c++)
    {
        int t=r;
        for(int i=r+1;i<n;i++)
        {
            if(fabs(a[c][i])>fabs(a[c][t]))
            {
                t=i;
                break;
            }
        }
寻找最大的数值是因为可以避免系数变得太大，精度较高

        if(fabs(a[t][c])<age) continue;
        
如果最大值是零,便可以直接跳过了
 之所以要小于1e-6，是因为c++浮点数的一种弊端，所以小于egs时，可以近似的看作是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];
        }
        
将当前第c位的值初始化为1，当然其他数也同样需要除以第c位，之所以从第n位开始，是为了不覆盖掉第c位的值
        for(int i=r+1;i<n;i++)
        {
            if(fabs(a[i][c])>age)
如果当前行数为零，说明不用进行削为零的操作
            {
                for(int j=n;j>=c;j--)
                {
                    a[i][j]-=a[i][c]*a[r][j];
                }
            }
        }
        r++;
    }
    if(r<n)
    {
        
此时c已经到了最后一列，所以当前行下面的所有数全是零，所以左边的系数矩阵全是零，
右边的只要有不为零的就是0=x无解，要都是零的话就是无穷解；
        for(int i=r;i<n;i++)
        {
            if(fabs(a[i][n])>age)
            {
                return 2;
            }
        }
        return 1;
    }
    for(int i=n-1;i>=0;i--)
    {
        for(int j=i+1;j<n;j++)
        {
这里是第i行第i+1个数往后的都要变成零，这样a[i][j]的值就是a[i][n]了；
因为a[j][j]的值都算过了是a[j][n],所以等式如下
            a[i][n]-=a[j][n]*a[i][j];
        }
    }
    return 0;
}
int main()
{ 
    cin >> n;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<=n;j++)
        {
            cin >> a[i][j];
        }
    }
    int t=gauss();

    if(t==0)
    {
        for(int i=0;i<n;i++)
        {
            if(fabs(a[i][n])<age) a[i][n]=0;
          printf("%.2f\n",a[i][n]);
        }
    }
    else if(t==1)
    {
        puts("Infinite group solutions");
    }
    else 
        puts("No solution");
    return 0;
}

第二种：异或高斯消元(同理)

#include<iostream>
#include <cmath>
using namespace std;
const int N=119;
int n;
int a[N][N];
int guase()
{
    int r,c;
    for(c=0,r=0;c<n;c++)
    {
        int t=r;
        for(int i=r;i<n;i++)
        {
            if(a[i][c])
            {
                t=i;
                break;
            }
        }
        if(a[t][c]==0) continue;
        for(int i=c;i<=n;i++)
        {
            swap(a[t][i],a[r][i]);
        }
//a^b^c=d;
//a=d^b^c;
        for(int i=r+1;i<n;i++)
        {
            if(a[i][c])
           {
               for(int j=c;j<=n;j++)
               {
                   a[i][j]^=a[r][j];
               }
           }
        }
        r++;
    }
    if(r<n)
    {
        for(int i=r;i<n;i++)
        {
            if(a[i][n])
            {
                return 2;
            }
            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()
{
    cin >> n;
    for(int i=0;i<n;i++)
      for(int j=0;j<=n;j++)
      {
          cin >> a[i][j];
      }
    int t=guase();
   
    if(t==0)
    {
        for(int i=0;i<n;i++) 
        { 
            cout << a[i][n] << endl;
        }
    }
    else if(t==1)
    {
        puts("Multiple sets of solutions");
    }
    else 
    {
        puts("No solution");
    }
    return 0;
}

加个小图片