# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a354901
Showing 1-1 of 1

%I A354901 #22 Apr 21 2024 22:11:47
%S A354901 1,2,3,4,6,7,5,8,12,10,14,13,11,15,9,16,24,20,28,22,30,18,26,25,21,29,
%T A354901 23,31,19,27,17,32,48,40,56,36,52,44,60,42,58,38,54,46,62,34,50,49,41,
%U A354901 57,37,53,45,61,43,59,39,55,47,63,35,51,33,64,96,80,112
%N A354901 a(n) = (b(2n) - 1)/2 - n for n > 0. To get b(n) start with A = n and then for i = 0..f(n) apply A := A + 2^i*T(A, f(n) - i) where T(n, k) = floor(n/2^k) mod 2 and f(n) = A000523(n).
%C A354901 Permutation of the natural numbers such that infinite iteration of g(n, m) converges to A059893(n) where g(n, m) = g(A004754(n),m-1)/2 for n > 0, m > 0 with g(n, 0) = a(n).
%C A354901 Subsequences from a(2^m) to a(2^(m+1) - 1) for m >= 0 contain all numbers k such that 2^m <= k < 2^(m+1).
%H A354901 <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>
%o A354901 (PARI) b(n)=my(L=logint(n,2),A=n); for(i=0,L, A+=2^i*bittest(A,L-i)); A;
%o A354901 a(n)=(b(2*n) - 1)/2 - n
%Y A354901 Cf. A000523, A004754, A059893.
%K A354901 nonn
%O A354901 1,2
%A A354901 _Mikhail Kurkov_, Jun 11 2022 [verification needed]

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE