矩阵快速幂的模板

struct matrix
{
    ll t[64][64];
}re,x;
int lim,n;
matrix operator * (matrix a,matrix b)
{
    int i,j,k;
    matrix c;
    for(i=0;i<=lim;i++)
        for(j=0;j<=lim;j++)
    {
        c.t[i][j]=0;
        for(k=0;k<=lim;k++)
        {
            c.t[i][j]+=a.t[i][k]*b.t[k][j]%mod;
            c.t[i][j]=c.t[i][j]%mod;
        }
    }
    return c;
}
void ksm()
{
    int bj=0;
    while(n)
    {
        if(n&1)
        {
            if(!bj)
            {
                re=x;
                bj=1;
            }
            else re=re*x;
        }
        x=x*x;
        n>>=1;
    }
}