51Nod 1363 最小公倍数之和

本文介绍了一种解决数学问题Σlcm(i,b)的方法,通过对表达式的逐步转换,将其简化为可计算的形式。利用了gcd、μ函数的性质以及S(n)的定义,最终通过枚举和递归的方式高效地解决了大规模数据输入的问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目大意

Σbi=1lcm(i,b)Σi=1blcm(i,b)
共50000组数据,每组数据b109b≤109

解题思路

划开式子
原式=bΣni=1 igcd(i,b)=b∗Σi=1n igcd(i,b)
通常涉及gcdgcd的,不妨设u=gcd(i,b)u=gcd(i,b),令符合条件的iigcd(i,b)=u的统计答案。
即原式=bΣbu=1 1uΣbi=1 i[gcd(i,b)=u]=b∗Σu=1b 1uΣi=1b i[gcd(i,b)=u]
=bΣu|b 1uΣbui=1iu[gcd(i,bu)=1]=b∗Σu|b 1uΣi=1buiu[gcd(i,bu)=1]
=bΣu|bΣbui=1i[gcd(i,bu)=1]=b∗Σu|bΣi=1bui[gcd(i,bu)=1]
通常涉及到[gcd(i,b)=1][gcd(i,b)=1]的,想到Σd|bμ(d)=[b=1]Σd|bμ(d)=[b=1]
所以
=bΣu|bΣbui=1iΣx|gcd(i,bu)μ(x)=b∗Σu|bΣi=1buiΣx|gcd(i,bu)μ(x)
现在是一个大难点。很容易将人卡死。式子gcd(i,bu)gcd(i,bu)好似是变量。
将它看成定的,再去确定x,ix,i等值的范围。
gcd(i,bu)=xgcd(i,bu)=x,则满足条件的ii一定是x的倍数。
=bΣu|bΣx|buμ(x)Σx|i,ibui=b∗Σu|bΣx|buμ(x)Σx|i,i≤bui
=bΣu|bΣx|buμ(x)xS(bux)=b∗Σu|bΣx|buμ(x)∗x∗S(bux)
其中S(n)S(n)表示1+2+...+n1+2+...+n
化动为定,枚举T=uxT=ux。可以发现,T出现的次数等于T的因数个数。进一步研究发现。。。。
则原式=ΣT|bS(bT)Σd|Tμ(d)d=ΣT|bS(bT)Σd|Tμ(d)∗d
对于50000个询问的数b,每一个用O()O(很小)就可以求出所有T,T|b的对答案的贡献。
直接递归就好了。
Σd|Tμ(d)d=Πpd     (1p)Σd|Tμ(d)∗d=Πp为d的质因数     (1−p)

总结

最重要的几个窍门吧

①化动为定,如果式子中某个量是变的,那么通常将它看成定的,再去确定相关的原来枚举的量的范围。
②划开[gcd(x,y)=a][gcd(x,y)=a]的式子,变成[gcd(xa,ya)=1][gcd(xa,ya)=1],再用Σd|gcd(xa,ya)μ(d)=[gcd(xa,ya)=1]Σd|gcd(xa,ya)μ(d)=[gcd(xa,ya)=1]
③如果实在没什么头绪,试推一下反演,设f,g表示什么。

代码

打得太丑了,于是就不贴出来了。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值