欧拉线性筛模板

inline void linear_shaker() {
bool vis[maxn];
memset(vis,false,sizeof(vis));
phi[1]=1;
for (int i=2;i<maxn;++i) {
if (!vis[i]) prime[++tot]=i,phi[i]=i-1;
for (int j=1;j<=tot&&i*prime[j]<maxn;++j) {
vis[i*prime[j]]=true;
if (i%prime[j]==0) {phi[i*prime[j]]=phi[i]*prime[j];break;}
else phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
}

int sum[maxn]={0},mu[maxn]={0,1},prime[maxn],num=0;
inline void get_mu() {
memset(mu,INF,sizeof(mu)),mu[1]=1;
for (register int i=2;i<maxn;++i) {
if (mu[i]==INF)
prime[++num]=i,mu[i]=-1;
for (int j=1;j<=num&&i*prime[j]<maxn;++j) {
if (i%prime[j]==0) {mu[i*prime[j]]=0;break;}
else mu[i*prime[j]]=-mu[i];
}
}
for (register int i=1;i<maxn;++i) sum[i]=sum[i-1]+mu[i];
}

inline void linear_shaker() {
memset(vis,0,sizeof(vis));
d[1]=1,t[1]=0;
for (register int i=2;i<MAXN;++i) {
if (!vis[i]) prime[++tot]=i,d[i]=2,t[i]=1;
for (int j=1;j<=tot&&i*prime[j]<MAXN;++j) {
vis[i*prime[j]]=true;
if (i%prime[j]==0) {
d[i*prime[j]]=d[i]/(t[i]+1)*(t[i]+2);
t[i*prime[j]]=t[i]+1;
break;
}
d[i*prime[j]]=d[i]<<1;
t[i*prime[j]]=1;
}
}
for (register int i=2;i<MAXN;++i) d[i]+=d[i-1];
}