OFFSET
1,4
COMMENTS
a(n)=1 for n prime = 2,3,5,17,257,...;
a(n)=2 if n = 4 or 9, or n is in A067466(k).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
a[n_] := Sum[Sum[If[(g = GCD[i, j]) + i*j/g == n, 1, 0], {j, 1, i}], {i, 1, n}]; Array[a, 100] (* Amiram Eldar, Jun 06 2022 *)
PROG
(PARI) a(n) = sum(i=1, n, sum(j=1, i, n == gcd(i, j)+lcm(i, j))); \\ Michel Marcus, Feb 16 2021
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, May 16 2002
STATUS
approved