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
, a
, n, ans
;
bool is_prime
;
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
);
}
return 0;
}