hiho一下 第九十三周 数论二·Eular质数筛法

快速得求解[1,N]这个区间内素数的个数

用到了Eular质数筛法

#include <stdio.h>
#include <string.h>
#define MAX_N 1000000
#define MAX_P 78498

char isPrime[MAX_N+1];
int primeList[MAX_P+1];

int main(){
	int N, i, j, primeCount;
	
	scanf("%d", &N);
	memset(isPrime, 1, sizeof(isPrime));
	primeCount = 0;
	for(i = 2; i <= N; i++){
		if(isPrime[i]){
			primeCount++;
			primeList[primeCount] = i;
		}
		for(j = 1; j <= primeCount && i*primeList[j] <= N; j++){
			isPrime[i*primeList[j]] = 0;
			if(i % primeList[j] == 0){
				break;
			}
		}
	}
	printf("%d\n", primeCount);
	
	return 0;
}