OFFSET
1,1
COMMENTS
A126715(n) = smallest odd k such that k*2^n-1 is prime, the primes are not always in increasing order.
Here the primes k*2^n-1 are always in increasing order.
The ratio sum_{1..N}a(n)/sum_{1..N}n is near 2*log(2) as N increases.
The ratio a(n)/n is always < 8 for n from 1 to 6000.
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..6000
EXAMPLE
3*2^1-1 = 5 is prime, a(1)=3 by definition.
3*2^2-1 = 11 is prime, 3 > 3/2 so a(2) = 3.
3*2^3-1 = 23 is prime, so a(3) = 3.
3*2^4-1 = 47 is prime, so a(4) = 3.
3*2^5-1 = 95 is composite.
5*2^5-1 = 159 is composite.
7*2^5-1 = 223 is prime so a(5) = 7.
PROG
(PFGW & SCRIPT)
SCRIPT
DIM j, -1
DIM n, 0
DIMS t
OPENFILEOUT myf, a(n).txt
LABEL loop1
SET n, n+1
IF n>6000 THEN END
LABEL loop2
SET j, j+2
SETS t, %d, %d\,; n; j
PRP j*2^n-1, t
IF ISPRP THEN GOTO a
GOTO loop2
LABEL a
WRITE myf, t
SET j, j/2
IF j%2==0 THEN SET j, j+1
GOTO loop1
(PARI) a=[3]; for(n=2, 100, k=floor(a[n-1]/2)+2; if(k%2==0, k++); t=2^n; while(!isprime(k*t-1), k+=2); a=concat(a, k)); a \\ Colin Barker, Jul 22 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Jul 22 2014
STATUS
approved