Maple中的快速模冪運算

Question

以一個給定的整數爲模來計算整數矩陣的冪的快速方法是什麼？

我想：

> M := Matrix([[1,1],[1,0]]); M^(10^12) mod 73; 

但是這是非常緩慢的，最有可能華普試圖先計算力（與龐大的數字），然後才取模73.我怎麼能說服它爲每個乘法做模數？

Answer 1

restart: 

M := Matrix([[1,1],[1,0]]): 

str:=time[real](): 

LinearAlgebra:-Modular:-MatrixPower(73, M, 10^12); 

           [46 46] 
           [  ] 
           [46 0] 

time[real]()-str; 

           0.040