快速幂

实际上是将幂次看为二进制下，

嘻嘻嘻

struct mat
{
    LL s[4][4];
};
void operator*=( mat &a,const mat &b)
        {
            static mat c;
            int i,j,k;
            memset(c.s,0,sizeof(c));
            for(k=0;k<3;k++)
                for(i=0;i<3;i++)
                    for(j=0;j<3;j++)
                        c.s[i][j]+=a.s[i][k]*b.s[k][j]%mod,c.s[i][j]%=mod;
            for(i=0;i<3;i++)
                for(j=0;j<3;j++)
                    a.s[i][j]=c.s[i][j];
        }
mat ans,t;

 for(LL i=k;i;i=i>>1,t*=t)
        if(i&1) ans*=t;