A191665 Dispersion of A042963 (numbers >1, congruent to 1 or 2 mod 4), by antidiagonals.
1, 2, 3, 5, 6, 4, 10, 13, 9, 7, 21, 26, 18, 14, 8, 42, 53, 37, 29, 17, 11, 85, 106, 74, 58, 34, 22, 12, 170, 213, 149, 117, 69, 45, 25, 15, 341, 426, 298, 234, 138, 90, 50, 30, 16, 682, 853, 597, 469, 277, 181, 101, 61, 33, 19, 1365, 1706, 1194, 938, 554
Offset: 1
Examples
Northwest corner: 1...2...5....10...21 3...6...13...26...53 4...9...18...37...74 7...14..29...58...117 8...17..34...69...138
Links
- Ivan Neretin, Table of n, a(n) for n = 1..5050 (first 100 antidiagonals, flattened)
Programs
-
Mathematica
(* Program generates the dispersion array T of the increasing sequence f[n] *) r = 40; r1 = 12; c = 40; c1 = 12; a = 2; b = 5; m[n_] := If[Mod[n, 2] == 0, 1, 0]; f[n_] := a*m[n + 1] + b*m[n] + 4*Floor[(n - 1)/2] Table[f[n], {n, 1, 30}] (* A042963: (2+4k,5+4k) *) 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}]] (* A191665 *) Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191665 *)
Comments