%I #8 Jul 07 2023 05:45:38
%S 9,12,16,21,28,37,50,66,88,118,157,209,279,372,496,661,882,1176,1568,
%T 2090,2787,3716,4955,6606,8808,11744,15659,20879,27838,37118,49490,
%U 65987,87983,117310,156414,208552,278069,370759,494345,659127
%N a(n) = 9 + floor( Sum_{j=1..n-1} a(j)/3 ).
%H G. C. Greubel, <a href="/A120154/b120154.txt">Table of n, a(n) for n = 1..1000</a>
%t A120154[n_]:= A120154[n]= 9 +Quotient[Sum[A120154[k], {k,n-1}], 3];
%t Table[A120154[n], {n,60}] (* _G. C. Greubel_, Jun 20 2023 *)
%o (Magma)
%o function f(n,a,b)
%o t:=0;
%o for k in [1..n-1] do
%o t+:= a+Floor((b+t)/3);
%o end for;
%o return t;
%o end function;
%o g:= func< n,a,b | f(n+1,a,b)-f(n,a,b) >;
%o A120154:= func< n | g(n,9,0) >;
%o [A120154(n): n in [1..60]]; // _G. C. Greubel_, Jun 20 2023
%o (SageMath)
%o @CachedFunction
%o def A120154(n): return 9 + (sum(A120154(k) for k in range(1,n)))//3
%o [A120154(n) for n in range(1,61)] # _G. C. Greubel_, Jun 20 2023
%Y Cf. A072493, A073941, A112088.
%K nonn
%O 1,1
%A _Graeme McRae_, Jun 10 2006
|