【C++】素数打表

前言：
本文素数打表的全局开法如下

#include<bits/stdc++.h>
#define MAXN 100
using namespace std;
//用bool类型作标志位数组完全是为了方便处理，因为bool为一个字节，用memset可以直接处理为true
bool isprime[MAXN+5];//isprime[i]=true表示数字i为质数 
int nxt[MAXN+5];//nxt[i]表示数字i的最小质因数在prime数组的位置
int prime[MAXN+5];//prime[i]表示为第i个质数的值，相当于一个质数表 

---

素数打表

素数打表，预处理素数是必备技巧之一，其中有个问题要注意：因为要开一个标志位的数组来对位记录是否为质数，所以空间复杂度上面有限制的，100000以内有9500+个质数 ，1000000以内有78000+个质数

这里列举两个常用素数打表的方法：

埃拉托斯特尼筛法

埃拉托斯特尼筛法，简称埃氏筛，时间复杂度达到O(n²)
理论依据：如果要得到一个自然数n以内的全部素数，就必须把不大于根号n的所有素数的倍数剔除掉，那么剩下来的数就都是素数了。
过程：先用2来筛，把其倍数标记成非素数；再用下一个标记不是非素数的数来筛，把其倍数标记成非素数……以此类推

代码示例

void Eratosthenes_Sieve()
{
	int i,j;
	memset(isprime,true,sizeof(isprime));
	for(i=2;i<=n;i++)
	{ 
		if(!isprime[i])
			continue;
		prime[++prime[0]]=i;
        nxt[i]=prime[0];
		for(j=i*i;j<=n;j+=i)
		{
			isprime[j]=false;
			if(!nxt[j])
				nxt[j]=prime[0];
		} 
	} 
}

欧拉筛法

欧拉筛法，又称线性筛，时间复杂度只有O(n)
理论依据：在埃氏筛法的基础上，让每一个合数都只被它的最小质因子筛选一次，以达到不重复筛选的目的，大大地节省了时间，从埃氏筛的O(n²)降到O(n)级别
方法：

1. 把isprime全部置true
2. 判断isprime是否为true，是质数直接加入prime数组，prime数组存放的就是我们处理出来的质数，
   prime[0]则用来记录加入的质数数量，同时把prime存放的位置赋值给nxt数组，nxt也是我们需要打的一个辅助表
3. 循环prime数组，把当前的数和质数表内的质数的得到的乘积置为false，相当于从质数表中剔除，
   同时为了达到只筛一次的目的，如果当前的数能被质数表的质数整除，就break出去

代码示例

void Euler_Sieve()
{
	int i,j;
    memset(isprime,true,sizeof(isprime));
	prime[0]=0;
    for(i=2;i<=n;i++) 
	{
        if(isprime[i])
		{
            prime[++prime[0]]=i;
            nxt[i]=prime[0];
        }
        for(j=1;j<=prime[0]&&i*prime[j]<=n;j++)
		{
            isprime[i*prime[j]]=false;
            nxt[i*prime[j]]=j;
            if(i%prime[j]==0)
                break;
        }
    }
}