OFFSET
1,2
LINKS
Miles Englezou, Table of n, a(n) for n = 1..7248
EXAMPLE
3 is a term since (2^2-1)*(2^2)^0 = 3.
8 is a term since (3^2-1)*(3^2)^0 = 8.
12 is a term since (2^2-1)*(2^2)^1 = 12.
144 is a term since 3 and 48 are terms, gcd(3, 48) = 1, and 3*48 = 144.
PROG
(PARI) lst(lim)= my(S = Set()); for(n = 1, lim, my(v = n^2); my(f = factor(n)); for(i = 1, #f~, v *= (1 - 1/f[i, 1]^2)); v = numerator(v); if(v <= lim, S = setunion(S, Set([v])))); vecsort(Vec(S))
CROSSREFS
KEYWORD
nonn
AUTHOR
Miles Englezou, Apr 15 2026
STATUS
approved
