login
Transpose of square array A285321.
3

%I #14 Nov 09 2024 23:45:04

%S 2,4,3,8,9,6,16,27,12,5,32,81,18,25,10,64,243,24,125,20,15,128,729,36,

%T 625,40,45,30,256,2187,48,3125,50,75,60,7,512,6561,54,15625,80,135,90,

%U 49,14,1024,19683,72,78125,100,225,120,343,28,21

%N Transpose of square array A285321.

%C See A285321.

%o (Scheme) (define (A285322 n) (A285321bi (A004736 n) (A002260 n))) ;; For A285321bi, see A285321.

%o (Python)

%o from functools import reduce

%o from operator import mul

%o from sympy import prime, primefactors

%o def a019565(n): return reduce(mul, (prime(i+1) for i, v in enumerate(bin(n)[:1:-1]) if v == '1')) if n > 0 else 1 # This function from Chai Wah Wu

%o def a007947(n): return 1 if n<2 else reduce(mul, primefactors(n))

%o def a065642(n):

%o if n==1: return 1

%o r=a007947(n)

%o n = n + r

%o while a007947(n)!=r:

%o n+=r

%o return n

%o def A(n, k): return a019565(k) if n==1 else a065642(A(n - 1, k))

%o for n in range(1, 11): print([A(n - k + 1, k) for k in range(1, n + 1)]) # _Indranil Ghosh_, Apr 19 2017

%Y Transpose: A285321.

%K nonn,tabl,changed

%O 1,1

%A _Antti Karttunen_, Apr 17 2017