cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A099837 Expansion of (1 - x^2) / (1 + x + x^2) in powers of x.

Original entry on oeis.org

1, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1
Offset: 0

Views

Author

Paul Barry, Oct 27 2004

Keywords

Comments

A transform of (-1)^n.
Row sums of Riordan array ((1-x)/(1+x), x/(1+x)^2), A110162.
Let b(n) = Sum_{k=0..floor(n/2)} binomial(n-k,k)(-1)^(n-2k). Then a(n) = b(n) - b(n-2) = A049347(n) - A049347(n-2) (n > 0). The g.f. 1/(1+x) of (-1)^n is transformed to (1-x^2)/(1+x+x^2) under the mapping G(x)->((1-x^2)/(1+x^2))G(x/(1+x^2)). Partial sums of A099838.
A(n) = a(n+3) (or a(n) if a(0) is replaced by 2) appears, together with B(n) = A049347(n) in the formula 2*exp(2*Pi*n*i/3) = A(n) + B(n)*sqrt(3)*i, n >= 0, with i = sqrt(-1). See A164116 for the case N=5. - Wolfdieter Lang, Feb 27 2014

Examples

			G.f. = 1 - x - x^2 + 2*x^3 - x^4 - x^5 + 2*x^6 - x^7 - x^8 + 2*x^9 - x^10 + ...
		

Crossrefs

Programs

  • Maple
    A099837 := proc(n)
        option remember;
        if n <=2 then
            op(n+1,[1,-1,-1]) ;
        else
            -procname(n-1)-procname(n-2) ;
        end if;
    end proc:
    seq(A099837(n),n=0..80) ; # R. J. Mathar, Apr 26 2022
  • Mathematica
    a[0] = 1; a[n_] := Mod[n+2, 3] - Mod[n, 3]; A099837 = Table[a[n], {n, 0, 71}](* Jean-François Alcover, Feb 15 2012, after Michael Somos *)
    LinearRecurrence[{-1, -1}, {1, -1, -1}, 50] (* G. C. Greubel, Aug 08 2017 *)
  • Maxima
    A099837(n) := block(
            if n = 0 then 1 else [2,-1,-1][1+mod(n,3)]
    )$ /* R. J. Mathar, Mar 19 2012 */
    
  • PARI
    {a(n) = [2, -1, -1][n%3 + 1] - (n == 0)}; /* Michael Somos, Jan 19 2012 */
    
  • PARI
    Vec((1-x^2)/(1+x+x^2) + O(x^20)) \\ Felix Fröhlich, Aug 08 2017

Formula

G.f.: (1-x^2)/(1+x+x^2).
Euler transform of length 3 sequence [-1, -1, 1]. - Michael Somos, Mar 21 2011
Moebius transform is length 3 sequence [-1, 0, 3]. - Michael Somos, Mar 22 2011
a(n) = -b(n) where b(n) = A061347(n) is multiplicative with b(3^e) = -2 if e > 0, b(p^e) = 1 otherwise. - Michael Somos, Jan 19 2012
a(n) = a(-n). a(n) = c_3(n) if n > 1 where c_k(n) is Ramanujan's sum. - Michael Somos, Mar 21 2011
G.f.: (1 - x) * (1 - x^2) / (1 - x^3). a(n) = -a(n-1) - a(n-2) unless n = 0, 1, 2. - Michael Somos, Jan 19 2012
Dirichlet g.f.: Sum_{n>=1} a(n)/n^s = zeta(s)*(3^(1-s)-1). - R. J. Mathar, Apr 11 2011
a(n+3) = R(n,-1) for n >= 0, with the monic Chebyshev T-polynomials R with coefficient table A127672. - Wolfdieter Lang, Feb 27 2014
For n > 0, a(n) = 2*cos(n*Pi/3)*cos(n*Pi). - Wesley Ivan Hurt, Sep 25 2017
From Peter Bala, Apr 20 2024: (Start)
a(n) is equal to the n-th order Taylor polynomial (centered at 0) of 1/c(x)^(2*n) evaluated at x = 1, where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the o.g.f. of the Catalan numbers A000108. Cf. A333093.
Row sums of the Riordan array A110162. (End)

A001350 Associated Mersenne numbers.

Original entry on oeis.org

