约数

约数

试除法求约数

Code

//获取一个数的所有约数，从小到大 
vector<int> get_divisors(int n) {
	vector<int> res;
	for(int i = 1; i <= n / i; i++) { //最多只要枚举到sqrt(n) 
		if(n % i == 0) {
			res.push_back(i);
			if(i != n / i) res.push_back(n / i); //把它成对的约数push进vector 
		}
	}
	sort(res.begin(), res.end());

	return res;
}

求一个数所有约数的个数

$N=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}...p_n^{{\alpha }_n}$ ，其中，$p$ 是 $N$ 的质因数，${\alpha }_i$ 是约数的指数，则一个数所有约数的个数为

$ans=(a_1+1)(a_2+1)(a_3+1)...(a_n+1)$

理由是，任何一个约数 $d$ 可以表示成 $p^{{\beta }_1}{p}^{{\beta }_2}{p}^{{\beta }_3}...p^{{\beta }_n},0\le {\beta }_i\le {\alpha }_i$ ，每一项的 ${\beta }_i$ 如果不同，那么约数 $d$ 就不相同（算数基本定理，每个数的因式分解是唯一的）所以，$n$ 的约数就跟 ${\beta }_i$ 的选法是一一对应的 ${\beta }_1$ 一共有 $0\sim {\alpha }_1$ 种选法，${\beta }_2$ 一共有 $0\sim {\alpha }_2$ 种选法 $\to$ ${\beta }_k$ 一共有 $0\sim {\alpha }_k$ 种选法，根据乘法原理，一共有 $ans$ 个约数

一个小结论：$int$ 范围内，约数最多的一个数的约数大概是1500个

Code

ll divnum(ll x) {
	unordered_map<ll, ll> primes;
	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 prime : primes) res = res * (prime.second + 1);
	
	return res;
}

练习

C. Common Divisors

约数之和

$N=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}...p_n^{{\alpha }_n}$ ，其中，$p$ 是 $N$ 的质因数，${\alpha }_i$ 是约数的指数，则一个数的所有约数之和为

$ans=(p_1^0+p_1^1+p_1^2+...p_1^{{\alpha }_1})\times (p_2^0+p_2^1+p_2^2+...p_2^{{\alpha }_2})\times (p_3^0+p_3^1+p_3^2+...p_3^{{\alpha }_3})\times ...\times (p_k^0+p_k^1+p_k^2+...p_k^{{\alpha }_k})$

使用乘法分配率展开，每一项都是约数

如何计算每一个项呢

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

通过对计算步骤的复现可以知道

当执行 1 次时，$t=(p+1)$

当执行 2 次时，$t=(p+1)\times p+1=p^2+p+1$

当执行 3 次时，$t=p^2+p+1\times p+1=p^3+p^2+p+1$

这样可以知道，执行 ${\alpha }_k$ 次之后，得到的式子就是

$t=p_1^0+p_1^1+p_1^2+...p_1^{{\alpha }_1}$

Code

给定 $n$ 个正整数 $a_i$ ，请你输出这些数的乘积的约数之和，答案对 $10^9+7$ 取模。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const long long mod = 1e9 + 7;

int main() {
    int T;
    scanf("%d", &T);
    unordered_map<int, int> primes;//存储分解完成的质数
    while(T--) {
        int x;
        scanf("%d", &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 prime : primes) {
        int p = prime.first, a = prime.second; //p表示约数 a表示约数的指数
        ll t = 1;
        while(a--) t = (t * p + 1) % mod;//根据公式计算每一项
        res = res * t % mod;
    }

    printf("%lld\n", res);

    return 0;
}

求最大公约数

使用

__gcd(a,b);