【题解】最小公倍数计数

给定 [ a , b ] [a,b] [a,b] ,求 lcm ( i , j ) \text{lcm}(i,j) lcm(i,j) [ a , b ] [a,b] [a,b] 中的有序点对 ( i , j ) (i,j) (i,j) 的数量。

a , b ≤ 1 0 11 a,b\leq 10^{11} a,b1011

sol :

考察 lcm ( i , j ) = n \text{lcm}(i,j)=n lcm(i,j)=n 时无序点对 ( i , j ) (i,j) (i,j) 的数量。

233

∑ i ∣ n ∑ j ∣ n [ i j ( i , j ) = n ] = ∑ k ∣ n ∑ i ∣ n k ∑ j ∣ n k [ ( i , j ) = 1 , i j = n k ] = ∑ k ∣ n ∑ i ∣ n k [ ( i , n i k ) = 1 ] = ∑ k ∣ n ∑ i ∣ n k ∑ j ∣ ( i , n i k ) μ ( j ) = ∑ j ∣ n μ ( j ) ∑ j ∣ i , i ∣ n j ∑ k ∣ n i j 1 = ∑ j 2 ∣ n μ ( j ) ∑ i ∣ n j 2 ∑ k ∣ n i j 2 1 = ∑ j 2 ∣ n μ ( j ) ∑ i ∣ n j 2 d ( n i j 2 ) \sum_{i|n}\sum_{j|n}[\frac{ij}{(i,j)}=n] \\ =\sum_{k|n}\sum_{i|\frac{n}{k}}\sum_{j|\frac{n}{k}}[(i,j)=1,ij=\frac{n}{k}] \\ =\sum_{k|n}\sum_{i|\frac{n}{k}}[(i,\frac{n}{ik})=1] \\ =\sum_{k|n}\sum_{i|\frac{n}{k}}\sum_{j|(i,\frac{n}{ik})}\mu(j) \\ =\sum_{j|n}\mu(j)\sum_{j|i,i|\frac{n}{j}}\sum_{k|\frac{n}{ij}}1 \\ =\sum_{j^2|n}\mu(j)\sum_{i|\frac{n}{j^2}}\sum_{k|\frac{n}{ij^2}}1 \\ =\sum_{j^2|n}\mu(j)\sum_{i|\frac{n}{j^2}}d(\frac{n}{ij^2}) injn[(i,j)ij=n]=kniknjkn[(i,j)=1,ij=kn]=knikn[(i,ikn)=1]=kniknj(i,ikn)μ(j)=jnμ(j)ji,ijnkijn1=j2nμ(j)ij2nkij2n1=j2nμ(j)ij2nd(ij2n)

∑ k = 1 n ∑ j 2 ∣ k μ ( j ) ∑ i ∣ k j 2 d ( k i j 2 ) = ∑ j 2 ≤ n μ ( j ) ∑ k = 1 [ n j 2 ] ∑ i ∣ k d ( i ) = ∑ j 2 ≤ n μ ( j ) ∑ k = 1 [ n j 2 ] d ( k ) [ n j 2 k ] \sum_{k=1}^n\sum_{j^2|k}\mu(j)\sum_{i|\frac{k}{j^2}}d(\frac{k}{ij^2}) \\ =\sum_{j^2\leq n}\mu(j)\sum_{k=1}^{[\frac{n}{j^2}]}\sum_{i|k}d(i) \\ =\sum_{j^2\leq n}\mu(j)\sum_{k=1}^{[\frac{n}{j^2}]}d(k)[\frac{n}{j^2k}] \\ k=1nj2kμ(j)ij2kd(ij2k)=j2nμ(j)k=1[j2n]ikd(i)=j2nμ(j)k=1[j2n]d(k)[j2kn]

事后来看,其实当时我做这道题的时候并没有抓住它的本质。

好在我们展开一下。

F ( n ) = ∑ i = 1 n d ( i ) [ n i ] = ∑ a ∑ b ∑ c [ a b c ≤ n ] F(n)=\sum_{i=1}^nd(i)[\frac{n}{i}] \\ =\sum_a\sum_{b}\sum_{c}[abc\leq n] F(n)=i=1nd(i)[in]=abc[abcn]

于是可以转化为经典的三元计数问题。

∑ j 2 ≤ n μ ( j ) F ( [ n j 2 ] ) \sum_{j^2\leq n}\mu(j)F([\frac{n}{j^2}]) j2nμ(j)F([j2n])


__EOF__

本文作者仰望星空的蚂蚁
本文链接https://www.cnblogs.com/cqbzly/p/17530183.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   仰望星空的蚂蚁  阅读(6)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」
点击右上角即可分享
微信分享提示