poj_3292

源地址：http://poj.org/problem?id=3292

题目大意不难懂。

这题其实就是个素数打表的变异，只要素数打表掌握好了，这道题目应该也不难。首先我们设一个数组,hnum[i]，如果其值等于0，我们就认为它是H-primes，那么题目中要求的H-semi-primes就是两个H-primes的乘积，那么我们令i=5,j=5开始，只要hnum[i]==0&&hnum[j]==0的话，hnum[i*j]就是所谓的H-semi-primes，我们把它记为1，否则-1.

最后统计下1-100w直接每个区间内的hnum[i]==1的情况就可以以O(1)的效率来实现计算了。

#include<stdio.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<time.h>
#include<queue>
#include<stack>
#include<iterator>
#include<math.h>
#include<stdlib.h>
#include<limits.h>
#include<set>
#include<map>
//#define ONLINE_JUDGE
#define eps 1e-8
#define INF 0x7fffffff
#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define sfs(a) scanf("%s",a)
#define sf(a) scanf("%d",&a)
#define sfI(a) scanf("%I64d",&a)
#define pf(a) printf("%d\n",a)
#define pfI(a) printf("%I64d\n",a)
#define pfs(a) printf("%s\n",a)
#define sfd(a,b) scanf("%d%d",&a,&b)
#define sft(a,b,c)scanf("%d%d%d",&a,&b,&c)
#define for1(i,a,b) for(int i=(a);i<b;i++)
#define for2(i,a,b) for(int i=(a);i<=b;i++)
#define for3(i,a,b)for(int i=(b);i>=a;i--)
#define MEM1(a) memset(a,0,sizeof(a))
#define MEM2(a) memset(a,-1,sizeof(a))
const double PI=acos(-1.0);
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
using namespace std;
#define ll long long
int n,m;
#define Mod 1000000007
#define N 510
#define M 1000100
const int size = 46340+4;
const int py = 10005;
int p;
int hnum[M];
int h;
void Init(){
	memset(hnum,0,sizeof hnum);
	for(ll i=5;i<M;i+=4){			//structure the table
		for(int ll j = 5;j<M;j+=4){
			if(i*j>=M) break;
				if(hnum[i] == 0 && hnum[j] == 0)
					hnum[i*j] = 1;
				else
					hnum[i*j] = -1;
		}
	}
	int res=0;
	for(int i=0;i<M;i++){	//record
		if(hnum[i] == 1)
			res++;
		hnum[i] = res;
	}

}
int main(){
#ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
//  freopen("out.txt","w",stdout);
#endif
    Init();
    while(scanf("%d",&h)!=EOF && h){
    	printf("%d %d\n",h,hnum[h]);
    }
return 0;
}