POJ 3292 Semi-prime H-numbers

最新推荐文章于 2021-07-26 12:46:55 发布

原创最新推荐文章于 2021-07-26 12:46:55 发布 · 454 阅读

0 ·

CC 4.0 BY-SA版权

组合数学数论专栏收录该内容

15 篇文章

订阅专栏

本文介绍了一种称为H-number的整数序列，并通过C++代码实现了判断这些整数是否为H-prime或H-semiprime的过程。文章通过筛法预处理了所有小于100万的H-semiprime数，并提供了一个查询接口来获取指定范围内H-semiprime数的数量。

H-number是指形如4n+1(n>=0)的整数，后又定义了好几种数。

水题。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
#define N 1000002
bool NotaHprime[N],IsaHsemiprime[N];
int ans[N];
int main()
{
    int i,j;
    memset(NotaHprime,false,sizeof(NotaHprime));
    memset(IsaHsemiprime,false,sizeof(IsaHsemiprime));
    for (i=5; i<1002; i+=4)
    {
        for (j=i; i*j<N; j+=4)
        {
            NotaHprime[i*j]=true;
        }
    }
    for (i=5; i<1002; i+=4)
    {
        for (j=i; i*j<N; j+=4)
        {
            if (!NotaHprime[i] && !NotaHprime[j])
                IsaHsemiprime[i*j]=true;
        }
    }
    ans[1]=0;
    for (i=2; i<N; i++)
    {
        if (IsaHsemiprime[i])
            ans[i]=ans[i-1]+1;
        else
            ans[i]=ans[i-1];
    }
    while (1)
    {
        scanf("%d",&i);
        if (i == 0)
            break;
        printf("%d %d\n",i,ans[i]);
    }
}