设A是数域上的一个n阶方阵，若在相同数域上存在另一个n阶矩阵B，使得：AB=BA=E（单位矩阵）。 则我们称B是A的逆矩阵，而A则被称为可逆矩阵。

现在告诉你一个n阶方阵A，求它的逆矩阵B。

输入数据有多组，第一行为数据组数T，接下来有T个矩阵。

每个矩阵的第一行为n(n<20)，表示n阶方阵，接下来的n行n列表示n*n的矩阵

矩阵元素为实数。

按照n行n列输出逆矩阵B，所有元素保留2位小数。数据保证一定有逆矩阵。

#include<stdio.h>

const int N=;

double a[N][N],b[N][N],c[N][N],bt[N][N],ct[N][N],ans[N][N];

int T,n;

void LU()

{

     for(int i=; i<n; i++)

        for(int j=; j<n; j++)b[i][j]=c[i][j]=bt[i][j]=ct[i][j]=ans[i][j]=;

    for(int i=; i<n; i++)b[i][i]=bt[i][i]=;

    double s;

    for(int i=; i<n; i++)

    {

        for(int j=i; j<n; j++)

        {

            s=;

            for(int k=; k<i; k++)

                s+=b[i][k]*c[k][j];

            c[i][j]=a[i][j]-s;

        }

        for(int j=i+; j<n; j++)

        {

            s=;

            for(int k=; k<i; k++)

                s+=b[j][k]*c[k][i];

            b[j][i]=(a[j][i]-s)/c[i][i];

        }

    }

    for(int i=; i<n; i++)

        for(int j=; j<i; j++)

        {

            s=;

            for(int k=; k<i; k++)

                s+=b[i][k]*bt[k][j];

            bt[i][j]=-s;

        }

    for(int i=; i<n; i++)

        ct[i][i]=/c[i][i];

    for(int i=; i<n; i++)

        for(int j=i-; j>=; j--)

        {

            s=;

            for(int k=j+; k<=i; k++)

                s+=c[j][k]*ct[k][i];

            ct[j][i]=-s/c[j][j];

        }

    for(int i=; i<n; i++)

        for(int j=; j<n; j++)

            for(int k=; k<n; k++)

                ans[i][j]+=ct[i][k]*bt[k][j];

    for(int i=; i<n; i++)

        for(int j=; j<n; j++)printf("%.2f%c",ans[i][j],j==n-?'\n':' ');

}

int main()

{

    scanf("%d",&T);

    while(T--)

    {

        scanf("%d",&n);

        for(int i=; i<n; i++)

            for(int j=; j<n; j++)

                scanf("%lf",&a[i][j]);

        LU();

    }

    return ;

}

#include<stdio.h>

#include<math.h>

#include<algorithm>

const double eps=1e-;

const int N=;

double a[N][N],b[N],x[N],s,t[N][N],ans[N][N];

int n;

void gauss()

{

    int i;

    for(int k=; k<=n; k++)

    {

        for(i=k; i<=n&&fabs(a[i][k])<eps; i++);

        if(i!=k)

        {

            for(int j=k; j<=n; j++)std::swap(a[i][j],a[k][j]);

            std::swap(b[i],b[k]);

        }

        for(i=k+; i<=n; i++)

        {

            s=a[i][k]/a[k][k];

            for(int j=k; j<=n; j++)a[i][j]-=a[k][j]*s;

            b[i]-=b[k]*s;

        }

    }

    for(i=n; i>=; --i)

    {

        s=b[i];

        for(int j=i+; j<=n; j++)s-=x[j]*a[i][j];

        x[i]=s/a[i][i];

        if(fabs(x[i])<eps)x[i]=;

    }

}

int main()

{

    int T;

    scanf("%d",&T);

    while(T--)

    {

        scanf("%d",&n);

        for(int i=; i<=n; i++)

        {

            for(int j=; j<=n; j++)

                scanf("%lf",&t[i][j]),a[i][j]=t[i][j];

            b[i]=;

        }

        for(int k=; k<=n; k++)

        {

            for(int i=; i<=n; i++)

            {

                for(int j=; j<=n; j++)a[i][j]=t[i][j];

                b[i]=;

            }

            b[k]=;

            gauss();

            for(int i=; i<=n; i++)ans[i][k]=x[i];

        }

        for(int i=; i<=n; i++)

            for(int j=; j<=n; j++)printf("%.2f%c",ans[i][j],j==n?'\n':' ');

    }

    return ;

}