关于素数问题相关算法模板

关于素数问题相关算法模板

原始筛法

• 问题背景：判断一个数是不是素数
• 原理：枚举所有可能的因数，有因数则不是素数
• 这里为什么用sqrt(n)，是因为没有必要枚举所有的数字，每一个小于根号n的因数i，一定有一个大于根号n的因数j与之相对应，也就是i*j=n,所以枚举小于等于根号n的因数即可。

void isPrimer(int n){
	for(int i=2;i<=sqrt(n);i++){
		if(n%i==0) return false;
	}
	return true;
}

埃氏筛法（普通筛法/埃拉托斯特尼(Eratosthenes)筛法）

• 问题背景：枚举n以内的所有素数
• 原理：先把所有整数列出来，从2开始遍历，判断每一个数，如果是素数，则将其所有范围内的倍数划掉；对于每一个数字i，如果它没有被划去，她一定是素数，因为他不是任何2到i-1数字的倍数，然后接着划去它的倍数就OK。

const int maxn = 100000000;
int prime[maxn+5];
int a[maxn+5];
int sumPrime(int n){
	int p = 0;//记录素数的个数
	for(int i=0; i<=n; i++) a[i] = 1;//将所有的数字状态置为1
	a[0] = 0; a[1] = 0;
	for(int i=2; i<=n; i++){
		if(a[i]){
			prime[p++] = i;//记录素数和个数
			for(int j = 2*i; j<=n; j+=i) a[j] = 0;//将其倍数划去 
		}
	}
}

线性筛（欧拉Euler筛）

• 问题背景：在O(n)的时间复杂度内求出1~n之间所有的质数
• 原理：prime[]数组中的素数是递增的,当i%prime[j]==0，那么iprime[j+1]这个合数肯定被prime[j]乘以某个数筛掉。原因：因为i中含有因子prime[j],prime[j]比prime[j+1]小，即i=kprime[j]，那么iprime[j+1]=（kprime[j]）prime[j+1]=k（prime[j]*prime[j+1]），接下去的素数同理。所以不用筛下去了。因此，在满足i%prime[j]==0这个条件之前以及第一次满足改条件时,prime[j]必定是prime[j]*i的最小因子。
• 通俗的理解：筛法很多数被处理了不止1遍，比如6，在素数为2的时候处理1次，为3时候又处理一次，因此又造成了不必要处理。此代码模板保证每个合数只会被它的最小质因子数筛去，因此每个数只会被标记一次，时间复杂度为O(n)

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) break;//保证每一个合数只会被它最小的质因子筛去，因此一个数只会被标记一次
        }
    }
}

具体算法定理可看：欧拉筛法求素数