login
A363455
The number of distinct primorial numbers (A002110) larger than 1 in the representation of A025487(n) as a product of primorial numbers.
3
0, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 1, 2, 1, 2, 3, 2, 2, 2, 2, 1, 2, 3, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 3, 2, 2, 2, 2, 2, 1, 3, 2, 1, 3, 2, 3, 2, 2, 2, 2, 2, 1, 3, 2, 2, 3, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 3, 2, 2, 3, 2, 3, 3
OFFSET
1,6
COMMENTS
The number of distinct exponents in the prime factorization of A025487(n).
The minimal number of powers of primorial numbers (A100778) in the representation of A025487(n) as a product of powers of primorial numbers.
The record values are all the nonnegative integers. The positions of the records are the positions of the terms of the Chernoff sequence (A006939) in A025487, i.e., the first position of k, for k = 0, 1, 2, ..., is A363456(k).
LINKS
FORMULA
a(n) = A071625(A025487(n)).
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
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 03 2023
STATUS
approved