【题型探究】公约数之和

整理的算法模板合集： ACM模板

---

母题：UVA11417 GCD

数据很小，我们只需要直接暴力枚举即可。

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
#include <iostream>

using namespace std;
typedef long long ll;
const int N = 5000007;

int gcd(int a, int b)
{
   
    if(b == 0) return a;
    return gcd(b, a % b);
}

ll solve(int n)
{
   
    ll res = 0;
    for(int i = 1; i <= n; ++ i){
   
        for(int j = i + 1; j <= n; ++ j){
   
            res += gcd(i, j);
        }
    }
    return res;
}
int n;
int main()
{
   
    while(scanf("%d", &n) != EOF && n){
   
        printf("%d\n", solve(n));
    }
    return 0;
}

拓展一、UVA11426 GCD - Extreme (II)

本题的数据开到了$4\ast 10^6$，所以暴力一定TLE。

啊，这，范围太大了，我可以莫比乌斯反演 + 整除分块，$O(\sqrt{n})$ 解决！

《算法竞赛入门经典训练指南》上提供了一个构造函数的方法：

我们设$f(n)=gcd(1,n)+gcd(2,n)+gcd(3,n)+...+gcd(n-1,n)$

答案很明显就是$S(n)=f(2)+f(3)+...+f(n)$

我们首先考虑如何求$f(n)$

我们可以先把 g c d