线性筛法(不)完全版
求出素数,欧拉函数,莫比乌斯函数,约数个数函数,约数和函数等
原理:
1.线性筛法保证每个数只被最小素因子更新(这也是复杂度为线性O(N)的保证)
2.大部分积性函数支持线性筛法求值 (这里包括完全积性函数和不完全积性函数)
3.相关公式
令 , 表示因子个数
有
对于下面两个函数,可以处理出每个数最小质因子的相关信息(次数/幂次值)等
代码
#include<bits/stdc++.h>
typedef long long ll;
using namespace std;
const int maxn = 1000009;
int prime[maxn],phi[maxn],mul[maxn],d[maxn],e[maxn];
ll p[maxn],s[maxn];
bool vis[maxn];
void init(int mr)
{
mul[1]=phi[1]=d[1]=p[1]=1;
for(int i=2;i<=mr;i++)
{
if(!vis[i])
{
prime[++prime[0]]=i;
phi[i]=i-1;mul[i]=-1;
d[i]=2;e[i]=1;
p[i]=i+1;s[i]=1ll*i*i;
}
for(int j=1;j<=prime[0]&&i*prime[j]<=mr;j++)
{
int num=i*prime[j];
vis[num]=1;
if(i%prime[j]==0)
{
phi[num]=phi[i]*prime[j];mul[num]=0;
d[num]=d[i]/(e[i]+1)*(e[i]+2);e[num]=e[i]+1;
p[num]=(s[i]*prime[j]-1)*p[i]/(s[i]-1);s[num]=s[i]*prime[j];
break;
}
else
{
phi[num]=phi[i]*(prime[j]-1);mul[num]=-mul[i];
d[num]=d[i]*d[prime[j]];e[num]=1;
p[num]=p[i]*p[prime[j]];s[num]=1ll*prime[j]*prime[j];
}
}
}
}
int main()
{
init(maxn-9);
return 0;
}
可能要注意的细节:积性函数必有
;乘法与除法的先后顺序,防止爆int和向下取证变成0等
欧拉函数取值在[2,n-1]之间,莫比乌斯函数只会是{-1,0,1},约数个数显然不会很大,最小素因子次数也不大
关于约数和函数,显然有
根据调和级数那一套理论,约数和取值小于
最小素因子幂次+1也不会大于n*n
-----------------------------------------------------------------分割线-----------------------------------------------------------
有了她,我们可以干什么?
线性时间复杂度内求出一个区间[1,x]内所有的函数值!
还有呢?
如果顺手处理出num数组表示i的最小素因子值
质因数分解!
然后就可以做题啦!=w=