Maple 中的快速模幂运算

Question

计算整数矩阵对给定整数取模的幂的快速方法是什么？

我试过了：

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

但这很慢，很可能 Maple 尝试先计算幂（有大量数字），然后才取模 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