A191672 Dispersion of A042965 (>1 and congruent to 0 or 1 or 3 mod 4), by antidiagonals.
1, 3, 2, 5, 4, 6, 8, 7, 9, 10, 12, 11, 13, 15, 14, 17, 16, 19, 21, 20, 18, 24, 23, 27, 29, 28, 25, 22, 33, 32, 37, 40, 39, 35, 31, 26, 45, 44, 51, 55, 53, 48, 43, 36, 30, 61, 60, 69, 75, 72, 65, 59, 49, 41, 34, 83, 81, 93, 101, 97, 88, 80, 67, 56, 47, 38
Offset: 1
Examples
Northwest corner: 1....3...5....8....12 2....4...7....11...16 6....9...13...19...27 10...15..21...29...40 14...20..28...39...53
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 = 3; b = 4; c2 = 5; m[n_] := If[Mod[n, 3] == 0, 1, 0]; f[n_] := a*m[n + 2] + b*m[n + 1] + c2*m[n] + 4*Floor[(n - 1)/3] Table[f[n], {n, 1, 30}] (* A042965 *) 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}]] (* A191672 *) Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191672 *)
Comments