login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A285631 a(n) = negative of exponent of s(1) 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) = -3, and s(1), s(2), s(3), s(4) are variables. 2
-3, -1, 0, 0, 0, 3, 5, 6, 8, 12, 15, 17, 21, 26, 30, 33, 39, 45, 50, 54, 62, 69, 75, 81, 90, 98, 105, 113, 123, 132, 140, 150, 161, 171, 180, 192, 204, 215, 225, 239, 252, 264, 276, 291, 305, 318, 332, 348, 363, 377, 393, 410, 426, 441, 459, 477, 494, 510, 530 (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+19) - 6*n - 57 = a(-n) for all n in Z.

0 = (a(n) - 2*a(n-1) + a(n-2) + 1) * (a(n) - 2*a(n+1) + a(n+2) + 1) for all n in Z.

G.f.: (-3 + 5*x - x^2 - x^3 + 3*x^5 - x^6 - x^7 + x^8 + 2*x^9 - x^10 - x^11 + 2*x^12 + x^13 - x^14 - x^15 + 3*x^16 - x^18 + 2*x^19 - x^20) / ((1 - x)^2 * (1 - x^19)).

EXAMPLE

G.f. = -3 - x + 3*x^5 + 5*x^6 + 6*x^7 + 8*x^8 + 12*x^9 + 15*x^10 + 17*x^11 + ...

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) = -3, s(6) = s3*(s4*s2^5 -s4^2*s1^4 -s1*s2^2*s3^3) / (s2*s1^5) gives a(6) = -5, and s(7) = (-s4^3*s1^3*s2 +s4*s2^3*s3^3 -s1*s3^6) / s1^6 gives a(7) = 6.

MATHEMATICA

a[ n_] := (3 n^2 + {35, 45, 30, 9, 39, 44, 24, 17, 42, 42, 17, 24, 44, 39, 9, 30, 45, 35, 0}[[Mod[n, 19, 1]]]) / 19 - 3;

PROG

(PARI) {a(n) = (3*n^2 + [0, 35, 45, 30, 9, 39, 44, 24, 17, 42, 42, 17, 24, 44, 39, 9, 30, 45, 35] [n%19+1]) / 19 - 3};

(PARI) {a(n) = my(m=abs(n)); polcoeff( (-3 + 5*x - x^2 - x^3 + 3*x^5 - x^6 - x^7 + x^8 + 2*x^9 - x^10 - x^11 + 2*x^12 + x^13 - x^14 - x^15 + 3*x^16 - x^18 + 2*x^19 - x^20) / ((1 - x)^2 * (1 - x^19)) + x * O(x^m), m)};

(MAGMA) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!((-3 +5*x-x^2-x^3+3*x^5-x^6-x^7+x^8+2*x^9-x^10-x^11+2*x^12+x^13-x^14-x^15 + 3*x^16-x^18+2*x^19-x^20)/((1-x)^2*(1-x^19)))); // G. C. Greubel, Jul 29 2018

CROSSREFS

Cf. A285640.

Sequence in context: A193525 A049828 A286131 * A316836 A058612 A099725

Adjacent sequences:  A285628 A285629 A285630 * A285632 A285633 A285634

KEYWORD

sign,easy

AUTHOR

Michael Somos, Apr 22 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 | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 23 16:52 EDT 2019. Contains 321432 sequences. (Running on oeis4.)