Semi-prime H-number(poj3292)

最新推荐文章于 2021-01-06 17:55:30 发布

原创最新推荐文章于 2021-01-06 17:55:30 发布 · 299 阅读

0 ·

CC 4.0 BY-SA版权

数论专栏收录该内容

25 篇文章

订阅专栏

本文介绍了一个基于素数筛选与同余定理的算法实现，通过C++代码详细展示了如何找出特定范围内的合成数，并统计这些合成数的数量。该算法首先通过筛选找到素数，再利用这些素数生成合成数。

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

#include <cstdio>
#include <iostream>
using namespace std;
const int H_MAX=1000001+16;
int f_prime[H_MAX],mix_prime[H_MAX];
int H_prime[H_MAX];
int cnt[H_MAX];
int tot;
void cl()
{
    int j,i;
    for(i=5;i<H_MAX;i+=4)
    {
        if(f_prime[i]) continue;
        H_prime[++tot]=i;
        for(j=i*5;j<H_MAX;j+=i*4)
         f_prime[j]=1;
    }

    for(i=1;i<=tot;i++)
     for(j=1;j<=i&&H_prime[i]*H_prime[j]<H_MAX;j++)
      mix_prime[H_prime[i]*H_prime[j]]=1;

    for(i=1;i<H_MAX;i++)
     cnt[i]=cnt[i-1]+mix_prime[i];
}
int main()
{
    cl();

    while(1)
    {
        int x;

        cin>>x;

        if(x) printf("%d %d\n",x,cnt[x]);
        else return 0;
    }
}