cp's OEIS Frontend

A153153 Permutation of natural numbers: A059893-conjugate of A003188.

%I A153153 #14 Feb 26 2020 03:58:32
%S A153153 0,1,3,2,5,6,7,4,9,10,15,12,13,14,11,8,17,18,23,20,29,30,27,24,25,26,
%T A153153 31,28,21,22,19,16,33,34,39,36,45,46,43,40,57,58,63,60,53,54,51,48,49,
%U A153153 50,55,52,61,62,59,56,41,42,47,44,37,38,35,32,65,66,71,68,77,78,75,72
%N A153153 Permutation of natural numbers: A059893-conjugate of A003188.
%H A153153 Antti Karttunen, <a href="/A153153/b153153.txt">Table of n, a(n) for n = 0..2047</a>
%H A153153 <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>
%F A153153 a(n) = A065190(A231550(n)). - _Yosu Yurramendi_, Jan 15 2020
%F A153153 a(1) = 1, a(2^(m+1)+2*k) = 2*a(2^(m+1)-1-k), a(2^(m+1)+2*k+1) = 2*a(2^m+k), m >= 0, 0 <= k < 2^m. - _Yosu Yurramendi_, Jan 25 2020
%o A153153 (R)
%o A153153 a <- 1
%o A153153 maxlevel <- 5 # by choice
%o A153153 #
%o A153153 for(m in 0:maxlevel) for(k in 0:(2^m-1)){
%o A153153   a[2^(m+1)+2*k  ] <- 2*a[2^(m+1)-1-k] + 1
%o A153153   a[2^(m+1)+2*k+1] <- 2*a[2^m+k]
%o A153153 }
%o A153153 a <- c(0,a)
%o A153153 # _Yosu Yurramendi_, Jan 25 2020
%Y A153153 Inverse: A153154. a(n) = A059893(A003188(A059893(n))).
%K A153153 nonn,base
%O A153153 0,3
%A A153153 _Antti Karttunen_, Dec 20 2008