普通快速幂:
ll quick_mod(ll a,int n)
{
ll sum=1;//注意是乘法,所以初始化为1
while(n)
{
if(n&1)//&运算符表示取这个数二进制的最后一位
{
sum=sum*a%mod;//若这位为1那么总和就应该乘当前的a
}
a=a*a%mod;//不论该位为0还是1,a都应该平方一波,具体自己理解
n>>=1;//把n的二进制右移一位,等价于n/=2;
}
return sum;
}
快速乘法:
ll cheng(ll a,ll n)
{
ll res=0;
a=a%mod;
n=n%mod;//注意先模一波,防止大数运算爆掉
while(n)//与快速幂基本一样,唯一的区别就是把*改成+
{
if(n&1)
{
res=(res+a)%mod;
}
a=(a+a)%mod;
n>>=1;
}
return res;
}
下面是矩阵的东西:
首先定一个矩阵结构体,我习惯这样因为写函数返回类型的时候直接写juzhen就可以
struct juzhen
{
ll data[maxn][maxn];
};
矩阵乘法:(优化版本)
juzhen mul(juzhen a,juzhen b)
{
juzhen c;
memset(c.data,0,sizeof(c.data));
for(int k=0;k<N;k++)
{
for(int i=0;i<N;i++)
{
if(!a.data[i][k]) continue;//优化在这个地方,如果为0说明它对这个位置贡献不了值,那么就直接跳过就可以
for(int j=0;j<N;j++)
{
if(!b.data[k][j]) continue;//同上
c.data[i][j]=(c.data[i][j]+a.data[i][k]*b.data[k][j]%mod)%mod;//这个地方自己推一下还是很好理解的,就是一点一点往上添东西的过程,而不是一次性填满,目的上面说了
}
}
}
return c;
}
矩阵快速幂:
juzhen quick_pow(juzhen a,ll b)
{
memset(res.data,0,sizeof(res.data));
for(int i=0,j=0;i<N&&j<N;i++,j++)
{
res.data[i][j]=1;//首先构造一个单位矩阵
}
while(b)
{
if(b&1)
{
res=mul(res,a);//与普通快速幂不同,矩阵快速幂就是把数的相乘变成矩阵与矩阵之间的乘法
}
a=mul(a,a);
b>>=1;
}
return res;
}