【数论】【模板】快速幂、快速乘法以及矩阵乘法、矩阵快速幂

普通快速幂：

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;
}