login
Symmetric square array A(n,k) = A122111(A122111(n) * A122111(k)), n >= 1, k >= 1, read by antidiagonals.
5

%I #16 Aug 04 2022 22:28:49

%S 1,2,2,3,3,3,4,5,5,4,5,6,7,6,5,6,7,10,10,7,6,7,10,11,9,11,10,7,8,11,

%T 14,14,14,14,11,8,9,12,13,15,13,15,13,12,9,10,15,20,22,22,22,22,20,15,

%U 10,11,14,21,18,17,21,17,18,21,14,11,12,13,22,25,28,26,26,28,25,22,13,12,13,20,17,21,33,30,19,30,33,21,17,20,13

%N Symmetric square array A(n,k) = A122111(A122111(n) * A122111(k)), n >= 1, k >= 1, read by antidiagonals.

%C A122111 is a self-inverse permutation, so this array represents a binary operation A(.,.) over the positive integers that is isomorphic to multiplication. Its primes are the positive powers of 2 (as defined by standard multiplication): 2, 4, 8, 16, 32, ... . The positive powers of 2, as defined by A(.,.), are the prime numbers as we usually understand them: 2, 3, 5, 7, 11, ... . - _Peter Munn_, Aug 04 2022

%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>

%F A(n, A061142(n)) = A003961(n). - _Peter Munn_, Aug 04 2022

%e The top left 15 X 15 corner of the array:

%e 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,

%e 2, 3, 5, 6, 7, 10, 11, 12, 15, 14, 13, 20, 17, 22, 21,

%e 3, 5, 7, 10, 11, 14, 13, 20, 21, 22, 17, 28, 19, 26, 33,

%e 4, 6, 10, 9, 14, 15, 22, 18, 25, 21, 26, 30, 34, 33, 35,

%e 5, 7, 11, 14, 13, 22, 17, 28, 33, 26, 19, 44, 23, 34, 39,

%e 6, 10, 14, 15, 22, 21, 26, 30, 35, 33, 34, 42, 38, 39, 55,

%e 7, 11, 13, 22, 17, 26, 19, 44, 39, 34, 23, 52, 29, 38, 51,

%e 8, 12, 20, 18, 28, 30, 44, 27, 50, 42, 52, 45, 68, 66, 70,

%e 9, 15, 21, 25, 33, 35, 39, 50, 49, 55, 51, 70, 57, 65, 77,

%e 10, 14, 22, 21, 26, 33, 34, 42, 55, 39, 38, 66, 46, 51, 65,

%e 11, 13, 17, 26, 19, 34, 23, 52, 51, 38, 29, 68, 31, 46, 57,

%e 12, 20, 28, 30, 44, 42, 52, 45, 70, 66, 68, 63, 76, 78, 110,

%e 13, 17, 19, 34, 23, 38, 29, 68, 57, 46, 31, 76, 37, 58, 69,

%e 14, 22, 26, 33, 34, 39, 38, 66, 65, 51, 46, 78, 58, 57, 85,

%e 15, 21, 33, 35, 39, 55, 51, 70, 77, 65, 57, 110, 69, 85, 91,

%o (PARI)

%o up_to = 105;

%o A122111(n) = if(1==n,n,my(f=factor(n), es=Vecrev(f[,2]),is=concat(apply(primepi,Vecrev(f[,1])),[0]),pri=0,m=1); for(i=1, #es, pri += es[i]; m *= prime(pri)^(is[i]-is[1+i])); (m));

%o A350066sq(n,k) = A122111(A122111(n)*A122111(k));

%o A350066list(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] = A350066sq(col,(a-(col-1))))); (v); };

%o v350066 = A350066list(up_to);

%o A350066(n) = v350066[n]; \\ _Antti Karttunen_, Dec 13 2021

%Y Cf. A122111, A297002 (main diagonal), A253550 (after its initial term, gives row 2 / column 2 from the second term onward).

%Y See the formula section for the relationships with A003961, A061142.

%Y Cf. also A003991, A129595, A331590.

%K nonn,tabl

%O 1,2

%A _Antti Karttunen_, Dec 13 2021