POJ 3292 Semi-prime H-numbers（另类筛选）

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

题意：给出一个数n，问从1~n，有多少个数，被拆分成除1以外的被4除余1的数的时候，只有两项。比如125=5*5*5，被拆成了三项，就不符合条件。

由于有多组数据，而且可能数据挺大，为了避免重复运算，所以选择打表，一次性算完最大的数据。对于一个除1以外被4除余1的数，如果它不能够被拆分，那么它乘以比它小的不能够被拆分的数，刚好就是满足题目条件的数。如果它可以被拆分，那么它乘以比它小的除1以外被4除余1的数，一定是不满足题目条件的数。

#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
#define N 1000009
using namespace std;
int cnt, sum, flag[N], a[N], n, ans[N];
bool is_prime[N];
int main()
{
    cnt = 0;
    for (int i = 1; i <= 250000; i++)
    {
        cnt++;
        a[cnt] = i*4+1;
    }
    memset(flag, 0, sizeof(flag));
    for (int i = 1; i <= cnt; i++)
    {
        if (flag[a[i]] == 0)
        {
            for (int j = 1; j <= i; j++)
                if (flag[a[j]] == 0)
            {
                if (a[i] * a[j] > 1000001) break;
                flag[a[i]*a[j]] = 1;
            }
        }
        else
            {
                for (int j = 1; j <= i; j++)
                {
                    if (a[i] * a[j] > 1000001) break;
                    flag[a[i]*a[j]] = -1;
                }
            }
    }
    ans[1] = 0;
    for (int i = 2; i <= 1000001; i++)
        if (flag[i] == 1)
        ans[i] = ans[i-1]+1;
    else ans[i] = ans[i-1];
    while (~scanf("%d", &n) && n != 0)
    {
        printf("%d %d\n", n, ans[n]);
    }
    return 0;
}