|
%I #7 Jun 17 2020 19:20:59
%S 1,2,3,4,5,4,7,8,9,9,11,8,13,16,16,16,17,27,19,27,25,25,23,16,25,25,
%T 27,27,29,27,31,32,36,36,36,16,37,36,36,81,41,27,43,64,64,49,47,32,49,
%U 64,49,64,53,81,49,81,64,64,59,81,61,64,64,64,64,64,67,64,64,64,71,32,73
%N a(1)=1 and for n>1: round(n^(1/Omega(n)))^Omega(n), where Omega(n) is the total number of prime factors of n (A001222).
%C A079867(n)<=a(n)<=A079869(n); A020639(n)<=a(n)<=A006530(n);
%C a(m)=m=A079867(m)=A079871(m) iff m is a prime power (A000961).
%t ron[n_]:=Module[{c=PrimeOmega[n]},Round[n^(1/c)]^c]; Join[{1},Array[ ron,80,2]] (* _Harvey P. Dale_, Jun 17 2020 *)
%Y a(n)=A079868(n)^A001222(n). Cf. A079872, A068794, A068795.
%K nonn
%O 1,2
%A _Reinhard Zumkeller_, Jan 13 2003
|
