约数——数学知识（C++）

一、定义

若整数n除以整数d的余数为0，即d能整除n，则称d是n的约数，n是d的倍数，记为d|n。
算数基本定理的推论
一个大于1的正整数N，如果它的标准分解式为：

那么它的正因数个数为

它的全体正因数之和为

二、约数的判定：试除法

AcWing 869. 试除法求约数
给定n个正整数ai，对于每个整数ai,请你按照从小到大的顺序输出它的所有约数。

输入格式
第一行包含整数n。

接下来n行，每行包含一个整数ai。

输出格式
输出共n行，其中第 i 行输出第 i 个整数ai的所有约数。

数据范围
1≤n≤100,
2≤ai≤2∗109
输入样例：

2
6
8

输出样例：

1 2 3 6 
1 2 4 8 

众所周知，除了可能有一对同样的数相乘等于x以外，其余的一对约数一个小于sqrt(x)、一个大于sqrt(x)

优化思路:
a|n a为n的约数, 那么 n/a|n ，n/a也是n的约数
并且一个数的因数的成对出现，大于sqrt(n)的因数只有一个。
反证：若大于sqrt(n)的因数有两个，sqrt(n)*sqrt(n)已经大于n，就矛盾了（因为一个数的因数不可能大于n）
所以只需1 ~ sqrt(n) 找约数即可

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

vector<int> get_divisors(int x)
{
    vector<int> res;
    for (int i = 1; i <= x / i; i ++ )//基本数学常识 只用除到sqrt(x)
        if (x % i == 0)// 此时i 为 n 的一个 因子（约数）
        {
            res.push_back(i);
            if (i != x / i) res.push_back(x / i);//防止i 与 n/ i 相等
        }
    sort(res.begin(), res.end());
    return res;
}

int main()
{
    int n;
    cin >> n;

    while (n -- )
    {
        int x;
        cin >> x;
        auto res = get_divisors(x);

        for (auto x : res) cout << x << ' ';作用就是迭代容器中所有的元素，每一个元素的临时名字就是x
        cout << endl;
    }

    return 0;
}

三、约数的应用1：求约数个数/求约数之和

定理：

1、对于一个数N，可以被分解质因数为（其中pi为质因数，ai为指数）N=(p1^a1) * (p2^a2) * (p3^a3) ……

2、那么N的约数个数为(a1+1) * (a2+1) * (a3+1)……

3、那么N的约数之和为(p1⁰+p1¹+……+p1^a1)*(p2⁰+p2¹+……+p2^a2) …… *(pk⁰+pk¹+……+pk^ak)

原因：

约数个数：只要是N的约数，分解出来的质因数就一定是与N的质 因 数的数字相同或者少几个、而且每个质因数的指数范围是0~ai。因为可以取到0也可以取到ai，所以每个质因数有(ai+1)种选法。

约数之和：把式子按乘法分配律展开可以发现是()+()+……+()约数个数个括号相加
所以约数个数即为各种选法相乘

经典例题

AcWing 870. 约数个数
给定n个正整数ai，请你输出这些数的乘积的约数个数，答案对10⁹+7取模。

输入格式
第一行包含整数n。

接下来n行，每行包含一个整数ai。

输出格式
输出一个整数，表示所给正整数的乘积的约数个数，答案需对109+7取模。

数据范围
1≤n≤100,
1≤ai≤2∗10⁹
输入样例：

3
2
6
8

输出样例：

12

#include <iostream>
#include <algorithm>
#include <unordered_map>
#include <vector>

using namespace std;

typedef long long LL;

const int N = 110, mod = 1e9 + 7;

int main()
{
    int n;
    cin >> n;

    unordered_map<int, int> primes;//f.first 存放底数,f.second存放f.first的指数

    while (n -- )
    {
        int x;
        cin >> x;

        for (int i = 2; i <= x / i; i ++ )
            while (x % i == 0)
            {
                x /= i;
                primes[i] ++ ;//用哈希表存因式分解中的指数
            }

        if (x > 1) primes[x] ++ ;//如果还没分解清,说明剩下的这个数是它的因数且只有一个
    }

    LL res = 1;
    for (auto p : primes) res = res * (p.second + 1) % mod;//约数个数=(a1+1)(a2+1)...(ak+1)

    cout << res << endl;

    return 0;
}

例题

AcWing 871. 约数之和
给定n个正整数ai，请你输出这些数的乘积的约数之和，答案对10⁹+7取模。

输入格式
第一行包含整数n。

接下来n行，每行包含一个整数ai。

输出格式
输出一个整数，表示所给正整数的乘积的约数之和，答案需对109+7取模。

数据范围
1≤n≤100,
1≤ai≤2∗10⁹
输入样例：

3
2
6
8

输出样例：

252

算法分析：
基本思想：
如果 N=p1^c1∗p2^c2∗…∗pk^ck
约数个数：(c1+1)∗(c2+1)∗…∗(ck+1)
约数之和： (p1⁰+p1¹+…+p1^c1)∗…∗(pk⁰+pk¹+…+pk^ck)

while (b -- ) t = (t * a + 1) % mod;

t=t∗p+1
t=1
t=p+1
t=p²+p+1
……
t=p^b+p^b−1+…+1

#include <iostream>
#include <algorithm>
#include <unordered_map>
#include <vector>

using namespace std;

typedef long long LL;

const int N = 110, mod = 1e9 + 7;

int main()
{
    int n;
    cin >> n;

    unordered_map<int, int> primes;

    while (n -- )
    {
        int x;
        cin >> x;

        for (int i = 2; i <= x / i; i ++ )
            while (x % i == 0)
            {
                x /= i;
                primes[i] ++ ;
            }

        if (x > 1) primes[x] ++ ; //注意是x>1
    }

    LL res = 1;
    for (auto p : primes)
    {
        LL a = p.first, b = p.second;//p为底数,a为指数
        LL t = 1;
        while (b -- ) t = (t * a + 1) % mod;//求出 p0一直加到p的k的次方 的和//先求出公式中的每个括号内的值
        res = res * t % mod;//把每个括号的值相乘
    }

    cout << res << endl;

    return 0;
}

四、约数的应用2：求最大公约数

求a和b的最大公约数:求a和b的最大公约数:欧几里得算法(辗转相除法)

gcd(a,b)=gcd(b,a%b);

证明:

设a和b的最大公约数为d.
∵d∣a且d∣b 且 a %b=a−⌊a/b⌋∗b
∴d∣(a−⌊a/b⌋∗b) 即d∣(a%b)
∴gcd(a,b)=gcd(b,a%b)得证

int gcd(int a, int b)
{
	return b ? gcd(b,a%b) : a;
} //当b被约成0的时候，此时对应的a即为最小公约数

经典例题

AcWing 872. 最大公约数
给定n对正整数ai,bi，请你求出每对数的最大公约数。

输入格式
第一行包含整数n。

接下来n行，每行包含一个整数对ai,bi。

输出格式
输出共n行，每行输出一个整数对的最大公约数。

数据范围
1≤n≤10⁵,
1≤ai,bi≤2∗10⁹
输入样例：

2
3 6
4 6

输出样例：

3
2

#include <iostream>
#include <algorithm>

using namespace std;


int gcd(int a, int b)
{
    return b ? gcd(b, a % b) : a;
}


int main()
{
    int n;
    cin >> n;
    while (n -- )
    {
        int a, b;
        scanf("%d%d", &a, &b);
        printf("%d\n", gcd(a, b));
    }

    return 0;
}