A191667 Dispersion of A016813 (4k+1, k>1), by antidiagonals.

1, 5, 2, 21, 9, 3, 85, 37, 13, 4, 341, 149, 53, 17, 6, 1365, 597, 213, 69, 25, 7, 5461, 2389, 853, 277, 101, 29, 8, 21845, 9557, 3413, 1109, 405, 117, 33, 10, 87381, 38229, 13653, 4437, 1621, 469, 133, 41, 11, 349525, 152917, 54613, 17749, 6485, 1877, 533

Offset

1,2

Comments

For a background discussion of dispersions, see A191426.

...

Each of the sequences (4n, n>2), (4n+1, n>0), (3n+2, n>=0), generates a dispersion. Each complement (beginning with its first term >1) also generates a dispersion. The six sequences and dispersions are listed here:

...

A191452=dispersion of A008586 (4k, k>=1)

A191667=dispersion of A016813 (4k+1, k>=1)

A191668=dispersion of A016825 (4k+2, k>=0)

A191669=dispersion of A004767 (4k+3, k>=0)

A191670=dispersion of A042968 (1 or 2 or 3 mod 4 and >=2)

A191671=dispersion of A004772 (0 or 1 or 3 mod 4 and >=2)

A191672=dispersion of A004773 (0 or 1 or 2 mod 4 and >=2)

A191673=dispersion of A004773 (0 or 2 or 3 mod 4 and >=2)

...

EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):

A191452 has 1st col A042968, all else A008486

A191667 has 1st col A004772, all else A016813

A191668 has 1st col A042965, all else A016825

A191669 has 1st col A004773, all else A004767

A191670 has 1st col A008486, all else A042968

A191671 has 1st col A016813, all else A004772

A191672 has 1st col A016825, all else A042965

A191673 has 1st col A004767, all else A004773

...

Regarding the dispersions A191670-A191673, there is a formula for sequences of the type "(a or b or c mod m)", (as in the Mathematica program below):

If f(n)=(n mod 3), then (a,b,c,a,b,c,a,b,c,...) is given by a*f(n+2)+b*f(n+1)+c*f(n), so that "(a or b or c mod m)" is given by a*f(n+2)+b*f(n+1)+c*f(n)+m*floor((n-1)/3)), for n>=1.

Links

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

Example

Northwest corner:
1....5....21....85....341
2....9....37....149...597
3....13...53....213...853
4....17...69....277...1109
6....25...101...405...1621

Mathematica

(* Program generates the dispersion array T of the increasing sequence f[n] *)
r = 40; r1 = 12;  c = 40; c1 = 12;
f[n_] := 4*n+1
Table[f[n], {n, 1, 30}]  (* A016813 *)
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}]] (* A191667 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191667  *)

Crossrefs

Row 1: A002450.

Cf. A004772, A016813, A191671, A191426.

Sequence in context: A276533 A303685 A189746 * A130329 A096035 A074769

Adjacent sequences: A191664 A191665 A191666 * A191668 A191669 A191670

Keyword

nonn,tabl

Author

Clark Kimberling, Jun 11 2011

Status

approved