快速幂+矩阵快速幂模板

快速幂

long long q_pow(long long a,long long k,long long b)//求a^b%b
{
    long long f=1,s=a;
    while(k)
    {
        if(k&1)
            f=f*s%b;
        s=s*s%b;
        k=k/2;
    }
    return f;
}

矩阵快速幂

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
using namespace std;
const int N=100;
int c[N][N];
node cmp(int a[][N],b[][N],int n)
{
    memset(c,0,sizeof(c));
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    for(int k=0;k<n;k++)
    c[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]=c[i][j];
}
node q_pow(int a[][N],int n)//计算a^b,n为矩阵的大小
{
    memset(res,0,sizeof(res));
    for(int i=0;i<n;i++)
        res[i][i]=1;
    while(n)
    {
        if(n&1)
            cmp(res,a,N);
        cmp(a,a,N);
        n>>=1;
    }
}