数论·素筛·POJ 3292·Semi-prime H-numbers

题目大意：

给定一个H-number，问有多少H-semi-primes：

定义一种数叫H-numbers，它是所有能除以四余一的数。

在H-numbers中分三种数：

1、H-primes，这种数只能被1和它本身整除，不能被其他的H-number整除，例如9是一个H-number，能被1，3，9整除，但3不是H-number，所以他是H-primes。

2、H-semi-primes是由两个H-primes相乘得出的。

3、剩下的是H-composite。

问给一个数，求1到这个数之间有多少个H-semi-primes。

解题思路：

像打素数表一样打出这个semi表；

AC代码;

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
#define   maxn          1010
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   ms(x,y)      memset(x,y,sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   PI           pair<int,int>
//#define   mp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10
typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;
typedef   pair<int ,int > P;
//#define N 100
const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;
const  ll    INF = (ll)1e18+300;
const int   maxd = 1000000 + 10;

int vis[maxd];
int pri[maxd];

void init() {
    int n = maxd;
    ms(vis, 0);
    ms(pri, 0);
    for (int i = 5; i < n; i += 4) {
        for (int j = 5; j < n; j += 4) {
            if(i * j > n) {
                break;
            }
            if(!vis[j] && !vis[i]) {
                vis[i*j] = 1;
            }
            else{
                vis[i*j] = -1;
            }
        }
    }
    int num = 0;
    for (int i = 0; i < n; i++) {
        if(vis[i] == 1) {
            num ++;
        }
        pri[i] = num;
    }
}

int main() {
    int n;
    init();
    while(cin >> n && n) {
        cout << n << " " << pri[n] << endl;
    }
}