首页 > 代码库 > POJ 3292 Semi-prime H-numbers (仿素数筛)
POJ 3292 Semi-prime H-numbers (仿素数筛)
题目地址:POJ 3292
先利用素数筛的原理把H_prime筛出来,然后打表。要预处理,否则TLE。
代码如下:
#include <iostream> #include <string.h> #include <math.h> #include <queue> #include <algorithm> #include <stdlib.h> #include <map> #include <set> #include <stdio.h> using namespace std; #define LL __int64 #define pi acos(-1.0) const int mod=100000000; const int INF=0x3f3f3f3f; const double eqs=1e-8; int h_prime[1000000], tot, check[1000100], vis[1000002], num[1000002]; void init() { int i, j; LL x, ans=0; tot=0; memset(check,0,sizeof(check)); for(i=5;i<=1000001;i+=4){ if(!check[i]){ h_prime[tot++]=i; } for(j=0;j<tot;j++){ if(i*h_prime[j]>1000001) break; check[i*h_prime[j]]=1; if(i%h_prime[j]==0) break; } } for(i=0;i<tot;i++){ for(j=0;j<tot;j++){ x=h_prime[i]*h_prime[j]; if(x>1000001) break; if(!vis[x]){ vis[x]=1; } } } for(i=0;i<1000001;i++){ if(vis[i]) ans++; num[i]=ans; } } int main() { int n, i, pos, j; LL ans, x; init(); while(scanf("%d",&n)!=EOF&&n){ printf("%d %d\n",n,num[n]); } return 0; }
POJ 3292 Semi-prime H-numbers (仿素数筛)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。