

A191663


Dispersion of A042948 (numbers >3, congruent to 0 or 1 mod 4), by antidiagonals.


34



1, 4, 2, 9, 5, 3, 20, 12, 8, 6, 41, 25, 17, 13, 7, 84, 52, 36, 28, 16, 10, 169, 105, 73, 57, 33, 21, 11, 340, 212, 148, 116, 68, 44, 24, 14, 681, 425, 297, 233, 137, 89, 49, 29, 15, 1364, 852, 596, 468, 276, 180, 100, 60, 32, 18, 2729, 1705, 1193, 937, 553
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

Row 1: A084639.
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:
...
A191663=dispersion of A042948 (0 or 1 mod 4 and >1)
A054582=dispersion of A005843 (0 or 2 mod 4 and >1; evens)
A191664=dispersion of A014601 (0 or 3 mod 4 and >1)
A191665=dispersion of A042963 (1 or 2 mod 4 and >1)
A191448=dispersion of A005408 (1 or 3 mod 4 and >1, odds)
A191666=dispersion of A042964 (2 or 3 mod 4)
...
EXCEPT for at most 2 initial terms (so that column 1 always starts with 1):
A191663 has 1st col A042964, all else A042948
A054582 has 1st col A005408, all else A005843
A191664 has 1st col A042963, all else A014601
A191665 has 1st col A014601, all else A042963
A191448 has 1st col A005843, all else A005408
A191666 has 1st col A042948, all else A042964
...
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

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


EXAMPLE

Northwest corner:
1...4...9....20...41
2...5...12...25...52
3...8...17...36...73
6...13..28...57...116
7...16..33...68...137


MATHEMATICA

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


CROSSREFS

Cf. A042948, A042964, A191426.
Sequence in context: A191739 A091450 A163253 * A206561 A008831 A289506
Adjacent sequences: A191660 A191661 A191662 * A191664 A191665 A191666


KEYWORD

nonn,tabl


AUTHOR

Clark Kimberling, Jun 11 2011


STATUS

approved



