码迷,mamicode.com
首页 > 其他好文 > 详细

POJ 3292 Semi-prime H-numbers

时间:2017-08-22 01:42:20      阅读:186      评论:0      收藏:0      [点我收藏+]

标签:iter   content   table   edit   display   stream   int   ons   semi   

题意:

H_Number 是一个比4的倍数多1的数,即4n + 1。H_Number 分为 H_Prime 和 H_Comosite。其中 H_Prime 仅能由1×h组成,而 H_Composite 除了1和h可有多个因子。H_Semi_Prime 表示仅有两个 H_Prime 因子(除了1和本身)。

问:给一个 H_Number,求出1到 H_Number 中有多少个 H_Semi_Prime。

思路:

打表,表分为三种情况,0表示H_Number,1表示H_Semi_Prime, - 1表示H_Composite。先将表都初始化为0,然后从小到大打。

#include <iostream>
#include <cstring>
using namespace std;

const int MAX_H = 1000001;
int H_number[MAX_H + 1];
int Count[MAX_H];
int H;
void init() {
	memset(H_number, 0, sizeof(H_number)); // H_number[x]==0表示x是H_Prime;
	for (int i = 5; i <= MAX_H; i += 4)
	for (int j = 5; j <= MAX_H; j += 4) {
		if (i*j > MAX_H) break;
		if (H_number[i] == 0 && H_number[j] == 0)
			H_number[i*j] = 1; // 表示i*j为H_semi_Prime
		else
			H_number[i*j] = -1; // 表示i*j为H_Composite
	}
	Count[1] = 0; // 统计
	for (int i = 5; i <= MAX_H; i += 4)
	if (H_number[i] == 1) Count[i] = Count[i - 4] + 1;
	else Count[i] = Count[i - 4];
}
void solve() {
	printf("%d %d\n", H, Count[H]);
}
int main()
{
	init();
	while (cin >> H && H) {
		solve();
	}
	return 0;
}

POJ 3292 Semi-prime H-numbers

标签:iter   content   table   edit   display   stream   int   ons   semi   

原文地址:http://www.cnblogs.com/demian/p/7407317.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!