

A191656


Dispersion of (2,4,5,7,8,10,...), by antidiagonals.


2



1, 2, 3, 4, 5, 6, 7, 8, 10, 9, 11, 13, 16, 14, 12, 17, 20, 25, 22, 19, 15, 26, 31, 38, 34, 29, 23, 18, 40, 47, 58, 52, 44, 35, 28, 21, 61, 71, 88, 79, 67, 53, 43, 32, 24, 92, 107, 133, 119, 101, 80, 65, 49, 37, 27, 139, 161, 200, 179, 152, 121, 98, 74, 56
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

Row 1: A006999.
For a background discussion of dispersions, see A191426.
...
Each of the sequences (3n, n>0), (3n+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:
...
A191449=dispersion of A008583 (0 mod 3)
A191451=dispersion of A016777 (1 mod 3)
A191450=dispersion of A016789 (2 mod 3)
A191656=dispersion of A001651 (1 or 2 mod 3)
A083044=dispersion of A007494 (0 or 2 mod 3)
A191655=dispersion of A032766 (0 or 1 mod 3)
...
EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):
A191449 has 1st col A001651, all else A008583
A191451 has 1st col A007494, all else A016777
A191450 has 1st col A032766, all else A016789
A191656 has 1st col A008583, all else A001651
A083044 has 1st col A016777, all else A083044
A191655 has 1st col A016789, all else A032766
...
There is a formula for sequences of the type "(a or b mod m)", (as in the Mathematica program below):
If f(n)=(n mod 2), then (a,b,a,b,a,b,...) is given by
a*f(n+1)+b*f(n), so that "(a or b mod m)" is given by
a*f(n+1)+b*f(n)+m*floor((n1)/2)), for n>=1.


LINKS

Table of n, a(n) for n=1..64.


EXAMPLE

Northwest corner:
1...2....4....7....11
3...5....8....13...20
6...10...16...25...38
9...14...22...34...52
12..19...29...44...67


MATHEMATICA

(* Program generates the dispersion array T of the increasing sequence f[n] *)
r = 40; r1 = 12; c = 40; c1 = 12;
a = 2; b = 4; m[n_] := If[Mod[n, 2] == 0, 1, 0];
f[n_] := a*m[n + 1] + b*m[n] + 3*Floor[(n  1)/2]
Table[f[n], {n, 1, 30}] (* A001651: (2+5k, 4+5k, k>=0) *)
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}]] (* A191656 array *)
Flatten[Table[t[k, n  k + 1], {n, 1, c1}, {k, 1, n}]] (* A191656 sequence *)


CROSSREFS

Cf. A001651, A008583, A191426.
Sequence in context: A195083 A073294 A073295 * A266190 A265572 A265556
Adjacent sequences: A191653 A191654 A191655 * A191657 A191658 A191659


KEYWORD

nonn,tabl


AUTHOR

Clark Kimberling, Jun 10 2011


STATUS

approved



