A191707 Dispersion of A016873, (numbers >1 and congruent to 1, 2, 3, or 4 mod 5), by antidiagonals.

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

Author

Clark Kimberling, Jun 12 2011

Keywords

• nonn
• tabl

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)).

Examples

			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
		

Links

• Ivan Neretin, Table of n, a(n) for n = 1..5050 (first 100 antidiagonals, flattened)

Crossrefs

Cf. A047201, A008587, A191702, A191426.

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=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  *)