矩阵快速幂模板

const int mod=1e9+7;

struct mat
{
    int m[15][15];
}a;                          //用结构体储存矩阵，方便记录

mat mul(mat x,mat y,int n)
{
    mat ret;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            ret.m[i][j]=0;
            for(int k=1;k<=n;k++)ret.m[i][j]=(ret.m[i][j]+x.m[i][k]*y.m[k][j]%mod)%mod;
        }
    }
    return ret;
}                                //矩阵乘法板子

mat km(mat a,int n,int p)
{
    mat ret;
    memset(ret.m,0,sizeof(ret.m));
    for(int i=1;i<=n;i++)ret.m[i][i]=1;
    while(p)
    {
        if(p&1)ret=mul(ret,a,n);
        a=mul(a,a,n);
        p>>=1;
    }
    return ret;
}                                        //快速幂的模板把乘法换成上面写的矩阵乘法就行