A110897 -id:A110897 - 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.

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, 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, 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

Offset: 0

Roger L. Bagula, Sep 20 2005

G. C. Greubel, Table of n, a(n) for n = 0..10000

Cf. A110897.

[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

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]];
Table[a[n], {n,0,100}]

@CachedFunction
def a(n): # a = A110898
    if n<2: return n
    else: return 1 +a(n-a(n-1)) -(n-a(n-1))%3 +(n-a(n-2))%4
[a(n) for n in range(101)] # G. C. Greubel, Mar 30 2024

Edited by G. C. Greubel, Mar 30 2024

A206012 Modular recursion: a(0)=a(1)=a(2)=a(3)=1, thereafter: a(n) equals a(n - 2) + a(n - 3) when n = 0 mod 5, a(n - 1) + a(n - 3) when n = 1 mod 5, a(n - 1) + a(n - 2) when n = 2 mod 5, a(n - 1) + a(n - 4) when n = 3 mod 5, and a(n - 1) + a(n - 2) + a(n - 3) otherwise.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 3, 5, 8, 16, 13, 21, 34, 50, 105, 84, 134, 218, 323, 675, 541, 864, 1405, 2080, 4349, 3485, 5565, 9050, 13399, 28014, 22449, 35848, 58297, 86311, 180456, 144608, 230919, 375527, 555983, 1162429, 931510, 1487493, 2419003, 3581432, 7487928

Offset: 0

Roger L. Bagula, Mar 19 2012

This sequence was inspired by the work of Paul Curtz on three part sequences. I did a three part version of this that gave a generating polynomial and got even more variance by adding two more modulo sequences.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,6,0,0,0,0,3,0,0,0,0,-1).

Cf. A194880, A195240, A185138, A206568, A110897, A110898.

m:=40; R:=PowerSeriesRing(Integers(), m); Coefficients(R!((x^15-x^13+x^12-2*x^10-2*x^9+2*x^8-x^7-3*x^6-4*x^5+3*x^4+x^3+x^2+x+1)/(x^15-3*x^10-6*x^5+1))); // Bruno Berselli, Mar 20 2012

a[0] = 1; a[1] = 1; a[2] = 1; a[3] = 1;a[n_Integer] := a[n]=If[Mod[n, 5] == 0, a[n - 2] + a[n - 3], If[Mod[n, 5] == 1, a[n - 1] + a[n - 3], If[Mod[n, 5] == 2, a[n - 1] + a[n - 2], If[Mod[n, 5] == 3, a[n - 1] + a[n - 4], a[n - 1] + a[n - 2] + a[n - 3]]]]];b = Table[a[n], {n, 0, 50}];(* FindSequenceFunction gives*);Table[c[n] = b[[n]], {n, 1, 16}];c[n_Integer] := c[n] = -c[-15 + n] + c[-10 + n] + 6 c[-5 + n];d = Table[c[n], {n, 1, Length[b]}]
CoefficientList[Series[(x^15-x^13+x^12-2*x^10-2*x^9+2*x^8-x^7-3*x^6-4*x^5+3*x^4+x^3+x^2+x+1)/(x^15-3*x^10-6*x^5+1),{x,0,1001}],x] (* Vincenzo Librandi, Apr 01 2012 *)

Vec((x^15-x^13+x^12-2*x^10-2*x^9+2*x^8-x^7-3*x^6-4*x^5+3*x^4+x^3+x^2+x+1)/(x^15-3*x^10-6*x^5+1)+O(x^99)) \\ Charles R Greathouse IV, Mar 19 2012

G.f.: (x^15 - x^13 + x^12 - 2x^10 - 2x^9 + 2x^8 - x^7 - 3x^6 - 4x^5 + 3x^4 + x^3 + x^2 + x + 1) / (x^15 - 3x^10 - 6x^5 + 1). - Alois P. Heinz, Mar 19 2012

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.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

SageMath

Extensions