约数相关算法

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

using namespace std;
typedef long long LL;




// 试除法求所有约数
vector<int> get_divisors(int n)
{
	vector<int> res;

	// 约数是成对出现的，一个大于sqrt(n),一个小于sqrt(n)
	for (int i = 1; i <= n / i; i++)
	{
		if (n % i == 0)
		{
			res.push_back(i);

			if (i != n / i)
			{
				res.push_back(n / i);
			}
		}
	}

	// 排序
	sort(res.begin(), res.end());

	return res;
}

// 约数个数
// 一个数的分解质因数为：p1^a1 * p2^a2 * p3^a3 … pk^ak
// 则该数的约数个数为：(a1 + 1) * (a2 + 1) * (a3 + 1) …… (ak + 1)
int get_divisors_num(int n)
{
	unordered_map<int, int> primes; // 哈希表存储质因子及其质数(pi 和 ai)

	// 分解质因数 
	for (int i = 2; i < n / i; i++)
	{
		if (n % i == 0)
		{
			while (n % i == 0)
			{
				n /= i;
				primes[i]++;
			}
		}
	}

	if (n > 1)
	{
		primes[n]++;
	}

	// 公式求约数个数
	int res = 1;

	for (auto prime : primes)
	{
		res *= (prime.second + 1);
	}

	return res;
}


// 约数之和
// 该数的约数之和为：(p1^0 + p1^1 + p1^2 + p1^3 + …… + p1^a1) * …… * (pk^0 + pk^1 + pk^2 + pk^3 + …… + pk^ak)
int get_divisors_sum(int n)
{
	unordered_map<int, int> primes; // 哈希表存储质因子及其质数(pi 和 ai)

	// 分解质因数 
	for (int i = 2; i < n / i; i++)
	{
		if (n % i == 0)
		{
			while (n % i == 0)
			{
				n /= i;
				primes[i]++;
			}
		}
	}

	if (n > 1)
	{
		primes[n]++;
	}

	// 公示求约数之和
	int res = 1;

	for (auto prime : primes)
	{
		int p = prime.first;
		int a = prime.second;

		LL t = 1;
		while (a--)
		{
			t = (t * p + 1);
		}

		res *= t;
	}

	return res;
}

// 最大公约数
int gcd(int a, int b)
{
	// 如果b为0，最大公约数为a
	// 如果b不为0，返回b和a%b的最小公约数
	return b ? gcd(b, a % b) : a;
}

int main()
{


	return 0;
}