cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A161342 Number of "ON" cubic cells at n-th stage in simple 3-dimensional cellular automaton: a(n) = A160428(n)/8.

This page as a plain text file.
%I A161342 #34 Aug 05 2023 13:18:07
%S A161342 0,1,8,15,64,71,120,169,512,519,568,617,960,1009,1352,1695,4096,4103,
%T A161342 4152,4201,4544,4593,4936,5279,7680,7729,8072,8415,10816,11159,13560,
%U A161342 15961,32768,32775,32824,32873,33216,33265,33608,33951,36352,36401,36744,37087,39488
%N A161342 Number of "ON" cubic cells at n-th stage in simple 3-dimensional cellular automaton: a(n) = A160428(n)/8.
%C A161342 First differences are in A161343. - _Omar E. Pol_, May 03 2015
%C A161342 From _Gary W. Adamson_, Aug 30 2016: (Start)
%C A161342 Let M =
%C A161342   1, 0, 0, 0, 0, ...
%C A161342   8, 0, 0, 0, 0, ...
%C A161342   7, 1, 0, 0, 0, ...
%C A161342   0, 8, 0, 0, 0, ...
%C A161342   0, 7, 1, 0, 0, ...
%C A161342   0, 0, 8, 0, 0, ...
%C A161342   0, 0, 7, 1, 0, ...
%C A161342   ...
%C A161342 Then M^k converges to a single nonzero column giving the sequence.
%C A161342 The sequence with offset 1 divided by its aerated variant is (1, 8, 7, 0, 0, 0, ...). (End)
%H A161342 Alois P. Heinz, <a href="/A161342/b161342.txt">Table of n, a(n) for n = 0..16383</a>
%H A161342 N. J. A. Sloane, <a href="/wiki/Catalog_of_Toothpick_and_CA_Sequences_in_OEIS">Catalog of Toothpick and Cellular Automata Sequences in the OEIS</a>
%H A161342 <a href="/index/Ce#cell">Index entries for sequences related to cellular automata</a>
%F A161342 From _Nathaniel Johnston_, Nov 13 2010: (Start)
%F A161342 a(n) = Sum_{k=0..n-1} 7^A000120(k).
%F A161342 a(n) = 1 + 7 * Sum_{k=1..n-1} A151785(k), for n >= 1.
%F A161342 a(2^n) = 2^(3n).
%F A161342 (End)
%F A161342 a(n) = Sum_{k=0..floor(log_2(n))} 7^k*A360189(n-1,k). - _Alois P. Heinz_, Mar 06 2023
%p A161342 b:= proc(n) option remember; `if`(n<0, 0,
%p A161342       b(n-1)+x^add(i, i=Bits[Split](n)))
%p A161342     end:
%p A161342 a:= n-> subs(x=7, b(n-1)):
%p A161342 seq(a(n), n=0..44);  # _Alois P. Heinz_, Mar 06 2023
%t A161342 A161342list[nmax_]:=Join[{0},Accumulate[7^DigitCount[Range[0,nmax-1],2,1]]];A161342list[100] (* _Paolo Xausa_, Aug 05 2023 *)
%Y A161342 Cf. A160410, A160428, A161343, A006046, A130665, A116520, A130667, A116522, A116526, A116525, A360189.
%K A161342 nonn
%O A161342 0,3
%A A161342 _Omar E. Pol_, Jun 14 2009
%E A161342 More terms from _Nathaniel Johnston_, Nov 13 2010