# poj 3292 Semi-prime H-numbers （打表 同余模运算）

H-prim为除了1和本身外不能被任何H-number整除的数

H-semi-prime为两个H-prim的乘积

#include <iostream>
#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <limits>
#include <stack>
#include <vector>
#include <map>

using namespace std;

#define N 1000020
#define INF 0xfffffff
#define PI acos (-1.0)
#define EPS 1e-8
#define met(a, b) memset (a, b, sizeof (a))

typedef long long LL;

int Hprim[N];

void Init ()
{
for (int i=5; i<N; i+=4)
{
for (int j=5; j<N; j+=4)
{
if (i*j>=N) break;

if (!Hprim[i] && !Hprim[j] && i*j%4==1)
Hprim[i*j] = 1;///<span style="font-family: KaiTi_GB2312;font-size:18px; line-height: 26px;">H-semi-prime</span>
else Hprim[i*j] = -1;
}
}

int cnt = 0;
for (int i=5; i<N; i++)
{
if (Hprim[i]==1) cnt++;
Hprim[i] = cnt;
}
}

int main ()
{
int n;
Init();
while (scanf ("%d", &n), n)
{
printf ("%d %d\n", n, Hprim[n]);
}
return 0;
}


#### POJ 3292 Semi-prime H-numbers

2013-12-28 13:45:58

#### POJ 3292 Semi-prime H-numbers（打表）

2013-10-18 19:58:41

#### Semi-prime H-numbers 筛选法

2013-09-28 09:04:30

#### POJ 3292 Semi-prime H-numbers

2014-07-21 11:15:31

#### Semi-prime H-numbers(POJ 3292)

2018-01-27 21:07:08

#### POJ 3292 Semi-prime H-numbers

2017-09-15 18:25:27

#### Semi-prime H-numbers POJ - 3292

2017-01-29 01:26:35

#### poj 3292 Semi-prime H-numbers

2014-07-23 21:31:08

#### 【POJ 3292】 Semi-prime H-numbers

2015-08-20 09:02:25

#### 【poj 3292】Semi-prime H-numbers

2016-03-31 15:43:30