luogu P3911

题目

修复了一些小bug,(。・_・。)ノI’m sorry~

给出集合 AA,设 n=∣An=A,求 ni=1nj=1lcm(Ai,Aj)ni=1nj=1lcm(Ai,Aj)

摸到这道题,我们马上就可以发现 AA 中的元素是无序的,题目中给出的式子没法化简,所以我们不得不考虑转化式子。

数据范围 1Ai5e41Ai5e4,这就很好办了……

开一个桶(这里用数组 tt 表示),t[i]t[i] 表示数 ii 在集合中出现的次数。

这样原式就转化为了:

ni=1nj=1lcm(i,j)×t[i]×t[j]ni=1nj=1lcm(i,j)×t[i]×t[j]

注意此后 nn 指集合 AA 中元素的最大值。

这样好像就可以转化 gcdgcd 求解了。

(为了此篇题解更通俗易懂,接下来的推导会有些冗余部分,见谅~)

愉快的推式子时间~:

ni=1nj=1lcm(i,j)×t[i]×t[j]ni=1nj=1i×jgcd(i,j)×t[i]×t[j]nd=1ni=1nj=1i×j[gcd(i,j)=d]×d×t[i]×t[j]nd=1ndi=1ndj=1id×jd[gcd(i,j)=1]×d×t[id]×t[jd]nd=1ndi=1ndj=1d×i×j×[gcd(i,j)=1]×t[id]×t[jd]nd=1ndi=1ndj=1d×i×j×ki,kjμ(k)×t[id]×t[jd]nd=1ndi=1ndj=1ki,kjμ(k)×d×i×j×t[id]×t[jd]nd=1ndk=1nkdi=1nkdj=1μ(k)×d×ik×jk×t[ikd]×t[jkd]ns=1s×ksμ(k)×k (nsi=1i×t[is])2ni=1nj=1lcm(i,j)×t[i]×t[j]ni=1nj=1i×jgcd(i,j)×t[i]×t[j]nd=1ni=1nj=1i×j[gcd(i,j)=d]×d×t[i]×t[j]nd=1ndi=1ndj=1id×jd[gcd(i,j)=1]×d×t[id]×t[jd]nd=1ndi=1ndj=1d×i×j×[gcd(i,j)=1]×t[id]×t[jd]nd=1ndi=1ndj=1d×i×j×ki,kjμ(k)×t[id]×t[jd]nd=1ndi=1ndj=1ki,kjμ(k)×d×i×j×t[id]×t[jd]nd=1ndk=1nkdi=1nkdj=1μ(k)×d×ik×jk×t[ikd]×t[jkd]ns=1s×ksμ(k)×k (nsi=1i×t[is])2

一些解释(大佬可直接跳过):

  1. s=kds=kd
  2. 结果中的平方为上面的式子中的 iijj,他俩其实是一样的。

最后,我们使用线性筛预处理 ksμ(k)×kksμ(k)×k,输入时直接处理 tt,然后暴力就完了。

代码按照上面的暴力即可,就不给了 QVQ。

posted @   aleph_blanc  阅读(28)  评论(0编辑  收藏  举报
编辑推荐:
· 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后端服务间歇性响应慢的问题排查记录
点击右上角即可分享
微信分享提示