快速幂/数论技巧专题

快速幂取模

long long pow(long long a,int k,long long mod)//快速幂取模
{
    long long b=1;
    while(k)
    {
        if(k&1) b=b*a%mod;
        a=a*a%mod;
        k>>=1;
    }
    return b;
}

矩阵快速幂

struct Mat
{
    long long m[3][3];
    Mat(){memset(m,0,sizeof(m));}
    void init(int c,int d,int x)
    {
        memset(m,0,sizeof(m));
        m[0][0]=d;m[0][1]=c;m[0][2]=x;m[1][0]=m[2][2]=1;
    }
};
Mat multi(const Mat &a,const Mat &b)
{
    Mat c;
    for(int i=0;i<3;i++)
    for(int j=0;j<3;j++)if(a.m[j][i]!=0)
    for(int k=0;k<3;k++)if(b.m[i][k]!=0)
    c.m[j][k]=(c.m[j][k]+a.m[j][i]*b.m[i][k]%MOD)%MOD;
    return c;
}
Mat pow(Mat &a,int k)
{
    Mat b;
    for(int i=0;i<3;i++) b.m[i][i]=1;
    while(k)
    {
        if(k&1) b=multi(b,a);
        a=multi(a,a);
        k>>=1;
    }
    return b;
}

POJ3233 Matrix Power Series 矩阵快速幂

HDU6395 Sequence 矩阵快速幂+分块

GCD/LCM

long long gcd(long long x,long long y)
{
    return y==0?x:gcd(y,x%y);
}
long long lcm(long long x,long long y)
{
    return x/gcd(x,y)*y;
}