OFFSET
1,1
COMMENTS
MATHEMATICA
z = 300;
A[n_, k_] := Module[{t, a, b}, t = (1 + Sqrt[5])/2;
a = Floor[n*(t + 1) + 1 + t/2]; b = Round[a*t]; ({b, a} . MatrixPower[{{1, 1}, {1, 0}}, k])[[2]]];
ts = Table[A[n, k], {n, 0, z - 1}, {k, 0, z - 1}]; (* A035506, Stolarsky array *)
W[n_, k_] := Fibonacci[k + 1] Floor[n*GoldenRatio] + (n - 1) Fibonacci[k];
tw = Table[W[n, k], {n, 1, z}, {k, 1, z}]; (* A035513, Wythoff array *)
diff = tw - ts;
u = Table[diff[[n]][[2]], {n, 1, z}]
Flatten[Position[u, 0]] (* A387349 *)
Flatten[Position[u, 1]] (* A387350 *)
Flatten[Position[u, -1]] (* A387351 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Aug 27 2025
STATUS
approved
