数论之素数

1.试除法判定质数

质数（素数）如果只包含一和它本身两个约数

bool check(int x)
{
    if(x<2)//要大于1
        return false;
    for(int i=2;i<x/i;i++)//这个地方最好写x/i
        if(x%i==0)
            return false;
        return true;
}

2分解质因子

void divide(int x)
{
    for(int i=2;i<=x/i;++i)
        if(x%i==0)
        {
            int s=0;
            while(x%i==0)x/=i,s++;
            cout<<i<<' '<<s;
        }
    if(x>1) cout<<x<<' '<<1<<endl;
}

埃氏筛

#define N 10000
int primes[N],cnt;
bool st[N];
void get_prime(int n)
{
    for(int i=2;i<=n;++i)
    {
        if(!st[i])//如果没有被筛掉
        {
            primes[cnt++]=i;//这个就是素数
        }
        for(int j=i+i;j<=n;j+=i)
            st[j]=true;
    }

}

欧拉筛 线性筛

核心：n只会被最小质因子筛掉
i%primes[j]==0 primes[j]一定是i的最小质因子 primes[j]也是primes[j]*i的最小质因子 这里之后就一定要break 后面就不是最小的质因子了

#define N 10000
int primes[N],cnt;
bool st[N];
void get_primes(int n)
{
    for(int i=2;i<=n;++i)
    {
        if(!st[i])
            primes[cnt++]=i;
        for(int j=0;primes[j]<=n/i;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)//primes[j]一定是i的最小质因子
                break;//这是最神奇的一步神来之笔！
        }

    }
}