OFFSET
1,2
COMMENTS
From Robert Israel, Jul 13 2016: (Start)
All terms except 1 are even.
Even terms are the numbers of the form n = (a+b)^2/8 such that for some d >= 1,
2^d <= n < 2^(d+1) and 2^(d+1)-1 = a*b. (End)
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
F:= proc(d) local v, R, X;
v:= 2^(d+1)-1;
R:= select(t-> t^2 < v, numtheory:-divisors(v));
op(select(t -> t >= (v+1)/2 and t < v+1, map(t -> (t+ v/t)^2/8, R)));
end proc:
sort(convert({1, seq(F(i), i=1..50)}, list)); # Robert Israel, Jul 13 2016
MATHEMATICA
Select[Range[10^6], Times @ Boole@ {IntegerQ@ Sqrt@ FromDigits[RotateRight@ #, 2], IntegerQ@ Sqrt@ FromDigits[RotateLeft@ #, 2]} &@ IntegerDigits[#, 2] == 1 &] (* Michael De Vlieger, Jun 29 2016 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Alex Ratushnyak, Jun 27 2016
STATUS
approved