最大公约数之和

技术最大公约数之和 最大公约数之和本文主要讲一下最大公约数的和的推导过程(因为其太过经典,其实是博主老忘)。
原式:
\[\sum_{i = 1}^n\sum_{j = 1}^n\gcd(i, j)
\

最大公约数之和

本文主要讲一下最大公约数的和的推导过程(因为其太过经典,其实是博主老忘)。

原式:

\[\ sum _ { I=1}^n\sum_{j=1}^n\gcd(i,j)

\]莫比乌斯反演经典入门题。

话不多说,进入正文。

\[\begin{aligned}

\ sum \ limits _ { I=1}^n\sum\limits_{j=1}^ngcd(i,j)\ \

=\ sum \ limits _ { k=1}^nk\sum\limits_{i=1}^n\sum\limits_{j=1}^n[gcd(i,j)=k]\ \

=\ sum \ limits _ { k=1}^nk\sum\limits_{i=1}^{ \左侧\左侧\地板\地板{ n } { k }右侧\地板} \ sum \ limits _ { j=1}^{ \左侧\地板\地板{ n } { k }右侧\地板}[gcd(i,j)=1] \\

=\ sum \ limits _ { k=1}^nk\sum\limits_{i=1}^{ \左侧\左侧\右侧\右侧{ k } \ sum \ limits _ { j=1}^{ \左侧\左侧\左侧\右侧{ k } \右侧\左侧} \(gcd(I,j))

\ end {对齐}

\]根据\(\=\* I \),即\(\(n)=\ sum \ limits _ { d | n } \ mu(d)\),得:

\[\ sum \ limits _ { k=1}^nk\sum\limits_{i=1}^{ \左侧\右侧\地板\地板{ n } { k }右侧\地板} \ sum \ limits _ { j=1}^{ \左侧\地板\地板{ n } { k }右侧\地板}\sum\limits_{d | (i,j)}\mu(d)

\]我们先考虑这样一个式子如何化简:

\[\ sum \ limites _ { I=1}^{ \左侧\左侧\右侧\右侧{ k } \ sum \ limites _ { d | I } \ mu(d)

\]把枚举\(i\)改成枚举\(d\),\(\左侧\地板\地板{ n } { k }右侧\地板\)以内是\(d\)的倍数的数有\(\左侧\地板\ dfrac { n } { dk } \右侧\地板\)个,得:

\[\ sum \ limits _ { d=1}^{ \左侧\左侧\地板\地板{ n } { k }右侧\地板} \左侧\地板\地板{ n } { dk }右侧\地板\亩(d)

\]我们先枚举\(d\),并把这个式子代入到刚才我们化简得那个式子中去:

\[\begin{aligned}

\ sum \ limits _ { k=1}^nk\sum\limits_{i=1}^{ \左侧\右侧\地板\地板{ n } { k }右侧\地板} \ sum \ limits _ { j=1}^{ \左侧\地板\地板\地板{ n } { k }右侧\地板}\sum\limits_{d | (i,j)}\mu(d) \

=\ sum \ limits _ { k=1}^nk\sum\limits_{d=1}^{ \left\lfloor\frac{n}{k}\right\rfloor}\left\lfloor\dfrac{n}{dk}\right\rfloor^2\mu(d)

\ end {对齐}

\]再令\(T=dk\),并枚举\(T\)(其实下面的式子和上面的式子里\(d\)和\(k\)反过来了,不过我懒得改了QwQ):

\[\ sum _ { t=1}^n\sum_{d \ mid t } d \ mu(\ frac TD)\ f floor \ frac nt\rfloor^2

\]至此,就已经是一般形式了,这个可以用整除分块快速求解。

但是,这道题还没有完,还可以进一步转化。

我们知道\(\varphi=\mu * id\),正好式子里存在!所以:

\[\ sum _ { t=1}^n\varphi(t)\lfloor\frac nt\rfloor^2

\]现在,这道题才算是真正结束了(感觉一下子式子里啥都没了QwQ)

\[\_EOF\_

\]

本文来自博客园,作者:xixike,转载请注明原文链接:https://www .cn博客。com/xixixike/p/15713088。超文本标记语言

内容来源网络,如有侵权,联系删除,本文地址:https://www.230890.com/zhan/154011.html

(0)

相关推荐

  • [NOIP2016普及组]魔法阵

    技术[NOIP2016普及组]魔法阵 [NOIP2016普及组]魔法阵不是枚举暴力,也不是推式子 $\mathcal{O(1)}$,而是两者的有机结合——通过数学推导减少枚举量,满足时间复杂度要求。很大

    礼包 2021年11月11日
  • 黄老,黄老道家与道家是什么关系

    技术黄老,黄老道家与道家是什么关系黄老道家是田氏取代姜子牙后裔齐国政权后,为证明其执政合法性而在思想领域新创造的学说,上联黄帝是为了便于修正老子道学理论,特别是田氏要找一个高于姜子牙祖宗炎帝的祖宗,当然最合适的就是黄帝了

    生活 2021年10月20日
  • MySQL优化的最基础操作是什么

    技术MySQL优化的最基础操作是什么这篇文章主要介绍“MySQL优化的最基础操作是什么”,在日常操作中,相信很多人在MySQL优化的最基础操作是什么问题上存在疑惑,小编查阅了各式资料,整理出简单好用的操作方法,希望对大家

    攻略 2021年11月15日
  • MySQL慢SQL采集方案分析

    技术MySQL慢SQL采集方案分析本篇内容主要讲解“MySQL慢SQL采集方案分析”,感兴趣的朋友不妨来看看。本文介绍的方法操作简单快捷,实用性强。下面就让小编来带大家学习“MySQL慢SQL采集方案分析”吧!作为一名M

    攻略 2021年11月5日
  • C#基础语法之-泛型

    技术C#基础语法之-泛型 C#基础语法之-泛型泛型:一共7个知识点
    1.引入泛型,延迟声明
    2.如何声明和使用泛型
    3.泛型的好处和原理
    4.泛型类,泛型方法,泛型接口,泛型委托
    5.泛型约束
    6.协变

    礼包 2021年12月18日
  • vue 中query和pramas中的那点事

    技术vue 中query和pramas中的那点事 vue 中query和pramas中的那点事//$router : 是路由操作对象,只写对象
    //$route : 路由信息对象,只读对象
    //操作 路

    礼包 2021年11月11日