OFFSET
0,2
COMMENTS
Possible upper bounds of twin primes pairs ending in 1: For a 30k + r "wheel", k > 0, r = 1, 13, 19 are the only possible values that can form an upper bound of a twin prime pair. The 30k+r wheel gives the sequence 1, 7, 11, 13, 17, 19, 23, 29 31, 37, 41, 43, 47, 49, 53, 59, ... which is frequently used in prime number sieves to skip multiples of 2, 3, 5. The fact that subtracting 2 from 30k+7, 11, 17, 23 will give us a multiple of 3 or 5 precludes these numbers from being an upper bound of a twin prime pair. This leaves us with r = 1, 13, 19 for k > 0 as the only possible cases to form an upper bound of a twin prime pair. 1, 13, 19 concludes the 6 numbers of the 8 number wheel that can form part of a twin prime pair.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..2000
Albert van der Horst, Counting Twin Primes
Index entries for linear recurrences with constant coefficients, signature (2,-1).
FORMULA
a(n) = 2*a(n-1) - a(n-2) for n > 1. - Vincenzo Librandi, Dec 30 2014
G.f.: (1 + 29*x)/(1 - x)^2. - Vincenzo Librandi, Dec 30 2014
E.g.f.: (1 + 30*x)*exp(x). - G. C. Greubel, Apr 04 2016
EXAMPLE
61 = 30 * 2 + 1, the upper part of the twin prime pair 59, 61.
MATHEMATICA
Range[1, 3001, 30] (* Vladimir Joseph Stephan Orlovsky, Jun 15 2011 *)
CoefficientList[Series[(1 + 29 x) / (1 - x)^2, {x, 0, 40}], x] (* Vincenzo Librandi, Dec 30 2014 *)
LinearRecurrence[{2, -1}, {1, 31}, 100] (* G. C. Greubel, Apr 04 2016 *)
PROG
(Magma) [30*n+1: n in [0..50]]; // Vincenzo Librandi, Jun 16 2011
(PARI) a(n)=30*n+1 \\ Charles R Greathouse IV, Oct 07 2015
(Scala) (0 to 49).map(30 * _ + 1) // Alonso del Arte, Jun 02 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Cino Hilliard, May 06 2007
STATUS
approved