日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

洛谷 - P4783 【模板】矩阵求逆(高斯消元求逆矩阵)

發布時間:2024/4/11 编程问答 45 豆豆
生活随笔 收集整理的這篇文章主要介紹了 洛谷 - P4783 【模板】矩阵求逆(高斯消元求逆矩阵) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

題目鏈接:點擊查看

題目大意:給出一個 n * n 的矩陣,求出其逆矩陣,mod 為 1e9 + 7,若不存在輸出?No Solution

題目分析:囤個模板,原理就是,初始時在原矩陣右側設置一個單位矩陣,隨著原矩陣進行高斯消元的過程,單位矩陣也跟隨進行一系列的變換,當原矩陣變為單位矩陣時,單位矩陣也就變成了逆矩陣,符號語言就是:

代碼:

#include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cmath> #include<cstring> #include<algorithm> #include<stack> #include<climits> #include<queue> #include<map> #include<set> #include<sstream> #include<cassert> #include<bitset> #include<unordered_map> using namespace std;typedef long long LL;typedef unsigned long long ull;const int inf=0x3f3f3f3f;const int N=410;const int mod=1e9+7;int n;LL a[N][N<<1];LL q_pow(LL a,LL b) {LL ans=1;while(b){if(b&1)ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans; }bool Gauss_j() {for(int i=1,r;i<=n;i++){r=i;for(int j=i+1;j<=n;j++)if(a[j][i]>a[r][i])r=j;if(r!=i)swap(a[i],a[r]);if(!a[i][i])return false;int kk=q_pow(a[i][i],mod-2);for(int k=1;k<=n;k++){if(k==i)continue;int p=a[k][i]*kk%mod;for( int j=i;j<=(n<<1);j++)a[k][j]=((a[k][j]-p*a[i][j])%mod+mod)%mod;}for(int j=1;j<=(n<<1);++j) a[i][j]=(a[i][j]*kk%mod);}return true; }int main() { #ifndef ONLINE_JUDGE // freopen("data.in.txt","r",stdin); // freopen("data.out.txt","w",stdout); #endif // ios::sync_with_stdio(false);scanf("%d",&n);for(int i=1;i<=n;i++){a[i][i+n]=1;for(int j=1;j<=n;j++)scanf("%lld",&a[i][j]);}if(Gauss_j()){for(int i=1;i<=n;i++){for(int j=n+1;j<=n<<1;j++)printf("%lld ",a[i][j]);puts("");}}elseputs("No Solution");return 0; }

?

總結

以上是生活随笔為你收集整理的洛谷 - P4783 【模板】矩阵求逆(高斯消元求逆矩阵)的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。