CAStep[rule_,a_]:=映射[rule[[10-#]]&,ListConvolve[{{0,2,0},{2,1,2},},a,2],{2}];
代码=401;阶段=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,阶段]]](*提取相关术语*)
|