$\gcd(a,b)=g$ とし,$\dfrac{a}{g}=x,\dfrac{b}{g}=y$ とすると,$x,y$ は正整数で,かつ $\mathrm{lcm}(a,b)=gxy$ である.このとき,$an^3+b=\mathrm{lcm}(a,b)n^2+\gcd(a,b)n$ は
$$gxn^3+gy=gxyn^2+gn$$
すなわち
$$g(n-y)(xn^2-1)=0$$
と変形でき,これを $n$ についての方程式と見たときその解は
$$n=y,\pm\sqrt{\frac{1}{x}}$$
である.$\pm\sqrt{\dfrac{1}{x}}$ は $x=1$ のとき,つまり $b$ が $a$ で割り切れるときのみ整数となる.これを踏まえて次のように場合分けをする.
以上より,
$$\sum_{a=1}^{100}\sum_{b=1}^{100}f(a,b)=2×100+3×382+1×9518=\mathbf{10864}$$
である.
この問題を解いた人はこんな問題も解いています