设 A A A为 n ∗ s n*s n∗s的矩阵 B B B为 s ∗ m s*m s∗m的矩阵,则 A ∗ B A*B A∗B的结果 C C C为 n ∗ m n*m n∗m的矩阵。
给出矩阵乘法的公式: C i j = ∑ k = 1 s A i k B k j C_{ij}=\sum_{k=1}^sA_{ik}B_{kj} Cij=∑k=1sAikBkj ,其中 i i i表示行, j j j表示列。
简单来说, C C C矩阵中的第 i i i行第 j j j列由 A A A矩阵的第 i i i行和 B B B矩阵的第 j j j列对应相乘得到。
因此,对于许多递推问题,我们可以构建转移矩阵来将递推转移到下一步的状态。
例题(洛谷P1939):
题目描述
a[1]=a[2]=a[3]=1
a[x]=a[x-3]+a[x-1] (x>3)
求a数列的第n项对1000000007(10^9+7)取余的值
输入格式
第一行一个整数T,表示询问个数。
以下T行,每行一个正整数n
输出格式
每行输出一个非负整数表示答案
说明
对于100%的数据 T<=100,n<=2*10^9
在这道题中,由于n的范围非常大,暴力递推无论是时间还是空间都不过关,这时我们就可以选择使用矩阵加速递推来完成题目
首先,题目中的递推公式为 a [ 1 ] = a [ 2 ] = a [ 3 ] = 1 , a [ x ] = a [ x − 3 ] + a [ x − 1 ] ( x > 3 ) a[1]=a[2]=a[3]=1,a[x]=a[x-3]+a[x-1] (x>3) a[1]=a[2]=a[3]=1,a[x]=a[x−3]+a[x−1](x>3)
则可见,当 n < = 3 n<=3