ランダムなメモランダム

雑草のような雑念と雑考

Ulam SpiralのMathematica code

k番目の数値のxy座標計算
coord[k_]:={0,0}+Apply[Plus,Transpose[{Flatten[Table[{Table[(-1)^(j-1),{i,1,j}],Table[0,{i,1,j}]},{j,1,300}]],Flatten[Table[{Table[0,{i,1,j}],Table[(-1)^(j-1),{i,1,j}]},{j,1,300}]]}]1;; k ]


グラフィック作成130×130 一般数はシアン、素数はシロ
cl = {};
Do[AppendTo[
cl, {If[PrimeQ[k], White, Cyan], Disk[coord[k], 1/2]}], {k, 0,
Length[Transpose[{Flatten[
Table[{Table[(-1)^(j - 1), {i, 1, j}], Table[0, {i, 1, j}]}, {j,
1, 130}]],
Flatten[Table[{Table[0, {i, 1, j}],
Table[(-1)^(j - 1), {i, 1, j}]}, {j, 1, 130}]]}]]}]