|
CAStep[rule_，a_]：=映射[rule[[10-#]]&，ListConvolve[{{0，2，0}，{2，1，2}，}，a，2]，{2}]；
代码=141；阶段=128；
规则=整数位数[code，2，10]；
g=2*级+1；（*网格最大尺寸*）
a=PadLeft[｛｛1｝｝，｛g，g｝，0，Floor[｛g，g｝/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}]；
映射[Function[Apply[Plus，Flatten[#1]]]，ca]（*每个阶段计数ON个细胞*）
| 