A246035 Number of odd terms in f^n, where f = (1/x+1+x)*(1/y+1+y).
1, 9, 9, 25, 9, 81, 25, 121, 9, 81, 81, 225, 25, 225, 121, 441, 9, 81, 81, 225, 81, 729, 225, 1089, 25, 225, 225, 625, 121, 1089, 441, 1849, 9, 81, 81, 225, 81, 729, 225, 1089, 81, 729, 729, 2025, 225, 2025, 1089, 3969, 25, 225, 225, 625, 225, 2025, 625, 3025, 121, 1089, 1089, 3025, 441, 3969, 1849, 7225, 9, 81, 81, 225, 81, 729, 225
Offset: 0
Keywords
Examples
Here is the neighborhood: [X, X, X] [X, X, X] [X, X, X] which contains a(1) = 9 ON cells. . From Omar E. Pol, Mar 17 2015: (Start) Apart from the initial 1, the sequence can be written also as an irregular tetrahedron T(s,r,k) = A139818(r+2) * a(k), s>=1, 1<=r<=s, 0<=k<=(A011782(s-r)-1) as shown below: .. 9; ... 9; 25; .......... 9, 81; 25; 121; .................... 9, 81, 81, 225; 25, 225; 121; 441; ........................................ 9, 81, 81, 225, 81, 729, 225, 1089; 25, 225, 225, 625; 121, 1089; 441; 1849; ... Note that every row r is equal to A139818(r+2) times the beginning of the sequence itself, thus in 3D every column contains the same number: T(s,r,k) = T(s+1,r,k). (End)
Links
- Alois P. Heinz, Table of n, a(n) for n = 0..8192
- Shalosh B. Ekhad, N. J. A. Sloane, and Doron Zeilberger, A Meta-Algorithm for Creating Fast Algorithms for Counting ON Cells in Odd-Rule Cellular Automata, arXiv:1503.01796 [math.CO], 2015; see also the Accompanying Maple Package.
- Shalosh B. Ekhad, N. J. A. Sloane, and Doron Zeilberger, Odd-Rule Cellular Automata on the Square Grid, arXiv:1503.04249 [math.CO], 2015.
- N. J. A. Sloane, On the No. of ON Cells in Cellular Automata, Video of talk in Doron Zeilberger's Experimental Math Seminar at Rutgers University, Feb. 05 2015: Part 1, Part 2
- N. J. A. Sloane, On the Number of ON Cells in Cellular Automata, arXiv:1503.01168 [math.CO], 2015.
- Index entries for sequences related to cellular automata
Crossrefs
Programs
-
Maple
C:=f->subs({x=1, y=1}, f); # Find number of ON cells in CA for generations 0 thru M defined by rule # that cell is ON iff number of ON cells in nbd at time n-1 was odd # where nbd is defined by a polynomial or Laurent series f(x, y). OddCA:=proc(f, M) global C; local n, a, i, f2, p; f2:=simplify(expand(f)) mod 2; a:=[]; p:=1; for n from 0 to M do a:=[op(a), C(p)]; p:=expand(p*f2) mod 2; od: lprint([seq(a[i], i=1..nops(a))]); end; f:=(1/x+1+x)*(1/y+1+y); OddCA(f, 70);
-
Mathematica
b[0] = 1; b[n_] := b[n] = Expand[b[n - 1]*(x^2 + x + 1)]; a[n_] := Count[CoefficientList[b[n], x], _?OddQ]^2; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Apr 30 2017 *)
Formula
a(n) = A071053(n)^2.
Comments