LightOJ1236 （分解质因数）

题意： 给定$n$，求出$lcm(i,j)=n$的数对个数。其中$1\le i\le j\le n$。
数据范围：$1\le n\le 10^{14}$

题解： 对于一个数$x$，分解质因数后得到：$x=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}...p_k^{{\alpha }_k}$
那么对于$lcm(A,B)=p_1^{max(a_1,b_1)}{p}_2^{max(a_2,b_2)}...p_k^{max(a_k,b_k)}$

• 对于$n$以内的数，$a_i,b_i\in [0,{\alpha }_i]$，所以对于$a_i,b_i$都有${\alpha }_i+1$种选择，不过注意当两者同取最大幂时，该幂只能被取一次，所以共$2{\alpha }_i+1$次。
• 其次，由于$A\le B$，而由于最大幂只能取一次，所以$A=n,B=n$只出现一次，其余均出现了$2$次。故最后答案是$\lfloor \frac{({\prod }_{i=1}^{k}2{\alpha }_i+1)}{2}\rfloor +1$

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N = 1e7 + 10;
int pri[N / 10], cnt;
bool st[N];

void xs(int n) {
	st[0] = st[1] = true;
	for(int i = 2; i <= n; ++i) {
		if(!st[i]) pri[++cnt] = i;
		for(int j = 1; j <= cnt && 1ll * i * pri[j] <= n; ++j) {
			st[i * pri[j]] = true;
			if(i % pri[j] == 0) break;
		}
	}
}

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

int lcm(int a, int b) {
	return a / gcd(a, b) * b;
}

long long pairsFormLCM( int n ) {
    long long res = 0;
    for( int i = 1; i <= n; i++ )
        for( int j = i; j <= n; j++ )
           if( lcm(i, j) == n ) printf("#(%d, %d\n)", i, j), res++; // lcm means least common multiple
    return res;
}

int main()
{
	xs(N - 1);
	int T; scanf("%d", &T);
	for(int k = 1; k <= T; ++k) {
		ll n; scanf("%lld", &n);
		
		ll res = 1;
		for(int j = 1; j <= cnt && 1ll * pri[j] * pri[j] <= n; ++j) {
			if(n % pri[j] == 0) {
				int temp = 0;
				while(n % pri[j] == 0) n /= pri[j], ++temp;
				res *= (2 * temp + 1);
			}
		}
		if(n > 1) res *= (2 * 1 + 1);
		res = res / 2 + 1;
		printf("Case %d: %lld\n", k, res);
	}
}