|
%I #14 Jan 06 2021 18:35:24
%S 1,2,3,2,2,2,2,4,4,5,5,6,6,7,7,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,8,8,8,
%T 8,9,9,9,9,10,10,10,10,11,11,11,11,12,12,12,12,13,13,13,13,14,14,14,
%U 14,15,15,15,15,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8
%N a(n) = floor(n/a(floor(n/2))); a(1) = 1.
%H Robert Israel, <a href="/A097051/b097051.txt">Table of n, a(n) for n = 1..10000</a>
%F If floor(log_2(n))=2k+1, then a(n) = floor(n/2^k). If floor(log_2(n))=2k, then a(n) = 2^k.
%e a(50)=floor(50/a(25))
%e ..... a(25)=floor(25/a(12))
%e ........... a(12)=floor(12/a(6))
%e ................. a(6)=floor(6/a(3))
%e ...................... a(3)=floor(3/a(1))
%e ........................... a(1)=1
%e ...................... a(3)=floor(3/a(1))=floor(3/1)=3
%e ................. a(6)=floor(6/a(3))=floor(6/3)=2
%e ........... a(12)=floor(12/a(6))=floor(12/2)=6
%e ..... a(25)=floor(25/a(12))=floor(25/6)=4
%e a(50)=floor(50/a(25))=floor(50/4)=12.
%p f:= proc(n) option remember; floor(n/procname(floor(n/2))) end proc:
%p f(1):= 1:
%p map(f, [$1..200]); # _Robert Israel_, Jan 06 2021
%t a[1] = 1; a[n_] := a[n] = Floor[n/a[Floor[n/2]]]; Table[ a[n], {n, 94}]
%Y Cf. A096036, A097053.
%K nonn
%O 1,2
%A _Reinhard Zumkeller_ and _Robert G. Wilson v_, Jul 21 2004
%E Formula added by _Max Alekseyev_, Mar 02 2011
|
