%I #22 Oct 04 2024 05:43:46
%S 5,11,11,17,29,29,37,53,67,83,127,127,131,173,251,257,293,347,367,521,
%T 541,631,733,853,967,1031,1361,1373,1693,1861,2053,2203,2203,2213,
%U 2411,2819,3137,3491,3727,4099,4493,4919,5051,5333,6247,6563,6863,6899,7927
%N Smallest prime > the n-th nontrivial power of a prime.
%H Amiram Eldar, <a href="/A060846/b060846.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = nextprime(A025475(n+1)) = A007918(A025475(n+1)) = Min{p| p>A025475(n+1)}. [corrected by _Michel Marcus_, Aug 24 2019]
%e 78125=5^7 is followed by 78137.
%t NextPrime[Select[Range[10^4], !PrimeQ[#] && PrimePowerQ[#] &]] (* _Amiram Eldar_, Oct 04 2024 *)
%o (PARI) ispp(x) = !isprime(x) && isprimepower(x);
%o lista(nn) = apply(x->nextprime(x), select(x->ispp(x), [1..nn])); \\ _Michel Marcus_, Aug 24 2019
%o (Python)
%o from sympy import primepi, integer_nthroot, nextprime
%o def A060846(n):
%o def bisection(f,kmin=0,kmax=1):
%o while f(kmax) > kmax: kmax <<= 1
%o while kmax-kmin > 1:
%o kmid = kmax+kmin>>1
%o if f(kmid) <= kmid:
%o kmax = kmid
%o else:
%o kmin = kmid
%o return kmax
%o def f(x): return int(n+x-sum(primepi(integer_nthroot(x,k)[0]) for k in range(2,x.bit_length())))
%o return nextprime(bisection(f,n,n)) # _Chai Wah Wu_, Sep 15 2024
%Y Cf. A025475, A000961, A001597, A001694, A007917, A007918, A013632, A013633, A049711, A060845, A068435, A246547.
%K nonn
%O 1,1
%A _Labos Elemer_, May 03 2001