随笔分类 - 莫比乌斯反演
摘要：胡策题，迫使懒得点技能树选手学算法。。。 至少杜教和min_25会了 不想写小结啊啊啊啊 那么对于这个题先上个莫反套路 然后用杜教乘个I把f*u给消掉，发现h就是f，这个可以用类似min_25的思想搞。。。
阅读全文
摘要：了然自闭。话说sdoi2018两道反演真的好吗 出门右拐%题解谢谢 终于见到一题不是gcd的反演了。。。 学了下大质数分解。。。注意pollard_rho随机的时候要选择一个常数c并且不再边，不然999999726000014413这种东西会让你怀疑人生 这种题暴露我的非酋属性了啊bzoj垫底loj
阅读全文
摘要：最终的柿子：sigema(min)T K(T)*S(T) K(T)=sigema(n/T)i sigema(m/T)j i*j S(T)=sigema(d|T)u(d)*d^d*T 收获：写柿子的时候不能省掉范围。。。要分块找向下整除
阅读全文
摘要：懵逼乌斯反演裸题 画完柿子只要预处理S(i)=sigema(d|i)d u(i/d)*d^k 线性筛它一手 i%prime[j]!=0的话，首先i的因子也是i*prime[j]的因子，但是由于多了一个质因子所以正负性反过来。其次多了这个prime[j]，i的所有因子都可以乘这个数，同样是它的因子，值
阅读全文
摘要：这题其实现在看来没那么难。 题解不写了，popoQQQ的论文讲的很清楚。 被卡常数至死。LOJACbzojTLE成SB 那个取mod啊，它是%2^31,所以说可以直接自然溢出。。。。最后and一个2147483647就可以了。 一个取mod这么恐怖，删了就3632ms，不删10sTLE
阅读全文
摘要：一个下午摸不到键盘系列。 看题解推柿子要爆炸。。 推到吐了，看到最后又变成约数个数心里一万只cnm飞过。。所以在电脑上就不重写了。 觉得这一辈子都不会忘记约数个数公式了。。 有(si)用(ji)的东(gong)西(shi)： d(nm)= sigema(x|n) sigema(y|m) [gcd(x
阅读全文
摘要：第三次复习反演。。。 感觉第二次学的时候还是有点用的 整理了下思路。可以去看一下我的blog 这题就是第二种类型的题。然而我不会化简K。。%PoPoQQQ
阅读全文
摘要：复习一波反演。 非常粗暴的就是F(1)+2。 注意范围应该是n-1，因为在坐标系上，只有在原点看，点(1,2)才会挡住(2,4)，最后加2就是上面和右边的两个点。 然而今天在caioj闲逛的时候发现了另外一种做法，就是用phi，抽象一下，设当前点为(x,y)那么应该xy互质才可以看到 （代码写得很丑
阅读全文
摘要：设F(i)=sigema(1~n) gcd(x,y)==i 的个数 ans=sigema(1~n,i为质数)i F(i); 设G(i)=sigema(1~n) gcd(x,y)==i 的倍数 的个数 那G(i)=(n/i)*(m/i)，反演。以上做法是和2818一样的，但是由于多组数据就会T 观察一
阅读全文
摘要：感觉莫反学的很不扎实。。很虚。。bs lmy 为什么这么说，是因为我连求F(i)要把他转换成F(1)再求都不知道。。。。 那么这题抽象点看ans就是等于这个 ans = n*m + sigema(1~n)x sigema(1~m)y (gcd(x,y)-1)*2 F(i)表示 gcd(x,y)==i
阅读全文
摘要：怎么感觉一直在做市选的说。。搞得在中山的我都没信心了。。。 好吧做这题主要是冲着莫比乌斯反演去的，然后实际上也是容斥原理的应用，跟反演没什么关系，但是莫比乌斯函数的一个应用。 首先将题目询问第k个无平方因子数，那么我立刻想到的是，这些数在莫比乌斯函数中的值是0（但是这个没什么用） 实际上的做法是二分
阅读全文
摘要：由于这个证明过程太过长。。。推荐大家看这个大佬的博文，我就是看这个学会的 http://www.cnblogs.com/chenyang920/p/4811995.html 原式 : G(n)=sigma(F(d)) (其中d|n，就是n%d==0) 反演公式: F(n)=sigma(U(n/d)*
阅读全文