A191728 Dispersion of A047209, (numbers >1 and congruent to 1 or 4 mod 5), by antidiagonals.
1, 4, 2, 11, 6, 3, 29, 16, 9, 5, 74, 41, 24, 14, 7, 186, 104, 61, 36, 19, 8, 466, 261, 154, 91, 49, 21, 10, 1166, 654, 386, 229, 124, 54, 26, 12, 2916, 1636, 966, 574, 311, 136, 66, 31, 13, 7291, 4091, 2416, 1436, 779, 341, 166, 79, 34, 15, 18229, 10229
Offset: 1
Examples
Northwest corner: 1....4....11...29....74 2....6....16...41....104 3....9....24...61....154 5....14...36...91....229 7....19...49...124...311 8....21...54...136...341
Links
- Ivan Neretin, Table of n, a(n) for n = 1..5050
Programs
-
Mathematica
(* Program generates the dispersion array t of the increasing sequence f[n] *) r = 40; r1 = 12; c = 40; c1 = 12; a=4; b=6; m[n_]:=If[Mod[n,2]==0,1,0]; f[n_]:=a*m[n+1]+b*m[n]+5*Floor[(n-1)/2] Table[f[n], {n, 1, 30}] (* A047209 *) 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}]] (* A191728 *) Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191728 *)
Comments