基础数论算法

• 最大公约数　

int gcd(int a,int b)
{
    int temp;
    while(b)
    {
        temp = b;
        b = a%b;
        a = b;
    }

    return a;
}

• 最小公倍数

int lcm(int a,int b)
{
    return a*b/gcd(a,b);
}

• 素数

　　　　判断n是否为素数

bool prime(int n)
{
    for(int i=0; i<=sqrt(n); i++)
    {
        if(n%i==0) return false;
    }
    return true;
}

　　　　判断n以内的素数：素数筛法

void prime()
{
    int prime[n+1];
    memset(prime,0,sizeof(prime));

    prime[1]=1;
    for(int i=2; i<=n; i+=2)
        prime[i]=1;

    for(int i=3; i<=n; i+=2)
        for(int j=i; j<=n; j+=i)
            prime[j]=1;
}