OFFSET
1,1
COMMENTS
Maple implementation: see A030513.
Numbers of the form p^25 (5th powers of A050997, subset of A010813) or p*q^12, where p and q are distinct primes. - R. J. Mathar, Mar 01 2010
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
FORMULA
A000005(a(n))=26.
MATHEMATICA
Select[Range[900000], DivisorSigma[0, #]==26&] (* Vladimir Joseph Stephan Orlovsky, May 05 2011 *)
PROG
(PARI) is(n)=numdiv(n)==26 \\ Charles R Greathouse IV, Jun 19 2016
(Python)
from sympy import integer_nthroot, primerange, primepi
from oeis_sequences.OEISsequences import bisection
def A137489(n):
def f(x): return int(n+x-sum(primepi(x//p**12) for p in primerange(integer_nthroot(x, 12)[0]+1))+primepi(integer_nthroot(x, 13)[0])-primepi(integer_nthroot(x, 25)[0]))
return bisection(f, n, n) # Chai Wah Wu, May 01 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
R. J. Mathar, Apr 22 2008
STATUS
approved
