%I #32 Mar 16 2021 10:22:56
%S 1,1,2,3,1,4,5,2,6,7,3,8,9,1,10,11,4,12,13,5,14,15,2,16,17,6,18,19,7,
%T 20,21,3,22,23,8,24,25,9,26,27,1,28,29,10,30,31,11,32,33,4,34,35,12,
%U 36,37,13,38,39,5,40,41,14,42,43,15,44,45,2,46,47,16,48,49,17,50,51,6,52,53,18,54,55,19,56,57,7
%N Row index of n in A191450: a(3n) = 2n, a(3n+1) = 2n+1, a(3n+2) = a(n+1).
%C a(n) gives the row index of n in square array A191450, or equally, the column index of n in A254051.
%H Antti Karttunen, <a href="/A253887/b253887.txt">Table of n, a(n) for n = 1..10000</a>
%F a(3n) = 2n, a(3n+1) = 2n+1, a(3n+2) = a(n+1).
%F a(n) = A126760(2n-1).
%F a(n) = A249746(A003602(A064216(n))). - _Antti Karttunen_, Feb 04 2015
%o (Scheme, with memoization-macro definec)
%o (definec (A253887 n) (if (= 2 (modulo n 3)) (A253887 (/ (+ 1 n) 3)) (inv_for_A032766 n)))
%o (define (inv_for_A032766 n) (+ (* 2 (floor->exact (/ n 3))) (modulo n 3)))
%o (Python)
%o def a(n):
%o if n%3==0: return 2*n//3
%o elif n%3==1: return 2*(n - 1)//3 + 1
%o else: return a((n - 2)//3 + 1)
%o print([a(n) for n in range(1, 101)]) # _Indranil Ghosh_, Jun 06 2017
%Y Odd bisection of A126760.
%Y Cf. A254046 (the corresponding column index).
%Y Cf. A003602, A032766, A064216, A253786, A253893, A249746.
%Y Cf. A191450, A254051, A254047, A254052, A254104.
%K nonn
%O 1,3
%A _Antti Karttunen_, Jan 22 2015
|