OFFSET
1,1
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
18 = 2*3*3 and all prime divisors are consecutive primes but the least prime signature is 12 = 2*2*3; so a(1) = 18.
MATHEMATICA
With[{nn = 4800}, Select[Range[2, nn], And[#1 != Times @@ MapIndexed[Prime[First@ #2]^#1 &, Sort[#3, Greater]], Last[#2] == Prime@ Length[#2]] & @@ Apply[Join, {{#1}, Transpose@ #2}] & @@ {#, FactorInteger[#]} &] ] (* Michael De Vlieger, Feb 06 2020 *)
PROG
(Python)
from itertools import count
from functools import lru_cache
from sympy import prime, integer_log, primorial
from oeis_sequences.OEISsequences import bisection
def A056808(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)
@lru_cache(maxsize=None)
def h(x, m): return sum(h(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 -= h(x, k)-g(x, k, 1)
return c
return bisection(f, n, n) # Chai Wah Wu, Mar 23 2026
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Alford Arnold, Aug 22 2000
EXTENSIONS
More terms from Larry Reeves (larryr(AT)acm.org), Nov 28 2000
STATUS
approved
