poj 3292 Semi-prime H-numbers

本文介绍了一种求解特定范围内H-合成数数量的算法。形如4n+1(n≥1)的数被称为H数,其中H-素数是指不能被其他H数整除的数,而H-合成数则是两个H-素数的乘积。文章通过筛选H-素数并枚举其乘积来找出所有H-合成数。

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

http://poj.org/problem?id=3292

 

题意:

形如4n+1(n>=1)的数称为H数

H数中,不能被除本身之外的H数整除的数 称为H-素数;能被除本身之外的H数整除的数 称为H-合数

两个H-素数(可以相同)的乘积 称为H-合成数

给定n,求n以内有多少H-合成数

 

猜想这样的数不多,考虑先用筛法求素数的方法筛出H-素数

然后两两枚举H-素数,从而求出H-合成数

筛 H-素数:

已知 i=4n+1(n>=1) 是一个H-素数

那么i+4*i*x (x>=1)一定是H数不是H-素数

i+4*i*x= 4n+1+16*x*n+4x=4*(n+4*x*n+x)+1

i+4*i*x=i*(1+4x)

 

#include<cstdio>

using namespace std;

#define N 1000002

bool is_H_prime[N];
int H_prime[N],m;

int Semi_prime[N];

int main()
{
    for(int i=5;i<N;i+=4)
    {
        if(is_H_prime[i]) continue;
        H_prime[++m]=i;
        for(int j=i*5;j<N;j+=4*i) is_H_prime[j]=true;
    }
    for(int i=1;i<=m;++i)
        for(int j=1;j<=i;++j)
        {
            if(H_prime[i]*H_prime[j]>=N) break;
            Semi_prime[H_prime[i]*H_prime[j]]=true;
        }
    for(int i=1;i<N;++i) Semi_prime[i]+=Semi_prime[i-1];
    int n;
    while(1)
    {
        scanf("%d",&n);
        if(!n) return 0;
        printf("%d %d\n",n,Semi_prime[n]);
    }
}

 

转载于:https://www.cnblogs.com/TheRoadToTheGold/p/8503700.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值