|
%I #8 Sep 23 2023 17:42:49
%S 1,1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,8,9,9,10,10,11,11,12,12,12,13,13,
%T 14,14,15,15,16,16,16,17,17,18,18,19,19,20,20,21,21,21,22,22,23,23,24,
%U 24,25,25,25,26,26,27,27,28,28,29,29,30,30,30,31,31,32,32,33,33,34,34,34
%N a(1) = 1, a(n) = ceiling((n-a(a(a(n-1))))/2).
%F a(n)=c*n+O(1) where c is the positive root of x^3+2*x-1=0, c=0.453397651516403767644746538...
%p a:= proc(n) a(n):= ceil((n-a(a(a(n-1))))/2) end: a(1):=1:
%p seq(a(n), n=1..100); # _Alois P. Heinz_, Sep 23 2023
%Y Cf. A005374.
%K nonn
%O 1,4
%A _Benoit Cloitre_, Nov 12 2003
|
