|
|
A066889
|
|
a(n) = g(P(n)+2) - P(n), where P(n) = Product_{k=1..n} Fibonacci(k) and g(i) is the smallest prime >= i.
|
|
2
|
|
|
2, 2, 3, 5, 7, 11, 17, 17, 37, 23, 47, 37, 29, 19, 47, 59, 19, 37, 71, 59, 31, 67, 239, 101, 739, 409, 43, 367, 167, 251, 73, 71, 419, 1567, 107, 83, 223, 191, 227, 449, 97, 173, 103, 523, 79, 137, 223, 1163, 661, 103, 103, 541, 227, 2383, 433, 71
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The first 169 terms are primes. Are all terms primes? See links for similar sequences.
Note that g is not the usual "nextprime" function. If the usual nextprime function is used, we get A286296.
|
|
LINKS
|
|
|
EXAMPLE
|
a(4) = 5 because Fibonacci(1)*Fibonacci(2)*Fibonacci(3)*Fibonacci(4) = 1*1*2*3 = 6, g(6+2) = 11, and 11 - 6 = 5.
|
|
MATHEMATICA
|
Join[{2, 2}, Drop[NextPrime[#+2]-#&/@FoldList[Times, Fibonacci[ Range[ 60]]], 2]] (* Harvey P. Dale, May 31 2017 *)
|
|
PROG
|
(PARI) { m=1; for (n=1, 1000, m*=fibonacci(n); write("b066889.txt", n, " ", nextprime(m+2) - m) ) } \\ Harry J. Smith, Apr 05 2010
(MuPAD) f := 1:for n from 1 to 100 do f := f*numlib::fibonacci(n):a := nextprime(f+2)-f:print(a) end_for
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Frank Buss (fb(AT)frank-buss.de), Jan 22 2002
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|