51nod1184 第N个质数

如题。$n \leq 1e9$。

方法零：二分，然后洲阁筛。要魔改一下的洲阁筛。跑得慢。卡卡能过。没意思。

 //#include<iostream>
 #include<cstring>
 #include<cstdlib>
 #include<cstdio>
 //#include<map>
 #include<math.h>
 //#include<time.h>
 //#include<complex>
 #include<algorithm>
 using namespace std;

 #define LL long long
 LL n,mzz=;
 #define maxn 640011
 #define maxm 10000011
 int prime[maxm/],sp[maxm],lp; bool notprime[maxm];
 void pre()
 {
     lp=;
     for (int i=;i<=mzz;i++)
     {
         if (!notprime[i]) {prime[++lp]=i;}
         sp[i]=sp[i-]+!notprime[i];
         for (int j=,tmp;j<=lp && 1ll*prime[j]*i<=mzz;j++)
         {
             notprime[tmp=i*prime[j]]=;
             if (i%prime[j]); else break;
         }
     }
 }

 LL G(int i,LL j)
 {
     if (i==) return (j+)>>;
     if (j<=prime[i]) return ;
     if (j<=mzz && j<=(LL)prime[i]*prime[i]) return sp[j]+-i;
     return (G(i-,j)-G(i-,j/prime[i]));
 }

 LL mg(LL n) {return G(lp,n);}

 LL calc(LL n)
 {
     if (n<=mzz) return sp[n];
     lp=; for (;(LL)prime[lp]*prime[lp]<=n;lp++);
     return mg(n)-+sp[prime[lp]];
 }

 int q;
 int main()
 {
     scanf("%d",&q); pre();
     LL L=,R=;
     if (q>=) L=; else R=;
     if (q>=) L=; else R=;
     if (q>=) L=; else R=;
     if (q>=) L=; else R=;
     if (q>=) L=; else R=;
     while (L<R)
     {
         LL mid=(L+R)>>,tmp;
         if ((tmp=calc(mid))>=q) R=mid;
         else L=mid+;
     }
     printf("%lld\n",L);
     return ;
 }