数论

参考刘汝佳《算法竞赛入门经典-训练指南》P119~P122

typedef long long LL; 
//最大公约数 
LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);   
} 
//最小公倍数 
LL lcm(LL a,LL b)
{
    return a/gcd(a,b)*b;
}
//素数筛 
void sieve(int n)
{
    int m=(int)sqrt(n+0.5);
    memset(vis,0,sizeof(vis));
    for(int i=2;i<=m;i++){
        if(!vis[i]){
        for(int j=i*i;j<=n;j+=i){
            vis[j]=1;
        }
    }
}
//生成素数表 ，放在prime中，返回素数个数
int gen_primes(int n)
{
    sieve(n);
    int c=0;
    for(int i=2;i<=n;i++){
        if(!vis[i])
            prime[c++]=i;}
    }   
    return c;
} 
//欧几里得算法 
void gcd(LL a,LL b,LL& d,LL& x,LL& y)
{
    if(!b)
        {d=a,x=1,y=0;}
    else{
        gcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}
//求逆元 
LL inv(LL a,Ll n)
{
    LL d,x,y;
    gcd(a,n,d,x,y);
    return d==1;?(x+n)%n:-1;
}