OFFSET
1,2
COMMENTS
Subsequence of A036436. - Michel Marcus, Jul 15 2014
LINKS
Jens Kruse Andersen, Table of n, a(n) for n = 1..2989
MAPLE
filter:= proc(n) local t;
t:= numtheory:-tau(n);
issqr(t) and issqr(n-t)
end proc;
select(filter, [$1..10^5]); # Robert Israel, Jul 20 2014
MATHEMATICA
okQ[n_] := With[{t = DivisorSigma[0, n]}, IntegerQ@Sqrt[t] && IntegerQ@Sqrt[n-t]];
Select[Range[10^5], okQ] (* Jean-François Alcover, Feb 08 2023 *)
PROG
(PARI) isok(n) = issquare(numdiv(n)) && issquare(n-numdiv(n)); \\ Michel Marcus, Jul 15 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Muehlfeld, Jul 13 2014
STATUS
approved