cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A113274 Record gaps between twin primes.

Original entry on oeis.org

2, 6, 12, 18, 30, 36, 72, 150, 168, 210, 282, 372, 498, 630, 924, 930, 1008, 1452, 1512, 1530, 1722, 1902, 2190, 2256, 2832, 2868, 3012, 3102, 3180, 3480, 3804, 4770, 5292, 6030, 6282, 6474, 6552, 6648, 7050, 7980, 8040, 8994, 9312, 9318, 10200, 10338, 10668
Offset: 1

Views

Author

Bernardo Boncompagni, Oct 21 2005

Keywords

Comments

a(n) mod 6 = 0 for each n>1.

Examples

			The first twin primes are 3,5 and 5,7 so a(0)=5-3=2. The following pair is 11,13 so a(1)=11-5=6. The following pair is 17,19 so 6 remains the record and no terms are added.
		

Crossrefs

The smallest primes originating the sequence are given in A113275. Cf. A008407, A005250, A002386.

Programs

  • Mathematica
    NextLowerTwinPrim[n_] := Block[{k = n + 6}, While[ !PrimeQ[k] || !PrimeQ[k + 2], k+=6]; k]; p = 5; r = 2; t = {2}; Do[ q = NextLowerTwinPrim[p]; If[q > r + p, AppendTo[t, q - p]; Print[{p, q - p}]; r = q - p]; p = q, {n, 10^9}]; t (* Robert G. Wilson v, Oct 22 2005 *)
    DeleteDuplicates[Differences[Select[Partition[Prime[Range[10^7]],2,1],#[[2]]-#[[1]] == 2&][[All,2]]],GreaterEqual] (* The program generates the first 27 terms of the sequence. *) (* Harvey P. Dale, Dec 31 2022 *)

Formula

a(n) = A036063(n) + 2.
a(n) = A036062(n) - A113275(n).
From Alexei Kourbatov, Dec 29 2011: (Start)
(1) Upper bound: gaps between twin primes are smaller than 0.76*(log p)^3, where p is the prime at the end of the gap.
(2) Estimate for the actual size of the maximal gap that ends at p: maximal gap = a(log(p/a)-1.2), where a = 0.76*(log p)^2 is the average gap between twin primes near p, as predicted by the Hardy-Littlewood k-tuple conjecture.
Formulas (1) and (2) are asymptotically equal as p tends to infinity. However, (1) yields values greater than all known gaps, while (2) yields "good guesses" that may be either above or below the actual size of known maximal gaps.
Both formulas (1) and (2) are derived from the Hardy-Littlewood k-tuple conjecture via probability-based heuristics relating the expected maximal gap size to the average gap. Neither of the formulas has a rigorous proof (the k-tuple conjecture itself has no formal proof either). In both formulas, the constant ~0.76 is reciprocal to the twin prime constant 1.32032...
(End)

Extensions

More terms from Robert G. Wilson v, Oct 22 2005
Corrected terms based on A036063, cross-checked with independent computations by Carlos Rivera and Richard Fischer (linked).
Terms up to a(72) are given in Kourbatov (2013), terms up to a(75) in Oliveira e Silva website.

A091592 Numbers n such that there are no twin primes between n^2 and (n+1)^2.

Original entry on oeis.org

1, 9, 19, 26, 27, 30, 34, 39, 49, 53, 77, 122
Offset: 1

Views

Author

Hugo Pfoertner, Jan 25 2004

Keywords

Comments

Numbers n such that there is no pair of twin primes P, P+2 with n^2 < P < P+2 < n^2+2*n.
The first 7 terms of this sequence were given by Ernst Jung in a discussion in the Newsgroup de.sci.mathematik entitled "Primzahlen zwischen (2x-1)^2 und (2x+1)^2" (primes between ...and...) with other significant contributions from Hermann Kremer and Rainer Rosenthal. It is conjectured that there are no further terms beyond a(12)=122. This has been tested to 50000 by Robert G. Wilson v.
Tested up to 10^7 and found no such numbers. - Arkadiusz Wesolowski, Jul 11 2011

Examples

			9 is a term because no twin primes are found in the interval [9^2,10^2].
		

Crossrefs

Programs

  • Maple
    isA091592 := proc(n) local p; p := nextprime(n^2) ; q := nextprime(p) ; while q < n^2+2*n do if q-p = 2 then RETURN(false) ; fi; p :=q ; q := nextprime(p) ; od: RETURN(true) ; end: for n from 1 do if isA091592(n) then printf("%d ",n) ; fi; od: # R. J. Mathar, Aug 26 2008
  • Mathematica
    fQ[n_] := StringCount[ ToString@ PrimeQ[ Range[n^2, (n + 1)^2]], "True, False, True"] == 0; lst = {}; Do[ If[ fQ@n, AppendTo[lst, n]], {n, 25000}]

Extensions

Edited by N. J. A. Sloane, Aug 31 2008 at the suggestion of Pierre CAMI

A036061 Increasing gaps among twin primes: the largest prime of the starting twin pair.

Original entry on oeis.org

5, 7, 19, 43, 73, 313, 349, 661, 2383, 5881, 13399, 18541, 24421, 62299, 187909, 687523, 688453, 850351, 2868961, 4869913, 9923989, 14656519, 17382481, 30752233, 32822371, 96894043, 136283431, 234966931, 248641039, 255949951, 390817729, 698542489, 2466641071
Offset: 1

Views

Author

Keywords

Comments

Has many terms in common with A054691, but neither of the two is a subsequence of the other one. - M. F. Hasler, May 07 2022

Crossrefs

Programs

  • Mathematica
    Block[{s = Select[Partition[Prime@ Range[10^7], 2, 1], Subtract @@ # == -2 &][[All, -1]], t}, t = Differences@ s; Map[s[[FirstPosition[t, #]]] &, Union@ FoldList[Max, t]][[All, 1]]] (* Michael De Vlieger, Jan 18 2019 *)

Formula

a(n) = A036062(n) - A036063(n).
a(n) = A113275(n) + 2.

Extensions

Terms 5, 7 prepended by Max Alekseyev, Nov 05 2015
a(17) corrected and a(31)-a(33) from Sean A. Irvine, Oct 21 2020

A036062 Increasing gaps among twin primes: the smallest prime of the second twin pair.

Original entry on oeis.org

5, 11, 29, 59, 101, 347, 419, 809, 2549, 6089, 13679, 18911, 24917, 62927, 188831, 688451, 689459, 851801, 2870471, 4871441, 9925709, 14658419, 17384669, 30754487, 32825201, 96896909, 136286441, 234970031, 248644217, 255953429
Offset: 1

Views

Author

Keywords

Crossrefs

Formula

a(n) = A036061(n) + A036063(n).

Extensions

Terms a(3)-a(41) are given by Rathbun (1998).
Corrected by Jud McCranie, Jan 04 2001
Terms up to a(72) are listed in Kourbatov (2013), terms up to a(75) on Oliveira e Silva's website, added by Max Alekseyev, Nov 06 2015
Showing 1-4 of 4 results.