O ( n 5 ) O(n^5) O(n5) 做法:
- 先求出
A
A
A 的伴随矩阵
A
∗
A^{*}
A∗,后利用
A
∗
A
∗
=
∣
A
∣
∗
E
⇒
A
−
1
=
A
∗
∣
A
∣
A*A^{*}=|A|*E\Rightarrow A^{-1}=\frac{A^{*}}{|A|}
A∗A∗=∣A∣∗E⇒A−1=∣A∣A∗ 求解
需要求 O ( n 2 ) O(n^2) O(n2) 次行列式
O ( n 4 ) O(n^4) O(n4) 做法:
- 对每一行来一波高斯消元
O
(
n
3
)
O(n^3)
O(n3) 做法:
首先介绍矩阵的初等变换(以下为初等行变换):
- 交换两行,记做 r i ↔ r j r_i\leftrightarrow r_j ri↔rj
- 将一行的所有元乘上数 k ≠ 0 k\neq 0 k=0
- 将一行的所有元的 k ≠ 0 k\neq 0 k=0 倍加到另一行上
- 若 A A A 能通过若干次初等行变换变为 B B B,则称 A A A 与 B B B 行等价, A ∼ B A\sim B A∼B(等价符号上有一个字母 r r r,没打出来,一下均用 ∼ \sim ∼ 表示行等价)
- 定理:设 A , B A,B A,B 为 m × n m\times n m×n 矩阵, A ∼ B A\sim B A∼B 的充分必要条件是存在 m m m 阶可逆矩阵 P P P 使得 P A = B PA=B PA=B
引入:初等矩阵: 把单位矩阵
E
E
E 经过一次初等变换得到的矩阵称为初等矩阵。
性质1:设
A
A
A 为
m
×
n
m\times n
m×n矩阵,对
A
A
A 进行一次初等行变换相当于在
A
A
A 的左边乘上对应的
m
m
m 阶初等矩阵
性质2:方阵
A
A
A 可逆的充分必要条件是存在有限个初等矩阵
P
1
,
P
2
,
.
.
.
,
P
l
P_1,P_2,...,P_l
P1,P2,...,Pl 使
A
=
P
1
P
2
.
.
.
P
l
A=P_1P_2...P_l
A=P1P2...Pl
充分性:因为初等矩阵均可逆,所以有限个初等矩阵的乘积仍可逆,故
A
A
A 可逆。
必要性:设
A
A
A 可逆,它经过有限次初等行变换变为行最简矩阵
B
B
B,那么存在初等矩阵
Q
1
,
Q
2
,
.
.
.
,
Q
l
Q_1,Q_2,...,Q_l
Q1,Q2,...,Ql 使
Q
1
Q
2
.
.
.
Q
l
A
=
B
Q_1Q_2...Q_lA=B
Q1Q2...QlA=B,因为
A
,
Q
i
A,Q_i
A,Qi 均可逆,故
B
B
B 可逆,故
∣
B
∣
≠
0
,
B
=
E
|B|\neq 0,B=E
∣B∣=0,B=E
此时
A
=
Q
1
−
1
Q
2
−
1
.
.
.
Q
l
−
1
B
=
Q
1
−
1
Q
2
−
1
.
.
.
Q
l
−
1
=
P
1
.
.
.
P
l
A=Q_1^{-1}Q_2^{-1}...Q_l^{-1}B=Q_1^{-1}Q_2^{-1}...Q_l^{-1}=P_1...P_l
A=Q1−1Q2−1...Ql−1B=Q1−1Q2−1...Ql−1=P1...Pl,此时
P
i
=
Q
i
−
1
P_i=Q_i^{-1}
Pi=Qi−1为初等矩阵
-
定理证明: A ∼ B ⇔ A\sim B \Leftrightarrow A∼B⇔ 存在有限个矩阵 P 1 , P 2 , . . . , P l P_1,P_2,...,P_l P1,P2,...,Pl 使 P l . . . P 2 P 1 A = B ⇔ P_l...P_2P_1A=B \Leftrightarrow Pl...P2P1A=B⇔ 存在可逆矩阵 P P P使 P A = B PA=B PA=B
-
推论:方阵 A A A 可逆的充分必要条件是 A ∼ E A\sim E A∼E
证明: A A A 可逆 ⇔ \Leftrightarrow ⇔ 存在可逆矩阵 P P P,使 P A = E ⇔ A ∼ E PA=E\Leftrightarrow A\sim E PA=E⇔A∼E -
本题解法:现在知道 A A A,求 P P P 使 P A = E PA=E PA=E
同时我们有 P E = E PE=E PE=E,那么我们有 P ( A , E ) = ( E , P ) P(A,E)=(E,P) P(A,E)=(E,P),于是我们对矩阵 ( A , E ) (A,E) (A,E) 做初等行变换,当 A A A 变为 E E E 时, E E E 就变成了 P P P,无解情况用推论判断即可,复杂度 O ( n 3 ) O(n^3) O(n3),非常好写
#include<bits/stdc++.h>
#define cs const
using namespace std;
int read(){
int cnt = 0, f = 1; char ch = 0;
while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
return cnt * f;
}
cs int Mod = 1e9 + 7;
typedef long long ll;
int add(int a, int b){ return a + b >= Mod ? a + b - Mod : a + b; }
int dec(int a, int b){ return a - b < 0 ? a - b + Mod : a - b; }
int mul(int a, int b){ ll r=(ll)a*b; return r>=Mod?r%Mod:r; }
int ksm(int a, int b){ int as=1; for(;b;b>>=1,a=mul(a,a)) if(b&1) as=mul(as,a); return as; }
void Add(int &a, int b){ a = add(a, b); }
void Mul(int &a, int b){ a = mul(a, b); }
void Dec(int &a, int b){ a = dec(a, b); }
cs int N = 405, M = N + N;
#define poly vector<int>
int n; poly a[N];
poly operator - (poly a, poly b){
for(int i=1; i<=n+n; i++) Dec(a[i],b[i]); return a;
}
poly operator * (poly a, int coe){
for(int i=1; i<=n+n; i++) Mul(a[i],coe); return a;
}
int main(){
n = read();
for(int i = 1; i <= n; i++){
a[i].resize(n+n+1);
for(int j = 1; j <= n; j++) a[i][j] = read();
}
for(int i = 1; i <= n; i++) a[i][i+n] = 1;
for(int i = 1; i <= n; i++){
int k = i;
for(int j = i+1; j <= n; j++) if(a[j][i]){ k = j; break; }
if(i ^ k) swap(a[i], a[k]);
int iv = ksm(a[i][i], Mod-2);
for(int j = i+1; j <= n; j++){
int coe = mul(iv, a[j][i]);
a[j] = a[j] - a[i] * coe;
}
}
for(int i = n; i >= 1; i--){
if(a[i][i] == 0){ puts("No Solution"); return 0; }
for(int j = i+1; j <= n; j++) if(a[i][j]) a[i] = a[i] - a[j] * a[i][j];
int iv = ksm(a[i][i],Mod-2); a[i] = a[i] * iv;
}
for(int i = 1; i <= n; i++){
for(int j = n+1; j <= n+n; j++) cout << a[i][j] << " ";
puts("");
} return 0;
}