[go: up one dir, main page]

login
Revision History for A065684 (Bold, blue-underlined text is an addition; faded, red-underlined text is a deletion.)

Showing entries 1-10 | older changes
Number of primes <= prime(n) which begin with a 5.
(history; published version)
#12 by Giovanni Resta at Fri Jun 15 08:09:44 EDT 2018
STATUS

reviewed

approved

#11 by Joerg Arndt at Fri Jun 15 08:09:25 EDT 2018
STATUS

proposed

reviewed

#10 by Michael De Vlieger at Thu Jun 14 21:35:16 EDT 2018
STATUS

editing

proposed

#9 by Michael De Vlieger at Thu Jun 14 21:35:14 EDT 2018
MATHEMATICA

Accumulate@ Array[Boole[First@ IntegerDigits@ Prime@ # == 5] &, 103] (* Michael De Vlieger, Jun 14 2018 *)

STATUS

proposed

editing

#8 by Jon E. Schoenfield at Thu Jun 14 21:29:19 EDT 2018
STATUS

editing

proposed

#7 by Jon E. Schoenfield at Thu Jun 14 21:29:17 EDT 2018
LINKS

Harry J. Smith, <a href="/A065684/b065684.txt">Table of n, a(n) for n = 1,...,1000</a>

EXAMPLE

a(i) = 1 for 2 < i < 16 and a(16) = 2 as 53 = A000040(16) is the second prime beginning with 5. a(664579) = 72951 (A000040(664579) = 9999991 is the largest prime < 10,000,000^7).

PROG

(PARI) digitsIn(x)= { local(d); if (x==0, return(1)); d=1 + log(x)\log(10); if (10^d == x, d++, if (10^(d-1) > x, d--)); return(d) } MSD(x)= { return(x\10^(digitsIn(x)-1)) } { a=0; p=2; for (n=1, 1000, q=prime(n); while (p <= q, if(MSD(p) == 5, a++); p=nextprime(p+1)); write("b065684.txt", n, " ", a) ) } [From _\\ _Harry J. Smith_, Oct 26 2009]

CROSSREFS
STATUS

approved

editing

#6 by Russ Cox at Fri Mar 30 18:50:19 EDT 2012
AUTHOR

_Reinhard Zumkeller (reinhard.zumkeller(AT)gmail.com), _, Nov 13 2001

Discussion
Fri Mar 30
18:50
OEIS Server: https://oeis.org/edit/global/246
#5 by Russ Cox at Fri Mar 30 17:24:33 EDT 2012
PROG

(PARI) digitsIn(x)= { local(d); if (x==0, return(1)); d=1 + log(x)\log(10); if (10^d == x, d++, if (10^(d-1) > x, d--)); return(d) } MSD(x)= { return(x\10^(digitsIn(x)-1)) } { a=0; p=2; for (n=1, 1000, q=prime(n); while (p <= q, if(MSD(p) == 5, a++); p=nextprime(p+1)); write("b065684.txt", n, " ", a) ) } [From _Harry J. Smith (hjsmithh(AT)sbcglobal.net), _, Oct 26 2009]

Discussion
Fri Mar 30
17:24
OEIS Server: https://oeis.org/edit/global/133
#4 by N. J. A. Sloane at Thu Nov 11 07:34:06 EST 2010
LINKS

Harry J. Smith, <a href="/A065684/b065684.txt">Table of n, a(n) for n=1,...,1000</a>

KEYWORD

base,nonn,new

#3 by N. J. A. Sloane at Tue Jun 01 03:00:00 EDT 2010
LINKS

Harry J. Smith, <a href="b065684.txt">Table of n, a(n) for n=1,...,1000</a>

PROG

(PARI) digitsIn(x)= { local(d); if (x==0, return(1)); d=1 + log(x)\log(10); if (10^d == x, d++, if (10^(d-1) > x, d--)); return(d) } MSD(x)= { return(x\10^(digitsIn(x)-1)) } { a=0; p=2; for (n=1, 1000, q=prime(n); while (p <= q, if(MSD(p) == 5, a++); p=nextprime(p+1)); write("b065684.txt", n, " ", a) ) } [From Harry J. Smith (hjsmithh(AT)sbcglobal.net), Oct 26 2009]

KEYWORD

base,nonn,new