|
%I #23 Jun 03 2019 09:14:57
%S 3,7,71,18089,10391023,781379079653017,2111421691000680031,
%T 1430286763442005122380663256416207
%N Primes which are the denominators of convergents of the continued fraction expansion of e.
%C The position of a(n) in A000040 (the prime numbers) is A102049(n) = A000720(a(n)). - _Jonathan Sondow_, Dec 27 2004
%C The next term has 166 digits. [_Harvey P. Dale_, Aug 23 2011]
%H Joerg Arndt, <a href="/A094008/b094008.txt">Table of n, a(n) for n = 1..10</a>
%H E. B. Burger, <a href="https://www.jstor.org/stable/2695737">Diophantine Olympics and World Champions: Polynomials and Primes Down Under</a>, Amer. Math. Monthly, 107 (Nov. 2000), 822-829.
%H J. Sondow, <a href="https://www.jstor.org/stable/27642006">A geometric proof that e is irrational and a new measure of its irrationality</a>, Amer. Math. Monthly 113 (2006) 637-641 (article), 114 (2007) 659 (addendum).
%H J. Sondow, <a href="http://arxiv.org/abs/0704.1282"> A geometric proof that e is irrational and a new measure of its irrationality</a>, arXiv:0704.1282 [math.HO], 2007-2010.
%H J. Sondow and K. Schalm, <a href="http://arxiv.org/abs/0709.0671">Which partial sums of the Taylor series for e are convergents to e? (and a link to the primes 2, 5, 13, 37, 463), II</a>, Gems in Experimental Mathematics (T. Amdeberhan, L. A. Medina, and V. H. Moll, eds.), Contemporary Mathematics, vol. 517, Amer. Math. Soc., Providence, RI, 2010; arXiv:0709.0671 [math.NT], 2007-2009.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/e.html">e</a>.
%F a(n) = A007677(A094007(n)) = A000040(A102049(n)).
%e a(1) = 3 because 3 is the first prime denominator of a convergent, 8/3, of the simple continued fraction for e
%t Block[{$MaxExtraPrecision=1000},Select[Denominator[Convergents[E,500]], PrimeQ]] (* _Harvey P. Dale_, Aug 23 2011 *)
%o (PARI)
%o default(realprecision,10^5);
%o cf=contfrac(exp(1));
%o n=0;
%o { for(k=1, #cf, \\ generate b-file
%o pq = contfracpnqn( vector(k,j, cf[j]) );
%o p = pq[1,1]; q = pq[2,1];
%o \\ if ( ispseudoprime(p), n+=1; print(n," ",p) ); \\ A086791
%o if ( ispseudoprime(q), n+=1; print(n," ",q) ); \\ A094008
%o ); }
%o /* _Joerg Arndt_, Apr 21 2013 */
%Y Cf. A094007.
%Y See also A000040, A000720, A007677, A102049.
%K nonn
%O 1,1
%A _Jonathan Sondow_, Apr 20 2004
| 