天天看點

ural 1086

題目大意:輸入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

繼續閱讀