poj3292 素数筛选

题目链接点击打开链接

题意：h数为4*n + 1 的数 ， h数可以分为unit (只有数字1) ， h素数 和 h合成数， h素数类似素数的定义。 h_semi_prime为两个h_素数的积。输入h，求1到h之间有多少个h_semi_prime。

思路： 主要是刷选出h素数，可模仿 埃氏筛法 。最开始素数筛选有问题导致wrong answer

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string>
#include <string.h>
#include <math.h>
using namespace std;
const int maxn = 1e6 + 5 ;

bool is_H_prime[maxn] , is_H_si[maxn] ;
int H_prime[maxn] , p = 0 , H[maxn];

void Prime()
{
    memset(is_H_prime , true , sizeof(is_H_prime)) ;
    for(int i = 5 ; i < maxn ; i += 4)
    {
        if(is_H_prime[i])
        {
            H_prime[p++] = i ;
            for(int j = 5*i ; j <  maxn ; j += 4*i ) //attention
            {
                is_H_prime[j] = false ;
            }
          /*
            for(int j = i + i ; j < maxn ; j += i)
            {
                if((j - 1)%4 == 0)
                {
                    is_H_prime[j] = false ;
                }
            }*/
        }
    }
}

void solve()
{
    memset(is_H_si , false , sizeof(is_H_si)) ;
    for(int i = 0 ; i < 166 ; i ++)
    {
        for(int j = i ;  H_prime[j] <= maxn/H_prime[i] ; j ++)
        {
            int t = H_prime[i] * H_prime[j] ;
            if( t < maxn)
                is_H_si[t] = true ;
        }

    }

    int Count = 0 ;
    for(int i = 1 ; i < maxn ; i ++ )
    {
        if (is_H_si[i]) Count ++ ;
        H[i] = Count ;
    }
}

int main()
{
    Prime() ; solve() ; // 打表，预处理
    int number ;
    while(scanf("%d" , &number) != EOF && number)
    {
        printf("%d %d\n" , number , H[number]) ;
    }
    return 0;
}