OFFSET
1,4
COMMENTS
LINKS
EXAMPLE
The top left corner of the array:
k= 1 2 3 4 5 6 7 8 9 10 11 12
2k= 2 4 6 8 10 12 14 16 18 20 22 24
--+-------------------------------------------------------------------------------
1 | -1, 1, 3, 11, 1, 21, 5, 49, 39, 23, -5, 87,
2 | -1, 7, 5, 71, 13, 85, -1, 463, 95, 149, 7, 605,
3 | -3, -1, 7, 93, -19, 189, -11, 1151, 357, 87, -37, 2023,
4 | -3, 23, -11, 645, 5, 495, -29, 9839, 165, 783, -13, 9757,
5 | -9, -73, -65, -465, -127, -273, -119, -721, 39, -903, -129, 2743,
6 | -9, -49, -119, 519, -103, -255, -105, 26399, -1377, 225, -227, 18649,
7 | -15, -217, -209, -2967, -231, -2679, -397, -36721, -2223, -2825, -351, -28937,
...
Terms of column 9: 39 (3*13), 95 (5*19), 357 (3*7*17), 165 (3*5*11), 39 (3*13), -1377 (- 3^4 * 17), -2223 (- 3^2 * 13 * 19), ..., show an ascending "prime thread" (3, 5, 7, 11, 13, 17, 19, ...) that is mentioned in comments.
PROG
(PARI)
up_to = 66;
A246278sq(row, col) = if(1==row, 2*col, my(f = factor(2*col)); for(i=1, #f~, f[i, 1] = prime(primepi(f[i, 1])+(row-1))); factorback(f));
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A372562sq(row, col) = A252748(A246278sq(row, col));
A372562list(up_to) = { my(v = vector(up_to), i=0); for(a=1, oo, for(col=1, a, i++; if(i > up_to, return(v)); v[i] = A372562sq(col, (a-(col-1))))); (v); };
v372562 = A372562list(up_to);
A372562(n) = v372562[n];
KEYWORD
sign,tabl
AUTHOR
Antti Karttunen, May 21 2024
STATUS
approved