NOIP 普及组初赛 · NOIP 2018 普及组初赛试题 第 22 题

完善程序

（最大公约数之和）下列程序想要求解整数 $n$ 的所有约数两两之间最大公约数的和对 $10007$ 求余后的值，试补全程序。（第一空 $2$ 分，其余 $3$ 分）

举例来说，$4$ 的所有约数是 $1, 2, 4$。$1$ 和 $2$ 的最大公约数为 $1$；$2$ 和 $4$ 的最大公约数为 $2$；$1$ 和 $4$ 的最大公约数为 $1$ 。于是答案为 $1 + 2 + 1 = 4$。

要求 getDivisor 函数的复杂度为 $O(\sqrt{n})$，gcd 函数的复杂度为$O(\log \max(a,b))$。

#include <iostream>
using namespace std;

const int N = 110000, P = 10007;
int n;
int a[N], len;
int ans;

void getDivisor() {
    len = 0;
    for (int i = 1; ① <= n; ++i)
        if (n % i == 0) {
          a[++len] = i;
          if ( ② != i) a[++len] = n / i;
        }
}

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

int main() {
    cin >> n;
    getDivisor();
    ans = 0;
    for (int i = 1; i <= len; ++i) {
        for (int j = i + 1; j <= len; ++j) {
        	ans = ( ⑤ ) % P;
        }
    }
    cout << ans << endl;
    return 0;
}

约数成对枚举至√n，GCD用辗转相除，累加需先加后取模。

循环条件勿用sqrt防浮点；n/i!=i防重复；GCD终止return a；递归参数a%b；累加ans+=gcd再取模。

考点：约数成对枚举｜欧几里得算法｜取模累加

程序在做什么

输入n，用O(√n)枚举所有约数存入数组a，再两两配对求GCD并累加对10007取余。