OFFSET
1,2
COMMENTS
LINKS
Daniel Suteu, Table of n, a(n) for n = 1..100
EXAMPLE
a(4) = 17 is a term because A124440(17) = 100 = 10^2.
MAPLE
N:= 40000: # for terms <= N
G:= add(numtheory:-mobius(n)*n*x^(2*n)/((1-x^n)*(1-x^(2*n))^2), n=1..N/2):
S:= series(G, x, N+1):
A66840:= [seq(coeff(S, x, j), j=1..N)]:
f:= proc(n) n*numtheory:-phi(n)/2 - A66840[n] end proc:
f(1):= 1: f(2):= 1:
select(t -> issqr(f(t)), [$1..N]);
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Feb 03 2021
EXTENSIONS
More terms from Daniel Suteu, Feb 03 2021
STATUS
approved