login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A257495 The number of iterations (x -> 2x+1) until a prime is found, starting with prime(n). 1
1, 1, 1, 2, 1, 4, 2, 2, 1, 1, 2, 2, 1, 24, 2, 1, 2, 4, 2, 4, 2552, 4, 1, 1, 4, 8, 4, 2, 2, 1, 6, 1, 3, 4, 2, 2, 2, 8, 4, 1, 1, 2, 1, 8, 3, 6, 4, 4, 2, 2, 1, 1, 2, 1, 2, 3, 8, 2, 4, 1, 12, 1, 2, 21, 4, 3, 2, 4, 6, 2, 11, 1, 2, 16, 4, 4, 2, 4, 2, 8, 1, 12, 1, 8 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,4

COMMENTS

The number of iterations is defined as in A050412 (always positive).

Sophie Germain primes correspond to values a(n)=1 (A156660).

The plot without largest outliers allows detail on lower bound trending. Such outliers begin beyond the 121st entry. The number of terminal primes (those from the terminating iteration) being Sophie Germain through the first 10000 seeds is approximately 910. The number of Sophie Germain primes expected below 10000 is approximately 156 (computationally the comparison is more complicated, obviously).

a(2552) corresponds to 12525084203....315016703 with 771 digits. See A171390. - Vincenzo Librandi, Apr 27 2015

LINKS

Bill McEachen, Table of n, a(n) for n = 1..7075, using ispseudoprime() in the Pari code.

FORMULA

a(n) = A050412(prime(n)). - Michel Marcus, Jun 08 2015

EXAMPLE

Starting from prime(6)=13, sequential values for evaluation are 2*13+1=27, 2*27+1=55, 2*55+1=111, 2*111+1=223. The first prime is encountered at the 4th iteration, thus a(6)=4.

MAPLE

A257395 := proc(n)

    A050412(ithprime(n)) ;

end proc: # R. J. Mathar, Jul 23 2015 reusing code from A050412

MATHEMATICA

Length@ NestWhileList[2 # + 1 &, Prime@ #, CompositeQ, {2, 1}] - 1 & /@ Range@ 120 (* Michael De Vlieger, Apr 26 2015 *)

PROG

(PARI) genit()={

my(maxx=122, istrt=1, opt=1); n=istrt; cnt=1; val=2*prime(n)+1;

prev=val; prcnt=0; while(n<=maxx, if( val%6!=1 && val%6!=5, cnt+=1; val=2*val+1 );

if(ispseudoprime(val), print1(cnt, ", "); if(opt>0&&ispseudoprime(2*val+1), prcnt+=1);

cnt=1; n+=1; val=2*prime(n)+1; prev=val ); if(!ispseudoprime(val), cnt+=1; val=2*val+1));

}

(PARI) a(n, k=prime(n))=my(t=1); while(!ispseudoprime(k=2*k+1), t++); t \\ Charles R Greathouse IV, May 22 2015

CROSSREFS

Cf. A050412 (Riesel problem).

Sequence in context: A239101 A145983 A245025 * A120025 A109090 A220780

Adjacent sequences:  A257492 A257493 A257494 * A257496 A257497 A257498

KEYWORD

nonn

AUTHOR

Bill McEachen, Apr 26 2015

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 18 14:35 EDT 2019. Contains 322209 sequences. (Running on oeis4.)