【学习笔记】[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 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力!
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」