|
|
A058233
|
|
Primes p such that p#+1 is divisible by the next prime after p.
|
|
8
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
No additional terms through the 100000th prime. - Harvey P. Dale, Mar 12 2014
a(5) > prime(1400000) = 22182343. - Robert Price, Apr 02 2018
|
|
LINKS
|
|
|
EXAMPLE
|
2*3*5*7*11*13*17+1 is divisible by 19.
|
|
MATHEMATICA
|
primorial[n_] := Product[ Prime[k], {k, 1, PrimePi[n]}]; Select[ Prime[ Range[1000]], Divisible[ primorial[#] + 1, NextPrime[#]] &] (* Jean-François Alcover, Aug 19 2013 *)
Module[{prs=Prime[Range[500]]}, Transpose[Select[Thread[{Rest[ FoldList[ Times, 1, prs]], prs}], Divisible[ First[#]+1, NextPrime[Last[#]]]&]][[2]]] (* Harvey P. Dale, Mar 12 2014 *)
|
|
PROG
|
(Python)
from sympy import nextprime
A058233_list, p, q, r = [], 2, 3, 2
for _ in range(10**3):
if (r+1) % q == 0:
r *= q
|
|
CROSSREFS
|
|
|
KEYWORD
|
nice,nonn,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|