A160727 -id:A160727 - cp's OEIS frontend

A160410 Number of "ON" cells at n-th stage in simple 2-dimensional cellular automaton (see Comments for precise definition).

0, 4, 16, 28, 64, 76, 112, 148, 256, 268, 304, 340, 448, 484, 592, 700, 1024, 1036, 1072, 1108, 1216, 1252, 1360, 1468, 1792, 1828, 1936, 2044, 2368, 2476, 2800, 3124, 4096, 4108, 4144, 4180, 4288, 4324, 4432, 4540, 4864, 4900, 5008, 5116, 5440, 5548, 5872, 6196

Offset: 0

Omar E. Pol, May 20 2009

On the infinite square grid, we consider cells to be the squares, and we start at round 0 with all cells in the OFF state, so a(0) = 0.
At round 1, we turn ON four cells, forming a square.
The rule for n > 1: A cell in turned ON iff exactly one of its four vertices is a corner vertex of the set of ON cells. So in each generation every exposed vertex turns on three new cells.
Therefore:
At Round 2, we turn ON twelve cells around the square.
At round 3, we turn ON twelve other cells. Three cells around of every corner of the square.
And so on.
For the first differences see the entry A161411.
Shows a fractal behavior similar to the toothpick sequence A139250.
A very similar sequence is A160414, which uses the same rule but with a(1) = 1, not 4.
When n=2^k then the polygon formed by ON cells is a square with side length 2^(k+1).
a(n) is also the area of the figure of A147562 after n generations if A147562 is drawn as overlapping squares. - Omar E. Pol, Nov 08 2009
From Omar E. Pol, Mar 28 2011: (Start)
Also, toothpick sequence starting with four toothpicks centered at (0,0) as a cross.
Rule: Each exposed endpoint of the toothpicks of the old generation must be touched by the endpoints of three toothpicks of new generation. (Note that these three toothpicks looks like a T-toothpick, see A160172.)
The sequence gives the number of toothpicks after n stages. A161411 gives the number of toothpicks added at the n-th stage.
(End)

			From _Omar E. Pol_, Sep 24 2015: (Start)
With the positive terms written as an irregular triangle in which the row lengths are the terms of A011782 the sequence begins:
    4;
   16;
   28,  64;
   76, 112, 148, 256;
  268, 304, 340, 448, 484, 592, 700, 1024;
  ...
Right border gives the elements of A000302 greater than 1.
This triangle T(n,k) shares with the triangle A256534 the terms of the column k, if k is a power of 2, for example, both triangles share the following terms: 4, 16, 28, 64, 76, 112, 256, 268, 304, 448, 1024, etc.
.
Illustration of initial terms, for n = 1..10:
.       _ _ _ _                         _ _ _ _
.      |  _ _  |                       |  _ _  |
.      | |  _|_|_ _ _ _ _ _ _ _ _ _ _ _|_|_  | |
.      | |_|  _ _     _ _     _ _     _ _  |_| |
.      |_ _| |  _|_ _|_  |   |  _|_ _|_  | |_ _|
.          | |_|  _ _  |_|   |_|  _ _  |_| |
.          |   | |  _|_|_ _ _ _|_|_  | |   |
.          |  _| |_|  _ _     _ _  |_| |_  |
.          | | |_ _| |  _|_ _|_  | |_ _| | |
.          | |_ _| | |_|  _ _  |_| | |_ _| |
.          |       |   | |   | |   |       |
.          |  _ _  |  _| |_ _| |_  |  _ _  |
.          | |  _|_| | |_ _ _ _| | |_|_  | |
.          | |_|  _| |_ _|   |_ _| |_  |_| |
.          |   | | |_ _ _ _ _ _ _ _| | |   |
.          |  _| |_ _| |_     _| |_ _| |_  |
.       _ _| | |_ _ _ _| |   | |_ _ _ _| | |_ _
.      |  _| |_ _|   |_ _|   |_ _|   |_ _| |_  |
.      | | |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _| | |
.      | |_ _| |                       | |_ _| |
.      |_ _ _ _|                       |_ _ _ _|
.
After 10 generations there are 304 ON cells, so a(10) = 304.
(End)

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. 31.
Omar E. Pol, Illustration of initial terms (2009)
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS
Index entries for sequences related to cellular automata

Cf. A000079, A000302, A048883, A139250, A147562, A160118, A160412, A160414, A161411, A160717, A160720, A160727, A256530, A256534.

