OFFSET
1,2
COMMENTS
One half times all singly even numbers (A016825) and all powers of them.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
1 = 2^1/2, 2 = 2^2/2, 3 = 6^1/2, 4 = 2^3/2, 5 = 10^1/2, 7 = 14^1/2, ...
MAPLE
filter:= proc(n) local v, m;
if n::odd then true
else v:= padic:-ordp(n, 2); m:= n/2^v; type(simplify(m^(1/(v+1))), integer) fi
end proc:
select(filter, [$1..200]); # Robert Israel, Oct 03 2017
MATHEMATICA
Block[{nn = 120, s}, s = Range[2, 2 nn, 4]; Union@ Flatten@ Map[#^Range@ Floor@ Log[#, 2 nn]/2 &, s]] (* Michael De Vlieger, Oct 02 2017 *)
PROG
(PARI) {a(n) = my(cnt, m, k); if( n<1, return(0), while( cnt<n, m++; k=2*m; ispower(k, , &k); if(k%4==2, cnt++)); m)};
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael Somos, Oct 02 2017
STATUS
approved
