|
|
A098920
|
|
Least k such that k*M#(n) + 1 is prime where M#(n) is the product of the first n Mersenne primes = Product_{j=1..n} A000668(j).
|
|
0
|
|
|
2, 2, 2, 6, 30, 10, 4, 26, 12, 10, 12, 436, 644, 424, 2164, 862, 408, 9558, 5226, 12350
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
EXAMPLE
|
2*(2^2-1)*(2^3-1)*(2^5-1) + 1 = 1303 prime so a(3)=2.
|
|
MATHEMATICA
|
With[{s = FoldList[Times, Array[2^MersennePrimeExponent@ # - 1 &, 16]]}, Array[Block[{k = 2}, While[! PrimeQ[k s[[#]] + 1], k++]; k] &, Length@ s]] (* Michael De Vlieger, Dec 27 2019 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
hard,more,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|