luogu P3911
修复了一些小bug,(。・_・。)ノI’m sorry~
给出集合 AA,设 n=∣A∣n=∣A∣,求 ∑ni=1∑nj=1lcm(Ai,Aj)∑ni=1∑nj=1lcm(Ai,Aj)。
摸到这道题,我们马上就可以发现 AA 中的元素是无序的,题目中给出的式子没法化简,所以我们不得不考虑转化式子。
数据范围 1≤Ai≤5e41≤Ai≤5e4,这就很好办了……
开一个桶(这里用数组 tt 表示),t[i]t[i] 表示数 ii 在集合中出现的次数。
这样原式就转化为了:
n∑i=1n∑j=1lcm(i,j)×t[i]×t[j]n∑i=1n∑j=1lcm(i,j)×t[i]×t[j]
注意此后 nn 指集合 AA 中元素的最大值。
这样好像就可以转化 gcdgcd 求解了。
(为了此篇题解更通俗易懂,接下来的推导会有些冗余部分,见谅~)
愉快的推式子时间~:
n∑i=1n∑j=1lcm(i,j)×t[i]×t[j]⇒n∑i=1n∑j=1i×jgcd(i,j)×t[i]×t[j]⇒n∑d=1n∑i=1n∑j=1i×j[gcd(i,j)=d]×d×t[i]×t[j]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1id×jd[gcd(i,j)=1]×d×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1d×i×j×[gcd(i,j)=1]×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1d×i×j×∑k∣i,k∣jμ(k)×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1∑k∣i,k∣jμ(k)×d×i×j×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑k=1⌊nkd⌋∑i=1⌊nkd⌋∑j=1μ(k)×d×ik×jk×t[ikd]×t[jkd]⇒n∑s=1s×∑k∣sμ(k)×k (⌊ns⌋∑i=1i×t[is])2n∑i=1n∑j=1lcm(i,j)×t[i]×t[j]⇒n∑i=1n∑j=1i×jgcd(i,j)×t[i]×t[j]⇒n∑d=1n∑i=1n∑j=1i×j[gcd(i,j)=d]×d×t[i]×t[j]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1id×jd[gcd(i,j)=1]×d×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1d×i×j×[gcd(i,j)=1]×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1d×i×j×∑k∣i,k∣jμ(k)×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑i=1⌊nd⌋∑j=1∑k∣i,k∣jμ(k)×d×i×j×t[id]×t[jd]⇒n∑d=1⌊nd⌋∑k=1⌊nkd⌋∑i=1⌊nkd⌋∑j=1μ(k)×d×ik×jk×t[ikd]×t[jkd]⇒n∑s=1s×∑k∣sμ(k)×k (⌊ns⌋∑i=1i×t[is])2
一些解释(大佬可直接跳过):
- s=kds=kd。
- 结果中的平方为上面的式子中的 ii 和 jj,他俩其实是一样的。
最后,我们使用线性筛预处理 ∑k∣sμ(k)×k∑k∣sμ(k)×k,输入时直接处理 tt,然后暴力就完了。
代码按照上面的暴力即可,就不给了 QVQ。
【推荐】还在用 ECharts 开发大屏?试试这款永久免费的开源 BI 工具!
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· dotnet 源代码生成器分析器入门
· ASP.NET Core 模型验证消息的本地化新姿势
· 对象命名为何需要避免'-er'和'-or'后缀
· SQL Server如何跟踪自动统计信息更新?
· AI与.NET技术实操系列:使用Catalyst进行自然语言处理
· dotnet 源代码生成器分析器入门
· 官方的 MCP C# SDK:csharp-sdk
· 一款 .NET 开源、功能强大的远程连接管理工具,支持 RDP、VNC、SSH 等多种主流协议!
· 一步一步教你部署ktransformers,大内存单显卡用上Deepseek-R1
· 一次Java后端服务间歇性响应慢的问题排查记录