高斯消元+完全主元素+若尔当方法


#include<bits/stdc++.h>

using namespace std;

const int maxn=50;

int n;
double a[maxn][maxn];
double ans[maxn];

void swap_r(int q,int p)
{
    for(int i=1;i<=n+1;i++){
        double t=a[p][i];
        a[p][i]=a[q][i];
        a[q][i]=t;
    }
}

void swap_c(int q,int p)
{
    for(int i=1;i<=n+1;i++){
        double t=a[i][p];
        a[i][p]=a[i][q];
        a[i][q]=t;
    }
}

void gs()
{
    for(int i=1;i<n;i++){
        double m=fabs(a[i][i]);
        int p=i,q=i;
        for(int j=i+1;j<=n;j++){
            for(int k=i;k<=n;k++){
                if(fabs(a[j][k])>m){
                    m=fabs(a[j][k]);
                    p=j;
                    q=k;
                }
            }
        }
        if(p!=i){
            swap_r(p,i);
        }
        if(q!=i){
            swap_c(q,i);
        }
        for(int j=i+1;j<=n;j++){
            if(a[i][j]==0.0){
                continue;
            }
            double t=a[j][i]/a[i][i];
            a[j][i]=0.0;
            for(int k=i+1;k<=n+1;k++){
                a[j][k]-=t*a[i][k];
            }
        }
    }
}

void red()
{
    for(int i=n;i>1;i--){
        for(int j=i-1;j>=1;j--){
            if(a[j][i]==0){
                continue;
            }
            double t=a[j][i]/a[i][i];
            a[j][i]=0.0;
            for(int k=i+1;k<=n+1;k++){
                a[j][k]-=t*a[i][k];
            }
        }
    }
}

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n+1;j++){
            scanf("%lf",&a[i][j]);
        }
    }
    for(int i=1;i<=n;i++){
        a[n+1][i]=i;
    }
    gs();
    red();
    for(int i=1;i<=n;i++){
        ans[(int)a[n+1][i]]=a[i][n+1]/a[i][i];
    }
    for(int i=1;i<=n;i++){
        printf("x%d = %f\n",i,ans[i]);
    }
    return 0;
}

 

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值