逆矩阵模板 #线性代数

最新推荐文章于 2024-08-27 18:55:41 发布

Samaritan_infi

最新推荐文章于 2024-08-27 18:55:41 发布

阅读量349

点赞数

分类专栏：线性代数

本文链接：https://blog.csdn.net/Samaritan_infi/article/details/77860188

版权

线性代数专栏收录该内容

1 篇文章 0 订阅

订阅专栏

#include<cstdio>
#include<algorithm>
#include<cmath>
#include <iostream>
using namespace std;
double a[105][205],b[105],ans;
int main()
{
    int n,q,m;
    while(scanf("%d",&n)!=EOF)
    {
        m=2*n;
        for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
       {
           scanf("%lf",&a[i][j]);
           if(i==j)
           a[i][n+j]=1;
           else
           a[i][n+j]=0;
       }
       for(int i=1;i<=n;i++)
       {
        int pos=i;
        for(;pos<=n;pos++)
            if(fabs(a[pos][i])>1e-6)break;
            if(pos!=i)
            for(int k=i;k<=m;k++)
                swap(a[i][k],a[pos][k]);
        double temp=a[i][i];
        for(int j=i;j<=m;j++)a[i][j]/=temp;
        for(int j=1;j<=n;j++)
        {
            if(j==i)continue;
            temp=a[j][i];
            for(int k=i;k<=m;k++)
            a[j][k]-=a[i][k]*temp;
        }
        }
        for(int i = 1 ; i <= n ; i ++ ){
            for(int j = 1 ; j <= 2*n ; j ++ )
                cout << a[i][j]<<" " ;
            cout << endl ;
        }

    }
}

can you hear me ?