用 Mathematica 搜索生命游戏中的静物

发表于 | 更新于

果壳网已死。临死前元胞自动机小组还回光返照了一段时间——虽然真正活跃的人大概不超过五个。

以下的代码就算是我给果壳网以及元胞自动机小组烧的一点纸钱吧。

a = Merge[
   k[#[[;; 2]], 
       CellularAutomaton["GameOfLife", #][[2, 2]]] -> #[[3]] & /@ 
    Tuples[Tuples[{0, 1}, 3], 3], Identity];
g[x_, y_, z_] := 
  FindPath[Flatten@
     MapIndexed[
      Outer[v[#2 - 1, Most[#1]] -> v[#2, Rest[#1]] &, ##, 1] &, 
      Lookup[a, 
       MapThread[
        k, {First[
          Partition[{x, y}, {2, 3}, {2, 1}, {{1, -2}, {-1, 2}}, 0]], 
         z}], {}]], v[0, {0, 0}], v[Length[z], {0, 0}], Length[z] + 2,
     All][[;; , 2 ;;, 2, 1]];
With[{w = 16, h = 16},
 o = Table[0, w + 2];
 i = g[o, o, o];
 f[l : {___, x_, y_}] := 
  If[Length@l == h + 2, 
   If[y == o && MemberQ[i, x], l, Missing["NotFound"]], 
   Catch[Do[
     If[MissingQ[#], , Throw[#]] &[f[Append[l, z]]], {z, 
      RandomSample[g[x, y, y]]}]; Missing["NotFound"]]];
 SparseArray[f[{o, RandomChoice[i]}]]]

返回的结果是一个SparseArray，可以用ArrayPlot来转换成图片。这是搜索16*16的静物。把With[{w = 16, h = 16},里的16换成别的数字可以搜索别的大小的静物，但不要把w改得太大。如果一次搜索的时间太久，可能是运气不好，可以用Alt+.来中断计算。

下面是5次搜索的结果：

x = 16, y = 16, rule = B3/S23
ob2o2bob2o5bo$2ob4obobo2b3o$9bo2bo$5o5b2ob3o$o2bobob3o2bo2bo$5bobo3bo2bo$4bo3b2ob3o$o2bob3o2bo$4obo2bo2bob3o$5bo2b2ob2o2bo$3obo4bo4bo$o2b2obo2bob2ob2o$bo3b2o3bobo2bo$2b2o10bo$obob3o4b3o$2o4bo4bo!
x = 16, y = 16, rule = B3/S23
bo11bobo$b3o6bobob2o$4bobob3obo$b2obob2o4bo$o2bo5b2obobo$2o8bo2bobo$6bobo2bobobo$3bobob5ob2o$b3obo$o4b3obob4o$b4o3b3o4bo$3bo2bo5bo2bo$5b2ob7o$b4o3bo$bo3bobo3bob2o$4b2ob2o2b2obo!
x = 16, y = 16, rule = B3/S23
2bobob2obobobo$bob2obob2ob2obo$bo4bo3bo3bo$2bob3o3bo3b2o$ob2o5b2obobo$2o3b2o4b2obo$3b2ob4o3bo$2obo5bob2o$2o3b2o2bobo$3bobobobobo$3obo2b2obob2o$o3bobo3bo3bo$4b2obobob2obo$2b2o3bobo3bo$bo2bob2o2bobo$2b2o5b2ob2o!
x = 16, y = 16, rule = B3/S23
2b2obo4bo3b2o$bob2obob3o2bobo$o5b2o3b3o$b5o3b2o4bo$2bo2bob3o2b4o$o5bo4b2o$2o6b3o3b2o$bo4bobo2b2obo$bobo2b2o3bo2bo$2ob3o6bob2o$bo4b2o3b2o2bo$bob2obo5bobo$2o2bo3bobo3b2o$2bobo2b2ob2o2bo$2bobobo2bo2bobo$3b2ob2o3b2ob2o!
x = 16, y = 16, rule = B3/S23
bobo2bob2o2bobo$ob6ob5obo$o14bo$b4ob5ob3o$4b2o4bobo$obo4bo2bo2bobo$2ob8o3b2o$3bo$2obo2bo2b2o$bob4o3bo2bobo$o9bobob2o$3obo4b2obo$3b2o3bo4b3o$2o5bobob2o2bo$ob3o3b2obo2bo$4bo7b2o!