|
CAStep[rule_，a_]：=映射[rule[[10-#]]&，ListConvolve[{{0，2，0}，{2，1，2}，}，a，2]，{2}]；
代码=966；阶段=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}]；
表[FromDigits[部分[ca[i]][[i]]，范围[i，2*i-1]]，10]，{i，1，阶段-1}]
| 