约数【试除法+倍数法】

定义：若整数 n 除以整数 d 的余数为 0，即 d 能整除 n ，则称 d 是 n 的约数，n 是 d 的倍数，记为 d | n。

 求 N 的正约数集合------试除法

         若 d >= sqrt(n) 是 n 的约数，则 n/d <= sqrt(n)也是 n 的约数。总之，约数是成对出现的（除了完全平方数，sqrt(n) 会单独出现）。因此，只需要扫描 d = 1~sqrt(n)，尝试 d neng能否整除 n，若能整除，则 n / d 也是 n 的约数。时间复杂度为 O(sqrt(n))。

#include <iostream>
#include <cstdio>
using namespace std;
const int maxn = 10000+7;
int main()
{
    int factor[maxn], m, n;
    m = 0;
    scanf("%d", &n);
    for(int i = 1; i*i <= n; i++) {
        if(n % i ==0){
            factor[++m] = i;
            if(i != n/i) factor[++m] = n/i;
        }
    }
    for(int i = 1; i <= m; i++)
        printf("%d\n", factor[i]);
    return 0;
}

试除法的推论：一个整数 n 的约数个数上界为 2*sqrt(n)。

求 1~n 每个数的正约数集合------倍数法：若用试除法分别求出 1~n 每个数的正约数集合，时间复杂度过高，为O(n*sqrt(n))。可以反过来考虑，对于每个数 d ，1~n 中以 d 为约数的数就是 d 的倍数 d,2d,3d......( n / d)*d。时间复杂度为 O(n*logn)。

#include <iostream>
#include <vector>
#include <cstdio>
#define ll long long
using namespace std;
const int maxn = 10000+7;
int main()
{
    int n;
    vector<int> factor[maxn];
    scanf("%d", &n);
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= n/i; j++){
            factor[i*j].push_back(i);
        }
    }
    for(int i = 1; i <= n; i++){
        for(int j = 0; j < factor[i].size(); j++){
            printf("%d ",factor[i][j]);
        }
        printf("\n");
    }
    return 0;
}

倍数法的推论：1~n 每个数的约数个数的总和大约为 nlogn。