LightOJ 1236 Pairs froming LCM 素因子分解

题意

给定一个数$n$，问在小于等于$n$的数中，有多少数对(a,b)满足$lcm(a,b)=n$，其中$a<b$

思路

$lcm=n$的数对中的两个数一定是n的因子
我们从因子的角度来看LCM：
对n进行素因子分解
$n=p_1^{a_1}p_2^{a_2}p_3^{a_3}...p_k^{a_k}$
那么它的因子$m_i$一定也可以表示为
$m_i=p_1^{b_{i1}}p_2^{b_{i2}}p_3^{b_{i3}}...p_k^{b_{ik}} ,0=<b_{ik}<=a_{ik}$
$lcm(m_i,m_j)=p_1^{max(b_{i1},b_{j1})}p_2^{max(b_{i2},b_{j2})}...p_k^{max(b_{ik},b_{jk})}$
$lcm(m_i,m_j)=n$
$max(b_{i1},b_{j1})=a_1,...,max(b_{ik},b_{jk})=a_k$
如果$b_{i1}$取$a_1$，那么$b_{j1}$可以取$0-a_1$，共有$a_1+1$种取法
如果$b_{i1}$取小于$a_1$的数，那么$b_{j1}$只能取$a_1$，共有$a_1$种取法，所以对于第一个素因子总共有$2a_1+1$种取法
所以对于所有的素因子来说总共有$\sum(2a_k+1)$种取法
由于本题目中要求数对是有序的，所以最终的结果要加一除以二