【学习笔记】[HDU 5407] CRB and Candies

题意:多组数据,每次给定一个 n ,求 l c m ( C ( n , 0 ) , C ( n , 1 ) , . . . , C ( n , n ) )   m o d   1 0 9 + 7 lcm(C(n,0),C(n,1),...,C(n,n))\bmod 10^9+7 lcm(C(n,0),C(n,1),...,C(n,n))mod109+7

solution:
直接算肯定不好算。

我们证明一个结论: l c m ( C ( n , 0 ) , C ( n , 1 ) , . . , C ( n , n ) ) = ( l c m ( 1 , 2 , . . . , n + 1 ) n + 1 ) lcm(C(n,0),C(n,1),..,C(n,n))=(\frac{lcm(1,2,...,n+1)}{n+1}) lcm(C(n,0),C(n,1),..,C(n,n))=(n+1lcm(1,2,...,n+1))

定义 vp(x) 表示 p 在 x 中的最大次幂。

看看 Kummer 定理:

vp(C(n+m,m)) = n + m 在 p 进制下进位的次数。

然后是一个推论:max(vp(C(n,x))) = n 在 p 进制下的位数 - 最低位不为 p-1 的位数。

可以贪心证明。

再回头看结论,我们需要证明:

v p ( l c m ( C ( n , 0 ) , C ( n , 1 ) , . . . , C ( n , n ) ) ) = v p ( l c m ( 1 , 2 , . . . , n + 1 ) n + 1 ) vp(lcm(C(n,0),C(n,1),...,C(n,n))) = vp(\frac{lcm(1,2,...,n+1)}{n+1}) vp(lcm(C(n,0),C(n,1),...,C(n,n)))=vp(n+1lcm(1,2,...,n+1))

先看右边式子。 vp(lcm(1,2,…,n+1)) = p 进制下 n+1 的位数

如果 n+1 = p^x , 那么 vp(lcm(C(n,0),C(n,1),…,C(n,n))) = 0, vp(\frac{n+1}{n+1} = 1) 。

如果 n+1 \ne p^x, 那么 vp(lcm(1,2,…,n+1)) = vp(lcm(1,2,…,n)) = n 在 p 进制下的位数 , vp(n+1) = n+1 在 p 进制下末尾 0 的个数 = n 在 p 进制下末尾 p-1 的个数 。

证毕哈哈。拿着这个式子 O(n) 预处理即可。


__EOF__

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