A191443 Dispersion of the sequence ([n*sqrt(3)+1]), where [ ]=floor, read by antidiagonals.
1, 2, 3, 4, 6, 5, 7, 11, 9, 8, 13, 20, 16, 14, 10, 23, 35, 28, 25, 18, 12, 40, 61, 49, 44, 32, 21, 15, 70, 106, 85, 77, 56, 37, 26, 17, 122, 184, 148, 134, 97, 65, 46, 30, 19, 212, 319, 257, 233, 169, 113, 80, 52, 33, 22, 368, 553, 446, 404, 293, 196, 139
Offset: 1
Examples
Northwest corner: 1....2....4....7....13 3....6....11...20...35 5....9....16...28...49 8....14...25...44...77 10...18...32...56...97
Programs
-
Mathematica
(* Program generates the dispersion array T of increasing sequence f[n] *) r=40; r1=12; c=40; c1=12; x = Sqrt[3]; f[n_] := Floor[n*x+1] (* complement of column 1 *) mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]] rows = {NestList[f, 1, c]}; Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}]; t[i_, j_] := rows[[i, j]]; TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191443 array *) Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191443 sequence *) (* Program by Peter J. C. Moses, Jun 01 2011 *)
Comments