积性函数、线性筛、莫比乌斯反演和一堆乱七八糟的题目

最新推荐文章于 2022-07-11 16:02:04 发布

EIKY

最新推荐文章于 2022-07-11 16:02:04 发布

阅读量633

点赞数 1

分类专栏：整理文章标签： ACM c++

整理专栏收录该内容

3 篇文章 0 订阅

订阅专栏

某菊苣的总结，比较全面= =。

转自：http://jcvb.is-programmer.com/posts/41846.html

·积性函数
定义在正整数集上的函数 $f(n)$ （称为算术函数），若 $\text g \text c \text d(a,b)=1$ 时有 $f(a)f(b)=f(ab)$ ，则 $f(x)$ 称为积性函数。
一个显然的性质：（非恒等于零的）积性函数 $f(n)$ 必然满足 $f(1)=1$ 。
定义逐点加法 $(f+g)(x)=f(x)+g(x),逐点乘法(f\cdot g)(x)=f(x)g(x)$ 。
一个比较显然的性质：若 $f,g$ 均为积性函数，则 $f\cdot g$ 也是积性函数。
积性函数的求值： $n=\prod p_i^{\alpha_i}$ ，则 $f(n)=\prod f(p_i^{\alpha_i})$ ，所以只要解决 $n=p^\alpha$ 时 $f(n)$ 的值即可。

例如：
恒为1的常函数 $1(n)=1$ ，
恒等函数 $\text{id}(n)=n$ ，
单位函数 $\epsilon(n)=[n==1]$ ，（这三个都是显然为积性）
欧拉函数 $\varphi(n)$ （只要证两个集合相等就能证明积性）
莫比乌斯函数 $\mu(n)$ （由定义也是显然的）

·Dirichlet卷积
对两个算术函数 $f,g$ ，定义其Dirichlet卷积为新函数 $f*g$ ，满足 $(f*g)(n)=\sum_{d|n}{f(d)g(n/d)}$ 。
一些性质：
交换律 $f*g=g*f$ ,
结合律 $(f*g)*h=f*(g*h)$ ,
单位元 $f*\epsilon=f$ ,
对逐点加法的分配律 $f*(g+h)=f*g+f*h$

重要性质：若 $f,g$ 均为积性函数，则 $f*g$ 也是积性函数。（展开式子即可证明）
$n$ 的约数个数 $d(n)$ 可以写成 $d(n)=(1*1)(n)$ ；约数和 $\sigma(n)$ 可以写成 $\sigma(n)=(1*\text{id})(n)$ ，由上面的性质可知这两个函数均是积性函数。

重要性质： $\sum_{d|n}{\mu(d)}=[n==1]$ ，即 $1*\mu=\epsilon$ 。（可用二项式定理证明）
重要性质： $\sum_{d|n}{\varphi(d)}=n$ ，即 $1*\varphi=\text{id}$ 。（n是质数时显然成立，再由积性得证）

·O(nlgn)预处理Dirichlet卷积
若已知 $f(i),g(i),i=1,2,\dots,n$ 的值，则可以在O(nlgn)时间内计算出 $(f*g)(i),i=1,2,\dots,n$ 。

intf[MAXN],g[MAXN],h[MAXN]={0};
voidcalc(intn)
{   
    for(inti=1;i*i<=n;i++)
        for(intj=i;i*j<=n;j++)
            if(j==i)h[i*j]+=f[i]*g[i];
            elseh[i*j]+=f[i]*g[j]+f[j]*g[i];
}

·线性筛O(n)预处理

intpr[MAXN],bo[MAXN]={0},tot=0;
voidsieve(intn){
    for(inti=2;i<=n;i++){
        if(!bo[i])pr[++tot]=i;
        for(intj=1;j<=tot && pr[j]*i<=n;j++){
            bo[i*pr[j]]=1;
            if(i%pr[j]==0)break;
        }
    }
}

可以证明每个合数只会被其最小质因子筛去，因此复杂度为线性。
如果能够充分挖掘 $f(n)$ 的性质，则可以用线性筛在O(n)时间内预处理 $f(i),i=1,2,\dots,n$ 。
关键在于对下面三种情况分别进行处理：
1. $x$ 是质数，求 $f(x)$
2. $p$ 是质数， $p\nmid x$ ，求 $f(px)$
3. $p$ 是质数， $p\mid x$ ，求 $f(px)$
第1种情况往往比较简单。如果能证明 $f$ 是积性函数，则容易知道第2种情况是 $f(px)=f(p)f(x)$ ，剩下只要解决第3种情况的递推就够了。
有时虽然 $f$ 不是积性函数，但也能够通过分析其性质分别解决这三种情况（后面会举一些例子）。

·一堆题目
以下均默认 $a\leq b$
1101: [POI2007]Zap
令 $a'=\left \lfloor a/d \right \rfloor,b'=\left \lfloor b/d \right \rfloor$
$\sum_{i=1}^a\sum_{j=1}^b[\text{gcd}(i,j)==d]$
$=\sum_{i=1}^{a'}\sum_{j=1}^{b'}[\text{gcd}(i,j)==1]$
用莫比乌斯函数的性质把求和的式子换掉，
$=\sum_{i=1}^{a'}\sum_{j=1}^{b'}\sum_{d|\text{gcd}(i,j)}\mu(d)$
其中 $d|\text{gcd}(i,j) \Leftrightarrow d|i \,\text{and} \,d|j$ ，更换求和指标，
$=\sum_{d=1}^{a'}\mu(d)\left \lfloor a'/d \right \rfloor \left \lfloor b'/d \right \rfloor$
容易知道 $\left\lfloor a'/d\right \rfloor$ 单调不上升，且最多有 $2\sqrt{a'}$ 种不同的取值。所以按取值分成 $O(\sqrt{n})$ 个段分别处理，一个连续段内的和可以用预处理出的莫比乌斯函数前缀和求出。

