素数筛选以及快速幂和慢速乘的实现

最大公约数

int gcd(int a,int b)//辗转相除法求最大公约数
{
    while(b)
    {
        int c=a%b;
        a=b;
        b=c;
    }
    return a;
}

最小公倍数
由于a*b有可能会溢出，因此可以先除再乘，因为GCD是它们的最大公约数，因此a/GCD一定是个整数

int lcm(int a,int b,int GCD)//最小公倍数=两个数的乘积/最大公约数
{
    return a/GCD*b;
}

快速幂和慢速乘

long long int fastPow(long long int a,long long int b,long long int p)//求(a^b)%p
{
    long long int ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ans;
}

long long int slowMul(long long int a,long long int b,long long int p) //a*b%p
{
    long long int ans=0;
    while(b)
    {
        if (b & 1) ans = (ans + a) % p;
        a = (a + a) % p;
        b >>= 1;
    }
    return ans;
}

素数的判定，时间复杂度是O(sqrt(n))

bool isPrime(int n)//判断是否是素数
{
    if(n<=1) return false;
    for(int i=2;i<=sqrt(n);i++)
    {
        if(n%i==0) return false;
    }
    return true;
}

埃氏筛选找素数 时间复杂度O(nloglogn)

const int maxn=1001;
bool p[maxn]{0};//如果i为素数，那么p[i]为false,否则为true
int prime[maxn];//存放素数
int pnum=0;//素数的个数
void findPrime()//获取2-1000的素数表
{
    for(int i=2;i<maxn;i++)
    {
        if(!p[i])//如果p[i]是素数，那么它的倍数都不是素数
        {
            prime[pnum++]=i;
            for(int j=i*i;j<maxn;j=j+i)
            {
                p[j]=true;
            }
        }
    }
}

欧拉筛法找素数 时间复杂度O(n)

const int maxn=1001;
bool p[maxn]{0};//如果i为素数，那么p[i]为false,否则为true
int prime[maxn];//存放素数
int pnum=0;//素数的个数
void findPrime()
{
    for(int i=2;i<maxn;i++)
    {
        if(!p[i]) prime[pnum++]=i;
        for(int j=0;j<pnum;j++)
        {
            if(i*prime[j]>maxn) break;
            p[i*prime[j]]=true;
            if(i%prime[j]==0) break;
        }
    }
}