数论2

数论-进阶#

扩展欧拉定理#

bφ(m) 则有 ababmodφ(m)+φ(m)(modm)

上帝与集合的正确用法#

思路

题目要你求:

2222modp

的值,多组询问

很显然,设 2222=x 则由于扩展欧拉定理:

2x2xmodφ(p)+φ(p)(modp)

然后出现了 2xmodφ(p) 这个式子,我们就可以不断递归求解了

狄利克雷卷积#

定义#

数论函数 f(n)g(n) 的狄利克雷卷积定义为:

h(n)=d|nf(d)g(nd)

记作:h=fg

定理#

两个积性函数的狄利克雷卷积也是积性函数

两个积性函数对应位置相乘所得到的函数也是积性函数

交换律结合律#

fg=gf

f(gh)=(fg)h

其中 f,g,h 为积性函数

I 函数#

I(n)=[x=1]

I 是单位元, 数论函数 f 都有 fI=If=f

狄利克雷逆#

fg=I,则 f,g 互为彼此的狄利克雷逆

常见积性函数及性质#

1 函数:1(n)=1

Id 函数:Id(n)=n

I 函数:I(n)=[n=1]

Id=φ1Idμ=φ

莫比乌斯反演#

莫比乌斯函数#

对于 μ(n)

n=1μ(n)=1

n=p1p2pk(其中 pi 为互异素数)μ(n)=(1)k

否则 μ(n)=0

由于 1μ=I,所以 1,μ 互为彼此的狄利克雷逆

定义#

设数论函数 f,称 F=f1 为莫比乌斯变换

设数论函数 F,称 f=Fμ 为莫比乌斯反演

优秀的性质#

对于约数个数函数 d(x) 有如下性质:

d(ij)=x|iy|j[gcd(x,y)=1]

d|nμ(d)d=φ(n)n

d|gcd(i,j)μ(d)=[gcd(i,j)=1]

Idμ=φ

P1447 [NOI2010] 能量采集#

此题即求:

2×i=1nj=1mgcd(i,j)n×m

后续我们假定 nm

显然我们需要求的就是

i=1nj=1mgcd(i,j)

按照常规套路,枚举 gcd(i,j) 等于多少

k=1nki=1nj=1m[gcd(i,j)=k]

常规莫比乌斯反演

k=1nki=1nkj=1mk[gcd(i,j)=1]

k=1nki=1nkj=1mkd|gcd(i,j)μ(d)

k=1nkd=1nkμ(d)nkdmkd

这样再使用整出分块足以通过此题

然而我们仍可以优化

我们令 T=kd

k=1nkT=1n[k|T]μ(Tk)nTmT

优先枚举 T

T=1nk|Tkμ(Tk)nTmT

由于有关狄利克雷卷积的常识 Idμ=φ,原式可以化简为

T=1nφ(T)nTmT

那么此题就可以在 O(min(n,m)) 的时间内解决

作者:Into_qwq

出处:https://www.cnblogs.com/into-qwq/p/16491444.html

版权:本作品采用「qwq」许可协议进行许可。

posted @   Into_qwq  阅读(33)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 使用C#创建一个MCP客户端
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· 按钮权限的设计及实现
more_horiz
keyboard_arrow_up dark_mode palette
选择主题
menu
点击右上角即可分享
微信分享提示