0, 1, 1, 4, 5, 11, 16, 29, 45, 76, 121, 199, 320, 521, 841, 1364, 2205, 3571, 5776, 9349, 15125, 24476, 39601, 64079, 103680, 167761, 271441, 439204, 710645, 1149851, 1860496, 3010349, 4870845, 7881196, 12752041, 20633239, 33385280, 54018521, 87403801, 141422324
Offset: 0

Views

Author

Keywords

Comments

a(n) is last term in the period of the continued fraction expansion of phi^n (phi being the golden number). E.g.: n=10, phi^10=[122,1,121,1,121,1,121,...] (and the period may only have 1 or 2 terms). Also, a(n) = floor(phi^n)-((n+1) mod 2), or a(n) = A014217(n)-((n+1) mod 2). - Thomas Baruchel, Nov 05 2002 [continued fraction value corrected by Jon E. Schoenfield, Jan 20 2019]
a(n) is the resultant of the polynomials x^2-x-1 and x^(n+1)-x^n-1 for n >= 1. - Richard Choulet, Aug 05 2007
This is a divisibility sequence; that is, if n divides m, then a(n) divides a(m). - Michael Somos, Feb 12 2012
Gives the number of arrangements of black and white beads on a necklace with a total of n beads satisfying (1) there is at least one black bead (2) between any two black beads the number of white beads is even and (3) rotations and flippings of a necklace are considered distinct (see Butler). - Peter Bala, Mar 06 2014
This is the case P1 = 1, P2 = 0, Q = -1 of the 3-parameter family of 4th-order linear divisibility sequences found by Williams and Guy. - Peter Bala, Mar 31 2014
The resultant of the (s_2, s_2+n) pair, where s_n(X) is X^n-X-1, is -a(n). See Rush link. - Michel Marcus, Sep 30 2019

Examples

			G.f. = x + x^2 + 4*x^3 + 5*x^4 + 11*x^5 + 16*x^6 + 29*x^7 + 45*x^8 + 76*x^9 + ...
n=1: a(9)/a(3) = 76/4 = 19; a(18)/a(6) = 5776/16 = 361 = 19^2. - _Bob Selcoe_, Jun 01 2014
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [Floor(-(1 - ((1 + Sqrt(5))/2)^n - (-(1 + Sqrt(5))/2)^(-n) + (-1)^n)): n in [0..40]]; // Vincenzo Librandi, Aug 15 2011
    
  • Maple
    A001350 := n -> add(binomial(k-1, 2*k-n)*n/(n-k), k=0..n-1);
    seq(A001350(n), n=0..39); # Peter Luschny, Sep 26 2014
  • Mathematica
    Clear[f, n]; f[n_] = -(1 - ((1 + Sqrt[5])/2)^n - (-(1 + Sqrt[5])/2)^(-n) + (-1)^n); Table[FullSimplify[ExpandAll[f[n]]], {n, 0, 30}] (* Roger L. Bagula and Gary W. Adamson, Nov 26 2008 *)
    a[ n_] := LucasL[n] - 1 - (-1)^n; (* Michael Somos, May 18 2015 *)
    a[ n_] := SeriesCoefficient[ x D[ Log[ 1 + x / (1 - x - x^2)], x], {x, 0, n}]; (* Michael Somos, May 18 2015 *)
    LinearRecurrence[{1, 2, -1, -1}, {0, 1, 1, 4}, 40] (* Jean-François Alcover, Jan 07 2019 *)
  • PARI
    {a(n) = fibonacci(n+1) + fibonacci(n-1) - 1 - (-1)^n};
    
  • PARI
    {a(n) = my(w = quadgen(5)); simplify( -(w^n - 1) * ((-1/w)^n - 1))}; /* Michael Somos, Feb 12 2012 */
    
  • Python
    from sympy import lucas
    def A001350(n): return lucas(n)-((n&1^1)<<1) # Chai Wah Wu, Sep 23 2023

Formula

