poj 3292 Semi-prime H-numbers

最新推荐文章于 2019-04-13 10:00:00 发布

算球？

最新推荐文章于 2019-04-13 10:00:00 发布

阅读量215

点赞数

CC 4.0 BY-SA版权

分类专栏： ----数论----

本文链接：https://blog.youkuaiyun.com/gyhguoge01234/article/details/78165799

----数论---- 专栏收录该内容

41 篇文章

订阅专栏

本文介绍了一种使用埃氏筛法改进版求解质数的方法，并通过C语言实现。该程序首先初始化一个足够大的整数数组来标记是否为质数，然后利用已知的质数去筛掉其倍数的合数。最后统计小于等于某个数的所有质数个数。

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

模仿埃氏筛法走一发。

#include <stdio.h>
#include <string.h>

const int MAXN = 1001000;
int isPrime[MAXN];
int prime[MAXN];
int mark[MAXN];
int cnt[MAXN];
int plen = 0;

void init()
{
    for(int i = 5; i < MAXN; i += 4)
    {
        if(!isPrime[i])
        {
            prime[plen++] = i;
            for(int j = i*5; j < MAXN; j += 4*i)
                isPrime[j] = 1;
        }
    }

    for(int i = 0; i < plen; ++i)
    {
        for(int j = 0; j <= i; ++j)
        {
            if(prime[i]*prime[j] > MAXN) break;
            mark[prime[i]*prime[j]] = 1;
        }
    }

    for(int i = 1; i < MAXN; ++i)
        cnt[i] = cnt[i-1]+mark[i];
}

int main()
{
    init();
    int num;
    while(scanf("%d",&num) && num)
    {
        printf("%d %d\n",num,cnt[num]);
    }
    return 0;
}