给定 n,k,求
∑i=1n∑j=1n(i+j)kμ2(gcd(i,j))gcd(i,j)
数据范围:1≤n≤5×106,1≤k≤1018
∑i=1n∑j=1n(i+j)kμ2(gcd(i,j))gcd(i,j)=∑d=1nμ2(d)dk+1∑i=1⌊nd⌋∑j=1⌊nd⌋(i+j)k[gcd(i,j)=1]=∑d=1nμ2(d)dk+1∑i=1⌊nd⌋∑j=1⌊nd⌋(i+j)k∑s|i∑s|jμ(s)=∑d=1nμ2(d)dk+1∑s=1⌊nd⌋μ(s)sk∑i=1⌊nsd⌋∑j=1⌊nsd⌋(i+j)k=∑q=1nqk∑d|qμ2(d)dμ(qd)∑i=1⌊nq⌋∑j=1⌊nq⌋(i+j)k=∑q=1n∑i=1⌊nq⌋∑j=1⌊nq⌋(i+j)kqk∑d|qμ2(d)dμ(qd) 令 S(n)=∑i=1n∑j=1n(i+j)k,g(n)=nk∑d|nμ2(d)dμ(nd)
∑i=1n∑j=1n(i+j)kμ2(gcd(i,j))gcd(i,j)=∑q=1nS(⌊nq⌋)g(q)
对于 g(n),线性筛出 ik,O(nlogn) 预处理即可。
考虑 S(n) 怎么求,令 f(n)=∑i=1nik,所以 S(n)=S(n−1)+(2n)k+2[f(2n−1)−f(n+1)],直接推即可。
左边整除分块,右边直接 O(nlogn) 暴力搞出来,时间复杂度 O(nlogn)。
sto
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· 25岁的心里话