struct inversematrix
void showb()
}ll qpow(ll x, ll n)
return res;
}int inverse()
for(int cur = 1, pit; cur <= n; ++cur)
if(pit > n)
return -1;
if(pit != cur)
int tmp = qpow(a[cur][cur], mod - 2);
for(int j = cur; j <= n; ++j)
a[cur][j] = 1ll * a[cur][j] * tmp % mod;
for(int j = 1; j <= n; ++j)
b[cur][j] = 1ll * b[cur][j] * tmp % mod;
for(int i = 1; i <= n; ++i)
}return n;
}} ma;
類似的方法可以找出使得a矩陣變換到b矩陣要左乘的矩陣。
分開a和b來寫可以減少矩陣第二維的長度,(在交換兩行時)對記憶體更友好。
矩陣的逆矩陣,伴隨矩陣
include include include include include include include include define n 100 using namespace std templateout type convert const in value t struct numb...
矩陣的求逆
最近做乙個加密演算法遇到需要計算矩陣的逆,閒著無聊,記錄一下,以後免得再麻煩。include include include define max 20 define e 0.000000001 計算矩陣src的模 double calculate a double src max int n fo...
逆矩陣 偽逆矩陣 資料的壓縮和復原
逆矩陣 偽逆矩陣 資料的壓縮和復原 這一塊知識雖然很簡單,但在光學各種實驗情況下經常用到,特此總結。矩陣的乘法 觀測矩陣 c a b c m 1 a m n b n 1 a可以理解為乙個轉化矩陣,或者說 觀測矩陣。矩陣b在觀測矩陣a上,觀察到的效果是矩陣c 這個思想在物理思想上非常重要,尤其是光學。...