OFFSET
1,1
COMMENTS
Heuristically, the density of the sequence in the primes should approach Artin's constant: 0.3739558136...
LINKS
Robert Israel, Table of n, a(n) for n = 1..500
EXAMPLE
For n = 2 one has a(2) = 3 since B(2) = 2 is a primitive root modulo 3.
MAPLE
filter:= proc(p) local b;
b:= combinat:-bell(p-1);
numtheory:-order(b, p) = p-1
end proc:
select(filter, [seq(ithprime(i), i=1..200)]); # Robert Israel, May 04 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Luis H. Gallardo, May 04 2022
EXTENSIONS
Corrected by Robert Israel, May 04 2023
STATUS
approved