数论·素筛·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
<< endl;
}
}