矩阵求逆模板（普通矩阵和01矩阵）

P4783 【模板】矩阵求逆

题意：

给定n*n的矩阵，求它的逆矩阵，对1e9+7取模。

数据范围：n<=400，a(i,j)<=1e9

解法：

线性代数学过.
流程基本就是:
将给定矩阵A利用初等变换化为单位阵,
变化过程中对另外一个单位矩阵B进行相同操作,
当矩阵A变成单位矩阵的时候,B就是A的逆矩阵了.

当矩阵A无法化为单位矩阵的时候当然就不存在逆矩阵了.

code：

#include<bits/stdc++.h>
using namespace std;
const int N=505;
const int mod=1e9+7;
int a[N][N<<1];
int n;
int ppow(int a,int b,int mod){
    int ans=1%mod;a%=mod;
    while(b){
        if(b&1)ans=1ll*ans*a%mod;
        a=1ll*a*a%mod;
        b>>=1;
    }
    return ans;
}
int Gauss_rev(int n){
    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j++){//找第i列非零的行换上来
            if(a[j][i]){
                swap(a[i],a[j]);
                break;
            }
        }
        if(!a[i][i])return 0;//无解
        int kk=ppow(a[i][i],mod-2,mod);//逆元
        for(int j=i;j<=n*2;j++){//当前行每一列都除以a[i][i]
            a[i][j]=1ll*a[i][j]*kk%mod;
        }
        for(int j=1;j<=n;j++){//其他行
            if(j!=i){
                kk=a[j][i];
                for(int k=i;k<=n*2;k++){
                    a[j][k]=(a[j][k]-1ll*kk*a[i][k]%mod+mod)%mod;
                }
            }
        }
    }
    return 1;
}
signed main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&a[i][j]);
        }
        a[i][i+n]=1;
    }
    if(!Gauss_rev(n)){
        puts("No Solution");
    }else{
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                printf("%d ",a[i][j+n]);
            }
            puts("");
        }
    }
    return 0;
}

---

01矩阵求逆（对bitset直接高斯消元就行了）：

#include<bits/stdc++.h>
using namespace std;
const int N=505;
bitset<N<<1>a[N];
int n;
int Gauss_rev(int n){
    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j++){//找第i列非零的行换上来
            if(a[j][i]){
                swap(a[i],a[j]);
                break;
            }
        }
        if(!a[i][i])return 0;//无解
        for(int j=1;j<=n;j++){//其他行
            if(a[j][i]&&j!=i){
                a[j]^=a[i];
            }
        }
    }
    return 1;
}
signed main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&a[i][j]);
        }
        a[i][i+n]=1;
    }
    if(!Gauss_rev(n)){
        puts("No Solution");
    }else{
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                printf("%d ",a[i][j+n]);
            }
            puts("");
        }
    }
    return 0;
}

---