POJ 3292 Semi-prime H-numbers（筛法）

Semi-prime H-numbers POJ - 3292

把所有的比4的倍数多1的数（称为H-number，为1，5，9，13，……）分成三类：1，H-prime，H-composite；

其中，定义H-composite为：不是1和H-prime的H-number；

定义H-prime为：只有1和这个H-prime自己可以整除自己的H-number。

比如9，虽然是合数，但是由于因数3不是H-number，所以9也是H-prime。

现在定义Semi-prime为H-number中由两个H-prime相乘所得的数，只能由两个！！

现在给你数个n（n<=1,000,001），找出n以内H-number中的Semi-prime，输出n和Semi-prime，用空格隔开

分析：

由于有多组数据，先用筛法把所有的H-prime找出来（注意筛的时候质数找的时候上界要用最大值的开方即i*i<=1000020/4，毕竟就算之后还有质数，也不会再筛出所给范围内的新H-composite了，也就是说筛的工作已经做完了），之后对于每一个输入的n再来找Semi-prime，找Semi-prime的思路是要分解每一个H-composite，这样有可能分解出来因子里面可能有不是H-number的数，比较不好操作；

考虑在筛H-prime的过程中，筛的时候不是乘1，2，3……，而是乘第一个H-number，第二个H-number……，可以在这个时候判断所乘的第j个H-number是不是H-prime，如果是的话，那么我们找到的这个H-composite就同时是一个Semi-prime，做上特殊记号就可以了；

但是我们在筛选的过程中，比如发现了第i个H-number是H-prime，就去筛H-composite，那么i之后的H-number中总有一些是我们还没有发现的H-composite，而且也不能通过这个H-prime筛出的，虽然这是暂时的，但是务必会影响我们对Semi-prime的判断，我们可以通过先将筛出的H-composite标记为H-composite再判断Semi-prime并标记为Semi-prime，这样的话，前面如果错误判断了一个Semi-prime，后面随着还未筛出的H-composite减少，错误的Semi-prime会被不断修正，从而保证最后的结果是正确的；

但是就这样交上去还是会te，大概是测试数据比较多吧，每一次都从头遍历的话还是受不了的，可以加入一个数组seminumber[n]，专门存1~n的semi-prime数量，这样只要遍历一遍就可以填满，以后每一次查找都会是常数时间，nice！！

还有就是为了在节约点空间，考虑到每一个数的格式都是4*i+1，所以数组的下标就设定为i好了，也是很方便的

最后还要补充一点就是，底下这个代码，一开始提交的时候是Runtime Eerror的，仔细找了一下原因，果然数组还是会越界，在下标等于1000020的时候越界，只是溢出了一个单位就RE，真是坑。。。虽然自己电脑上没啥事。。。

再说一句，你会发现，如果将这里的定义平行为一般的素数，合数定义，方法是惊人一致的。。。

#include<iostream>
#include<cstring>
using namespace std;
int no_hprime[1000020/4+1],seminumber[1000020/4+1];
int main(){
	memset(no_hprime,0,sizeof(no_hprime));
	memset(seminumber,0,sizeof(seminumber));
	for(int i=1;i*i<=1000020/4;i++){
		if(no_hprime[i]==0){
			for(int j=1;4*i*j+i+j<=1000020/4;j++){
				no_hprime[4*i*j+i+j]=1;
				if(no_hprime[j]==0){
					no_hprime[4*i*j+i+j]=2;
				}
			}
		}
	}
	int sum=0;
	for(int i=1;i<=1000020/4;i++){
		if(no_hprime[i]==2)sum++;
		seminumber[i]=sum;
	}
	int n;
	cin>>n;
	while(n!=0){
		cout<<n<<' '<<seminumber[(n-1)/4]<<endl;
		cin>>n;
	}
	return 0;
}