OFFSET
1,2
COMMENTS
Conjecture: Apart from 0, and the further exclusions noted below, the sequence gives the values of c/6 such that an infinite number of primes, p, result in both p^3+c and p^3-c being positive primes. Taking the complement we say: the excluded c/6 values are {1,6} mod 7. See A005097 for a conjecture on the modulo patterns of excluded c/6 values for the general case of p^q + c and p^q - c both prime, for any q > 0, and see A047222 for q=2. Note that polynomial factorization also excludes a few c/6 values. This occurs here when c is an even cube (A016743), which requires a further exclusion of certain c/6 values in this sequence when (6c)^3/6 == 0 (mod 7), or c/6 = {0, 12348, 98784, ...}. - Richard R. Forberg, Jun 28 2016
LINKS
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,1,-1).
FORMULA
G.f.: x^2*(2*x^2 + 3*x + 2)*(x^2 - x + 1) / ( (x^4 + x^3 + x^2 + x + 1)*(x-1)^2 ). - R. J. Mathar, Dec 04 2011
a(n) = a(n-1) + a(n-5) - a(n-6). - Wesley Ivan Hurt, Sep 03 2022
MATHEMATICA
Table[7 n + {0, 2, 3, 4, 5}, {n, 0, 12}] // Flatten (* or *)
Select[Range[0, 79], ! MemberQ[{1, 6}, Mod[#, 7]] &] (* or *)
Rest@ CoefficientList[Series[x^2 (2 x^2 + 3 x + 2) (x^2 - x + 1)/((x^4 + x^3 + x^2 + x + 1) (x - 1)^2), {x, 0, 57}], x] (* Michael De Vlieger, Jul 25 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved