%I #18 Sep 15 2024 22:00:16
%S 3,7,7,13,23,23,31,47,61,79,113,113,127,167,241,251,283,337,359,509,
%T 523,619,727,839,953,1021,1327,1367,1669,1847,2039,2179,2179,2207,
%U 2399,2803,3121,3469,3719,4093,4483,4909,5039,5323,6229,6553,6857,6883,7919
%N Largest prime < a nontrivial power of a prime.
%H Harry J. Smith, <a href="/A060845/b060845.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) = prevprime[A025475(n)] = A007917[A025475(n)] = Max{p| p < A025475(n)}
%e 78125=5^7 follows 78121
%t Take[NextPrime[#,-1]&/@Union[Flatten[Table[Prime[p]^n,{n,2,20},{p,25}]]], 50] (* _Harvey P. Dale_, Mar 26 2012 *)
%o (PARI) { m=1; for (n=1, 1000, m++; while(sigma(m)*eulerphi(m)*(1 - isprime(m)) <= (m - 1)^2, m++); write("b060845.txt", n, " ", precprime(m - 1)); ) } \\ _Harry J. Smith_, Jul 19 2009
%o (Python)
%o from sympy import primepi, integer_nthroot, prevprime
%o def A060845(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 prevprime(bisection(f,n,n)) # _Chai Wah Wu_, Sep 15 2024
%Y Cf. A025475, A000961, A001597, A001694, A007917, A007918, A013632, A013633, A049711, A060846.
%K nonn
%O 1,1
%A _Labos Elemer_, May 03 2001