RasterGraphics[state_?MatrixQ,colors_Integer:2,opts___]:=
Graphics[Raster[Reverse[1-state/(colors -1)]],
AspectRatio ->(AspectRatio/.{opts}/.AspectRatio ->Automatic),
Frame ->True, FrameTicks ->None, GridLines ->None];
rule=1340761804646523638425234105559798690663900360577570370705802859623\
705267234688669629039040624964794287326910250673678735142700520276191850\
5902735959769690
Show[GraphicsArray[Map[RasterGraphics,CellularAutomaton[{rule, {2,
{{4,2,1}, {32,16,8}, {256,128,64}}}, {1,1}}, {{{1,1}, {1,1}}, 0}, 9,-10]]]];
ca=CellularAutomaton[{rule,{2,{{4,2,1},{32,16,8},{256,128,64}}},{1,
1}},{{{1,1},{1,1}},0},99,-100];
Table[Total[ca[[i]],2],{i,1,Length[ca]}]
(* John W. Layman, Sep 01 2009; Sep 02 2009 *)
a[n_] := 4*Sum[3^DigitCount[k, 2, 1], {k, 0, n-1}];
Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Nov 17 2017, after N. J. A. Sloane *)

A160410(n)=sum(i=0,n-1,3^norml2(binary(i)))<<2 \\ M. F. Hasler, Dec 04 2012

Equals 4*A130665. This provides an explicit formula for a(n). - N. J. A. Sloane, Jul 13 2009

a(2^k) = (2*(2^k))^2 for k>=0.

Edited by David Applegate and N. J. A. Sloane, Jul 13 2009

A160414 Number of "ON" cells at n-th stage in simple 2-dimensional cellular automaton (same as A160410, but a(1) = 1, not 4).

Original entry on oeis.org

0, 1, 9, 21, 49, 61, 97, 133, 225, 237, 273, 309, 417, 453, 561, 669, 961, 973, 1009, 1045, 1153, 1189, 1297, 1405, 1729, 1765, 1873, 1981, 2305, 2413, 2737, 3061, 3969, 3981, 4017, 4053, 4161, 4197, 4305, 4413, 4737, 4773, 4881, 4989, 5313, 5421, 5745

Offset: 0

Omar E. Pol, May 20 2009

The structure has a fractal behavior similar to the toothpick sequence A139250.
First differences: A161415, where there is an explicit formula for the n-th term.
For the illustration of a(24) = 1729 (the Hardy-Ramanujan number) see the Links section.

			From _Omar E. Pol_, Sep 24 2015: (Start)
With the positive terms written as an irregular triangle in which the row lengths are the terms of A011782 the sequence begins:
1;
9;
21,    49;
61,    97,  133,  225;
237,  273,  309,  417,  453, 561,  669,  961;
...
Right border gives A060867.
This triangle T(n,k) shares with the triangle A256530 the terms of the column k, if k is a power of 2, for example both triangles share the following terms: 1, 9, 21, 49, 61, 97, 225, 237, 273, 417, 961, etc.
.
Illustration of initial terms, for n = 1..10:
.       _ _ _ _                       _ _ _ _
.      |  _ _  |                     |  _ _  |
.      | |  _|_|_ _ _ _ _ _ _ _ _ _ _|_|_  | |
.      | |_|  _ _     _ _   _ _     _ _  |_| |
.      |_ _| |  _|_ _|_  | |  _|_ _|_  | |_ _|
.          | |_|  _ _  |_| |_|  _ _  |_| |
.          |   | |  _|_|_ _ _|_|_  | |   |
.          |  _| |_|  _ _   _ _  |_| |_  |
.          | | |_ _| |  _|_|_  | |_ _| | |
.          | |_ _| | |_|  _  |_| | |_ _| |
.          |  _ _  |  _| |_| |_  |  _ _  |
.          | |  _|_| | |_ _ _| | |_|_  | |
.          | |_|  _| |_ _| |_ _| |_  |_| |
.          |   | | |_ _ _ _ _ _ _| | |   |
.          |  _| |_ _| |_   _| |_ _| |_  |
.       _ _| | |_ _ _ _| | | |_ _ _ _| | |_ _
.      |  _| |_ _|   |_ _| |_ _|   |_ _| |_  |
.      | | |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _| | |
.      | |_ _| |                     | |_ _| |
.      |_ _ _ _|                     |_ _ _ _|
.
After 10 generations there are 273 ON cells, so a(10) = 273.
(End)

Paolo Xausa, Table of n, a(n) for n = 0..10000
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.]
Omar E. Pol, Illustration of initial terms
Omar E. Pol, Illustration of the structure after 24th stage (contains 1729 ON cells)
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS
Index entries for sequences related to cellular automata

Cf. A001235, A011541, A011782, A000225, A060867, A139250, A147562, A160117, A160118, A160410, A160412, A161415, A160720, A160727, A151725, A256530, A256534.

