%I #9 Sep 21 2025 20:40:04
%S 3,5,5,6,14,9,9,21,44,13,12,41,57,90,21,8,31,219,135,230,25,18,57,93,
%T 629,261,324,33,17,51,277,161,2529,411,560,37,15,122,111,933,345,4211,
%U 609,702,45,18,101,1094,213,2857,461,9519,831,1034,57,30,85,393,4402,387,5325,741,13337,1281,1652,61
%N Square array A(n, k) = 3*A246278(n, k) - sigma(A246278(n, k)), read by falling antidiagonals. A033885(n) = 3*n-sigma(n) as applied to the prime shift array.
%H Antti Karttunen, <a href="/A388285/b388285.txt">Table of n, a(n) for n = 1..22155; the first 210 antidiagonals</a>
%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>.
%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>.
%F A(n, k) = A033885(A246278(n, k)) = 3*A246278(n, k) - A355927(n, k).
%F A(n, k) = A246278(n,k) + A378979(n, k).
%F A(n, k) = A388284(n,k) * A355925(n, k).
%e The top left corner of the array:
%e k= | 1 2 3 4 5 6 7 8 9 10 11 12
%e 2k= | 2 4 6 8 10 12 14 16 18 20 22 24
%e ----+-------------------------------------------------------------------------
%e 1 | 3, 5, 6, 9, 12, 8, 18, 17, 15, 18, 30, 12,
%e 2 | 5, 14, 21, 41, 31, 57, 51, 122, 101, 85, 61, 165,
%e 3 | 9, 44, 57, 219, 93, 277, 111, 1094, 393, 453, 147, 1377,
%e 4 | 13, 90, 135, 629, 161, 933, 213, 4402, 1477, 1113, 239, 6519,
%e 5 | 21, 230, 261, 2529, 345, 2857, 387, 27818, 3381, 3777, 471, 31413,
%e 6 | 25, 324, 411, 4211, 461, 5325, 561, 54742, 6973, 5973, 711, 69207,
%e 7 | 33, 560, 609, 9519, 741, 10333, 939, 161822, 11553, 12573, 1005, 175641,
%e Columns 60 and 90:
%e k= | 60 90
%e 2k= | 120 180
%e ----+------------------------------
%e 1 | 0, -6,
%e 2 | 915, 1501,
%e 3 | 13899, 19221,
%e 4 | 79947, 125097,
%e 5 | 513525, 604797,
%e 6 | 1272093, 1660317,
%e 7 | 3935343, 4391493,
%e 8 | 8512059, 10293513,
%e 9 | 20603199, 25961461,
%o (PARI)
%o up_to = 66;
%o A033885(n) = (3*n-sigma(n));
%o 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));
%o A388285sq(row,col) = A033885(A246278sq(row,col));
%o A388285list(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] = A388285sq(col,(a-(col-1))))); (v); };
%o v388285 = A388285list(up_to);
%o A388285(n) = v388285[n];
%Y Cf. A033885, A246278, A355925, A355927, A378979, A388284.
%K sign,tabl
%O 1,1
%A _Antti Karttunen_, Sep 21 2025