【刷题】数学知识——约数：求约数个数

假设一个数n可以拆成$n=p_1^{a_1}\ast p_2^{a_2}\ast ...\ast p_k^{a_k}$，
其中p是质数，如$36=2^2\ast 3^2$,
那么约数的个数就有$(a_1+1)\ast (a_2+1)\ast ...\ast (a_k+1)$个，即将是质数的约数个数+1后相乘
如36的约数个数就有$(2+1)\ast (2+1)=9$个：1 2 3 4 6 9 12 18 36。

这是因为每个约数必然是由几个质数组合而成，不同的约数所选的质数组合不一样。

而可选的质数只有$a_1$个$p_1$，$a_2$个$p_2$，…，$a_k$个$p_k$。
问题转化成：从这些数中，共有多少种不同的选法。

因为第i个质数能选0~$a_i$次，有$a_i+1$种选择，那么所有k个质数的选法就是$(a_1+1)\ast (a_2+1)\ast ...\ast (a_k+1)$种

所以只要用试除法统计所有数的所有质数个数，将各个数的质数个数相加得到$a_1,...,a_k$，再用$(a_1+1)\ast (a_2+1)\ast ...\ast (a_k+1)$相乘，就能得到答案。

PS. int范围内约数个数最多的数，约数有1500个左右。

#include <iostream>
#include <unordered_map> 
using namespace std;

long long ans = 1;
const int mod = 1e9 + 7;

int n, a;

unordered_map<int, int> primes;

int main() {
	scanf("%d", &n);
	while (n -- ) {
		scanf("%d", &a);
		
		for (int i = 2; i * i <= a; i ++ ) {
			while (a % i == 0) {
				a /= i;
				primes[i] ++;  // 各个数的约数个数相加 
			}
		}
		if (a > 1) primes[a] ++; // a是质数时 
	}
	for (auto prime: primes) {
		ans = ans * (prime.second + 1) % mod;
	}
	printf("%lld\n", ans); 
	return 0;
}

这边和试除法有些不同：
1、没有同时把$i$和$a/i$统计进去，这是因为这题只要统计是质数的约数。
2、为了确保把最大的质数统计进去，在结尾加了if (a > 1) primes[a] ++。
如$105=3\ast 5\ast 7$，只会在前面的循环中统计3和5，此时a=7，$\sqrt{7}<5$，退出循环，不会记录下7，因此要加上这句。
至于为什么至多只有一个质数没在循环中统计进去，这是因为枚举到$a=i\ast i$时，$i+1$到$a$至多只有1个质数是a的约数。如果有两个质数是a的约数，那么这两个质数相乘>a，矛盾。