高斯消元法（n^3）顺序，主元，约当

模板题目

顺序消去法

//顺序高斯消元
//拿列主元随便改的
int n;
double a[111][111],ans[111];
signed main()
{
    cin>>n;
    rpp(i,n) rpp(j,n+1) cin>>a[i][j];

    for(int i=1;i<=n;++i)//枚举当前项
    {
        for(int j=i+1;j<=n;++j)
        {
            double tmp=a[j][i]/a[i][i];
            for(int k=i;k<=n+1;++k) // 第i次迭代的时候前i-1项肯定是0，直接从i开始
                a[j][k]-=a[i][k]*tmp;
        }
        // rpp(i,n) 
        // {
        //     rpp(j,n+1) 
        //     {
        //         cout<<a[i][j]<<" ";
        //     }
        //     cout<<endl;
        // }
        // cout<<endl;
    }
    for(int i=n;i>=1;--i)
    {
        ans[i]=a[i][n+1];
        for(int j=i+1;j<=n;++j)//他下面的都剪掉
            ans[i]-=a[i][j]*ans[j];
        ans[i]/=a[i][i];
    }
    rpp(i,n) printf("%.2f\n",ans[i]);
    cout<<endl;
    stop;
    return 0;
}

列主元消除法

//列主元消除法（加减消元+代入消元）
//和顺序消除法相比，只多了一个选择每项最大主元
//和全主元消除法相比，差别在于全主元是在右下角的矩阵内选择最大值作为主元
//列主元只在该列寻找
//判断解的情况
//对于无解的判断：某一行前 nnn 个数均为 0，最后的结果却不为 000。
//对于无数解的判断：某一行 n+1n+1n+1 个数均为 000.
int n;
double a[111][111],ans[111];
signed main()
{
    cin>>n;
    rpp(i,n) rpp(j,n+1) cin>>a[i][j];

    for(int i=1;i<=n;++i)//枚举当前项
    {
        //(1)找到主元所在行，换到当前行
        int mx=i;
        for(int j=i+1;j<=n;++j)//找到该项最大的那一个方程
            if(fabs(a[j][i])>fabs(a[mx][i])) mx=j;
        if(mx!=i) swap(a[mx],a[i]);//交换这两行
        if(!a[i][i]) 
        {
            cout<<"No Solution"<<endl;
            return 0;
        }
        //(2)将该行主元系数置为1
        double tmp=a[i][i];
        rpp(j,n+1) a[i][j]/=tmp;
        //(3)将该行以下的该项系数全部变成0
        for(int j=i+1;j<=n;++j)
        {
            tmp=a[j][i];
            for(int k=i;k<=n+1;++k) // 第i次迭代的时候前i-1项肯定是0，直接从i开始
                a[j][k]-=a[i][k]*tmp;
        }
    }


    //(4)回代
    for(int i=n;i>=1;--i)
    {
        ans[i]=a[i][n+1];
        for(int j=i+1;j<=n;++j)//他下面的都剪掉
            ans[i]-=a[i][j]*ans[j];
    }
    rpp(i,n) printf("%.2f\n",ans[i]);
    cout<<endl;
    stop;
    return 0;
}

全主元消除法

//全主元消除法
// 和列主元相比只改了（1），然后加一个找pos的就行
int n;
int pos[111];//i列为第pos[i]项
double a[111][111],ans[111];
signed main()
{
    cin>>n;
    rpp(i,n) pos[i]=i;
    rpp(i,n) rpp(j,n+1) cin>>a[i][j];

    for(int i=1;i<=n;++i)//枚举当前项
    {
        //(1)在右下角中查找主元
        int r=i,s=i;
        for(int _=i+1;_<=n;++_)
            for(int __=i+1;__<=n;++__)
                if(fabs(a[_][__])>fabs(a[r][s])) r=_,s=__;
        if(r!=i) swap(a[i],a[r]);//先交换行
        if(s!=i)//然后交换列，同时更新pos数组
        {
            swap(pos[s],pos[i]);
            rpp(_,n) swap(a[_][i],a[_][s]);
        }
        if(!a[i][i])
        {
            cout<<"No Solution"<<endl;
            return 0;
        }
        //(2)将该行主元系数置为1
        double tmp=a[i][i];
        rpp(j,n+1) a[i][j]/=tmp;
        //(3)将该行以下的该项系数全部变成0
        for(int j=i+1;j<=n;++j)
        {
            tmp=a[j][i];
            for(int k=i;k<=n+1;++k) // 第i次迭代的时候前i-1项肯定是0，直接从i开始
                a[j][k]-=a[i][k]*tmp;
        }
    }
    //(4)回代
    for(int i=n;i>=1;--i)
    {
        ans[i]=a[i][n+1];
        for(int j=i+1;j<=n;++j)//他下面的都剪掉
            ans[i]-=a[i][j]*ans[j];
    }
    map<int,int>mp;
    rpp(i,n) mp[pos[i]]=i;//第i项在那一列
    rpp(i,n) printf("%.2f\n",ans[mp[i]]);
    cout<<endl;
    return 0;
}

高斯–约当消元

//列主元高斯--约当消元（精度最高的（相对地）
//缺点：只能判断是否有唯一解，无法判断究竟是无解还是无数解。
// 1.选择一个尚未被选过的未知数作为主元，选择包含这个主元的方程。
// 2.将这个方程主元的系数化为1。
// 3.通过加减消元，消掉其它方程中的这个未知数。
// 4.重复以上步骤，直到把所有式子变成形如: a1+b0+c*0……=d 
double a[111][111];
signed main()
{
    int n;cin>>n;
    rpp(i,n) rpp(j,n+1) cin>>a[i][j];
    rpp(i,n)//枚举当前要消除的项（列）
    {
        //(1)找到一个主元
        int mx=i;
        for(int j=i+1;j<=n;++j)//找到该项最大的那行方程//每个方程只能有一个主元/小于i的不用比了
            if(fabs(a[j][i])>fabs(a[mx][i])) mx=j;
        rpp(j,n+1) swap(a[mx][j],a[i][j]);//交换这两行
        if(!a[i][i])//最大值等于0则说明该列都为0，肯定无解 
        {
            puts("No Solution");
            return 0;
        }
        //(2)将其他行中的该项消除
        rpp(j,n)
        {
            if(j==i) continue;
            double tmp=a[j][i]/a[i][i];
            for(int k=i;k<=n+1;++k)//i列即之前肯定已经消完了
                a[j][k]-=tmp*a[i][k];
        }
        // rpp(_,n) 
        // {
        //     rpp(__,n+1) printf("%6.2f",a[_][__]);
        //     cout<<endl;
        // }
        // cout<<endl;//输出中间值
    } 
    //(3)输出最终结果
    rpp(i,n)  printf("%.2f\n",a[i][n+1]/a[i][i]);
    stop;
    return 0;
}