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

%I #6 Aug 11 2012 09:38:41

%S 1,2,3,2,5,2,7,2,3,3,11,2,13,4,4,2,17,3,19,3,5,5,23,2,5,5,3,3,29,3,31,

%T 2,6,6,6,2,37,6,6,3,41,3,43,4,4,7,47,2,7,4,7,4,53,3,7,3,8,8,59,3,61,8,

%U 4,2,8,4,67,4,8,4,71,2,73,9,4,4,9,4,79,2,3,9,83,3,9,9,9,3,89,3,10,5,10,10

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

%C A079866(n)<=a(n)<=A079870(n); A020639(n)<=a(n)<=A006530(n);

%C a(m)=A079866(m)=A079870(m) iff m is a prime power (A000961).

%t Join[{1},Table[Floor[n^(1/PrimeOmega[n])+1/2],{n,2,100}]] (* _Harvey P. Dale_, Aug 11 2012 *)

%Y Cf. A079869(n)=a(n)^A001222(n), A079881.

%K nonn

%O 1,2

%A _Reinhard Zumkeller_, Jan 13 2003