矩阵快速幂的初步理解

关于矩阵快速幂，首先来了解一下什么叫快速幂。
上代码：

#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)
我们可以推出

根据矩阵乘法的定义可推出下一位。
所以要求多少位，就可以乘以多少次方的矩阵。

类似的递推问题也可以用这个来解决。