OFFSET
1,1
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
MATHEMATICA
With[{P = FoldList[Times, Prime@ Range@ Max@ #]}, Map[P[[#]] &, #]] &@ Map[PrimeNu@ # &, Select[Range[10^4], Last[#] == Length[#] &@ PrimePi@ FactorInteger[#][[All, 1]] &]] (* Michael De Vlieger, Feb 06 2020 *)
PROG
(Python)
from math import prod
from itertools import count
from functools import lru_cache
from sympy import prime, integer_log, primorial, primefactors
from oeis_sequences.OEISsequences import bisection
def A080404(n):
@lru_cache(maxsize=None)
def g(x, m): return sum(g(x//(prime(m)**i), m-1) for i in range(1, integer_log(x, prime(m))[0]+1)) if m-1 else x.bit_length()-1
def f(x):
c = n+x
for k in count(1):
if primorial(k)>x:
break
c -= g(x, k)
return c
return prod(primefactors(bisection(f, n+1, n+1))) # Chai Wah Wu, Mar 18 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Mar 19 2003
STATUS
approved
