A161342 -id:A161342 - cp's OEIS frontend

A360189 Triangle T(n,k), n>=0, 0<=k<=floor(log_2(n+1)), read by rows: T(n,k) = number of nonnegative integers <= n having binary weight k.

1, 1, 1, 1, 2, 1, 2, 1, 1, 3, 1, 1, 3, 2, 1, 3, 3, 1, 3, 3, 1, 1, 4, 3, 1, 1, 4, 4, 1, 1, 4, 5, 1, 1, 4, 5, 2, 1, 4, 6, 2, 1, 4, 6, 3, 1, 4, 6, 4, 1, 4, 6, 4, 1, 1, 5, 6, 4, 1, 1, 5, 7, 4, 1, 1, 5, 8, 4, 1, 1, 5, 8, 5, 1, 1, 5, 9, 5, 1, 1, 5, 9, 6, 1, 1, 5, 9, 7, 1

Offset: 0

Alois P. Heinz, Mar 04 2023

T(n,k) is defined for all n >= 0 and k >= 0. Terms that are not in the triangle are zero.

			T(6,2) = 3: 3, 5, 6, or in binary: 11_2, 101_2, 110_2.
T(15,3) = 4: 7, 11, 13, 14, or in binary: 111_2, 1011_2, 1101_2, 1110_2.
Triangle T(n,k) begins:
  1;
  1, 1;
  1, 2;
  1, 2, 1;
  1, 3, 1;
  1, 3, 2;
  1, 3, 3;
  1, 3, 3, 1;
  1, 4, 3, 1;
  1, 4, 4, 1;
  1, 4, 5, 1;
  1, 4, 5, 2;
  1, 4, 6, 2;
  1, 4, 6, 3;
  1, 4, 6, 4;
  1, 4, 6, 4, 1;
  ...

Alois P. Heinz, Rows n = 0..2^12-1, flattened
Peter J. Taylor, Closed form for Sum_{k=0..n} [wt(k) = m] where wt(n) is the binary weight of n, answer to question on MathOverflow (2024).
Wikipedia, Iverson bracket

Columns k=0-2 give: A000012, A029837(n+1) = A113473(n) for n>0, A340068(n+1).
Last elements of rows give A090996(n+1).
Cf. A000027, A000120, A000225, A000788, A006046, A007318, A116520, A116522, A116525, A116526, A130665, A130667, A161342, A231500, A231501, A231502, A272020, A361257.

b:= proc(n) option remember; `if`(n<0, 0,
      b(n-1)+x^add(i, i=Bits[Split](n)))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n)):
seq(T(n), n=0..23);

