cp's OEIS Frontend

A110898 a(n) = 1 + a(n - a(n-1)) (n - a(n-1)) mod 3 + (n - a(n-2)) mod 4, with a(0) = 0, a(1) = 1.

%I A110898 #15 Mar 31 2024 08:48:03
%S A110898 0,1,3,3,2,6,1,3,8,3,5,2,8,8,4,4,9,8,6,11,6,7,5,7,10,7,11,9,12,7,8,6,
%T A110898 10,9,7,14,6,9,6,12,14,11,6,9,15,9,12,15,10,15,13,5,17,7,15,14,8,18,
%U A110898 16,10,12,18,16,15,15,14,20,15,6,18,17,17,17,7,16,9,15,15,19,13,16,13,21,17
%N A110898 a(n) = 1 + a(n - a(n-1)) (n - a(n-1)) mod 3 + (n - a(n-2)) mod 4, with a(0) = 0, a(1) = 1.
%H A110898 G. C. Greubel, <a href="/A110898/b110898.txt">Table of n, a(n) for n = 0..10000</a>
%t A110898 a[n_]:= a[n]= If[n<2, n, 1 +a[n-a[n-1]] +Mod[n-a[n-1],3] -Mod[n-a[n-2],3]];
%t A110898 Table[a[n], {n,0,100}]
%o A110898 (Magma) [n le 2 select n-1 else 1 + Self(n-Self(n-1)) - ((n+2-Self(n-1)) mod 3) + ((n+3-Self(n-2)) mod 4) : n in [1..101]]; // _G. C. Greubel_, Mar 30 2024
%o A110898 (SageMath)
%o A110898 @CachedFunction
%o A110898 def a(n): # a = A110898
%o A110898     if n<2: return n
%o A110898     else: return 1 +a(n-a(n-1)) -(n-a(n-1))%3 +(n-a(n-2))%4
%o A110898 [a(n) for n in range(101)] # _G. C. Greubel_, Mar 30 2024
%Y A110898 Cf. A110897.
%K A110898 nonn,easy,less
%O A110898 0,3
%A A110898 _Roger L. Bagula_, Sep 20 2005
%E A110898 Edited by _G. C. Greubel_, Mar 30 2024