OFFSET
1,2
COMMENTS
Perrin-like prime-divisibility sequence, but based upon template 7=5+2 in place of 5=3+2.
1. Apparently identical to A007387 but for latter's third term 3. 2. Attention directed to remainder upon division of a term by its (composite) argument, when latter =1 or 5 (mod 6). Possible factorization tool for impostor candidate primes? 3. Recurrence period, any length-five string of term values (mod 6) found in the sequence: 2^3*13*31, to Perrin's three-term period of 7*13. Note 13= 2*6+1, 31 = 5*6+1. 4. Query: Smallest pseudoprime >9. 5. Query: Closed form for n-th term.
Semiprimes a= 9, 14, 34, 57, 91 etc. are at the indices n=9, 12, 16, 17, 19, 21, 24, 25, 26, 31, 32, 40, 44, 45, 51, 53, 59, 66, 72, 76, 80, 110 etc. - R. J. Mathar, Nov 24 2007
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (0,1,0,0,1).
FORMULA
O.g.f.: -x*(2+5*x^3)/(-1+x^2+x^5). - R. J. Mathar, Nov 24 2007
Rewritten, Mathar's o.g.f. resembles a logarithmic derivative: -(5*x^4 + 2*x) / (x^5 +x^2-1). Any significance? - G. Reed Jameson (Reedjameson(AT)yahoo.com), Dec 13 2007, Dec 16 2007
a(-n) = A136598(n).
MATHEMATICA
LinearRecurrence[{0, 1, 0, 0, 1}, {0, 2, 0, 2, 5}, 60] (* Harvey P. Dale, Oct 21 2015 *)
PROG
(PARI) {a(n) = if( n<0, n = 1 - n; polsym(x^5 + x^2 - 1, n)[n], n++; polsym(x^5 - x^3 - 1, n)[n])} /* Michael Somos, Feb 12 2012 */
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
G. Reed Jameson (Reedjameson(AT)yahoo.com), Nov 23 2007
EXTENSIONS
More terms from R. J. Mathar, Nov 24 2007
STATUS
approved