当我有些理解线性筛

最近在学习数论，重新看了莫比乌斯反演，欧拉函数等。。。

发现线性筛有很大的用处。。。。

1.筛素数

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5+5;

bool not_prime[N];
int prime[N];
int tot;
void get_prime(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!not_prime[i])prime[++tot]=i;
        for(int j=1;j<=tot&&prime[j]*i<=n;j++)
        {
            not_prime[prime[j]*i]=1;
            if(i%prime[j]==0)break;
        }
    }
}
int main()
{
    get_prime(1000);
}

时间复杂度O（n）;

其中i%prime[j]==0，break;

这一句很重要，因为这样才使我们的复杂度降低到了o（n）；

根据唯一分解定理：

我们i%prime[j]==0退出即保证了每个合数只能被他的最小质因数筛掉，即每个数访问了一次。

设

x1=i*prime[j];

x2=i*prime[j+1];

又因为i%prime[j]==0;

所以prime[j]为i的最小质因子。

所以prime[j]也为x2的最小质因子。所以prime[j]后边的就不用继续筛了。

 

如果明白了这一点就能理解，为什么在线性筛的过程中能够求出来莫比乌斯函数的值。

#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int N = 1e5+5;
bool not_prime[N];
int prime[N];
int mo[N];
int tot;

void init(int n)
{
    mo[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!not_prime[i])
        {
            prime[++tot]=i;
            mo[i]=-1;
        }
        for(int j=1;j<=tot&&prime[j]*i<=n;j++)
        {
            not_prime[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                mo[i*prime[j]]=0;
                break;
            }
            else mo[i*prime[j]]=-mo[i];
        }
    }
}

int main()
{
    //freopen("D:\ou2.txt","w",stdout);
    init(1000);
}

首先如果一个数是素数，那么mo[i]=-1；

根据莫比乌斯函数的定义：

当i%prime[j]==0的时候，i*prime[j]的出现了次方为两次的质因数（k*prime[j]^2）。所以当前值为0；

否则的话说明i与prime[j]互质，由于莫比乌斯函数是积性函数，所以mo[i*prime[j]]=mo[i]*mo[prime[j]]。