login
A273283(n)^Omega(n), where Omega = A001222.
5

%I #23 Jun 03 2016 15:02:51

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

%T 27,125,29,125,31,32,49,49,49,81,37,49,49,81,41,125,43,125,125,49,47,

%U 243,49,125,121,125,53,81,121,81,121,121,59,81,61,121,125,64,121,125,67,125,121,125,71

%N A273283(n)^Omega(n), where Omega = A001222.

%C a(n) >= A079871(n) for any n>=2.

%C a(n) = n iff n is the power of a prime (A246655).

%H Giuseppe Coppoletta, <a href="/A273285/b273285.txt">Table of n, a(n) for n = 2..10000</a>

%F a(n) = A273283(n)^A001222(n).

%e a(33) = 49 because Omega(33)=2 and 5^2 < 33 < 7^2.

%e a(2403) = 11^4 (and A273284(2403)=7^4) because A001222(3^3*89)=4 and 7^4 < n < 11^4.

%t Table[NextPrime[(Ceiling[n^(1/PrimeOmega[n])] - 1)]^PrimeOmega[n], {n,2,50}] (* _G. C. Greubel_, May 26 2016 *)

%o (Sage) s=sloane.A001222; [next_prime(ceil(n^(1/s(n)))-1)^s(n) for n in (2..82)]

%Y Cf. A273283, A273284, A273291, A079871, A001222, A246655.

%K nonn

%O 2,1

%A _Giuseppe Coppoletta_, May 20 2016