login
A191707
Dispersion of A016873, (numbers >1 and congruent to 1, 2, 3, or 4 mod 5), by antidiagonals.
10
1, 2, 5, 3, 7, 10, 4, 9, 13, 15, 6, 12, 17, 19, 20, 8, 16, 22, 24, 26, 25, 11, 21, 28, 31, 33, 32, 30, 14, 27, 36, 39, 42, 41, 38, 35, 18, 34, 46, 49, 53, 52, 48, 44, 40, 23, 43, 58, 62, 67, 66, 61, 56, 51, 45, 29, 54, 73, 78, 84, 83, 77, 71, 64, 57, 50, 37
OFFSET
1,2
COMMENTS
For a background discussion of dispersions and their fractal sequences, see A191426. For dispersions of congruence sequences mod 3 or mod 4, see A191655, A191663, A191667.
...
Each of the sequences (5n, n>1), (5n+1, n>1), (5n+2, n>=0), (5n+3, n>=0), (5n+4, n>=0), generates a dispersion. Each complement (beginning with its first term >1) also generates a dispersion. The ten sequences and dispersions are listed here:
...
A191702=dispersion of A008587 (5k, k>=1)
A191703=dispersion of A016861 (5k+1, k>=1)
A191704=dispersion of A016873 (5k+2, k>=0)
A191705=dispersion of A016885 (5k+3, k>=0)
A191706=dispersion of A016897 (5k+4, k>=0)
A191707=dispersion of A047201 (1, 2, 3, 4 mod 5 and >1)
A191708=dispersion of A047202 (0, 2, 3, 4 mod 5 and >1)
A191709=dispersion of A047207 (0, 1, 3, 4 mod 5 and >1)
A191710=dispersion of A032763 (0, 1, 2, 4 mod 5 and >1)
A191711=dispersion of A001068 (0, 1, 2, 3 mod 5 and >1)
...
EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):
A191702 has 1st col A047201, all else A008587
A191703 has 1st col A047202, all else A016861
A191704 has 1st col A047207, all else A016873
A191705 has 1st col A032763, all else A016885
A191706 has 1st col A001068, all else A016897
A191707 has 1st col A008587, all else A047201
A191708 has 1st col A042968, all else A047203
A191709 has 1st col A042968, all else A047207
A191710 has 1st col A042968, all else A032763
A191711 has 1st col A042968, all else A001068
...
Regarding the dispersions A191670-A191673, there is a formula for sequences of the type "(a or b or c or d mod m)", (as in the relevant Mathematica programs):
...
If f(n)=(n mod 3), then (a,b,c,d,a,b,c,d,a,b,c,d,...) is given by a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n); so that for n>=1, "(a, b, c, d mod m)" is given by
a*f(n+3)+b*f(n+2)+c*f(n+1)+d*f(n)+m*floor((n-1)/4)).
LINKS
Ivan Neretin, Table of n, a(n) for n = 1..5050 (first 100 antidiagonals, flattened)
EXAMPLE
Northwest corner:
1....2....3....4....6
5....7....9....12...16
10...13...17...22...28
15...19...24...31...39
20...26...33...42...53
25...32...41...52...66
MATHEMATICA
(* Program generates the dispersion array T of the increasing sequence f[n] *)
r = 40; r1 = 12; c = 40; c1 = 12;
a=2; b=3; c2=4; d=6; m[n_]:=If[Mod[n, 4]==0, 1, 0];
f[n_]:=a*m[n+3]+b*m[n+2]+c2*m[n+1]+d*m[n]+5*Floor[(n-1)/4]
Table[f[n], {n, 1, 30}] (* A047201 *)
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}]] (* A191707 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191707 *)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Clark Kimberling, Jun 12 2011
STATUS
approved