Morass

Sieve

Mar 11th, 2016
118
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.16 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define PB push_back
  4. #define ZERO (1e-10)
  5. #define INF (1<<30)
  6. #define LINF (1LL<<60)
  7. #define NINF -INF
  8. #define CL(A,I) (memset(A,I,sizeof(A)))
  9. #define add(A,B) (g[A].PB(B),g[B].PB(A))
  10. #define ADD(A,B,V) (add(A,B),v[A].PB(V),v[B].PB(V))
  11. #define DEB printf("DEB!\n");
  12. #define D(X) cout<<"  "<<#X": "<<X<<endl;
  13. #define EQ(A,B) (A+ZERO>B&&A-ZERO<B)
  14. typedef long long ll;
  15. typedef pair<ll,ll> pll;
  16. typedef vector<int> vi;
  17. typedef pair<int,int> ii;
  18. #define IN(n) int n; cin >> n;
  19. #define FOR(i, m, n) for (int i(m); i < n; i++)
  20. #define REP(i, n) FOR(i, 0, n)
  21. #define F(n) REP(i, n)
  22. #define FF(n) REP(j, n)
  23. #define FT(m, n) FOR(k, m, n)
  24. #define PR(A,X) (X?A.second:A.first)
  25. #define MX (86600000LL)
  26. #define IS(A) ((cn[A>>6]>>(A&63))&1)
  27. #define ST(A) (cn[A>>6]|=1ULL<<((A&63)))
  28. int P[1<<23]={2},L(1);
  29. ll cn[1<<23];
  30. void gen(){
  31.     ST(0),ST(1);
  32.     int Q(1+sqrt(MX));
  33.     for(int k(3);k<=Q;k+=2)if(!IS(k))
  34.         for(ll h(k<<1),j(k*k);j<MX;j+=h)ST(j);
  35.     for(int i(3);i<MX;i+=2)if(!IS(i))P[L++]=i;
  36. }
  37. int tt,N;
  38. int main(void){
  39.     gen();
  40.     IN(tt)F(tt)scanf("%d",&N),printf("%d\n",P[N-1]);
  41.     return 0;
  42. }
Advertisement
Add Comment
Please, Sign In to add comment