A290450 Primes with property that the next prime has the same last digit.

139, 181, 241, 283, 337, 409, 421, 547, 577, 631, 691, 709, 787, 811, 829, 887, 919, 1021, 1039, 1051, 1153, 1171, 1249, 1399, 1471, 1627, 1637, 1699, 1723, 1801, 1879, 2017, 2029, 2053, 2089, 2143, 2521, 2647, 2719, 2731, 2767, 2887, 2917, 3001, 3089, 3109, 3361, 3413, 3517, 3547, 3571

Offset

1,1

Comments

Starts off the same as A031928, primes p such that the next prime is p + 10. First term that differs is 887, since 897 = 3 * 13 * 23 and the next prime is 907.

As the primes get larger and more sparsely distributed, the difference between successive primes is less likely to be less than 10.

One might expect that a prime is 1/4 as likely to be followed by a prime with the same least significant digit in base 10 (since the possibilities are 1, 3, 7, 9).

One might also expect this sequence to consist of a quarter of the primes. And yet pi(a(50)) = pi(3547) = 497; the 200th prime is 1223.

Links

Marius A. Burtea, Table of n, a(n) for n = 1..5005

Evelyn Lamb, "Peculiar Pattern Found in 'Random' Prime Numbers", Nature, March 14, 2016, republished by Scientific American.

Robert J. Lemke Oliver and Kannan Soundararajan, Unexpected biases in the distribution of consecutive primes, arXiv preprint arXiv:1603.03720 [math.NT], submitted on March 11, 2016.

Formula

A031928 UNION A031938 UNION A124596 UNION A126721 UNION ... - R. J. Mathar, Jan 23 2022

Example

139 is in the sequence because the immediately following prime is 149, which also ends in 9.
But 149 is not in the sequence because the next prime after that one is 151, which ends in 1, not 9.

Mathematica

Select[Partition[Prime[Range[500000]], 2, 1], Mod[#[[1]], 10] == Mod[#[[2]], 10] &][[All, 1]] ( * Harvey P. Dale, Aug 21 2017 *)
Module[{nn=500000, prs, p}, prs=Prime[Range[nn]]; p=Divisible[#, 10]&/@ Differences[prs]; Pick[Most[prs], p]] (* Harvey P. Dale, Aug 22 2017 *)

Prog

(Magma) f:=func<p, m|p mod 10 eq m and NextPrime(p) mod 10 eq m>; a:=[]; for p in PrimesUpTo(4000) do if f(p, 1) or f(p, 3) or f(p, 7) or f(p, 9) then Append(~a, p); end if; end for; a; // Marius A. Burtea, Oct 16 2019

Crossrefs

Cf. A031928 (subset), A050434 (with 2 digits).

Sequence in context: A020357 A050967 A071382 * A209618 A031928 A095649

Adjacent sequences: A290447 A290448 A290449 * A290451 A290452 A290453

Keyword

nonn,base

Author

Alonso del Arte, Aug 06 2017

Status

approved