二階矩陣整數冪(II)
程式由網友roviury提供。
程式可以計算一個2×2矩陣的冪。 這個程式的速度比第(I)版快,特別是計算較大正整數冪時,速度明顯快很多。
網更新日期: 2011年11月8日
程式 (217 bytes)
?→A: ?→B: ?→C: ?→D: ?→M: C=0 => Goto 3:
A+D→X: . 5( X - √( X2 - 4(AD - BC→Y: X - Y→X:
(X - D)÷C→A: (Y - D)÷C→B: 1→C: 1→D: Goto 1:
Lbl 3: BA≠BD=> Goto 4: A^M◢ BMA^(M-1◢ 0◢ D^M: Goto 2:
Lbl 4: A→X: D→Y: (D-A)÷B→D: 1→A: 1→B:
Lbl 1: X^M→X: Y^M→Y: 1÷(AD - BC→M: ADXM - BCYM◢
MAB(Y - X◢ MCD(X - Y◢ ADYM - BCXM: Lbl 2
例題: 計算以下矩陣及其逆矩陣:
按 Prog 1 再按 1 EXE 2 EXE 3 EXE 4 EXE 3 EXE (顯示37)
EXE (顯示54) EXE (顯示81) EXE (顯示118) EXE (顯示 -59/4)
EXE (顯示 27/4) EXE (顯示 81/8) EXE (顯示 -37/8)
所以
例題2: 計算
按 Prog 1 再按 1 EXE 0 EXE 0 EXE 1EXE 2009 EXE (顯示1)
EXE (顯示0) EXE (顯示0) EXE (顯示1)
所以