AcWing算法基础课笔记——约数

约数

（1）试除法求约数

时间复杂度 O(sqrt(n))

题目

题目描述 ：

给定 n 个正整数 ai，对于每个整数 ai，请你按照从小到大的顺序输出它的所有约数。

输入输出格式 ：

输入

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

输出

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

输入输出样例 ：

输入

2
6
8

输出

1 2 3 6
1 2 4 8

代码

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

using namespace std;

vector<int> get_divisors(int n) {
	vector<int> res;
	
	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;
} 

int main() {
	int n;
	cin >> n;
	
	while(n -- ) {
		int x;
		cin >> x;
		auto res = get_divisors(x);
		for(auto t : res) cout << t << ' ';
		cout << endl;
	}
	
	return 0;
}

（2）约数个数

如果一个数N可以被写成
$$N=p_1^{{\alpha }_1}\times p_2^{{\alpha }_2}\cdots \times p_k^{{\alpha }_k}$$
那么它的约数个数为

$({\alpha }_1+1)\times ({\alpha }_2+1)\cdots \times ({\alpha }_k+1)$

int范围内约数个数最多大概1500

约数个数为
$$(p_1^0+p_1^1+\cdots +p_1^{{\alpha }_1})\times \cdots \times (p_k^0+p_k^1+\cdots +p_k^{{\alpha }_k})$$

题目——约数个数

题目描述

给定 𝑛 个正整数 𝑎_𝑖，请你输出这些数的乘积的约数个数，答案对 10^9+7 取模。

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

接下来 𝑛 行，每行包含一个整数 𝑎_𝑖。

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

数据范围
1≤n≤100,1≤ 𝑎_𝑖≤2×109

输入样例：

3
2
6
8

输出样例：

12

代码

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

using namespace std;

typedef long long LL;

const int 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] ++; 
	}
	
	LL res = 1;
	for(auto prime : primes) res = res * (prime.second + 1) % mod;
	
	cout << res << endl;
	
	return 0;
} 

（3）约数之和

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

using namespace std;

typedef long long LL;

const int 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] ++; 
	}
	
	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;
	}
	
	cout << res << endl;
	
	return 0;
} 

（4）欧几里得算法

也叫辗转相除法

题目——最大公约数

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

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

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

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

数据范围
1≤n≤105,
1≤ai,bi≤2∗109
输入样例
2
3 6
4 6
输出样例
3
2

代码

#include<iostream>

using namespace std;

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

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