Count prime hoj 区间求素数个数（大范围）

/*区间筛素数
简述：有的时候，我们需要知道某个特定区间的素数(区间大小较小，但数可能很大)。
那么数组就开不下，这时候我们仍然可以使用筛法，只是所有的下标都进行了偏移。
大家理解下面这段代码可以先用普通筛法写，然后数组下标集体移动即可。*/

#include <stdio.h>
#include <cmath>
#include <cstring>
const int maxn = 1000001;

int PrimeList[maxn];
int PrimeNum;
bool IsNotPrime[maxn]; // IsNotPrime[i] = 1表示i + L这个数是素数.

void SegmentPrime(int L, int U)
{
// 求区间[L, U]中的素数.
int i, j;
int SU = sqrt(1.0 * U);
int d = U - L + 1;
for (i = 0; i < d; i++) IsNotPrime[i] = 0; // 一开始全是素数.
for (i = (L % 2 != 0); i < d; i += 2) IsNotPrime[i] = 1; // 把偶数的直接去掉.
for (i = 3; i <= SU; i += 2)
{
if (i > L && IsNotPrime[i - L]) continue; // IsNotPrime[i - L] == 1说明i不是素数.
j = (L / i) * i; // j为i的倍数，且最接近L的数.
if (j < L) j += i;
if (j == i) j += i; // i为素数，j = i说明j也是素数，所以直接 + i.
j = j - L;
for (; j < d; j += i) IsNotPrime[j] = 1; // 说明j不是素数(IsNotPrime[j - L] = 1).
}
if (L <= 1) IsNotPrime[1 - L] = 1;
if (L <= 2) IsNotPrime[2 - L] = 0;
PrimeNum = 0;
for (i = 0; i < d; i++) if (!IsNotPrime[i]) PrimeList[PrimeNum++] = i + L;
}
int main()
{
int l,r;
while(scanf("%d%d",&l,&r)==2)
{
SegmentPrime(l,r);
printf("%d\n",PrimeNum);
}
return 0;
}