矩阵乘法是数学线性代数里面的东西,可以用在方程式的推导,使递推变成一个矩阵乘法,就比如:f(n)=f(n-1)-2*f(n-2)+1;让你求f(n)的值,如果一个一个去递推有时候会有点费时间,但是可以转换成矩阵的乘法,然后再加上快速幂就行了。
矩阵乘法的计算
在计算两个矩阵的时候可以在纸上画,矩阵相乘的结果矩阵中,就比如下图左边那个两个红框重合位置的结果就相当于这一行第一个乘以这一列第一个加上这一行第二个乘以这一列第二个加上这一行第三个乘以这一列第三个,即1*1+2*2+3*3=14.其它位置也是同样规则。
struct matrix// 矩阵
{
ll x[3][3];
};
matrix multiply(matrix &a, matrix &b)//矩阵乘法的代码
{
matrix c;
memset(c.x,0,sizeof(c.x));
for(int i = 0; i <3; i++)
{
for(int j = 0; j <3; j++)
{
for(int k = 0; k <3; k++)
{
c.x[i][j] += a.x[i][k]*b.x[k][j];
}
c.x[i][j] %= m;//如果题上需要取模的话
}
}
return c;
}
快速幂
就不说了直接套模板就行,矩阵乘法经常和快速幂一块儿用。
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
矩阵推导
那么问题是你怎么将一个方程的问题对应的矩阵推导出来,就举一个简单的例子。
例子:f(n)=f(n-1)-2*f(n-2)+1。
由这个例子可以看出这相当于是一个递推,所以,由等号右边可知由右边三个元素可推出左边,所以可制造一个矩阵包含这三个元素,然后根据方程式推导出下一个对应的矩阵。如图。
矩阵题只需要自己推出来矩阵就行了,剩下的只需要再写一个矩阵快速幂模板就行。
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
样例代码:
struct matrix//本题所需的只是3阶矩阵 可以根据题意自己设置
{
ll x[3][3];
};
matrix mpow(matrix a, ll m)
{
matrix b;
memset(b.x,0,sizeof(b.x));
for(int i=0; i<3; i++)//首先需要有一个初始化矩阵存放结果,只需要把i==j的位置放1就行,那么任何矩阵乘以该矩阵都不会变。
b.x[i][i]=1;
while(m) //快速幂模板
{
if(m & 1)b = multiply(b, a);
a = multiply(a, a);
m >>= 1;
}
return b; //返回结果矩阵
}
matrix a= {1,2,1,1,0,0,0,0,1};//首先是自己推导矩阵
int main()
{
mod=m;
if(n <= 2)//意思是问的就只是前两个值 直接输出就行
{
printf("%lld\n", n % m);
return 0;
}
matrix s=a;
s = mpow(s, n-2);//矩阵快速幂 使s矩阵乘以n-2次,然后返回结果矩阵s
ll ans=((2ll*s.x[0][0]%mod+s.x[0][1])%mod+s.x[0][2])%mod;//这个位置就是根据矩阵乘法规则自己去算的结果,当然也可以直接在矩阵快速幂的时候就把初始矩阵加进去也行,但是要注意乘法规则。
printf("%lld\n", ans);
return 0;
}
以上就是关于一些矩阵快速幂的浅解。
总结
总而言之,这一类题最重要的就是方程式对应矩阵的推导,只要矩阵推导出来就行,然后就是关于矩阵乘法的计算,注意顺序,否则容易出错。