OFFSET
1,1
LINKS
EXAMPLE
The top left corner of the array:
k= | 1 2 3 4 5 6 7 8 9 10 11
----+------------------------------------------------------------------------
1 | 13, 19, 17, 21, 29, 73, 41, 33, 41, 121, 45,
2 | 28, 46, 34, 46, 52, 316, 78, 64, 70, 562, 82,
3 | 66, 92, 86, 96, 116, 1866, 188, 126, 146, 2762, 176,
4 | 134, 210, 148, 176, 190, 7190, 274, 218, 260, 11690, 274,
5 | 276, 326, 320, 342, 386, 35124, 494, 452, 474, 42962, 540,
6 | 404, 528, 430, 482, 560, 71384, 628, 586, 664, 95714, 716,
7 | 630, 704, 698, 800, 834, 187902, 944, 936, 1004, 215210, 1038,
8 | 818, 990, 932, 970, 1084, 304058, 1422, 1160, 1198, 375290, 1274,
9 | 1220, 1538, 1266, 1404, 1496, 661412, 1706, 1542, 1634, 846650, 1772,
10 | 1770, 1892, 1944, 2060, 2118, 1515570, 2480, 2234, 2408, 1645322, 2582,
PROG
(PARI)
up_to_cols = 210;
A326127(n) = (sigma(n)-core(n)-n);
A388984list(up_to) = { my(v=vector(up_to), i=0); forstep(n=2, oo, 2, if(isprime(core(n)) && !isprimepower(n), i++; v[i] = n; if(i==up_to, return(v)))); };
v388984 = A388984list(up_to_cols);
A388984(n) = v388984[n];
up_to = binomial(up_to_cols+1, 2)
A388981sq(row, col) = if(1==row, A388984(col), my(f = factor(A388984(col))); for(i=1, #f~, f[i, 1] = prime(primepi(f[i, 1])+(row-1))); factorback(f));
A388987sq(row, col) = A326127(A388981sq(row, col));
A388987list(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] = A388987sq(col, (a-(col-1))))); (v); };
v388987 = A388987list(up_to);
A388987(n) = v388987[n];
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Antti Karttunen, Sep 22 2025
STATUS
approved
