【题解】P2303 [SDOI2012] Longge 的问题

题目链接

\(\displaystyle\sum_{i=1}^n \gcd(i,n)\)

将这个柿子展开变复杂,得到

\[\displaystyle\sum_{i=1}^{n}\sum_{d \mid n,d \mid i}d[\gcd(i,n)=d] \]

因为仅有 \(d\mid n\)\(d \mid i\) 的时候可能会累加 \(d\) ,则有效的 \(i\) 一定是 \(d\) 的倍数,即 \(i = kd\)

那么原式就可以化为

\[\displaystyle\sum_{d \mid n}\sum_{k=1}^{kd\leq n} d [\gcd(kd,n)=d] \]

在这个柿子中,\(d\) 是常量,将其提出来,并将 \(kd\leq n\) 移项得 \(k \leq \dfrac{n}{d}\) 得到

\[\displaystyle\sum_{d\mid n}d\sum_k^{\frac{n}{d}}[\gcd(kd,n)=d] \]

艾弗森括号里面的也能同时除个 \(d\) 得到 \(\displaystyle\sum_{d\mid n}d\sum_k^{\frac{n}{d}}[\gcd(k,\dfrac{n}{d})=1]\)

由于 \(\displaystyle\varphi(n)=\sum_{i=1}^{n}[\gcd(i,n)]\)

所以原式就化成了 \(\displaystyle\sum_{d\mid n}d\times \varphi(\dfrac{n}{d})\)

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;

ll getPhi(ll n){
    ll res = n;
    for(ll i = 2;i*i<=n;i++){
        if(n % i == 0) res = res / i * (i - 1);
        while(n % i == 0) n /= i;
    }
    if(n > 1) res = res / n * (n - 1);
    return res;
}

ll N;

int main(){
    scanf("%lld",&N);
    ll res = 0;
    ll i = 1;
    for(;i*i<N;i++){
        if(N % i == 0) res += i * getPhi(N / i) + (N / i) * getPhi(i);
    }
    if(i * i == N) res += i * getPhi(i);
    printf("%ll\n",res);
    return 0;
}

原文地址:http://www.cnblogs.com/burnling/p/16917575.html

1. 本站所有资源来源于用户上传和网络,如有侵权请邮件联系站长! 2. 分享目的仅供大家学习和交流,请务用于商业用途! 3. 如果你也有好源码或者教程,可以到用户中心发布,分享有积分奖励和额外收入! 4. 本站提供的源码、模板、插件等等其他资源,都不包含技术服务请大家谅解! 5. 如有链接无法下载、失效或广告,请联系管理员处理! 6. 本站资源售价只是赞助,收取费用仅维持本站的日常运营所需! 7. 如遇到加密压缩包,默认解压密码为"gltf",如遇到无法解压的请联系管理员! 8. 因为资源和程序源码均为可复制品,所以不支持任何理由的退款兑现,请斟酌后支付下载 声明:如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性