logo

Matt Blog

春节

数学学习笔记

2024-07-10 Views2295字13 min read

狄利克雷卷积 Dirichlet

定义

其本质是一种运算,可以用 * 表示。

h(x)=f(x)g(x)=dxf(d)g(xd)=ab=xf(a)g(b)h(x)=f(x)*g(x)=\sum_{d|x}f(d)g(\frac{x}{d})=\sum_{ab=x}f(a)g(b)

积性函数 f(x)f(x) 满足 (a,b)=1,f(ab)=f(a)f(b)\forall (a,b)=1,f(ab)=f(a)f(b)

完全积性函数 f(x)f(x) 满足 f(ab)=f(a)f(b)f(ab)=f(a)f(b)

可以使用线性筛算出积性函数。

性质

交换律:fg=gff*g=g*f

结合律:(fg)h=f(gh)(f*g)*h=f*(g*h)

分配律(f+g)h=fh+gh(f+g)*h=f*h+g*h

ffgg 是积性函数,则 fgf*g 也是积性函数。

ff 是积性函数,则 ff 的逆元也是积性函数。

逆元

单位函数 ε\varepsilonfε=ff*\varepsilon =f

gg 满足 fg=εf*g=\varepsilon,则称 g(x)g(x)f(x)f(x) 的逆元。逆元唯一。

ps:目前还不太懂这是什么,但好像可以理解为 ε(n)={1n=10n=1\varepsilon(n)=\begin{cases}1&n=1\\0&n\not=1\end{cases}

g(x)=ε(x)dx,d=1f(d)f(xd)1g(x)= \dfrac{\varepsilon(x)-\sum_{d|x,d\not=1} f(d)f(\dfrac{x}{d})}{1}

莫比乌斯函数 Möbius

定义式

n=p1a1p2a2...pkakμ(n)={1n=10ik,ai>1(1)kelsen=p_1^{a_1}p_2^{a_2}...p_k^{a_k}\\ \mu(n)=\begin{cases} 1 &n=1 \\ 0 & \exist i \le k,a_i>1 \\ (-1)^k& else \end{cases}

性质

μ(x)\mu(x) 是积性函数。

dnμ(d)={1n=10n=1\sum_{d|n}\mu(d)=\begin{cases}1&n=1\\0&n\not=1\end{cases}

dnμ(d)=ε(n),μ1=ε\sum_{d|n}\mu(d)=\varepsilon(n),\mu*1=\varepsilon。在狄利克雷卷积中,莫比乌斯函数是常数函数 11 的逆元。

证明

n=p1a1p2a2...pkak,n=p1p2...pkμ(k×a2)=0dnμ(d)=dnμ(d)n=p_1^{a_1}p_2^{a_2}...p_k^{a_k},n'=p_1p_2...p_k\\ \because \mu(k \times a^2)=0 \\ \therefore \sum_{d|n}\mu(d)=\sum_{d|n'}\mu(d)

则现在就是尝试在 kk 个互异质数中选择一个集合,若集合大小为奇数,答案 -1,否则答案 +1。k>1k>1 时,dnμ(d)=i=0k(ki)×(1)i\sum_{d|n'}\mu(d)=\sum_{i=0}^k\binom{k}{i}\times(-1)^i。二项式定理是 (a+b)n=i=0n(ni)aibni(a+b)^n=\sum_{i=0}^n\binom{n}{i}a^ib^{n-i},将 a=1,b=1a=-1,b=1 带入则是 (1+(1))n=i=0n(ni)×(1)i=0(1+(-1))^n=\sum_{i=0}^n\binom{n}{i}\times(-1)^i=0。故 k>1k>1 时,dnμ(d)=0\sum_{d|n}\mu(d)=0。剩下情况不难自证。

ps:以上证明只是本人在初学时的结论,不保证正确,等待未来的我修改或确认。

反演结论

[gcd(i,j)=1]=dgcd(i,j)μ(d)[\gcd(i,j)=1]=\sum_{d|\gcd(i,j)}\mu(d)

例题

「HAOI 2011」Problem b

i=xnj=ym[gcd(i,j)=k],(1T,x,y,n,m,k105)\displaystyle{\sum_{i=x}^n\sum_{j=y}^m[\gcd(i,j)=k]} ,(1 \le T,x,y,n,m,k \le 10^5)

(1)f(a,b,k)=i=1aj=1b[gcd(i,j)=k]f(a,b,k)=\sum_{i=1}^a\sum_{j=1}^b[\gcd(i,j)=k] \tag{1}

(2)i=xnj=ym[gcd(i,j)=k]=f(n,m,k)f(x1,m,k)f(n,y1,k)+f(x1,y1,k)\sum_{i=x}^n\sum_{j=y}^m[\gcd(i,j)=k]=f(n,m,k)-f(x-1,m,k)-f(n,y-1,k)+f(x-1,y-1,k) \tag{2}

(3)f(a,b,k)=i=1akj=1bk[gcd(i,j)=1]f(a,b,k)=\sum_{i=1}^{\left \lfloor \frac{a}{k} \right \rfloor}\sum_{j=1}^{\left \lfloor \frac{b}{k} \right \rfloor}[\gcd(i,j)=1] \tag{3}

(4)=i=1akj=1bkdgcd(i,j)μ(d)=\sum_{i=1}^{\left \lfloor \frac{a}{k} \right \rfloor}\sum_{j=1}^{\left \lfloor \frac{b}{k} \right \rfloor}\sum_{d|\gcd(i,j)}\mu(d) \tag{4}

(5)=dμ(d)i=1ak[di]j=1bk[dj]=\sum_{d}\mu(d)\sum_{i=1}^{\left \lfloor \frac{a}{k} \right \rfloor}[d|i]\sum_{j=1}^{\left \lfloor \frac{b}{k} \right \rfloor}[d|j] \tag{5}

(6)=d=1min(ak,bk)μ(d)akdbkd=\sum_{d=1}^{\min(\left \lfloor \frac{a}{k} \right \rfloor,\left \lfloor \frac{b}{k} \right \rfloor)}\mu(d)\left \lfloor \frac{\left \lfloor \frac{a}{k}\right \rfloor}d \right \rfloor\left \lfloor \frac{\left \lfloor \frac{b}{k}\right \rfloor}d \right \rfloor \tag{6}

(1)(1)(2)(2) 的过程是二维差分,(3)(3)(4)(4) 的过程就是反演结论的使用,(5)(5)(6)(6) 是因为 i=1ak[di]\displaystyle\sum_{i=1}^{\left \lfloor \frac{a}{k} \right \rfloor}[d|i] 就是 [1,ak][1,\left \lfloor \frac{a}{k} \right \rfloor] 中是 dd 的倍数的个数,另一个同理。

使用 O(V)O(V) 线性筛预处理出 μ(i)\mu(i)i=1nμ(i)\sum_{i=1}^n\mu(i),使用整除分块在 O(V)O(\sqrt V) 算出 (6)(6)。总时间复杂度为 O(V+V)O(V+\sum\sqrt V)

定义

由欧拉定理可知,(a,m)=1(a,m)=1aφ(m)1(modm)a^{\varphi(m)} \equiv 1(\bmod m)。故存在最小正整数满足 an1(modm)a^n \equiv 1(\bmod m),称 nnaamm 的阶,记作 ordm(a)\operatorname{ord}_m(a)

性质

a,a2...aordm(a)a,a^2...a^{\operatorname{ord}_m(a)}mm 两两不同余。

apaq(modm)a^p \equiv a^q(\bmod m),则 pq(modordm(a))p\equiv q(\bmod \operatorname{ord}_m(a))

(a,m)=(b,m)=1ordm(ab)=ordm(a)ordm(b)gcd(ordm(a),ordm(b))=1(a,m)=(b,m)=1 \\ \operatorname{ord}_m(ab)=\operatorname{ord}_m(a)\operatorname{ord}_m(b) \Leftrightarrow gcd(\operatorname{ord}_m(a),\operatorname{ord}_m(b))=1

原根 primitive-root

定义

(g,m)=1(g,m)=1,且 ordm(g)=φ(m)\operatorname{ord}_m(g)=\varphi(m),则称 gg 为模 mm 的原根。

性质

对于 m3,(g,m)=1m \ge 3,(g,m)=1,则 gg 是模 mm 的原根 \Leftrightarrow 对于 φ(m)\varphi(m) 的每个质因子 pp,都有 gφ(m)p1(modm)g^{\frac{\varphi(m)}{p}} \not\equiv 1(\bmod m)

mm 有原根,则原根个数为 φ(φ(m))\varphi(\varphi(m))

原根存在定理:一个数 mm 存在原根当且仅当 m=2,4,pa,2pam=2,4,p^a,2p^app 为奇质数。

aa 的最小原根大约在 a4\sqrt[4]a,long long 内的都可以暴力枚举。

ggnn 的最小原根,则 nn 的原根集合 S={gs(s,φ(n))=1}S=\{g^s|(s,\varphi(n))=1\}

部分证明

由第 1 条得,p is prime and pφ(n),gφ(n)p=1(modn)\forall \text{p is prime and }p| \varphi(n),g^{\frac{\varphi(n)}{p}}\not=1(\bmod n)

gsφ(n)p=(gφ(n))spg^{s^{\frac{\varphi(n)}{p}}}=(g^{\varphi(n)})^\frac{s}{p},若 gcd(s,φ(n))>1\gcd(s,\varphi(n))>1,则有一个 pp 使 sp\frac{s}{p} 为正整数,由欧拉定理得 gφ(n)1(modn)g^{\varphi(n)}\equiv1(\bmod n),矛盾,故 (s,φ(n))=1(s,\varphi(n))=1,证明其必要性。

Burnside&Pólya 定理

(1)[X/G]=1GgGXg[X/G]=\dfrac{1}{|G|}\sum_{g\in G}|X^g| \tag{1}

即本质不同的等价类个数是不动点个数的平均值。Xg|X^g| 表示 XX 在置换 gg 下的不动点个数。

(2)[X/G]=1GgGmc(g)[X/G]=\dfrac1{|G|}\sum_{g\in G}m^{c(g)} \tag{2}

c(g)c(g) 表示置换 gg 形成了多少个环。

因为一个点在置换 gg 上成为不定点当且仅当这个点所在的环是同一种颜色,所以对每个环单独染色即可。

ps:(1)(1) 中的 GG 是定义在染色方案上的,即有 nmn^m 种染色方案,(2)(2) 中的 GG 是定义在 nn 个点的轮换上,即有 nn 种旋转方式。

EOF