G.f.: x*(1+x^2)/((1-x^2)*(1-x-x^2)). - Simon Plouffe in his 1992 dissertation
a(n) = a(n-1) + a(n-2) + 1 -(-1)^n. a(-n) = (-1)^n * a(n).
a(n) = A050140(Fibonacci(n)). - Thomas Baruchel, Nov 05 2002
Convolution of F(n) and {1, 0, 2, 0, 2, ...}. a(n) = Sum_{k=0..n} ((1+(-1)^k)-0^k)*F(n-k) = Sum_{k=0..n} F(k)*((1+(-1)^(n-k))-0^(n-k)). - Paul Barry, Jul 19 2004
a(n) = 2*A074331(n) - A000045(n). - Paul Barry, Jul 19 2004
a(n) = Lucas_number(n) - 1 - (-1)^n = A000032(n) - 1 - (-1)^n. - Hieronymus Fischer, Feb 18 2006
a(n) = -(1 - ((1 + sqrt(5))/2)^n - (-(1 + sqrt(5))/2)^(-n) + (-1)^n). - Roger L. Bagula and Gary W. Adamson, Nov 26 2008
a(n) = n * Sum_{k=1..n} (Sum_{i=ceiling((n-k)/2)..(n-k)} (binomial(i,n-k-i)*binomial(k+i-1,k-1))/k*(-1)^(k+1)), n>0. - Vladimir Kruchinin, Sep 03 2010
a(n) = a(n-1) + 2*a(n-2) - a(n-3) - a(n-4). - Colin Barker, Apr 11 2014
a(n) = sqrt(A152152(n)). - Colin Barker, Apr 11 2014
a(n) = a(2*n)/A000032(n) when n is odd; a(n) = a(2*n)/(A000032(n+2)) when n is even. - Bob Selcoe, Jun 01 2014
a(12n+6)/a(4n+2) = (a(6n+3)/a(2n+1))^2. - Bob Selcoe, Jun 01 2014
a(n) = Sum_{k=0..n-1} binomial(k-1, 2*k-n)*n/(n-k). - Peter Luschny, Sep 26 2014
From Peter Bala, Mar 19 2015: (Start)
a(n) = -(alpha^n - 1)*(beta^n - 1), where alpha = 1/2*(1 + sqrt(5)) and beta = (1/2)*(1 - sqrt(5)).
a(n) = -det(I - M^n) where I is the 2 X 2 identity matrix and M = [ 1, 1; 1, 0 ]. Cf. A129744.
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + Sum_{n >= 1} Fibonacci(n)*x^n. Cf. A004146. (End)
a(n) = A052952(n-1) + A052952(n-3). - R. J. Mathar, Jul 02 2018
a(n) = (L(2*n+1) - L(n+1)) mod (L(n+1)-1) for n > 0 where L(k)=A000032(k). - Art Baker, Jan 17 2019
a(n) = Sum_{j=n..2*n-1} L(j) mod Sum_{j=0..n-1} L(j) where L(j)=A000032(j). - Art Baker, Jan 20 2019
Convolution of (1, 0, 3, 0, 5, 0, 7, ...) and (1, 1, 1, 2, 3, 5, 8, 13, ...). - Gary W. Adamson, Jul 08 2019
a(n) = Sum_{d|n} d*A060280(d) = Sum_{d|n} A031367(d). [Baake, Roberts, Weiss, eq(2)]. - R. J. Mathar, Oct 19 2021

Extensions

Additional comments from Michael Somos, Aug 01 2002

A117179 Riordan array ((1-x^2)/(1+x^2)^2,x/(1+x^2)).

Original entry on oeis.org

1, 0, 1, -3, 0, 1, 0, -4, 0, 1, 5, 0, -5, 0, 1, 0, 9, 0, -6, 0, 1, -7, 0, 14, 0, -7, 0, 1, 0, -16, 0, 20, 0, -8, 0, 1, 9, 0, -30, 0, 27, 0, -9, 0, 1, 0, 25, 0, -50, 0, 35, 0, -10, 0, 1, -11, 0, 55, 0, -77, 0, 44, 0, -11, 0, 1
Offset: 0

Views

Author

Paul Barry, Mar 01 2006

Keywords

Comments

Inverse is A117178. Row sums are (-1)^n*A098554(n+1). Diagonal sums are 1,0,-2,0,2,0,-2,... with g.f. (1-x^2)/(1+x^2).
Apparently, with the rows de-aerated and then reversed, this matrix becomes signed A034807 with the twos on the diagonal removed. Apparently, |D(2n,k+1)| = |D(2(n-1),k+1)| + |D(2n,k)| where D(n,k) is the k-th element on the n-th diagonal. - Tom Copeland, May 30 2017

Examples

			Triangle begins
1,
0, 1,
-3, 0, 1,
0, -4, 0, 1,
5, 0, -5, 0, 1,
0, 9, 0, -6, 0, 1,
-7, 0, 14, 0, -7, 0, 1,
0, -16, 0, 20, 0, -8, 0, 1
		
Showing 1-3 of 3 results.