read("transforms") ; isA000079 := proc(n) if type(n,'even') then nops(numtheory[factorset](n)) = 1 ; else false ; fi ; end proc:
A048883 := proc(n) 3^wt(n) ; end proc:
A161415 := proc(n) if n = 1 then 1; elif isA000079(n) then 4*A048883(n-1)-2*n ; else 4*A048883(n-1) ; end if; end proc:
A160414 := proc(n) add( A161415(k),k=1..n) ; end proc: seq(A160414(n),n=0..90) ; # R. J. Mathar, Oct 16 2010

A160414list[nmax_]:=Accumulate[Table[If[n<2,n,4*3^DigitCount[n-1,2,1]-If[IntegerQ[Log2[n]],2n,0]],{n,0,nmax}]];A160414list[100] (* Paolo Xausa, Sep 01 2023, after R. J. Mathar *)

my(s=-1, t(n)=3^norml2(binary(n-1))-if(n==(1<Altug Alkan, Sep 25 2015

a(n) = 1 + 4*A219954(n), n >= 1. - M. F. Hasler, Dec 02 2012

a(2^k) = (2^(k+1) - 1)^2. - Omar E. Pol, Jan 05 2013

Edited by N. J. A. Sloane, Jun 15 2009 and Jul 13 2009

More terms from R. J. Mathar, Oct 16 2010

A161415 First differences of A160414.

Original entry on oeis.org

1, 8, 12, 28, 12, 36, 36, 92, 12, 36, 36, 108, 36, 108, 108, 292, 12, 36, 36, 108, 36, 108, 108, 324, 36, 108, 108, 324, 108, 324, 324, 908, 12, 36, 36, 108, 36, 108, 108, 324, 36, 108, 108, 324, 108, 324, 324, 972, 36, 108, 108, 324, 108, 324, 324, 972, 108, 324, 324

Offset: 1

Omar E. Pol, May 20 2009, Jun 13 2009

nonn

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
Omar E. Pol, Illustration of initial terms [From _Omar E. Pol_, Nov 11 2009]
D. Applegate, Omar E. Pol, N. J. A. Sloane, The toothpick sequence and other sequences from cellular automata, arXiv:1004.3036 [math.CO] [From _R. J. Mathar_, Oct 16 2010]

Cf. A139250, A139251, A160411, A160413, A160414, A160417.
Cf. A160727.
Cf. A048883, A161411, A162349. [From Omar E. Pol, Nov 11 2009]

Contribution from R. J. Mathar, Oct 16 2010: (Start)
isA000079 := proc(n) if type(n,'even') then nops(numtheory[factorset](n)) = 1 ; else false ; fi ; end proc:
A048883 := proc(n) 3^wt(n) ; end proc:
A161415 := proc(n) if n = 1 then 1; elif isA000079(n) then 4*A048883(n-1)-2*n ; else 4*A048883(n-1) ; end if; end proc: seq(A161415(n),n=1..90) ; (End)

a[1] = 1; a[n_] := 4*3^DigitCount[n-1, 2, 1] - If[IntegerQ[Log[2, n]], 2n, 0];
Array[a, 60] (* Jean-François Alcover, Nov 17 2017, after N. J. A. Sloane *)

For n > 1, a(n) = 4*A048883(n-1), except a(n) = 4*A048883(n-1) - 2n if n is a power of 2. - N. J. A. Sloane, Jul 13 2009

More terms from R. J. Mathar, Oct 16 2010

A219954 (A160414(n)-1)/4, n >= 1.

Original entry on oeis.org

0, 2, 5, 12, 15, 24, 33, 56, 59, 68, 77, 104, 113, 140, 167, 240, 243, 252, 261, 288, 297, 324, 351, 432, 441, 468, 495, 576, 603, 684, 765, 992, 995, 1004, 1013, 1040, 1049, 1076, 1103, 1184, 1193, 1220, 1247, 1328, 1355, 1436

Offset: 1

M. F. Hasler, Dec 01 2012

nonn

As is obvious from symmetry, for all n>=1 one has A160414(n) = 1 + a multiple of 4.
Partial sums of 3^A000120(n-1)-(n/2 if n=2^k), n>1.
Numbers whose first differences are in A160727. - Omar E. Pol, Dec 02 2012

Paolo Xausa, Table of n, a(n) for n = 1..10000

Cf. A160727.

A219954list[nmax_]:=Accumulate[Table[If[n==1,0,3^DigitCount[n-1,2,1]-If[IntegerQ[Log2[n]],n/2,0]],{n,nmax}]];A219954list[100] (* Paolo Xausa, Sep 01 2023 *)

my(s=-1,t(n)=3^norml2(binary(n-1))-if(n==(1<

cp's OEIS Frontend

A160410 Number of "ON" cells at n-th stage in simple 2-dimensional cellular automaton (see Comments for precise definition).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A160414 Number of "ON" cells at n-th stage in simple 2-dimensional cellular automaton (same as A160410, but a(1) = 1, not 4).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A161415 First differences of A160414.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A219954 (A160414(n)-1)/4, n >= 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI