引理 111:
d(i⋅j)=∑x∣i∑y∣j[gcd(x,y)=1]
d (i \cdot j) = \sum_{x \mid i} \sum_{y \mid j} [\gcd (x,y) = 1]
d(i⋅j)=x∣i∑y∣j∑[gcd(x,y)=1]
证明:
一对满足要求的 x,yx, yx,y 可以构造出 xjyx \frac{j}{y}xyj 作为 ijijij 的因数。
考虑 ijijij 的质因子 ppp。
gcd(x,y)=1\gcd (x, y) = 1gcd(x,y)=1 表示都在 jjj 里取出 ppp 这个因子,只有在 jjj 的 ppp 这个因子被取空后才能拿 iii 里的 ppp 这个因子,所以不会出现类似于从 iii 里取出一个 ppp 再从 jjj 取出一个 ppp 和从 jjj 里取出两个 ppp 这两种等价方案都算了一次的情况。
引理 222:
d(i⋅j⋅k)=∑x∣i∑y∣j∑z∣k[gcd(x,y)=gcd(y,z)=gcd(x,z)=1]d (i \cdot j \cdot k) = \sum_{x \mid i} \sum_{y \mid j} \sum_{z \mid k} [gcd (x, y) = gcd (y, z) = gcd (x, z) = 1]d(i⋅j⋅k)=∑x∣i∑y∣j∑z∣k[gcd(x,y)=gcd(y,z)=gcd(x,z)=1]
证明:类似于引理 111。
化式子
∑i=1a∑j=1b∑k=1cd(i⋅j⋅k)∑i=1a∑j=1b∑k=1c∑u∣i∑v∣j∑w∣k[gcd(u,v)=gcd(u,w)=gcd(v,w)=1]∑u=1a∑v=1b∑w=1c⌊au⌋⌊bv⌋⌊cw⌋[gcd(u,v)=gcd(u,w)=gcd(v,w)=1]∑u=1a∑v=1b∑w=1c⌊au⌋⌊bv⌋⌊cw⌋∑d∣gcd(u,v)μ(d)[gcd(u,w)=gcd(v,w)=1]∑u=1a∑v=1b∑w=1c⌊au⌋⌊bv⌋⌊cw⌋∑d∣u,d∣vμ(d)[gcd(u,w)=gcd(v,w)=1]∑dμ(d)∑u=1⌊ad⌋∑v=1⌊bd⌋∑w=1c⌊aud⌋⌊bvd⌋⌊cw⌋[gcd(ud,w)=gcd(vd,w)=1]∑dμ(d)∑w=1c⌊cw⌋∑u=1⌊ad⌋⌊aud⌋[gcd(ud,w)=1]∑v=1⌊bd⌋⌊bvd⌋[gcd(vd,w)=1] \begin{aligned} \sum_{i = 1}^{a} \sum_{j = 1}^{b} \sum_{k = 1}^{c} d (i \cdot j \cdot k) \\ \sum_{i = 1}^{a} \sum_{j = 1}^{b} \sum_{k = 1}^{c} \sum_{u \mid i} \sum_{v \mid j} \sum_{w \mid k} [gcd (u, v) = gcd (u, w) = gcd (v, w ) = 1] \\ \sum_{u = 1}^{a} \sum_{v = 1}^{b} \sum_{w = 1}^{c} \lfloor \frac{a}{u} \rfloor \lfloor \frac{b}{v} \rfloor \lfloor \frac{c}{w} \rfloor [gcd (u, v) = gcd (u, w) = gcd (v, w ) = 1] \\ \sum_{u = 1}^{a} \sum_{v = 1}^{b} \sum_{w = 1}^{c} \lfloor \frac{a}{u} \rfloor \lfloor \frac{b}{v} \rfloor \lfloor \frac{c}{w} \rfloor \sum_{d \mid gcd (u, v)} \mu (d) [gcd (u, w) = gcd (v, w ) = 1] \\ \sum_{u = 1}^{a} \sum_{v = 1}^{b} \sum_{w = 1}^{c} \lfloor \frac{a}{u} \rfloor \lfloor \frac{b}{v} \rfloor \lfloor \frac{c}{w} \rfloor \sum_{d \mid u, d \mid v} \mu (d) [gcd (u, w) = gcd (v, w) = 1] \\ \sum_{d} \mu (d) \sum_{u = 1}^{\lfloor \frac{a}{d} \rfloor} \sum_{v = 1}^{\lfloor \frac{b}{d} \rfloor} \sum_{w = 1}^{c} \lfloor \frac{a}{ud} \rfloor \lfloor \frac{b}{vd} \rfloor \lfloor \frac{c}{w} \rfloor [gcd (ud, w) = gcd (vd, w) = 1] \\ \sum_{d} \mu (d) \sum_{w = 1}^{c} \lfloor \frac{c}{w} \rfloor \sum_{u = 1}^{\lfloor \frac{a}{d} \rfloor} \lfloor \frac{a}{ud} \rfloor [gcd (ud, w) = 1] \sum_{v = 1}^{\lfloor \frac{b}{d} \rfloor} \lfloor \frac{b}{vd} \rfloor [gcd (vd, w) = 1] \end{aligned} i=1∑aj=1∑bk=1∑cd(i⋅j⋅k)i=1∑aj=1∑bk=1∑cu∣i∑v∣j∑w∣k∑[gcd(u,v)=gcd(u,w)=gcd(v,w)=1]u=1∑av=1∑bw=1∑c⌊ua⌋⌊vb⌋⌊wc⌋[gcd(u,v)=gcd(u,w)=gcd(v,w)=1]u=1∑av=1∑bw=1∑c⌊ua⌋⌊vb⌋⌊wc⌋d∣gcd(u,v)∑μ(d)[gcd(u,w)=gcd(v,w)=1]u=1∑av=1∑bw=1∑c⌊ua⌋⌊vb⌋⌊wc⌋d∣u,d∣v∑μ(d)[gcd(u,w)=gcd(v,w)=1]d∑μ(d)u=1∑⌊da⌋v=1∑⌊db⌋w=1∑c⌊uda⌋⌊vdb⌋⌊wc⌋[gcd(ud,w)=gcd(vd,w)=1]d∑μ(d)w=1∑c⌊wc⌋u=1∑⌊da⌋⌊uda⌋[gcd(ud,w)=1]v=1∑⌊db⌋⌊vdb⌋[gcd(vd,w)=1]
令
qry1(d,w)=∑u=1⌊ad⌋⌊aud⌋[gcd(ud,w)=1],qry2(d,w)=∑v=1⌊bd⌋⌊bvd⌋[gcd(vd,w)=1]
qry1 (d, w) = \sum_{u = 1}^{\lfloor \frac{a}{d} \rfloor} \lfloor \frac{a}{ud} \rfloor [gcd (ud, w) = 1], qry2 (d, w) = \sum_{v = 1}^{\lfloor \frac{b}{d} \rfloor} \lfloor \frac{b}{vd} \rfloor [gcd (vd, w) = 1]
qry1(d,w)=u=1∑⌊da⌋⌊uda⌋[gcd(ud,w)=1],qry2(d,w)=v=1∑⌊db⌋⌊vdb⌋[gcd(vd,w)=1]
直接暴力求 qry1,qry2qry1, qry2qry1,qry2 的时间复杂度为 O(n52)O (n ^ {\frac{5}{2}})O(n25) (令 a,b,c,na, b, c, na,b,c,n 同阶)
证明:
当 d≤ad \leq \sqrt {a}d≤a 时,1≤w,u≤a1 \leq w, u \leq a1≤w,u≤a,这部分时间复杂度为 O(n52)O (n ^ {\frac{5}{2}})O(n25)
当 d>ad > \sqrt {a}d>a 时,1≤w≤a,1≤⌊ad⌋≤sqrt(a)1 \leq w \leq a, 1\leq \lfloor \frac{a}{d} \rfloor \leq sqrt (a)1≤w≤a,1≤⌊da⌋≤sqrt(a),这部分时间复杂度为 O(n52)O (n ^ {\frac{5}{2}})O(n25)。
又是推错的一些没用的式子
∑i=1a∑j=1b∑k=1cd(ijk)∑i=1a∑j=1b∑k=1c∑d∣ijk1∑i=1a∑j=1b∑d∑k=1,d∣ijkc1∑i=1a∑j=1b∑d⌊cdgcd(d,ij)⌋∑i=1a∑j=1b∑d⌊c⋅gcd(d,ij)d⌋∑i=1a∑g∑d⌊c⋅gd⌋∑j=1b[gcd(d,ij)=g]∑i=1a∑g∑d⌊c⋅gd⌋∑j=1b[gcd(dg,ijg)=1]∑g∑d⌊c⋅gd⌋∑t∑i=1,i∣ta[gcd(dg,tg)=1]
\begin{aligned}
\sum_{i = 1}^{a}\sum_{j = 1}^{b} \sum_{k = 1}^{c} d (ijk)
\\
\sum_{i = 1}^{a}\sum_{j = 1}^{b} \sum_{k = 1}^{c} \sum_{d \mid ijk} 1
\\
\sum_{i = 1}^{a}\sum_{j = 1}^{b} \sum_{d} \sum_{k = 1, d \mid ijk}^{c} 1
\\
\sum_{i = 1}^{a}\sum_{j = 1}^{b} \sum_{d} \lfloor \frac{c}{\frac{d}{\gcd (d, ij)}} \rfloor
\\
\sum_{i = 1}^{a}\sum_{j = 1}^{b} \sum_{d} \lfloor \frac{c \cdot \gcd (d, ij)}{d} \rfloor
\\
\sum_{i = 1}^{a} \sum_{g} \sum_{d} \lfloor \frac{c \cdot g}{d} \rfloor \sum_{j = 1}^{b}[\gcd (d, ij) = g]
\\
\sum_{i = 1}^{a} \sum_{g} \sum_{d} \lfloor \frac{c \cdot g}{d} \rfloor \sum_{j = 1}^{b}[\gcd (\frac{d}{g}, \frac{ij}{g}) = 1]
\\
\sum_{g} \sum_{d} \lfloor \frac{c \cdot g}{d} \rfloor \sum_{t} \sum_{i = 1, i \mid t}^{a} [\gcd (\frac{d}{g}, \frac{t}{g}) = 1]
\end{aligned}
i=1∑aj=1∑bk=1∑cd(ijk)i=1∑aj=1∑bk=1∑cd∣ijk∑1i=1∑aj=1∑bd∑k=1,d∣ijk∑c1i=1∑aj=1∑bd∑⌊gcd(d,ij)dc⌋i=1∑aj=1∑bd∑⌊dc⋅gcd(d,ij)⌋i=1∑ag∑d∑⌊dc⋅g⌋j=1∑b[gcd(d,ij)=g]i=1∑ag∑d∑⌊dc⋅g⌋j=1∑b[gcd(gd,gij)=1]g∑d∑⌊dc⋅g⌋t∑i=1,i∣t∑a[gcd(gd,gt)=1]