OFFSET
1,6
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
e[1] = 0; e[n_] := Length[Union[FactorInteger[n][[;; , 2]]]]; s = {0}; Do[If[GreaterEqual @@ (f = FactorInteger[n])[[;; , 2]] && PrimePi[f[[-1, 1]]] == Length[f], AppendTo[s, e[n]]], {n, 2, 10000}]; s
PROG
(Python)
from itertools import count
from functools import lru_cache
from sympy import prime, integer_log, primorial, factorint
from oeis_sequences.OEISsequences import bisection
def A363455(n):
@lru_cache(maxsize=None)
def g(x, m, j): return sum(g(x//(prime(m)**i), m-1, i) for i in range(j, integer_log(x, prime(m))[0]+1)) if m-1 else max(0, x.bit_length()-j)
def f(x):
c = n-1+x
for k in count(1):
if primorial(k)>x:
break
c -= g(x, k, 1)
return c
return len(set(factorint(bisection(f, n, n)).values())) # Chai Wah Wu, Mar 23 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 03 2023
STATUS
approved
