|
CAStep[rule_，a_]：=映射[rule[[10-#]]&，ListConvolve[{{0，2，0}，{2，1，2}，}，a，2]，{2}]；
代码=589；阶段=128；
规则=整数位数[code，2，10]；
g=2*级+1；（*网格最大尺寸*）
a=PadLeft[{{1}}，{g，g}，0，Floor[{g，c}/2]]；（*电网上的初始ON电池*）
ca=a；
ca=表[ca=CAStep[规则，ca]，{n，1，阶段+1}]；
PrependTo[ca，a]；
（*修剪整个网格以反映每个阶段一个单元格的增长*）
k=（长度[ca[[1]]]+1）/2；
ca=表[表[部分[ca[[n]][[j]]，范围[k+1-n，k-1+n]]，{j，k+1-n，k-1+n}]，{n，1，k}]；
on=映射[Function[Apply[Plus，Flatten[#1]]]，ca]（*在每个阶段计数on单元格*）
第[on，2^Range[0，Log[2，stages]]]部分（*提取相关术语*）
| 