Semi-prime H-number(poj3292)

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

wspl654321

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

阅读量273

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/qq_36799943/article/details/72677434

版权

数论专栏收录该内容

25 篇文章 0 订阅

订阅专栏

#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;
    }
}