|
CAStep[rule_，a_]：=映射[rule[[10-#]]&，ListConvolve[{{0，2，0}，{2，1，2}，}，a，2]，{2}]；
代码=659；阶段=128；
规则=整数位数[code，2，10]；
g=2*级+1；（*网格最大尺寸*）
a=PadLeft[{{1}}，{g，g}，0，Floor[{g，c}/2]]；（*电网上的初始ON电池*）
ca=a；
ca=表[ca=CAStep[rule，ca]，{n，1，stages+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^范围[0，日志[2，阶段]]]（*提取相关术语*）
| 