A162957 -id:A162957 - cp's OEIS frontend

A162956 a(0) = 0, a(1) = 1; a(2^i + j) = 3a(j) + a(j + 1) for 0 <= j < 2^i.

0, 1, 1, 4, 1, 4, 7, 13, 1, 4, 7, 13, 7, 19, 34, 40, 1, 4, 7, 13, 7, 19, 34, 40, 7, 19, 34, 46, 40, 91, 142, 121, 1, 4, 7, 13, 7, 19, 34, 40, 7, 19, 34, 46, 40, 91, 142, 121, 7, 19, 34, 46, 40, 91, 142, 127, 40, 91, 148, 178, 211, 415, 547, 364, 1, 4, 7, 13, 7, 19, 34, 40, 7, 19, 34, 46

Offset: 0

Gary W. Adamson, Jul 18 2009

nonn

2^n term triangle by rows, analogous to A160552 but multiplier is "3" instead of "2"
Row sums = powers of 5: (1, 5, 25, 125, 625,...).
Rows tend to A162957, obtained by taking (1, 3, 0, 0, 0,...) * A162956.

			The triangle begins:
0;
1;
1, 4;
1, 4, 7, 13;
1, 4, 7, 13, 7, 19, 34, 40;
1, 4, 7, 13, 7, 19, 34, 40, 7, 19, 34, 46, 40, 91, 142, 121;
...
Row 4 = (1, 4, 7, 13, 7, 19, 34, 40): brings down (1, 4, 7, 13) then 7 = 3*1 + 4, 19 = 3*4 + 7, 34 = 3*7 + 13, 40 = 3*13 + 1.

Alois P. Heinz, Table of n, a(n) for n = 0..16383
David Applegate, Omar E. Pol and N. J. A. Sloane, The Toothpick Sequence and Other Sequences from Cellular Automata, Congressus Numerantium, Vol. 206 (2010), 157-191. [There is a typo in Theorem 6: (13) should read u(n) = 4.3^(wt(n-1)-1) for n >= 2.]
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS

Cf. A162957, A162958, A160552.

Cf. A170838-A170852, A170854-A170872.

a:= proc(n) option remember; `if`(n<2, n,
      (j-> 3*a(j)+a(j+1))(n-2^ilog2(n)))
    end:
seq(a(n), n=0..100);  # Alois P. Heinz, Jan 28 2017

row[0] = {0}; row[1] = {1}; row[n_] := row[n] = Join[row[n-1], 3 row[n-1] + Append[Rest[row[n-1]], 1]]; Table[row[n], {n, 0, 7}] // Flatten (* Jean-François Alcover, Mar 13 2017 *)

Follows the same analogous procedure as A160552 but multiplier is 3 instead of 2. (n+1)-th row brings down n-th row and appends to the right and equal number of terms following the rules: from left to right,let a = last term, b = current term, c = next term. Then c = 3*a + b except for the rightmost term = 3*a + 1.

Edited with more terms by N. J. A. Sloane, Jan 02 2010

A162958 Equals A162956 convolved with (1, 3, 3, 3, ...).

Original entry on oeis.org

1, 4, 10, 19, 25, 40, 67, 94, 100, 115, 142, 175, 208, 280, 388, 469, 475, 490, 517, 550, 583, 655, 763, 850, 883, 955, 1069, 1201, 1372, 1696, 2101, 2344, 2350, 2365, 2392, 2425, 2458, 2530, 2638, 2725, 2758, 2830, 2944, 3076, 3247, 3571, 3976, 4225, 4258

Offset: 1

Gary W. Adamson, Jul 18 2009

nonn

Can be considered a toothpick sequence for N=3, following rules analogous to those in A160552 (= special case of "A"), A151548 = special case "B", and the toothpick sequence A139250 (N=2) = special case "C".
To obtain the infinite set of toothpick sequences, (N = 2, 3, 4, ...), replace the multiplier "2" in A160552 with any N, getting a triangle with 2^n terms. Convolve this A sequence with (1, N, 0, 0, 0, ...) = B such that row terms of A triangles converge to B.
Then generalized toothpick sequences (C) = A convolved with (1, N, N, N, ...).
Examples: A160552 * (1, 2, 0, 0, 0,...) = a B-type sequence A151548.
A160552 * (1, 2, 2, 2, 2,...) = the toothpick sequence A139250 for N=2.
A162956 is analogous to A160552 but replaces "2" with the multiplier "3".
Row terms of A162956 tend to A162957 = (1, 3, 0, 0, 0, ...) * A162956.
Toothpick sequence for N = 3 = A162958 = A162956 * (1, 3, 3, 3, ...).
Row sums of "A"-type triangles = powers of (N+2); since row sums of A160552 = (1, 4, 16, 64, ...), while row sums of A162956 = (1, 5, 25, 125, ...).
Is there an illustration of this sequence using toothpicks? - Omar E. Pol, Dec 13 2016

Alois P. Heinz, Table of n, a(n) for n = 1..16384
David Applegate, Omar E. Pol and N. J. A. Sloane, The Toothpick Sequence and Other Sequences from Cellular Automata, Congressus Numerantium, Vol. 206 (2010), 157-191. [There is a typo in Theorem 6: (13) should read u(n) = 4.3^(wt(n-1)-1) for n >= 2.]
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS

Cf. A139250, A152548, A160552, A162956, A163267.

Third diagonal of A163311.

b:= proc(n) option remember; `if`(n<2, n,
      (j-> 3*b(j)+b(j+1))(n-2^ilog2(n)))
    end:
a:= proc(n) option remember;
      `if`(n=0, 0, a(n-1)+2*b(n-1)+b(n))
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, Jan 28 2017

b[n_] := b[n] = If[n<2, n, Function[j, 3*b[j]+b[j+1]][n-2^Floor[Log[2, n]] ]];
a[n_] := a[n] = If[n == 0, 0, a[n-1] + 2*b[n-1] + b[n]];
Array[a, 100] (* Jean-François Alcover, Jun 11 2018, after Alois P. Heinz *)

Clarified definition by Omar E. Pol, Feb 06 2017

cp's OEIS Frontend

A162956 a(0) = 0, a(1) = 1; a(2^i + j) = 3a(j) + a(j + 1) for 0 <= j < 2^i.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions