数论初步(更新)

1.GCD与LCD

int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}

白书上的神奇代码。。。

2.扩展欧几里得

void egcd(LL a,LL b,LL &d,LL &x,LL &y)
{
    if(!b)
    {
        d=a;x=1;y=0;
    }
    else
    {
        egcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}

3.幂取模

LL pow_mod(LL a,LL n,LL m)
{
    if(!n)return 1;
    LL x=pow_mod(a,n/2,m);
    LL ans=x*x%m;
    if(n%2==1)ans=ans*a%m;
    return ans;
}