LightOJ1197 Help Hanzo 区间筛

LightOJ1197 Help Hanzo

标签

• 区间筛法

前言

• 我的csdn和博客园是同步的，欢迎来访danzh-博客园~

简明题意

• 求区间质数个数（区间长度<=1e6）

思路

• 区间[l,r]中的和数，质因子一定有一个<=sqrt( r)。因此提前筛出[1,sqrt( r)]中的质数，然后枚举这些质数的倍数就可以了。

注意事项

• 注意在区间中移动的指针应该设成long long
• 特判一下l=1的情况

总结

• 无

AC代码

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;

const int maxn = 1e5 + 10;

int mod;

bool no_prime[maxn];
int prime[maxn];
int shai(int n)
{
    int cnt = 0;
    no_prime[1] = 1;

    for (int i = 2; i <= n; i++)
    {
        if (!no_prime[i])
            prime[++cnt] = i;

        for (int j = 1; j <= cnt && prime[j] * i <= n; j++)
        {
            no_prime[prime[j] * i] = 1;
            if (i % prime[j] == 0) break;
        }
    }
    return cnt;
}

bool no_prime1[(int)1e6 + 10];
int cal(int l, int r)
{
    memset(no_prime1, 0, sizeof no_prime1);
    int cnt = 0;
    for (int i = 1; prime[i] <= sqrt(r); i++)
    {
        for (long long j = max(l % prime[i] == 0 ? l : (l / prime[i] + 1) * prime[i], 2 * prime[i]); j <= r; j += prime[i])
            if (no_prime1[j - l] == 0)
                no_prime1[j - l] = 1, cnt++;
    }
    return r - l + 1 - cnt + (l == 1 ? -1 : 0);
}

void solve()
{
    shai(maxn - 10);

    int t;
    scanf("%d", &t);
    for (int i = 1; i <= t; i++)
    {
        int l, r;
        scanf("%d%d", &l, &r);

        printf("Case %d: %d\n", i, cal(l, r));
    }
}

int main()
{
    freopen("Testin.txt", "r", stdin);
    solve();
    return 0;
}