|
%I #8 Dec 07 2019 12:18:29
%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 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
%O 1,1
%A _Antti Karttunen_, Apr 17 2017
|
