|
| |
|
|
A120862
|
|
Fixed-j dispersion for Q = 13: array D(g,h) (g, h >= 1), read by ascending antidiagonals.
|
|
7
|
|
|
|
1, 2, 10, 3, 20, 109, 4, 30, 218, 1189, 5, 43, 327, 2378, 12970, 6, 53, 469, 3567, 25940, 141481, 7, 63, 578, 5116, 38910, 282962, 1543321, 8, 76, 687, 6305, 55807, 424443, 3086642, 16835050, 9, 86, 829, 7494, 68777, 608761, 4629963, 33670100, 183642229
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,2
|
|
|
COMMENTS
|
For each positive integer n, there exists a unique pair (j,k) of positive integers such that (j+k+1)^2 - 4*k = 13*n^2; in fact, j(n) = A120869(n) and k(n) = A120870(n).
Suppose g >= 1 and let j = j(g). The numbers in row g of array D are among those n for which (j+k+1)^2 - 4*k = 13*n^2 for some k; that is, j stays fixed and k and n vary - hence the name "fixed-j dispersion". (The fixed-k dispersion for Q=13 is A120863.)
Every positive integer occurs exactly once in array D and every pair of rows are mutually interspersed. That is, beginning at the first term of any row having greater initial term than that of another row, all the following terms individually separate the individual terms of the other row.
|
|
|
LINKS
|
|
|
|
FORMULA
|
Define f(n) = floor(r*n) - floor(3*F(n)), where r = (11 + 3*sqrt(13))/2 and F(n) is the fractional part of (1 + sqrt(13))*n/2. Let D(g,h) be the term in row g and column h of the array to be defined:
D(1,1) = 1; D(1,2) = f(1); and D(1,h) = 11*D(1,h-1) - D(1,h-2) for h >= 3.
For arbitrary g >= 1, once row g is defined, define D(g+1,1) = least positive integer not in rows 1,2,...,g; D(g+1,2) = f(D(g+1,1)); and D(g+1,h) = 11*D(g+1,h-1) - D(g+1,h-2) for h >= 3. All rows after row 1 are thus inductively defined. [Corrected using Conjecture 1 in Kimberling (2007) by Petros Hadjicostas, Jul 07 2020]
|
|
|
EXAMPLE
|
Northwest corner:
1, 10, 109, 1189, ...
2, 20, 218, 2378, ...
3, 30, 327, 3567, ...
4, 43, 469, 5116, ...
5, 53, 578, 6305, ...
6, 63, 687, 7494, ...
...
|
|
|
PROG
|
(PARI) f(n) = floor((11 + 3*sqrt(13))/2*n) - floor(3*frac((1 + sqrt(13))*n/2));
unused(listus) = {my(v=vecsort(Vec(listus))); for (i=1, vecmax(v), if (!vecsearch(v, i), return (i)); ); };
D(nb) = {my(m = matrix(nb, nb), t); my(listus = List); for (g=1, nb, if (g==1, t = 1, t = unused(listus)); m[g, 1]=t; listput(listus, t); t = f(t); m[g, 2]=t; listput(listus, t); for (h=3, nb, t = 11*m[g, h-1] - m[g, h-2]; m[g, h] = t; listput(listus, t); ); ); m; };
lista(nb) = {my(m=D(nb)); for (n=1, nb, for (j=1, n, print1(m[n-j+1, j], ", "); ); ); } \\ Michel Marcus, Jul 09 2020
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
|
|
|
AUTHOR
|
|
|
|
EXTENSIONS
|
|
|
|
STATUS
|
approved
|
| |
|
|
