A191724 Dispersion of A047218, (numbers >1 and congruent to 0 or 3 mod 5), by antidiagonals.
1, 3, 2, 8, 5, 4, 20, 13, 10, 6, 50, 33, 25, 15, 7, 125, 83, 63, 38, 18, 9, 313, 208, 158, 95, 45, 23, 11, 783, 520, 395, 238, 113, 58, 28, 12, 1958, 1300, 988, 595, 283, 145, 70, 30, 14, 4895, 3250, 2470, 1488, 708, 363, 175, 75, 35, 16, 12238, 8125, 6175
Offset: 1
Examples
Northwest corner: 1....3....8....20....50 2....5....13...33....83 4....10...25...63....158 6....15...38...95....238 7....18...45...113...283 9....23...58...145...363
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=3; b=5; 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}] (* A047218 *) 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}]] (* A191724 *) Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191724 *)
Comments