The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A285640 a(n) = negative of exponent of s(2) in factorization of s(n) where s(1)^2 * s(n+2) * s(n-2) = s(2)^2 * s(n+1) * s(n-1) - s(3) * s(1) * s(n)^2, s(n) = -s(-n) for all n in Z,  with a(0) = -2, and s(1), s(2), s(3), s(4) are variables. 2
 -2, 0, -1, 0, 0, 0, 1, 0, 2, 1, 3, 3, 4, 5, 5, 7, 6, 9, 9, 11, 12, 13, 15, 15, 18, 18, 21, 22, 24, 26, 27, 30, 30, 34, 35, 38, 40, 42, 45, 46, 50, 51, 55, 57, 60, 63, 65, 69, 70, 75, 77, 81, 84, 87, 91, 93, 98, 100, 105, 108, 112, 116, 119, 124, 126, 132, 135 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,1 COMMENTS The sequence s(n) is a generalized Somos-4 sequence. Also, s(n) is a Laurent polynomial in the variables s(1), s(2), s(3), s(4) with only powers of s(1), s(2) appearing in the denominator. LINKS G. C. Greubel, Table of n, a(n) for n = 0..5000 FORMULA a(n) = a(n+16) - n - 8 = a(-n) for all n in Z. G.f.: (-2 + 4*x - 3*x^2 + 2*x^3 - x^4 + x^6 - 2*x^7 + 3*x^8 - 3*x^9 + 3*x^10 - 2*x^11 + x^12 - x^14 + 2*x^15 - x^16) / ((1 - x)^2 * (1 - x^16)). EXAMPLE G.f. = -2 - x^2 + x^6 + 2*x^8 + x^9 + 3*x^10 + 3*x^11 + 4*x^12 + 5*x^13 + ... Let s(1) = s1, s(2) = s2, s(3) = s3, s(4) = s4. Then s(5) = (s4*s2^3 -s3^3*s1) / s1^3 gives a(5) = 0, s(6) = s3*(s4*s2^5 -s4^2*s1^4 -s1*s2^2*s3^3) / (s2*s1^5) gives a(6) = 1, s(7) = (-s4^3*s1^3*s2 +s4*s2^3*s3^3 -s1*s3^6) / s1^6 gives a(7) = 0. MATHEMATICA a[ n_] := (n^2 - {1, 6, 3, 4, 5, 2, 7, 0, 7, 2, 5, 4, 3, 6, 1, 8}[[Mod[n, 16, 1]]]^2 ) / 32; PROG (PARI) {a(n) = (n^2 - [8, 1, 6, 3, 4, 5, 2, 7, 0, 7, 2, 5, 4, 3, 6, 1][n%16+1]^2) / 32}; (PARI) {a(n) = my(m=abs(n)); polcoeff( (-2 + 4*x - 3*x^2 + 2*x^3 - x^4 + x^6 - 2*x^7 + 3*x^8 - 3*x^9 + 3*x^10 - 2*x^11 + x^12 - x^14 + 2*x^15 - x^16) / ((1 - x)^2 * (1 - x^16)) + x * O(x^m), m)}; (MAGMA) m:=60; R:=PowerSeriesRing(Integers(), m); Coefficients(R!((-2 +4*x-3*x^2+2*x^3-x^4+x^6-2*x^7+3*x^8-3*x^9+3*x^10-2*x^11+x^12-x^14 + 2*x^15-x^16)/((1-x)^2*(1-x^16)))); // G. C. Greubel, Jul 29 2018 CROSSREFS Cf. A285631. Sequence in context: A215935 A270573 A096271 * A231189 A219558 A279210 Adjacent sequences:  A285637 A285638 A285639 * A285641 A285642 A285643 KEYWORD sign,easy AUTHOR Michael Somos, Apr 23 2017 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified May 21 06:36 EDT 2022. Contains 353889 sequences. (Running on oeis4.)