OFFSET
1,2
MAPLE
sqrf:= n-> mul(i[1]^irem(i[2], 2), i=ifactors(n)[2]):
f:=proc(n) global sqrf; option remember;
if n = 1 then 1 elif n=2 then 4 else sqrf(f(n-1)+f(n-2)); fi; end;
[seq(f(n), n=1..140)];
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, May 05 2016
STATUS
approved