HOJ 10148 Counting Prime（筛素数模板）

最新推荐文章于 2020-12-03 23:59:39 发布

麻由由哒哟

最新推荐文章于 2020-12-03 23:59:39 发布

阅读量794

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM题解

本文链接：https://blog.youkuaiyun.com/luovilonia/article/details/38822835

ACM题解专栏收录该内容

135 篇文章

订阅专栏

本文介绍了一种高效的素数筛选算法，并提供了完整的C++代码实现。该算法利用标记数组存储素数状态，通过迭代筛选得到指定范围内的所有素数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

记录筛素数的板子。

AC代码：

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <cstring>

using namespace std;

const int maxn = 4194304;
const int prin = 300000;

int prime[prin];
///仅包含奇数，isprime[k]代表2*k+3
bool isprime[maxn>>1];
int ptop = 0;

void sieve()
{
        memset(isprime, true, sizeof(isprime));
        prime[0] = 2;
        int sq = sqrt(maxn);
        for(int i = 3; i < maxn; i += 2)
        {
                if(isprime[(i-3)>>1])
                {
                        if(i <= sq+3)
                        {
                                for(int j = i*i; j < maxn; j += 2*i)
                                {
                                        isprime[(j-3)>>1] = false;
                                }
                        }
                        prime[++ptop] = i;
                }
        }
        return ;
}

int bin(int x)
{
        return (lower_bound(prime, prime+ptop+1, x) - prime);
}

int main()
{
        sieve();
        int n, m;
        while(~scanf("%d%d", &n, &m))
        {
        
                int tapn = bin(n);
                int tapm = bin(m);
                //cout << tapn << ' ' << tapm << endl;
                if(m == 2 || ((m&1) && isprime[(m-3)>>1]))
            tapm++;
        printf("%d\n", tapm-tapn);
        }
        return 0;
}