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;
}
}
矩阵快速幂的模板
猜你喜欢
转载自blog.csdn.net/becky_w/article/details/80059849
今日推荐
周排行