#pragma GCC optimize("Ofast") #include using namespace std; using ll=long long; ll counting_primes(ll x){ if(x<=70){ vector p_table(x+1,1); p_table[0]=0; p_table[1]=0; for(int i=2;i*i<=x;i++){ if(p_table[i]==0){ continue; } for(int j=2;i*j<=x;j++){ p_table[i*j]=0; } } ll ret=0; for(int i=1;i<=x;i++){ if(p_table[i]){ ret++; } } return ret; } //P2(x,a)を計算 int N=1; while(((ll)N+1)*((ll)N+1)*((ll)N+1)<=x){ N++; } int NN=(int)pow(x,2.0/3.0);//誤差が怖い vector p_table(N+1,1); vector P; p_table[0]=0; p_table[1]=0; for(int i=2;i*i<=N;i++){ if(p_table[i]==1){ for(int j=2;i*j<=N;j++){ p_table[i*j]=0; } } } for(int i=1;i<=N;i++){ if(p_table[i]==1){ P.push_back(i); } } int n=P.size(); int pi1=n; vector cnt(N+1); int pi2=-1; int x_half=1; while(((ll)x_half+1)*((ll)x_half+1)<=x){ x_half++; } ll S=0; vector I(N+1); for(int j=2;(j-1)*N+1<=NN;j++){ for(int i=1;i<=N;i++){ p_table[i]=1; } for(int i=0;i(x/((ll)min(j*N+1,NN+1)),(ll)N); ll right_tmp=min(x/(ll)((j-1)*N+1),(ll)x_half); if(left_tmp mu(N+1,1); mu[0]=0; for(int i=0;i f(N+1,0); f[1]=0; for(int i=0;i0;i++,tmp_N>>=1){ m++; } vector phi(n+1,0); vector> a(m+1,vector(N+2,0)); for(int k=1;(k-1)*N+1<=NN;k++){ for(int i=0;i<=m;i++){ for(int j=1;j<=(N/(1LL<(0,x/(((ll)min(k*N,NN)+1)*(ll)P[b])); int right=(int)min((ll)N,x/((ll)((k-1)*N+1)*(ll)P[b])); for(int i=left+1;i<=right;i++){ if(f[i]>P[b]&&mu[i]!=0){ int y=(int)(x/((ll)i*(ll)P[b])); int l=y-(k-1)*N; ll sum=phi[b]; for(int i=0;l>0;i++,l>>=1){ if((l&1)==1){ sum+=(ll)a[i][l]; } } //cout<0;i++,l>>=1){ if((l&1)==1){ phi[b]+=(ll)a[i][l]; } } int index=((k-1)*N/P[b]+1)*P[b]; while(index<=k*N&&index<=NN){ if(a[0][index-(k-1)*N]==1){ int l=index-(k-1)*N; for(int i=0;i<=m;i++){ a[i][(l+(1<>x; //cout<>L>>R; cout<