POJ
来源:互联网 发布:中国微观数据库 编辑:程序博客网 时间:2024/06/07 08:41
点击打开题目链接
Description
This problem is based on an exercise of David Hilbert, who pedagogically suggested that one study the theory of 4n+1 numbers. Here, we do only a bit of that.
An H-number is a positive number which is one more than a multiple of four: 1, 5, 9, 13, 17, 21,... are the H-numbers. For this problem we pretend that these are the only numbers. The H-numbers are closed under multiplication.
As with regular integers, we partition the H-numbers into units, H-primes, and H-composites. 1 is the only unit. An H-number h is H-prime if it is not the unit, and is the product of two H-numbers in only one way: 1 × h. The rest of the numbers are H-composite.
For examples, the first few H-composites are: 5 × 5 = 25, 5 × 9 = 45, 5 × 13 = 65, 9 × 9 = 81, 5 × 17 = 85.
Your task is to count the number of H-semi-primes. An H-semi-prime is an H-number which is the product of exactly two H-primes. The two H-primes may be equal or different. In the example above, all five numbers are H-semi-primes. 125 = 5 × 5 × 5 is not an H-semi-prime, because it's the product of three H-primes.
Input
Each line of input contains an H-number ≤ 1,000,001. The last line of input contains 0 and this line should not be processed.
Output
For each inputted H-number h, print a line stating h and the number of H-semi-primes between 1 and h inclusive, separated by one space in the format shown in the sample.
Sample Input
21 857890
Sample Output
21 085 5789 62
题目不难,只是借这个题写一下艾式筛法。
艾式筛法求素数代码(对于每一个素数,它的N倍一定不是素数):
for(int i=2;i<maxn;i++){ if(!vis[i]) { vis[i]=1; prime[cnt++]=i; for(int j=i*i;j<maxn;j*=i) vis[j]=1; }}
附上AC代码:
#include<iostream>using namespace std;const int maxn=1000000+5;int cnt[maxn];int a;void is_Hprime(){ for(int i=5;i<maxn;i+=4) { for(int j=5;j<maxn;j+=4) { if(i*j>maxn)break; if(!cnt[i]&&!cnt[j]) { cnt[i*j]=1; } else cnt[i*j]=-1; } } int num=0; for(int i=25;i<maxn;i++) { if(cnt[i]==1) num++; cnt[i]=num; }}int main(){ is_Hprime(); while(cin>>a,a) cout<<a<<' '<<cnt[a]<<endl;}
- POJ
- poj
- POJ
- POJ
- poj
- poj
- POJ
- POJ
- poj
- POJ
- POJ
- POJ
- POJ
- POJ
- POJ
- POJ
- POJ
- POJ
- android网页打印,安卓网页打印,h5页面打印,浏览器打印,js打印工具
- 排序小结(二)
- Nginx反向代理与负载均衡
- 进程间通信系列(3)匿名管道范例代码
- 今日头条2017春招研发岗笔试题——Paragraph
- POJ
- Android之冷启动和热启动,以及代码
- MySql基础
- C# 消息处理学习总结
- 实现24点游戏-穷举法
- 抽象类、抽象方法和内部类
- 移动端开发-WebApp-淘宝商城(响应式布局:任何手机屏幕尺寸自适应)
- android 蓝牙4.0开--发关于BLE和iBeacon的一些笔记
- linux中nm、ldd、readelf命令