POJ 3292 埃氏筛法 + 二分

题意

传送门 POJ 3292

仿埃氏筛法筛 $H-prime$，再求小于 $MAX\_H$ 的 $H-semi-prime$ ，二分答案即可。需要注意的是，$H-prime$ 非真正的素数，可能会存在不同的组合得到相同的乘积，二分答案前要去重。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
#define min(a,b)    (((a) < (b)) ? (a) : (b))
#define max(a,b)    (((a) > (b)) ? (a) : (b))
#define abs(x)    ((x) < 0 ? -(x) : (x))
#define INF 0x3f3f3f3f
#define delta 0.85
#define eps 1e-3
#define PI 3.14159265358979323846
#define MAX_N 1000005
using namespace std;
int H, cnt, h_cnt;
bool is_prime[MAX_N];
int prime[MAX_N / 2], h_prime[MAX_N];

void sieve(int n){
	int limit = sqrt((double)n);
	for(int i = 5; i <= n; i++) is_prime[i] = 1;
	for(int i = 5; i <= n; i += 4){
		if(is_prime[i]){
			prime[cnt++] = i;
			if(i <= limit){
				for(int j = i * i; j <= n; j += i) is_prime[j] = 0;
			}
		}
	}
}

int main(){
	sieve(MAX_N);
	for(int i = 0; i < cnt; i++){
		for(int j = i; j < cnt; j++){
			if(MAX_N / prime[j] >= prime[i]) h_prime[h_cnt++] = prime[i] * prime[j];
			else break;
		}
	}
	sort(h_prime, h_prime + h_cnt);
	int sz = unique(h_prime, h_prime + h_cnt) - h_prime;
	while(~scanf("%d", &H) && H){
		printf("%d %d\n", H, upper_bound(h_prime, h_prime + sz, H) - h_prime);
	}
	return 0;
}