poj 3292 Semi-prime H-numbers

题意： 
一个H-number是所有的模四余一的数。(x=4*k+1) 
如果一个H-number是H-primes 当且仅当它的因数只有1和它本身（除1外）。
一个H-number是H-semi-prime当且仅当它只由两个H-primes的乘积表示。
H-number剩下其他的数均为H-composite。
给你一个数h,问1到h有多少个H-semi-prime数。

 

Sample Input

 

21 
85
789
0

 

Sample Output

 

21 0
85 5
789 62
刷水题练手
首先一个数i是H数，那么i(4x+1)也是H数，且是H合数
显然x不能为0
因为i%4=1，(4x+1)%4=1
复杂度是o（N/4*(lnN/4)/4）≈800000

#include<iostream>
#include<cstdio>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long lol;
bool b[1000002];
int semi[1000002];
lol ans[1000002]; 
int main()
{lol i,j;
int x;
  for (i=5;i<=1000001;i+=4)
    {
      for (j=5;j*i<=1000001;j+=4)
    {
      if (semi[i]==0&&semi[j]==0)
        semi[i*j]=1;
      else semi[i*j]=-1;
    }
    } 
  for (i=1;i<=1000001;i++)
    if (semi[i]==1)
    ans[i]=ans[i-1]+semi[i];
    else ans[i]=ans[i-1];
  while (cin>>x && x)
    {
      printf("%d %lld\n",x,ans[x]);
    }
}