login
a(n) = floor(n^(1/omega(n))).
2

%I #27 May 15 2024 22:30:26

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

%T 31,32,5,5,5,6,37,6,6,6,41,3,43,6,6,6,47,6,49,7,7,7,53,7,7,7,7,7,59,3,

%U 61,7,7,64,8,4,67,8,8,4,71,8,73,8,8,8,8,4,79,8,81,9,83,4,9,9,9,9,89,4,9,9

%N a(n) = floor(n^(1/omega(n))).

%H Harry J. Smith, <a href="/A066296/b066296.txt">Table of n, a(n) for n = 2..1000</a> [Corrected by Sean A. Irvine]

%F a(n) = floor(n^(1/omega(n))) = floor(n^(1/A001221(n))).

%F a(n) = n if n is a prime power (A246655).

%t Table[Floor[n^(1/PrimeNu[n])], {n, 2, 50}] (* _G. C. Greubel_, May 08 2017 *)

%o (PARI) { for (n=2, 1000, a=floor(n^(1/omega(n))); write("b066296.txt", n, " ", a) ) } \\ _Harry J. Smith_, Feb 08 2010

%Y Cf. A001221, A000040, A002808.

%K nonn

%O 2,1

%A _Labos Elemer_, Dec 12 2001