素数筛选模板

埃式筛法

首先将2到n范围内的整数写下来,其中2是最小的素数。将表中所有的2的倍数划去,表中剩下的最小的数字就是3,他不能被更小的数整除,所以3是素数。再将表中所有的3的倍数划去……以此类推,如果表中剩余的最小的数是m,那么m就是素数。然后将表中所有m的倍数划去,像这样反复操作,就能依次枚举n以内的素数,这样的时间复杂度是O(nloglogn)。

代码:

memset(check, 0, sizeof(check)); //check为0代表素数,1代表非素数
int tot = 0;
for (int i = 2; i <= n; ++i)
{
  if (!check[i])
  {
    prime[tot++] = i;
  }
  for (int j = i+i; j <= n; j += i)
  {
    check[j] = 1;
  }
}

线性筛(欧拉筛法)

每个合数只会被它的最小质因数筛去,因此每个数只会被标记一次,所以时间复杂度是O(n)

代码:

const int MAXN = 10000;
int prime[MAXN+1];
bool is_pri[MAXN*100];
int tot;
void getprime() {
	memset(is_pri, false, sizeof(is_pri));
	for(int i = 2; i <= MAXN; i++) {
		if(!is_pri[i])
			prime[tot++] = i;
		for(int j = 0; j < tot; j++) {
			if(i * prime[j] > MAXN)
				break;
			is_pri[i*prime[j]] = true;
			if(i % prime[j] == 0)
				break;	
		}	
	}
}

猜你喜欢

转载自blog.csdn.net/Sclong0218/article/details/83382218