题意
传送门 POJ 3292
仿埃氏筛法筛 ,再求小于 的 ,二分答案即可。需要注意的是, 非真正的素数,可能会存在不同的组合得到相同的乘积,二分答案前要去重。
#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;
}