ural 1086

来源:互联网 发布:30岁程序员找不到工作 编辑:程序博客网 时间:2024/06/05 00:38

题目大意:输入T个正整数k,输出第k个素数。

Time Limit:2000MS     Memory Limit:65536KB     64bit IO Format:%I64d & %I64u

规模: T为正整数,0<ki<=15000。

理论基础:仅能被自身和1整除的正整数叫做素数,特别的,1既不是素数,也不是质数。

题目分析:水题,预处理出前15000个素数即可。

代码如下:

#include<iostream>#include<cstring>#include<string>#include<cstdlib>#include<cstdio>#include<cmath>#include<algorithm>#include<queue>using namespace std;typedef double db;#define DBG 1#define maa (1<<31)#define mii ((1<<31)-1)#define sl(c) ((int)(c).size())    //取字符串长度;#define forl(i, a, b) for(int i = (a); i <  (b); ++i)    //不带边界值循环,正序#define forle(i, a, b) for(int i = (a); i <= (b); ++i)   //带边界值循环,正序#define forh(i, a, b) for(int i = (a); i >  (b); --i)     //不带边界值,逆序#define forhe(i, a, b) for(int i = (a); i >= (b); --i)        //带边界值,逆序#define forlc(i, a, b) for(int i##_b = (b), i = (a); i <  i##_b; ++i)  //带别名的循环,用于不可改变值#define forlec(i, a, b) for(int i##_b = (b), i = (a); i <= i##_b; ++i)#define forgc(i, a, b) for(int i##_b = (b), i = (a); i >  i##_b; --i)#define forgec(i, a, b) for(int i##_b = (b), i = (a); i >= i##_b; --i)#define forall(i, v   )  forl(i, 0, sz(v))   //循环所有#define forallc(i, v   ) forlc(i, 0, sz(v))#define forlla(i, v   ) forhe(i, sz(v)-1, 0)#define forls(i, n, a, b) for(int i = a; i != b; i = n[i])   //搜表用#define rep(n)  for(int               repp = 0; repp <    (n); ++repp)#define repc(n) for(int repp_b = (n), repp = 0; repp < repp_b; ++repp)#define rst(a, v) memset(a, v, sizeof a)   //把字符v填充到a  reset 重置#define cpy(a, b) memcpy(a, b, sizeof a)   //copy b 的sizeof(a)个字符到a#define rstn(a, v, n) memset(a, v, (n)*sizeof((a)[0]))  //把字符v填充到a[n]之前的字节#define cpyn(a, b, n) memcpy(a, b, (n)*sizeof((a)[0]))    //copy b 的 n 个字符到a#define ast(b) if(DBG && !(b)) { printf("%d!!|\n", __LINE__); while(1) getchar(); }  //调试#define dout DBG && cout << __LINE__ << ">>| "#define pr(x) #x"=" << (x) << " | "#define mk(x) DBG && cout << __LINE__ << "**| "#x << endl#define pra(arr, a, b)  if(DBG) {\    dout<<#arr"[] |" <<endl; \    forlec(i, a, b) cout<<"["<<i<<"]="<<arr[i]<<" |"<<((i-(a)+1)%8?" ":"\n"); \    if((b-a+1)%8) puts("");\}                                                             //数列查看#define rd(type, x) type x; cin >> x   //读数inline int     rdi() { int d; scanf("%d", &d); return d; }inline char    rdc() { scanf(" "); return getchar(); }inline string  rds() { rd(string, s); return s; }inline db rddb() { db d; scanf("%lf", &d); return d; }template<class T> inline bool updateMin(T& a, T b) { return a>b? a=b, true: false; }template<class T> inline bool updateMax(T& a, T b) { return a<b? a=b, true: false; }bool isprime[165001]={true,true};int nprime[15000];void psieve(int n){    int t=(int)sqrt((float)n),temp,i,j;    for(i=2;i<=t;i++)    {        if(!isprime[i])        {            for(j=i;(temp=(j*i))<=n;j++)isprime[temp]=true;        }    }}int main(){    int cnt=0;    int n=165001,number;    psieve(n);    for(int i=2;cnt<15000;i++)if(!isprime[i])nprime[cnt++]=i;    scanf("%d",&number);    while(number--)    {        int i=rdi();        printf("%d\n",nprime[i-1]);    }return 0;}

其中,用到了素数筛法,效率较快,但这种筛法存在重复操作,例如:18=9*2=3*6,会被筛两次。另有一种筛法,随后会见到。
by : Jsun_moon  http://blog.csdn.net/Jsun_moon