关于矩阵快速幂,首先来了解一下什么叫快速幂。
上代码:
#include <iostream>
#include <cstdio>
using namespace std;
long long f(int n,int k)
{
int ans=1;
while(k)
{
if(k&1)
ans*=n;
k>>=1;
n=n*n;
}
}
int main(){
int n,k;
scanf("%d %d",&n,&k);
printf("%lld",f(n,k-1));
return 0;
}
这个代码就是一个简单的快速幂的模板了,这个代码可以快速的计算比较大的幂次。原理就是把幂次改成二进制,然后从地位开始扫描,若扫描到二进制中的1,就乘这个位数的次数,例如:999次方改成二进制就是1111100111,然后如果第一位是1,就乘一次方,第二位是1,就乘二次方以此类推,就可以把要计算999次的算法改成只需要少数计算的算法了,大大减少了计算量。
而矩阵快速幂就是把快速幂中的乘法改成矩阵的乘法,然后就可以快速的计算出高次方矩阵的值了。
然后矩阵乘法的定义就是:
由上图可知道矩阵乘法的定义,所以我们可以很简单的写出矩阵乘法的函数,再用快速幂的办法写出矩阵快速幂的函数。
代码如下
#include <iostream>
#include <cstdio>
using namespace std;
const int N=10;
int tmp[N][N];
void multi(int a[][N],int b[][N],int n) // 矩阵乘法
{
memset(tmp,0,sizeof tmp);
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
for(int k=0;k<n;k++)
tmp[i][j]+=a[i][k]*b[k][j];
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
a[i][j]=tmp[i][j];
}
int res[N][N];
void Pow(int a[][N],int n) //快速幂函数
{
memset(res,0,sizeof res);//n是幂,N是矩阵大小
for(int i=0;i<N;i++) res[i][i]=1;
while(n)
{
if(n&1)
multi(res,a,N);//res=res*a;复制直接在multi里面实现了;
multi(a,a,N);//a=a*a
n>>=1;
}
}
这两个函数就可以解决矩阵快速幂的计算了。
但是矩阵快速幂有什么用呢?
举个例子:
斐波那契数列大家都很熟悉了,如果要计算上亿个斐波那契函数,普通的循环或者递归肯定会超时。那这个时候就能使用矩阵快速幂了
有定义可知道,f(n)=f(n-1)+f(n-2)
我们可以推出
根据矩阵乘法的定义可推出下一位。
所以要求多少位,就可以乘以多少次方的矩阵。
类似的递推问题也可以用这个来解决。