login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A091859
a(1) = 2, a(n) is the least prime of the form 2k*a(n-1)+1 where 2k is the least even number not used earlier.
0
2, 5, 31, 311, 3733, 59729, 1194581, 57339889, 1605516893, 22477236503, 1933042339259, 143045133105167, 2574812395893007, 87543621460362239, 700348971682897913, 58829313621363424693, 1411903526912722192633
OFFSET
1,1
EXAMPLE
a(2) = 5 = 2*2+1, But a(3) = 31 though 11= 5*2 +1 is also prime.
PROG
Contribution from Rick L. Shepherd, Sep 08 2009: (Start)
(PARI) {a=[2]; terms=1; used=Set();
while(terms<30, k2=0; until(isprime(k2*a[terms]+1),
until(setsearch(used, k2)==0, k2=k2+2));
a=concat(a, k2*a[terms]+1); terms++; used=setunion(used, Set([k2]))); a} (End)
CROSSREFS
Sequence in context: A192397 A097396 A056788 * A085873 A303289 A051048
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Mar 13 2004
EXTENSIONS
a(6) onwards by Rick L. Shepherd, Sep 08 2009
STATUS
approved