|
%I #31 Sep 04 2022 14:46:00
%S 1,3,5,17,29,139,249,64570209,129140169,34315253252541,68630377364913,
%T 1044297913696328396542704032390321722034449074468444246791788357605,
%U 2088595827392656793085408064780643444068898148936888424953199350297
%N Let M(1)=0 and for n>1, B(n)=(M(ceiling(n/2))+M(floor(n/2))+2)/2, M(n)=3^B(n)+M(floor(n/2))+1. This sequence gives B(n).
%C The largest power of 3 in M(n) = A230640(n).
%H Max Alekseyev, <a href="/A230639/a230639.txt">Table n, expression for a(n) for n=2..100</a>
%H Max A. Alekseyev and N. J. A. Sloane, <a href="https://arxiv.org/abs/2112.14365">On Kaprekar's Junction Numbers</a>, arXiv:2112.14365, 2021; Journal of Combinatorics and Number Theory, 2022 (to appear).
%H <a href="/index/Coi#Colombian">Index entries for Colombian or self numbers and related sequences</a>
%p f:=proc(n) option remember; local B, M;
%p if n<=1 then RETURN([0, 0]);
%p else
%p B:=(f(ceil(n/2))[2] + f(floor(n/2))[2] + 2)/2;
%p M:=3^B+f(floor(n/2))[2]+1; RETURN([B, M]); fi;
%p end proc;
%p [seq(f(n)[1], n=1..9)];
%Y Cf. A230093, A230640 (for M(n)).
%Y Related base-3 sequences: A053735, A134451, A230641, A230642, A230643, A230853, A230854, A230855, A230856, A230639, A230640, A010063 (trajectory of 1)
%K nonn,base
%O 2,2
%A _N. J. A. Sloane_, Oct 31 2013
%E Terms a(10) onward from _Max Alekseyev_, Nov 02 2013
| 