if(a>b)swap(a,b);
intnex;
ll ans=0;
for(inti=1;i<=a;i=nex+1){
    nex=min(a/(a/i),b/(b/i));
    ans+=1ll*(s[nex]-s[i-1])*(a/i)*(b/i);
}

整道题目预处理O(n)，每次询问O(sqrt(n))。

2005: [Noi2010]能量采集
$\sum_{i=1}^a\sum_{j=1}^b\text{gcd}(i,j)$
$=\sum_{d=1}^a \sum_{i=1}^{\left \lfloor a/d \right \rfloor}\sum_{j=1}^{\left \lfloor b/d \right \rfloor}d[\text{gcd}(i,j)==1]$
套用上题，
$=\sum_{d=1}^a \sum_{d'=1}^{\left \lfloor a/d \right \rfloor}d\mu(d')\left \lfloor a/dd' \right \rfloor \left \lfloor b/dd' \right \rfloor$
关键一步，作代换 $D=dd'$ ，把求和指标d,d'改为D,d，（容易证明这是等价的）
$=\sum_{D=1}^a \sum_{d|D}d\mu(D/d)\left \lfloor a/D\right \rfloor \left \lfloor b/D\right \rfloor$
$=\sum_{D=1}^a \left \lfloor a/D\right \rfloor \left \lfloor b/D\right \rfloor (\text{id}*\mu)(D)$
由 $\varphi*1=\text{id}$ ，有 $\text{id}*\mu=\varphi$ ，
$=\sum_{D=1}^a \left \lfloor a/D\right \rfloor \left \lfloor b/D\right \rfloor \varphi(D)$
最后处理方式同上题。（不过原题只有一个询问，不分段也不要紧）

2693: jzptab
$\sum_{i=1}^a\sum_{j=1}^b ij/\text{gcd}(i,j)$
$=\sum_{d=1}^a\sum_{i=1}^{\left \lfloor a/d\right \rfloor}\sum_{j=1}^{\left \lfloor b/d\right \rfloor}ijd[\text{gcd}(i,j)==1]$
类似的方法将方括号换掉，化简为
$=\sum_{d=1}^a\sum_{d'=1}^{\left \lfloor a/d\right \rfloor}\sum_{i=1}^{\left \lfloor a/dd' \right \rfloor}\sum_{j=1}^{\left \lfloor b/dd' \right \rfloor}d\mu(d')d'^2ij$
由等差数列求和，
$=\sum_{d=1}^a\sum_{d'=1}^{\left \lfloor a/d \right \rfloor}d\mu(d')d'^2\text{sum}(\left \lfloor a/dd' \right \rfloor,\left \lfloor b/dd' \right \rfloor)$
其中 $\text{sum}(a,b)=a(a+1)b(b+1)/4$ ，
像上题一样代换求和指标，
$=\sum_{D=1}^a D \text{sum}(\left \lfloor a/D \right \rfloor,\left \lfloor b/D \right \rfloor)\sum_{d'|D}d'\mu(d')$
$=\sum_{D=1}^a D\cdot((\text{id}\cdot \mu)*1)(D) \text{sum}(\left \lfloor a/D \right \rfloor,\left \lfloor b/D \right \rfloor)$
处理方式同上。

3309: DZY Loves Math
类似的方法化简得
$\sum_{D=1}^a (\mu * f)(D) \left \lfloor a/D \right \rfloor \left \lfloor b/D \right \rfloor$
记 $g=\mu * f$ ， $g$ 并不是积性函数。
观察 $g(n)$ 的取值可以发现，
$g(1)=0$ ；
否则，当 $n=(\prod_{i=1}^k p_i)^\alpha=\prod_{i=1}^k p_i^\alpha$ 时， $g(n)=(-1)^{k+1}$ ；
否则， $n$ 中存在两个质因子的次数不相同， $g(n)=0$ 。
上述性质容易通过展开 $g(\prod p_i^{\alpha_i})$ 得到证明。
这样，在线性筛的同时记录 $n$ 的最小质因子 $p_1$ 的次数 $\alpha_1$ ，以及 $\frac{n}{p_1^{\alpha_1}}$ 即可。

EIKY

关注

1
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
积性函数、线性筛、莫比乌斯反演和一堆乱七八糟的题目

某菊苣的而总结，比较全面= =。转自：http://jcvb.is-programmer.com/posts/41846.html·积性函数定义在正整数集上的函数（称为算术函数），若时有，则称为积性函数。一个显然的性质：（非恒等于零的）积性函数必然满足。定义逐点加法。一个比较显然的性质：若均为积性函数，则也是积性函数。积性函数的求值：，则，所以只要解决时的值即可
复制链接

扫一扫

专栏目录