题目大意
求Σbi=1lcm(i,b)Σi=1blcm(i,b)
共50000组数据,每组数据b≤109b≤109。
解题思路
划开式子
原式=b∗Σni=1 igcd(i,b)=b∗Σi=1n igcd(i,b)
通常涉及gcdgcd的,不妨设u=gcd(i,b)u=gcd(i,b),令符合条件的ii即的统计答案。
即原式=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一定是的倍数。
=b∗Σu|bΣx|buμ(x)Σx|i,i≤bui=b∗Σu|bΣx|buμ(x)Σx|i,i≤bui
=b∗Σu|bΣx|buμ(x)∗x∗S(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=Πp为d的质因数 (1−p)Σ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表示什么。
代码
打得太丑了,于是就不贴出来了。