最初的应用:线性筛素数。
http://blog.csdn.net/leolin_/article/details/6642126
http://www.cnblogs.com/grubbyskyer/p/3852421.html
void get_p(){
memset(vis,1,sizeof(vis));
for(int i=2;i<=N;i++){
if(vis[i]) p[++p[0]]=i;
for(int j=1;j<=p[0]&&i*p[j]<=N;j++){
vis[i*p[j]]=false;
if(i%p[j]==0) break;
}
}
}
如果线性筛只能用在搞素数上就没什么卵用了,因为有取模运算,所以比最简单的筛法差不了多少,范围小的话甚至更慢,10^7以上才会快点。
其实线性筛还可以求一些满足积性的数论函数,效率还是O(n)。
只要在线性筛素数代码上加点东西就行了。
甚至是一坨奇奇怪怪的东西,只要是积性的一般都能求。需要多维护几个辅助数组。
求
ϕ
:
void get_phi(){
memset(vis,1,sizeof(vis));
phi[1]=1;
for(int i=2;i<=N;i++){
if(vis[i]) p[++p[0]]=i, phi[i]=i-1;
for(int j=1;j<=p[0]&&i*p[j]<=N;j++){
vis[i*p[j]]=false;
if(i%p[j]==0){ phi[i*p[j]]=phi[i]*p[j]; break; }
phi[i*p[j]]=phi[i]*phi[p[j]];
}
}
}
求 μ :
void get_mu(){
memset(vis,1,sizeof(vis));
mu[1]=1;
for(int i=2;i<=N;i++){
if(vis[i]) p[++p[0]]=i, mu[i]=-1;
for(int j=1;j<=p[0]&&i*p[j]<=N;j++){
vis[i*p[j]]=false;
if(i%p[j]==0){ mu[i*p[j]]=0; break; }
mu[i*p[j]]=-mu[i];
}
}
}
求F(i)表示i的约数和:
void get_F(){
memset(vis,1,sizeof(vis)); p[0]=0;
F[1].Fir=1;
for(int i=2;i<=N;i++){
if(vis[i]) p[++p[0]]=i, F[i].Fir=1+i, t1[i]=1+i, t2[i]=i;
for(int j=1,k;j<=p[0]&&(k=p[j]*i)<=N;j++){
vis[k]=false;
if(i%p[j]==0){
t2[k]=t2[i]*p[j]; t1[k]=t1[i]+t2[k];
F[k].Fir=F[i].Fir/t1[i]*t1[k];
break;
}
F[k].Fir=F[p[j]].Fir*F[i].Fir;
t2[k]=p[j]; t1[k]=1+p[j];
}
}
for(int i=1;i<=N;i++) F[i].Sec=i;
sort(F+1,F+1+N);