OFFSET
0,1
COMMENTS
Old name was: A four-term repeating sequence for constructing a summation sequence from negative to positive infinity containing all primes except 2 and 5.
a(n) allows for the creation of an infinite summation sequence, s(n), extending from negative to positive infinity. (See Formula section below.) With appropriate initialization, letting "s(n+)" be the set positive s(n) values, and "s(n-)" be the absolute value of the set of negative s(n) values, the following applies:
s(n+) includes all primes of the form 4*m+1 with m>=2. Thus excluding 5.
s(n-) includes all primes of the form 4*m+3 with m>=0.
Together these include all primes (except 2 and 5) without duplication.
The primes "p(+)" within s(n+) "appear" in the form 3*p(+) within s(n-).
The primes "p(-)" within s(n-) "appear" in the form 3*p(-) within s(n+).
By using this simple repeating pattern, rather than the two well known linear formulas above, all primes (except 2 and 5) are included via a single construction mechanism, and all integers ending in the digit 5 are excluded mathematically, resulting in fewer nonprimes among the values of s(n) than there are in the combination of 4*m+1 and 4*m+3.
(NOTE: In the above "m" is not that same index as "n").
This is one of only two such repeating sequences with the property of generating a summation sequence that includes all integers ending in 1,3,7 or 9, and thus all primes except 2 and 5 (for the other see A226294). Both have the same density of primes in s(n), because both generate only 40% of the integers (in absolute value). And both presumably have the same average density of primes in positive vs. negative values of s(n).
Also, continued fraction expansion of 4 + sqrt(646)/6. - Bruno Berselli, Jun 20 2013
LINKS
Index entries for linear recurrences with constant coefficients, signature (0,0,0,1).
FORMULA
For generating the summation sequence s, start with s(0) = 1, and a(0) = 8.
For positive values of s(n): s(n+1) = s(n) + a(n).
For negative values of s(n): s(n-1) = s(n) - a(n-1). Here, n is negative.
All values of a(n) are positive regardless of index. For example: a(-1) = a(-2) = a(-3) = 4; a(-4) = 8. Thus the simple pattern of a(n) and the simple arithmetic for generating s(n), are maintained across the n=0 boundary, in a manner similar to extending Fibonacci numbers to negative indices.
From Bruno Berselli, Jun 20 2013: (Start)
G.f.: 4*(2+x+x^2+x^3)/((1-x)*(1+x)*(1+x^2)).
a(n) = 4 + (1 + (-1)^n)*(1 + I^(n*(n+1))). (End)
From Wesley Ivan Hurt, Jul 09 2016: (Start)
a(n) = a(n-4) for n>3.
a(n) = 5 + I^(2*n) + I^(-n) + I^n.
a(n) = 5 + cos(n*Pi) + 2*cos(n*Pi/2) + I*sin(n*Pi). (End)
EXAMPLE
s(1) = 9, s(2) = 13, s(3) = 17, s(4) = 21, s(5) = 29, s(6) = 33, s(7) = 37.
s(-1) = -3, s(-2) = -7, s(-3) = -11, s(-4) = -19, s(-5) = -23, s(-6) = -27, s(-7) = -31.
MAPLE
seq(op([8, 4, 4, 4]), n=0..40); # Wesley Ivan Hurt, Jul 09 2016
MATHEMATICA
Flatten[Table[{8, 4, 4, 4}, {20}]] (* Bruno Berselli, Jun 20 2013 *)
PROG
(Magma) &cat [[8, 4, 4, 4]^^30]; // Wesley Ivan Hurt, Jul 09 2016
(PARI) a(n)=if(n%4, 4, 8) \\ Charles R Greathouse IV, Jul 17 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Richard R. Forberg, Jun 01 2013
EXTENSIONS
Simpler name from Joerg Arndt, Jun 16 2013
STATUS
approved