埃式筛法

埃式筛法

• 关于单个素数的判定一般时间复杂度为O（n^1/2）
• 但要枚举n以内的素数时，采用循环判断时，的时间复杂度较高
• 因此对于枚举2到n的素数时要用到埃式筛法

1.求n以内的素数

我们以20以内的数为例

int solve(int n) {
	int ans=0;
	memset(prime, true, n+1);
	prime[0] = prime[1] = false;
	for (int i = 2; i <= n; i++) {
		if (prime[i]) {
			ans++;
			for (int j = 2 * i; j <= n; j *= 2) prime[j] = false;
		}
	}
	return ans;
}

• 这里有两个初始化数组的方法

#include"algorithm"
fill(begin,end,value);

#include"string.h"
memset(name,value,sizeof());

---

2.求某区间内的素数

给定整数a与b,请问区间[a,b)内有多少个素数？

• 若采用1的方法，则会造成无用的计算
• 因为一个数b的最大的约束是 根号b,那么我们利用 根号b以内的数来 筛选[a,b)内的数

如图：示例为[a,b)的筛选
利用prime1来筛选prime2

#include"stdio.h"
#include"stdlib.h"
#include"string.h"
#include"math.h"
#define MAX 100
bool prime1[MAX];
bool prime2[MAX];

int solve(int a,int b) {
	int ans = 0;
	memset(prime1, true, b-a);
	memset(prime2, true, (int)sqrt(b)+1);
	prime2[0] = prime2[1] = false;
	for (int i = 2; i*i < b; i++) {
		if (prime2[i]) {
			for (int j = 2 * i; j*j < b; j *= 2) prime2[j] = false;
			for (int j = a; j < b; j++)	 {
				if (j%i == 0)  prime1[j - a] = false;
				}
			}
		}
	for (int i = 0; i < b - a; i++) {
		if (prime1[i] == true) {
			ans++;
		}
   }
	return ans;
}
int main() {
	int a,b;
	scanf_s("%d%d", &a,&b);
	printf("%d", solve(a,b));
	system("pause");
}