利用了矩阵结合律，先算出构造递推矩阵自乘的结果，再与初始矩阵相乘。 #include using namespace std; typedef long long ll; #define ls (o<<1) #define rs (o<<1|1) #define pb push_back //#define a(i,j) a[(i)*(m+2)+(j)] //m是矩阵的列数 const double PI= acos(-1.0); const int M = 1e5+