T(n,k) = my(v1); v1 = Vecrev(binary(n+1)); v1 = Vecrev(select(x->(x>0),v1,1)); sum(j=0, min(k,#v1-1), binomial(v1[j+1]-1,k-j)) \\ Mikhail Kurkov, Nov 27 2024

T(n,k) = T(n-1,k) + [A000120(n) = k] where [] is the Iverson bracket and T(n,k) = 0 for n<0.
T(2^n-1,k) = A007318(n,k) = binomial(n,k).
T(n,floor(log_2(n+1))) = A090996(n+1).
Sum_{k>=0} T(n,k) = n+1.
Sum_{k>=0} k * T(n,k) = A000788(n).
Sum_{k>=0} k^2 * T(n,k) = A231500(n).
Sum_{k>=0} k^3 * T(n,k) = A231501(n).
Sum_{k>=0} k^4 * T(n,k) = A231502(n).
Sum_{k>=0} 2^k * T(n,k) = A006046(n+1).
Sum_{k>=0} 3^k * T(n,k) = A130665(n).
Sum_{k>=0} 4^k * T(n,k) = A116520(n+1).
Sum_{k>=0} 5^k * T(n,k) = A130667(n+1).
Sum_{k>=0} 6^k * T(n,k) = A116522(n+1).
Sum_{k>=0} 7^k * T(n,k) = A161342(n+1).
Sum_{k>=0} 8^k * T(n,k) = A116526(n+1).
Sum_{k>=0} 10^k * T(n,k) = A116525(n+1).
Sum_{k>=0} n^k * T(n,k) = A361257(n).
T(n,k) = Sum_{j=0..min(k, A000120(n+1)-1)} binomial(A272020(n+1,j+1)-1,k-j) for n >= 0, k >= 0 (see Peter J. Taylor link). - Mikhail Kurkov, Nov 27 2024

A116522 a(0)=1, a(1)=1, a(n)=7a(n/2) for n=2,4,6,..., a(n)=6a((n-1)/2)+a((n+1)/2) for n=3,5,7,....

Original entry on oeis.org

0, 1, 7, 13, 49, 55, 91, 127, 343, 349, 385, 421, 637, 673, 889, 1105, 2401, 2407, 2443, 2479, 2695, 2731, 2947, 3163, 4459, 4495, 4711, 4927, 6223, 6439, 7735, 9031, 16807, 16813, 16849, 16885, 17101, 17137, 17353, 17569, 18865, 18901, 19117, 19333

Offset: 0

Roger L. Bagula, Mar 15 2006

nonn

A 7-divide version of A084230.
The Harborth: f(2^k) = 3^k suggests that a family of sequences of the form: f(2^k) = prime(n)^k.
From Gary W. Adamson, Aug 27 2016: (Start)
Let M = the production matrix below. Then lim_{k->infinity} M^k generates the sequence with offset 1 by extracting the left-shifted vector.
  1, 0, 0, 0, 0, ...
  7, 0, 0, 0, 0, ...
  6, 1, 0, 0, 0, ...
  0, 7, 0, 0, 0, ...
  0, 6, 1, 0, 0, ...
  0, 0, 7, 0, 0, ...
  0, 0, 6, 1, 0, ...
  ...
The sequence divided by its aerated variant is (1, 7, 6, 0, 0, 0, ...). (End)

Alois P. Heinz, Table of n, a(n) for n = 0..16383 (first 2501 terms from G. C. Greubel)
H. Harborth, Number of Odd Binomial Coefficients, Proc. Amer. Math. Soc. 62, 19-22, 1977.
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Identities and periodic oscillations of divide-and-conquer recurrences splitting at half, arXiv:2210.10968 [cs.DS], 2022, pp. 27, 33.
Eric Weisstein's World of Mathematics, Stolarsky-Harborth Constant

Cf. A000120, A006046, A077465, A130665, A116520, A130667, A161342, A360189.

a:=proc(n) if n=0 then 0 elif n=1 then 1 elif n mod 2 = 0 then 7*a(n/2) else 6*a((n-1)/2)+a((n+1)/2) fi end: seq(a(n),n=0..47);
# second Maple program:
b:= proc(n) option remember; `if`(n<0, 0,
      b(n-1)+x^add(i, i=Bits[Split](n)))
    end:
a:= n-> subs(x=6, b(n-1)):
seq(a(n), n=0..44);  # Alois P. Heinz, Mar 06 2023

b[0] := 0; b[1] := 1; b[n_?EvenQ] := b[n] = 7*b[n/2]; b[n_?OddQ] := b[n] = 6*b[(n - 1)/2] + b[(n + 1)/2]; a = Table[b[n], {n, 1, 25}]

G.f.: (r(x) * r(x^2) * r(x^4) * r(x^8) * ...), where r(x) = (1 + 7x + 6x^2).
a(n) = Sum_{k=0..n-1} 6^wt(k), where wt = A000120. - Mike Warburton, Mar 14 2019
a(n) = Sum_{k=0..floor(log_2(n))} 6^k*A360189(n-1,k). - Alois P. Heinz, Mar 06 2023

Edited by N. J. A. Sloane, Apr 16 2005

A160428 Number of ON cells at n-th stage of three-dimensional version of the cellular automaton A160410, using cubes.

Original entry on oeis.org

0, 8, 64, 120, 512, 568, 960, 1352, 4096, 4152, 4544, 4936, 7680, 8072, 10816, 13560, 32768, 32824, 33216, 33608, 36352, 36744, 39488, 42232, 61440, 61832, 64576, 67320, 86528, 89272, 108480, 127688, 262144, 262200, 262592, 262984, 265728, 266120, 268864, 271608

Offset: 0

Omar E. Pol, Jun 01 2009

nonn

Michael De Vlieger, Table of n, a(n) for n = 0..1000
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.]
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Identities and periodic oscillations of divide-and-conquer recurrences splitting at half, arXiv:2210.10968 [cs.DS], 2022, p. 33.
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS
Index entries for sequences related to cellular automata

Cf. A139250, A160119, 160379, A160410, A160429, A161340, A161342.

a[n_] := 8*Sum[7^DigitCount[k, 2, 1], {k, 0, n - 1}]; Array[a, 40, 0] (* Michael De Vlieger, Nov 01 2022 *)

a(n) = 8 * Sum_{k=0..n-1} 7^A000120(k)

a(n) = 8 + 56 * Sum_{k=1..n-1} A151785(k) for n >= 1

Formulas and more terms from Nathaniel Johnston, Nov 13 2010

More terms from Michael De Vlieger, Nov 01 2022

A116525 a(0)=1, a(1)=1, a(n) = 11a(n/2) for even n, and a(n) = 10a((n-1)/2) + a((n+1)/2) for odd n >= 3.

Original entry on oeis.org

0, 1, 11, 21, 121, 131, 231, 331, 1331, 1341, 1441, 1541, 2541, 2641, 3641, 4641, 14641, 14651, 14751, 14851, 15851, 15951, 16951, 17951, 27951, 28051, 29051, 30051, 40051, 41051, 51051, 61051, 161051, 161061, 161161, 161261, 162261, 162361, 163361, 164361

Offset: 0

Roger L. Bagula, Mar 15 2006

nonn

From Gary W. Adamson, Aug 30 2016: (Start)
Let M =
    1,  0,  0, 0, 0, ...
   11,  0,  0, 0, 0, ...
   10,  1,  0, 0, 0, ...
    0, 11,  0, 0, 0, ...
    0, 10,  1, 0, 0, ...
    0,  0, 11, 0, 0, ...
    0,  0, 10, 1, 0, ...
  ...
Then lim_{k->infinity} M^k converges to a single nonzero column giving the sequence.
The sequence divided by its aerated variant is (1, 11, 10, 0, 0, 0, ...). (End)

Alois P. Heinz, Table of n, a(n) for n = 0..16383 (first 2501 terms from G. C. Greubel)
H. Harborth, Number of Odd Binomial Coefficients, Proc. Amer. Math. Soc. 62, 19-22, 1977.
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Identities and periodic oscillations of divide-and-conquer recurrences splitting at half, arXiv:2210.10968 [cs.DS], 2022, pp. 27, 33.
Eric Weisstein's World of Mathematics, Stolarsky-Harborth Constant

Cf. A000120, A006046, A077465, A084230.

Cf. A130665, A116520, A130667, A116522, A161342, A116526, A360189.

a:=proc(n) if n=0 then 0 elif n=1 then 1 elif n mod 2 = 0 then 11*a(n/2) else 10*a((n-1)/2)+a((n+1)/2) fi end: seq(a(n),n=0..42);

b[0] := 0; b[1] := 1; b[n_?EvenQ] := b[n] = 11*b[n/2]; b[n_?OddQ] := b[n] = 10*b[(n - 1)/2] + b[(n + 1)/2]; a = Table[b[n], {n, 1, 25}]

Let r(x) = (1 + 11x + 10x^2).  The sequence is r(x) * r(x^2) * r(x^4) * r(x^8) * ... - Gary W. Adamson, Aug 30 2016
a(n) = Sum_{k=0..n-1} 10^wt(k), where wt = A000120. - Mike Warburton, Mar 14 2019
a(n) = Sum_{k=0..floor(log_2(n))} 10^k*A360189(n-1,k). - Alois P. Heinz, Mar 06 2023

Edited by N. J. A. Sloane, Apr 16 2005

A116526 a(0)=1, a(1)=1, a(n) = 9a(n/2) for even n >= 2, and a(n) = 8a((n-1)/2) + a((n+1)/2) for odd n >= 3.

Original entry on oeis.org

0, 1, 9, 17, 81, 89, 153, 217, 729, 737, 801, 865, 1377, 1441, 1953, 2465, 6561, 6569, 6633, 6697, 7209, 7273, 7785, 8297, 12393, 12457, 12969, 13481, 17577, 18089, 22185, 26281, 59049, 59057, 59121, 59185, 59697, 59761, 60273, 60785, 64881, 64945, 65457, 65969

Offset: 0

Roger L. Bagula, Mar 15 2006

nonn

A 9-divide version of A084230.
The interest this one has is in the prime form of even odd 2^n+1, 2^n.
From Gary W. Adamson, Aug 30 2016: (Start)
Let M =
  1, 0, 0, 0, 0, ...
  9, 0, 0, 0, 0, ...
  8, 1, 0, 0, 0, ...
  0, 9, 0, 0, 0, ...
  0, 8, 1, 0, 0, ...
  0, 0, 9, 0, 0, ...
  0, 0, 8, 1, 0, ...
  ...
Then M^k converges to a single nonzero column giving the sequence.
The sequence divided by its aerated variant is (1, 9, 8, 0, 0, 0, ...). (End)

Alois P. Heinz, Table of n, a(n) for n = 0..16383 (first 2501 terms from G. C. Greubel)
H. Harborth, Number of Odd Binomial Coefficients, Proc. Amer. Math. Soc. 62, 19-22, 1977.
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Identities and periodic oscillations of divide-and-conquer recurrences splitting at half, arXiv:2210.10968 [cs.DS], 2022, pp. 27, 33.
Eric Weisstein's World of Mathematics, Stolarsky-Harborth Constant

Cf. A000120, A006046, A077465, A130665, A116520, A130667, A116522, A161342, A116525, A360189.

a:=proc(n) if n=0 then 0 elif n=1 then 1 elif n mod 2 = 0 then 9*a(n/2) else 8*a((n-1)/2)+a((n+1)/2) fi end: seq(a(n),n=0..45);

b[0] := 0; b[1] := 1; b[n_?EvenQ] := b[n] = 9*b[n/2]; b[n_?OddQ] := b[n] = 8*b[(n - 1)/2] + b[(n + 1)/2]; a = Table[b[n], {n, 1, 25}]

a(n) = Sum_{k=0..n-1} 8^wt(k), where wt = A000120. - Mike Warburton, Mar 14 2019

a(n) = Sum_{k=0..floor(log_2(n))} 8^k*A360189(n-1,k). - Alois P. Heinz, Mar 06 2023

Edited by N. J. A. Sloane, Apr 16 2006

A161430 Fractal triangle whose virtual skeleton is a polyedge as the Y-toothpick structure of A160120.

Original entry on oeis.org

0, 9, 33, 57, 123, 147

Offset: 0

Omar E. Pol, Jun 17 2009

Number of ON states after n generations of cellular automaton based on the infinite triangular grid.

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, A160120, A161341, A161342, A161343.

A256141 Square array read by antidiagonals upwards: T(n,k), n>=0, k>=0, in which row n lists the partial sums of the n-th row of the square array of A256140.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 4, 1, 1, 5, 7, 9, 5, 1, 1, 6, 9, 16, 11, 6, 1, 1, 7, 11, 25, 19, 15, 7, 1, 1, 8, 13, 36, 29, 28, 19, 8, 1, 1, 9, 15, 49, 41, 45, 37, 27, 9, 1, 1, 10, 17, 64, 55, 66, 61, 64, 29, 10, 1, 1, 11, 19, 81, 71, 91, 91, 125, 67, 33, 11, 1, 1, 12, 21, 100, 89, 120, 127, 216, 129, 76, 37, 12, 1

Offset: 0

Omar E. Pol, Mar 16 2015

Questions:
Is also A130667 a row of this square array?
Is also A116522 a row of this square array?
Is also A116526 a row of this square array?
Is also A116525 a row of this square array?
Is also A116524 a row of this square array?

			The corner of the square array with the first 15 terms of the first 12 rows looks like this:
--------------------------------------------------------------------------
A000012: 1, 1, 1,  1,  1,  1,  1,   1,   1,   1,   1,   1,   1,   1,   1
A000027: 1, 2, 3,  4,  5,  6,  7,   8,   9,  10,  11,  12,  13,  14,  15
A006046: 1, 3, 5,  9, 11, 15, 19,  27,  29,  33,  37,  45,  49,  57,  65
A130665: 1, 4, 7, 16, 19, 28, 37,  64,  67,  76,  85, 112, 121, 148, 175
A116520: 1, 5, 9, 25, 29, 45, 61, 125, 129, 145, 161, 225, 241, 305, 369
A130667? 1, 6,11, 36, 41, 66, 91, 216, 221, 246, 271, 396, 421, 546, 671
A116522? 1, 7,13, 49, 55, 91,127, 343, 349, 385, 421, 637, 673, 889,1105
A161342: 1, 8,15, 64, 71,120,169, 512, 519, 568, 617, 960,1009,1352,1695
A116526? 1, 9,17, 81, 89,153,217, 729, 737, 801, 865,1377,1441,1953,2465
.......: 1,10,19,100,109,190,271,1000,1009,1090,1171,1900,1981,2710,3439
A116525? 1,11,21,121,131,231,331,1331,1341,1441,1541,2541,2641,3641,4641
.......: 1,12,23,144,155,276,397,1728,1739,1860,1981,3312,3422,4764,6095

Index entries for sequences related to cellular automata

Cf. A000120, A255740, A255741, A256140.
First five rows are A000012, A000027, A006046, A130665, A116520. Row 7 is A161342.
First eight columns are A000012, A000027, A005408, A000290, A028387, A000384, A003215, A000578. Column 9 is A081437. Column 11 is A015237. Columns 13-15 are A005915, A005917, A000583.

A161343 a(n) = 7^A000120(n).

Original entry on oeis.org

1, 7, 7, 49, 7, 49, 49, 343, 7, 49, 49, 343, 49, 343, 343, 2401, 7, 49, 49, 343, 49, 343, 343, 2401, 49, 343, 343, 2401, 343, 2401, 2401, 16807, 7, 49, 49, 343, 49, 343, 343, 2401, 49, 343, 343, 2401, 343, 2401, 2401, 16807, 49, 343, 343, 2401, 343, 2401, 2401, 16807, 343, 2401, 2401, 16807, 2401, 16807, 16807, 117649

Offset: 0

Omar E. Pol, Jun 14 2009

nonn

Also first differences of A161342.
From Omar E. Pol, May 03 2015: (Start)
It appears that when A151785 is regarded as a triangle in which the row lengths are the powers of 2, this is what the rows converge to.
Also this is also a row of the square array A256140.
(End)

			From _Omar E. Pol_, May 03 2015: (Start)
Also, written as an irregular triangle in which the row lengths are the terms of A011782, the sequence begins:
1;
7;
7, 49;
7, 49, 49, 343;
7, 49, 49, 343, 49, 343, 343, 2401;
7, 49, 49, 343, 49, 343, 343, 2401, 49, 343, 343, 2401, 343, 2401, 2401, 16807;
...
Row sums give A055274.
Right border gives A000420.
(End)

Cf. A000420, A011782, A055274, A160410, A151785, A161411, A160428, A160429, A161342, A256140, A256141.

Rows of the array A256140 are: A000007, A000012, A001316, A048883, A102376, A256135, A256136, this sequence.

a(n) = 7^hammingweight(n); \\ Omar E. Pol, May 03 2015

a(n) = A000420(A000120(n)). - Omar E. Pol, May 03 2015

G.f.: Product_{k>=0} (1 + 7*x^(2^k)). - Ilya Gutkovskiy, Mar 02 2017

More terms from Sean A. Irvine, Mar 08 2011
New name from Omar E. Pol, May 03 2015
a(52)-a(63) from Omar E. Pol, May 16 2015

cp's OEIS Frontend

A360189 Triangle T(n,k), n>=0, 0<=k<=floor(log_2(n+1)), read by rows: T(n,k) = number of nonnegative integers <= n having binary weight k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

PARI

Formula

A116522 a(0)=1, a(1)=1, a(n)=7*a(n/2) for n=2,4,6,..., a(n)=6*a((n-1)/2)+a((n+1)/2) for n=3,5,7,....

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A160428 Number of ON cells at n-th stage of three-dimensional version of the cellular automaton A160410, using cubes.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A116525 a(0)=1, a(1)=1, a(n) = 11*a(n/2) for even n, and a(n) = 10*a((n-1)/2) + a((n+1)/2) for odd n >= 3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A116526 a(0)=1, a(1)=1, a(n) = 9*a(n/2) for even n >= 2, and a(n) = 8*a((n-1)/2) + a((n+1)/2) for odd n >= 3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A161430 Fractal triangle whose virtual skeleton is a polyedge as the Y-toothpick structure of A160120.

Views

Author

Keywords

Comments

Links

Crossrefs

A256141 Square array read by antidiagonals upwards: T(n,k), n>=0, k>=0, in which row n lists the partial sums of the n-th row of the square array of A256140.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

A161343 a(n) = 7^A000120(n).

Views

Author

Keywords

Comments

Examples

A116522 a(0)=1, a(1)=1, a(n)=7a(n/2) for n=2,4,6,..., a(n)=6a((n-1)/2)+a((n+1)/2) for n=3,5,7,....

A116525 a(0)=1, a(1)=1, a(n) = 11a(n/2) for even n, and a(n) = 10a((n-1)/2) + a((n+1)/2) for odd n >= 3.

A116526 a(0)=1, a(1)=1, a(n) = 9a(n/2) for even n >= 2, and a(n) = 8a((n-1)/2) + a((n+1)/2) for odd n >= 3.