素數篩法

生成MAX_N內的素數表

#define MAX_N XXXXXXX
#define MAX_PRIME MAX_N
vector primes;
vector is_prime;
void init_primes()
{
	is_prime = vector(MAX_PRIME + 1, true); //初始化爲true
	is_prime[0] = is_prime[1] = false;
	for (int i = 2; i <= MAX_PRIME; ++i)
	{
		if (is_prime[i])   //優化
		{
			primes.push_back(i); 
			for (int j = i * 2; j <= MAX_PRIME; j += i)
			{
				is_prime[j] = false;   //將素數的倍數2倍、3倍等全都賦值爲false
			}
		}
	}
}