AKOJ-1695-找素数
题意:
给定区间L,R。
计算区间中素数个数。
2 <= L,R <= 2147483647, R-L <= 1000000。
思路:
素数区间筛
先筛(2-sqrt(r))。
再用(2-sqrt(r))中的素数筛(l-r)。
代码:
1.自己写的区间筛,将筛2-sqrt(r) 分开了。
#include <iostream> #include <string> #include <queue> #include <memory.h> #include <math.h> using namespace std; typedef long long LL; const int MAXN = 1000010; int isPrimersmall[MAXN]; int isPrimer[MAXN]; int main() { LL l,r; scanf("%lld%lld",&l,&r); isPrimersmall[0] = 1; isPrimersmall[1] = 1; int x = sqrt(r); for (int i = 2;i<=x;i++) { if (isPrimersmall[i] == 0) { for (int j = i * i; j <= x; j += i) isPrimersmall[j] = 1; } } for (int i = 2;i<=x;i++) { if (isPrimersmall[i] == 0) { for (int j = (int)((l + (i-1))/ i) * i;j <= r; j += i) { if (j != i) isPrimer[j-l] = 1; } } } int sum = 0; for (int i = 0;i<=r-l;i++) { if (isPrimer[i] == 0) sum++; } printf("%d\n",sum); return 0; }
2.将筛(2-sqrt(r))和(l-r)放在一起筛。
#include <iostream> #include <string> #include <queue> #include <memory.h> #include <math.h> using namespace std; typedef long long LL; const int MAXN = 1000010; int isPrimersmall[MAXN]; int isPrimer[MAXN]; int segement_sieve(LL a,LL b) { for (LL i = 0;i*i <= b;i++) isPrimersmall[i] = 1; for (LL i = 0;i <= b-a;i++) isPrimer[i] = 1; for (LL i = 2;i*i <= b;i++) { if (isPrimersmall[i]) { for (LL j = i*2;j*j <= b;j+=i) isPrimersmall[j] = 0; for (LL j = max(2LL,(a+i-1)/i) * i;j <= b ;j+=i) isPrimer[j-a] = 0; } } int sum = 0; for (int i = 0;i <= b-a;i++) if (isPrimer[i]) sum++; return sum; } int main() { LL a,b; scanf("%lld%lld",&a,&b); printf("%d\n",segement_sieve(a,b)); return 0; }