OFFSET
1,2
COMMENTS
a(1) = 1 and all other terms are composite, of course.
LINKS
John Tyler Rascoe, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
Select[DivisorSigma[0, Range[200]], !PrimeQ[#]&] (* Harvey P. Dale, Mar 20 2015 *)
PROG
(PARI) for(i=1, 200, if(!isprime(numdiv(i)), print1(numdiv(i)", "))) \\ Franklin T. Adams-Watters, Apr 09 2009
(Python)
from sympy import primepi, integer_nthroot, primerange, divisor_count
def A141241(n):
def f(x): return int(n+sum(primepi(integer_nthroot(x, k-1)[0]) for k in primerange(x.bit_length()+1)))
m, k = n, f(n)
while m != k: m, k = k, f(k)
return divisor_count(m) # Chai Wah Wu, Feb 22 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Jun 16 2008
EXTENSIONS
More terms from Franklin T. Adams-Watters, Apr 09 2009
STATUS
approved