线性筛三合一，强大O(n)

校内CJOJ2395by Jesse Liu

筛法三合一 Euler、Möbius、Prime函数

基于数论的积性函数

gcd(a,b)=1  则  ƒ(ab)=ƒ(a)ƒ(b)

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <map>
#include <set>
using namespace std;
#define file(x) freopen(x".in","r",stdin),freopen(x".out","w",stdout)
inline void read(int &ans){
  ans=0;char x=getchar();int f=0;
  while(x<'0'||x>'9'){if(x=='-')f=1;x=getchar();}
  while(x>='0'&&x<='9')ans=ans*10+x-'0',x=getchar();
  if(f)ans=-ans;
}
 
const int maxn=(int)1e7+10;
int phi[maxn],mu[maxn],p[maxn],flag[maxn],cnt;
void sieve(int n){
  mu[1]=phi[1]=1;
  for(int i=2;i<=n;i++){
    if(!flag[i])p[++cnt]=i,mu[i]=-1,phi[i]=i-1;
    for(int j=1;j<=cnt && i*p[j]<=n;j++){
      flag[i*p[j]]=1;
      if(i%p[j]==0){mu[i*p[j]]=0,phi[i*p[j]]=phi[i]*p[j];break;}
      mu[i*p[j]]=-mu[i],phi[i*p[j]]=phi[i]*(p[j]-1);
    }
  }
}
 
int main(){
  sieve((int)1e7);
  int T,op,n;
  for(read(T);T--;){
    read(op),read(n);
    printf("%d\n",op==1?p[n]:op==2?mu[n]:phi[n]);
  }
  return 0;
}

进一步学习的建议，Jesse Liu