OFFSET
0,5
COMMENTS
a(n-2) is the total degree of the irreducible factor F(n) of the n-th Somos polynomial. - Michael Somos, Jul 06 2011
LINKS
Matthew House, Table of n, a(n) for n = 0..10000
Michael Somos, Somos Polynomials.
Index entries for linear recurrences with constant coefficients, signature (2,-1,0,0,1,-2,1).
FORMULA
G.f.: x^3*(x^2+1)/ ((1-x)^3 * (1+x+x^2+x^3+x^4)). a(n) = +2*a(n-1) -a(n-2) +a(n-5) -2*a(n-6) +a(n-7). - R. J. Mathar, Apr 15 2010
Euler transform of length 5 sequence [2, 1, 0, -1, 1]. - Michael Somos, Jul 04 2011
a(1-n) = a(n). a(n) = a(n-5) + 2*n - 6 for all n in Z. - Michael Somos, Jul 04 2011
a(n) = a(n-1) + a(n-5) - a(n-6) + 2 for all n in Z. - Michael Somos, Jul 06 2011
a(n) = (1/5) * ( n^2 - n + [0,0,-2,-1,-2](mod 5) ). - Ralf Stephan, Aug 11 2013
a(n) - 2*a(n+1) + a(n+2) = (n == 1 (mod 5)) + (n == 3 (mod 5)) for all n in Z. - Michael Somos, Oct 19 2014
Sum_{n>=3} 1/a(n) = 50/9 - sqrt(2*(5+sqrt(5)))*Pi/3 + tan(Pi/(2*sqrt(5)))*Pi/sqrt(5). - Amiram Eldar, Oct 01 2022
EXAMPLE
G.f. = x^3 + 2*x^4 + 4*x^5 + 6*x^6 + 8*x^7 + 11*x^8 + 14*x^9 + 18*x^10 + 22*x^11 + ...
F(5) = y + 1 is of degree a(3) = 1, F(6) = y*z + y + z is of degree a(4) = 2.
MATHEMATICA
a[ n_] := Quotient[ n (n - 1), 5]; (* Michael Somos, Oct 19 2014 *)
PROG
(PARI) {a(n) = n * (n - 1) \ 5}; /* Michael Somos, Jul 04 2011 */
(Magma) [Floor(n*(n-1)/5): n in [0..50]]; // G. C. Greubel, Oct 28 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved