筛法求素数

求素数是程序设计比赛中经常遇到的问题，最基本的方法是通过素数的定义直接判断，只能被1和它本身整除的数就是素数了。这种方法适合判断单个数是否为素数，当要求一个范围内素数而这个范围又比较大时，这种方法就不太使用了，甚至程序要运行几分钟才能算出结果。

筛法的思想是去除要求范围内所有的合数，剩下的就是素数了，而任何合数都可以表示为素数的乘积，因此如果已知一个数为素数，则它的倍数都为合数。
普通的线性筛法：

int su[N];//这里面存放的是素数
int book[N+5];//su:素数集合,book[]标记是否为素数,注意一定要加x否则数组溢出特别难查出来
int cnt;//cnt代表的是在su[]这个数组中，因为要往里塞题给范围内的素数嘛，所以这个cnt就是塞素数时的下角标，当然了，还充当输出时的输出上界
void xianxingshai()
{
    cnt=1;
    memset(book,-1,sizeof(book));
    for(int i=2;i<maxx;i++)
    {
        if(book[i]==-1)
        {
            su[cnt++]=i;
        }
        for(int j=i*2;j<maxx;j+=i)
            book[j]=0;
    }
}

普通的线性筛法虽然大大缩短了求素数的时间，但是实际上还是做了许多重复运算，比如2*3=6，在素数2的时候筛选了一遍，在素数为3时又筛选了一遍。如果只筛选小于等于素数i的素数与i的乘积，既不会造成重复筛选，又不会遗漏。时间复杂度几乎是线性的。
优化后的线性筛法：

long long prime[N+5];//这个集合中存的是质数的集合
long long num_prime;//num_pirme记录素数个数
int book[N];//book[]为-1是素数，为0是合数
void su()
{
    memset(cnt,0,sizeof(cnt));
    cnt=0;
    cnt[0]=cnt[1]=1;
    for(int i=2; i<100010; i++)
    {
        if(!cnt[i])
            prime[cnt++]=i;
        for(int j=0; j<cnt&&prime[j]*i<100010; j++)
        {
            cnt[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
}

埃氏筛：

const int N = 100000 + 5;
bool prime[N];
void init()
{
    for(int i = 2; i < N; i ++)
        prime[i] = true;
    for(int i = 2; i*i < N; i ++) //判断改成i*i<N
    {
        if(prime[i])
        {
            for(int j = i*i; j < N; j += i) //从i*i开始就可以了
            {
                prime[j] = false;
            }
        }
    }
}
int main()
{
    init();
}

线性筛$O(n)$

const int N = 100000 + 5;
bool prime[N];//prime[i]表示i是不是质数
int p[N], tot;//p[N]用来存质数
void init()
{
    for(int i = 2; i < N; i ++)
        prime[i] = true;//初始化为质数
    for(int i = 2; i < N; i++)
    {
        if(prime[i])
            p[tot ++] = i;//把质数存起来
        for(int j = 0; j < tot && i * p[j] < N; j++)
        {
            prime[i * p[j]] = false;
            if(i % p[j] == 0)
                break;//保证每个合数被它最小的质因数筛去
        }
    }
}
int main()
{
    init();
}

单纯判断一个数是否为素数高效做法：

bool judge(ll num)
{
	if (num == 2 || num == 3 || num == 5)
		return true;
	if (num % 2 == 0 || num % 3 == 0 || num % 5 == 0 || num == 1)
		return false;
	ll c = 7;
	int maxc = (int)(sqrt(num));
	while (c <= maxc)
	{
		if (num % c == 0)
			return false;
		c += 4;
		if (num % c == 0)
			return false;
		c += 2;
		if (num % c == 0)
			return false;
		c += 4;
		if (num % c == 0)
			return false;
		c += 2;
		if (num % c == 0)
			return false;
		c += 4;
		if (num % c == 0)
			return false;
		c += 6;
		if (num % c == 0)
			return false;
		c += 2;
		if (num % c == 0)
			return false;
		c += 6;
	}
	return true;
}