POJ 3293（素数筛变形+思维）

题意

对于一个正整数，如果它模4得1，则称这个数为H数；如果H数是质数，则称为质H数；如果一个H数是由两个质H数相乘得到的，就称之为H-semi-prime数。现在给定一个n，求1-n有少个H-semi-prime数。

分析

利用素数筛的思想，枚举两个H数，让其从小到大相乘，一开始从5*5,5*9..取的两个数一定是质H数，得到的就是要求的数，然后打个标记。在代码中，v数组保存的信息，如果v[i]=0表示i没有被筛到，那么他一定是素数；或者v[i]要么能被两个素数相乘得到，这时候就标记1，是题目要求的数；要么是被两个数a,b相乘得到，且a,b中至少一个是合数，这样的我们给标记为-1，表示它既不是素数，也不是题目要求的合数。也就是说我们把筛出来的合数打上两个标记，一个是满足题目要求的合数，另一个是不满足的合数

代码

#include <iostream>
#include <cmath>
#include <cstring>
using namespace std;
const int MAXN = 1000001;
int v[MAXN + 5], h[MAXN + 5];
void is_H()
{
    memset(v, 0, sizeof(v));
    for (int i = 5; i <= MAXN; i += 4)
        for (int j = 5; j <= MAXN; j += 4)
        {
            if (j*i > MAXN) break;
            else if (!v[i] && !v[j]) v[i*j] = 1;
            else v[i*j] = -1;
        }
    for (int j = 5; j <= MAXN; j++)
    {
        if (v[j] == 1) h[j] = h[j - 1] + 1;
        else h[j] = h[j - 1];
    }
}
int main()
{
    int n;
    is_H();
    while (cin >> n && n)
    {
        cout << n << " " << h[n] << endl;
    }
    return 0;
}