login
a(1)=1 and for n>1: floor(n^(1/Omega(n)))^Omega(n), where Omega(n) is the total number of prime factors of n (A001222).
4

%I #6 May 19 2018 12:08:52

%S 1,2,3,4,5,4,7,8,9,9,11,8,13,9,9,16,17,8,19,8,16,16,23,16,25,25,27,27,

%T 29,27,31,32,25,25,25,16,37,36,36,16,41,27,43,27,27,36,47,32,49,27,49,

%U 27,53,16,49,16,49,49,59,16,61,49,27,64,64,64,67,64,64,64,71,32,73,64

%N a(1)=1 and for n>1: floor(n^(1/Omega(n)))^Omega(n), where Omega(n) is the total number of prime factors of n (A001222).

%C a(n)<=A079869(n); A020639(n)<=a(n)<=A006530(n);

%C a(m)=m=A079869(m)=A079871(m) iff m is a prime power (A000961).

%t Join[{1},Table[Floor[n^(1/PrimeOmega[n])]^PrimeOmega[n],{n,2,80}]] (* _Harvey P. Dale_, May 19 2018 *)

%Y a(n)=A079866(n)^A001222(n), cf. A068794, A068795.

%K nonn

%O 1,2

%A _Reinhard Zumkeller_, Jan 13 2003