\(a={p_1} ^ {a_1} *{p_1} ^ {a_1} *..........*{p_n} ^ {a_n}\)
\(b={p_1} ^ {b_1} *{p_1} ^ {b_1} *..........*{p_n} ^ {b_n}\)
\(lcm(a,b)={p_1} ^ {max(a_1,b_1)} *{p_2} ^ {max(a_2,b_2)} *..........*{p_n} ^ {max(a_n,b_n)}=n\)
假定\(a<=b\)
所以对n进行质因数分解，计算出每个质因数的指数部分，比如其中一部分\({p_n}^k\)则说明\(max(a_n,b_n)=k\)，那么如果\(a_n=k\),那么\(b_n\)有\(k+1\)种取值方法，同理如果\(b_n=k\),那么\(a_n\)有\(k+1\)种取值方法,那么对于这个质因数我们有\(2*(k+1)-1\)种取值方法,一开始\(ans=1\),对于每个质因数乘以其贡献,那么除了\(a=b=n\)的情况，其他都计算了两次,由于最后我们要的是\((a<=b)\)的方案数，那么\(ans=ans/2+1\)即可
来源：https://www.cnblogs.com/graytido/p/11897008.html