线性筛及其扩展——学习笔记

最初的应用：线性筛素数。
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)。
只要在线性筛素数代码上加点东西就行了。
甚至是一坨奇奇怪怪的东西，只要是积性的一般都能求。需要多维护几个辅助数组。
求$\phi$:

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]];
        }
    }
}

求$\mu$:

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);