OFFSET
0,2
COMMENTS
Start with cell (0,0) active; at each succeeding stage the cells that share exactly one edge with an active cell change their state.
The pattern has 4-fold symmetry; sequence just counts cells in one quadrant.
This is not the CA discussed by Singmaster in the reference given in A079314. That was an error based on my misreading of the paper. - N. J. A. Sloane, Aug 05 2009
REFERENCES
D. Singmaster, On the cellular automaton of Ulam and Warburton, M500 Magazine of the Open University, #195 (December 2003), pp. 2-7.
LINKS
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.]
D. Singmaster, On the cellular automaton of Ulam and Warburton, 2003 [Cached copy, included with permission]
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS
PROG
(PARI) M=matrix(101, 101); M[1, 1]=1; for(s=1, 100, c=[]; a=M[1, 1]; for(x=2, 100, for(y=2, 100, a+=M[x, y]; if(M[x-1, y]+M[x+1, y]+M[x, y-1]+M[x, y+1]==1, c=concat(c, [[x, y]]) )); a+=M[x, 1]+M[1, x]; if(M[x, 2]==0 && M[x-1, 1]+M[x+1, 1]==1, c=concat(c, [[x, 1]]) ); if(M[2, x]==0 && M[1, x-1]+M[1, x+1]==1, c=concat(c, [[1, x]]) )); print1(a, ", "); for(i=1, length(c), M[c[i][1], c[i][2]]=1-M[c[i][1], c[i][2]]) ) \\ Max Alekseyev, Feb 02 2007
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Feb 12 2003
EXTENSIONS
More terms from Max Alekseyev, Feb 02 2007
Edited by N. J. A. Sloane, Aug 05 2009
STATUS
approved