LightOJ 1236 Pairs Forming LCM (唯一分解定理)

题意：求 (i, j) 的组数，使得 lcm(i, j) = n (i ≤ j)。

题解：唯一分解定理
对于因子公式 n = p1^e1* p2^e2*…*pm^em ，要使得lcm(i, j) = n，那么a与b中，对于pi的指数，至少应有一个是等于ei的。那么对于一个因子，个数就是2 *ei + 1，遍历所有因子，乘起来即可。
注意i ≤ j，那么除了(n, n)之外其他组都出现两次。

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<fstream>
#include<map>
#define ll long long
using namespace std;

int t;
ll n;
//素数筛选和合数分解
const int MAXN = 1e7 + 10;
bool notprime[MAXN];
int prime[MAXN / 10];
void getPrime() {
	memset(notprime, false, sizeof(notprime));
	notprime[0] = notprime[1] = true;
	prime[0] = 0;
	for (int i = 2; i < MAXN; i++) {
		if (!notprime[i])prime[++prime[0]] = i;
		for (int j = 1; j <= prime[0] && prime[j] < MAXN / i; j++) {
			notprime[prime[j] * i] = true;
			if (i % prime[j] == 0) break;
		}
	}
}
long long factor[1000][2];
int fatCnt;
int getFactors(long long x) {
	fatCnt = 0;
	long long tmp = x;
	for (int i = 1; prime[i] <= tmp / prime[i]; i++) {
		factor[fatCnt][1] = 0;
		if (tmp % prime[i] == 0) {
			factor[fatCnt][0] = prime[i];
			while (tmp % prime[i] == 0) {
				factor[fatCnt][1]++;
				tmp /= prime[i];
			}
			fatCnt++;
		}
	}
	if (tmp != 1) {
		factor[fatCnt][0] = tmp;   //存放分解的素数
		factor[fatCnt++][1] = 1;   //存放此素数出现次数
	}
	return fatCnt;   //分解的素数个数（不相同）
}
int main() {
	getPrime();
	scanf("%d", &t);
	int cas = 0;
	while (t--) {
		int ans = 1;
		scanf("%lld", &n);
		int num = getFactors(n);
		for (int i = 0; i < num; i++) {
			ans *= 2 * factor[i][1] + 1;
		}
		ans = (ans + 1) / 2;
		printf("Case %d: %d\n", ++cas, ans);
	}
	return 0;
}