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.

Previous Showing 31-40 of 133 results. Next

A054878 Number of closed walks of length n along the edges of a tetrahedron based at a vertex.

Original entry on oeis.org

1, 0, 3, 6, 21, 60, 183, 546, 1641, 4920, 14763, 44286, 132861, 398580, 1195743, 3587226, 10761681, 32285040, 96855123, 290565366, 871696101, 2615088300, 7845264903, 23535794706, 70607384121, 211822152360, 635466457083
Offset: 0

Views

Author

Paolo Dominici (pl.dm(AT)libero.it), May 23 2000

Keywords

Comments

Number of closed walks of length n at a vertex of C_4, the cyclic graph on 4 nodes. 3*A015518(n) + a(n) = 3^n. - Paul Barry, Feb 03 2004
Form the digraph with matrix A = [0,1,1,1; 1,0,1,1; 1,1,0,1; 1,0,1,1]; a(n) corresponds to the (1,1) term of A^n. - Paul Barry, Oct 02 2004
Absolute values of A084567 (compare generating functions).
For n > 1, 4*a(n)=A218034(n)= the trace of the n-th power of the adjacency matrix for a complete 4-graph, a 4 X 4 matrix with a null diagonal and all ones for off-diagonal elements. The diagonal elements for the n-th power are a(n) and the off-diagonal are a(n)+1 for an odd power and a(n)-1 for an even (cf. A001045). - Tom Copeland, Nov 06 2012

Crossrefs

Row n=4 of A109502. A084567 (signed version).
{a(n)/3} for n>0 is A015518, non-closed walks.
Cf. A001045, A078008, A097073, A115341, A015518 (sequences where a(n)=3^n-a(n-1)). - Vladimir Joseph Stephan Orlovsky, Dec 11 2008

Programs

  • Magma
    [(3^n+(-1)^n*3)/4: n in [0..35]]; // Vincenzo Librandi, Jun 30 2011
    
  • Maple
    A054878:=n->(3^n + (-1)^n*3)/4: seq(A054878(n), n=0..50); # Wesley Ivan Hurt, Sep 16 2017
  • Mathematica
    Table[(3^n + (-1)^n*3)/4, {n, 0, 26}] (* or *)
    CoefficientList[Series[1/4*(3/(1 + x) + 1/(1 - 3 x)), {x, 0, 26}], x] (* Michael De Vlieger, Sep 15 2017 *)
  • PARI
    a(n) = (3^n + 3*(-1)^n)/4; \\ Altug Alkan, Sep 17 2017

Formula

a(n) = (3^n + (-1)^n*3)/4.
G.f.: 1/4*(3/(1+x) + 1/(1-3*x)).
E.g.f.: (exp(3*x) + 3*exp(-x))/4. - Paul Barry, Apr 20 2003
a(n) = 3^n - a(n-1) with a(0)=0. - Labos Elemer, Apr 26 2003
G.f.: (1 - 3*x^2 - 2*x^3)/(1 - 6*x^2 - 8*x^3 - 3*x^4) = (1 - 3*x^2 - 2*x^3)/charpoly(adj(C_4)). - Paul Barry, Feb 03 2004
From Paul Barry, Oct 02 2004: (Start)
G.f.: (1-2*x)/(1 - 2*x - 3*x^2).
a(n) = 2*a(n-1) + 3*a(n-2). (End)
G.f.: 1 - x + x/Q(0), where Q(k) = 1 + 3*x^2 - (3*k+4)*x + x*(3*k+1 - 3*x)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Oct 07 2013
a(n+m) = a(n)*a(m) + a(n+1)*a(m+1)/3. - Yuchun Ji, Sep 12 2017
a(n) = 3*a(n-1) + 3*(-1)^n. - Yuchun Ji, Sep 13 2017
From Peter Bala, May 28 2024: (Start)
a(n) = (-1)^n + Sum_{k = 1..n} (-1)^(n-k)*binomial(n, k)*4^(k-1).
G.f.: A(x) = 1/(1 - x^2) o 1/(1 - x^2), where o denotes the black diamond product of power series as defined by Dukes and White. Cf. A015575.
The black diamond product A(x) o A(x) is the g.f. for the number of closed walks of length n at a vertex along the edges of the 15-simplex. (End)

A108561 Triangle read by rows: T(n,0)=1, T(n,n)=(-1)^n, T(n+1,k)=T(n,k-1)+T(n,k) for 0 < k < n.

Original entry on oeis.org

1, 1, -1, 1, 0, 1, 1, 1, 1, -1, 1, 2, 2, 0, 1, 1, 3, 4, 2, 1, -1, 1, 4, 7, 6, 3, 0, 1, 1, 5, 11, 13, 9, 3, 1, -1, 1, 6, 16, 24, 22, 12, 4, 0, 1, 1, 7, 22, 40, 46, 34, 16, 4, 1, -1, 1, 8, 29, 62, 86, 80, 50, 20, 5, 0, 1, 1, 9, 37, 91, 148, 166, 130, 70, 25, 5, 1, -1, 1, 10, 46, 128, 239, 314, 296, 200, 95, 30, 6, 0, 1, 1, 11, 56, 174, 367
Offset: 0

Views

Author

Reinhard Zumkeller, Jun 10 2005

Keywords

Comments

Sum_{k=0..n} T(n,k) = A078008(n);
Sum_{k=0..n} abs(T(n,k)) = A052953(n-1) for n > 0;
T(n,1) = n - 2 for n > 1;
T(n,2) = A000124(n-3) for n > 2;
T(n,3) = A003600(n-4) for n > 4;
T(n,n-6) = A001753(n-6) for n > 6;
T(n,n-5) = A001752(n-5) for n > 5;
T(n,n-4) = A002623(n-4) for n > 4;
T(n,n-3) = A002620(n-1) for n > 3;
T(n,n-2) = A008619(n-2) for n > 2;
T(n,n-1) = n mod 2 for n > 0;
T(2*n,n) = A072547(n+1).
Sum_{k=0..n} T(n,k)*x^k = A232015(n), A078008(n), A000012(n), A040000(n), A001045(n+2), A140725(n+1) for x = 2, 1, 0, -1, -2, -3 respectively. - Philippe Deléham, Nov 17 2013, Nov 19 2013
(1,a^n) Pascal triangle with a = -1. - Philippe Deléham, Dec 27 2013
T(n,k) = A112465(n,n-k). - Reinhard Zumkeller, Jan 03 2014

Examples

			From _Philippe Deléham_, Nov 17 2013: (Start)
Triangle begins:
  1;
  1, -1;
  1,  0,  1;
  1,  1,  1, -1;
  1,  2,  2,  0,  1;
  1,  3,  4,  2,  1, -1;
  1,  4,  7,  6,  3,  0,  1; (End)
		

Crossrefs

Cf. A007318 (a=1), A008949(a=2), A164844(a=10).
Similar to the triangles A035317, A059259, A080242, A112555.
Cf. A072547 (central terms).

Programs

  • GAP
    Flat(List([0..13],n->List([0..n],k->Sum([0..k],i->Binomial(n,i)*(-2)^(k-i))))); # Muniru A Asiru, Feb 19 2018
  • Haskell
    a108561 n k = a108561_tabl !! n !! k
    a108561_row n = a108561_tabl !! n
    a108561_tabl = map reverse a112465_tabl
    -- Reinhard Zumkeller, Jan 03 2014
    
  • Maple
    A108561 := (n, k) -> add(binomial(n, i)*(-2)^(k-i), i = 0..k):
    seq(seq(A108561(n,k), k = 0..n), n = 0..12); # Peter Bala, Feb 18 2018
  • Mathematica
    Clear[t]; t[n_, 0] = 1; t[n_, n_] := t[n, n] = (-1)^Mod[n, 2]; t[n_, k_] := t[n, k] = t[n-1, k] + t[n-1, k-1]; Table[t[n, k], {n, 0, 13}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 06 2013 *)
  • Sage
    def A108561_row(n):
        @cached_function
        def prec(n, k):
            if k==n: return 1
            if k==0: return 0
            return -prec(n-1,k-1)-sum(prec(n,k+i-1) for i in (2..n-k+1))
        return [(-1)^k*prec(n, k) for k in (1..n-1)]+[(-1)^(n+1)]
    for n in (1..12): print(A108561_row(n)) # Peter Luschny, Mar 16 2016
    

Formula

G.f.: (1-y*x)/(1-x-(y+y^2)*x). - Philippe Deléham, Nov 17 2013
T(n,k) = T(n-1,k) + T(n-2,k-1) + T(n-2,k-2), T(0,0)=T(1,0)=1, T(1,1)=-1, T(n,k)=0 if k < 0 or if k > n. - Philippe Deléham, Nov 17 2013
From Peter Bala, Feb 18 2018: (Start)
T(n,k) = Sum_{i = 0..k} binomial(n,i)*(-2)^(k-i), 0 <= k <= n.
The n-th row polynomial is the n-th degree Taylor polynomial of the rational function (1 + x)^n/(1 + 2*x) about 0. For example, for n = 4, (1 + x)^4/(1 + 2*x) = 1 + 2*x + 2*x^2 + x^4 + O(x^5). (End)

Extensions

Definition corrected by Philippe Deléham, Dec 26 2013

A007613 a(n) = (8^n + 2*(-1)^n)/3.

Original entry on oeis.org

1, 2, 22, 170, 1366, 10922, 87382, 699050, 5592406, 44739242, 357913942, 2863311530, 22906492246, 183251937962, 1466015503702, 11728124029610, 93824992236886, 750599937895082, 6004799503160662, 48038396025285290, 384307168202282326, 3074457345618258602
Offset: 0

Views

Author

Keywords

Comments

Also, the cogrowth sequence of C3 X C3 = ; that is, the number of words of length 3n that reduce to the identity. - Sean A. Irvine, Nov 04 2024

References

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

Crossrefs

Programs

Formula

a(n) = A078008(3*n). - Paul Barry, Nov 29 2003
From Paul Barry, Mar 24 2004: (Start)
a(n) = (A082311(n) + (-1)^n)/2.
a(n) = (A001045(3*n+1) + (-1)^n)/2. (End)
a(n) = Sum_{k=0..n} binomial(3*n, 3*k). - Paul Barry, Jan 13 2005
a(n) = 8*a(n-1) + 6*(-1)^n. - Paul Curtz, Nov 19 2007
From Colin Barker, Sep 29 2014: (Start)
a(n) = 7*a(n-1) + 8*a(n-2).
G.f.: (1-5*x) / ((1+x)*(1-8*x)). (End)
E.g.f.: (1/3)*(exp(8*x) + 2*exp(-x)). - G. C. Greubel, Apr 23 2023

Extensions

More terms from Colin Barker, Sep 29 2014

A084247 a(n) = -a(n-1) + 2*a(n-2), a(0)=1, a(1)=2.

Original entry on oeis.org

1, 2, 0, 4, -4, 12, -20, 44, -84, 172, -340, 684, -1364, 2732, -5460, 10924, -21844, 43692, -87380, 174764, -349524, 699052, -1398100, 2796204, -5592404, 11184812, -22369620, 44739244, -89478484, 178956972, -357913940, 715827884, -1431655764, 2863311532
Offset: 0

Views

Author

Paul Barry, May 23 2003

Keywords

Comments

Row sums of triangle in A112555. - Philippe Deléham, Sep 21 2009

Crossrefs

Programs

  • Magma
    [(4-(-2)^n)/3: n in [0..40]]; // Vincenzo Librandi, Aug 13 2011
    
  • Mathematica
    LinearRecurrence[{-1,2},{1,2},40] (* Harvey P. Dale, Nov 05 2016 *)
  • PARI
    Vec((1+3*x)/((1-x)*(1+2*x)) + O(x^40)) \\ Michel Marcus, Feb 25 2016
    
  • SageMath
    [(4-(-2)^n)/3 for n in range(41)] # G. C. Greubel, Apr 24 2023

Formula

Binomial transform of A084246.
a(n) = A077925(n-1) + 1.
a(n) = (4 - (-2)^n)/3.
G.f.: (1+3*x)/((1-x)*(1+2*x)).
E.g.f.: (4*exp(x) - exp(-2*x))/3.
a(n) = (-1)^(n+1)*(A078008(n+1) - A078008(n)). - Paul Curtz, Jun 30 2008, Feb 24 2021
a(n) = -2*a(n-1) + 4 . - Philippe Deléham, Sep 15 2009
a(n+1) = 2*A151575(n). - Philippe Deléham, Sep 21 2009
a(n) = A077925(n) + 3*A077925(n-1). - R. J. Mathar, Feb 24 2021

A130249 Maximal index k of a Jacobsthal number such that A001045(k)<=n (the 'lower' Jacobsthal inverse).

Original entry on oeis.org

0, 2, 2, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8
Offset: 0

Views

Author

Hieronymus Fischer, May 20 2007

Keywords

Comments

Inverse of the Jacobsthal sequence (A001045), nearly, since a(A001045(n))=n except for n=1 (see A130250 for another version). a(n)+1 is equal to the partial sum of the Jacobsthal indicator sequence (see A105348).

Examples

			a(12)=5, since A001045(5)=11<=12, but A001045(6)=21>12.
		

Crossrefs

For partial sums see A130251.
Other related sequences A130250, A130253, A105348, A001045, A130233, A130241.
Cf. A000523, A078008 (runlengths).

Programs

  • Magma
    [Floor(Log(3*n+1)/Log(2)): n in [0..30]]; // G. C. Greubel, Jan 08 2018
    
  • Mathematica
    Table[Floor[Log[2, 3*n + 1]], {n, 0, 50}] (* G. C. Greubel, Jan 08 2018 *)
  • PARI
    for(n=0, 30, print1(floor(log(3*n+1)/log(2)), ", ")) \\ G. C. Greubel, Jan 08 2018
    
  • PARI
    a(n) = logint(3*n+1, 2); \\ Ruud H.G. van Tol, May 12 2024
    
  • Python
    def A130249(n): return (3*n+1).bit_length()-1 # Chai Wah Wu, Jun 08 2022

Formula

a(n) = floor(log_2(3n+1)).
a(n) = A130250(n+1) - 1 = A130253(n) - 1.
G.f.: 1/(1-x)*(Sum_{k>=1} x^A001045(k)).
a(n) = A000523(3*n+1). - Ruud H.G. van Tol, May 12 2024

A131531 Period 6: repeat [0, 0, 1, 0, 0, -1].

Original entry on oeis.org

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

Views

Author

Paul Curtz, Aug 26 2007

Keywords

Comments

Also: partial sums of A092220 shifted by two indices. - R. J. Mathar, Feb 08 2008
From Paul Curtz, Jun 05 2011: (Start)
The square array of this sequence in the top row and further rows defined as first differences of preceding rows starts (see A167613):
. 0, 0, 1, 0, 0, -1, ...
. 0, 1, -1, 0, -1, 1, ... = A092220,
. 1, -2, 1, -1, 2, -1, ... = A131556,
. -3, 3, -2, 3, -3 2, ...
. 6, -5, 5, -6, 5, -5, ...
. -11, 10, -11, 11, -10, 11, ...
. 21, -21, 22, -21, 21, -22, ...
. -42, 43, -43, 42, -43, 43, ...
The main diagonal in this array is A001045; the first superdiagonal is the negated elements of A001045, the second superdiagonal is A078008.
The left column of the array is basically the inverse binomial transform, (-1)^n * A024495(n), assuming offset 0.
The second column of the array is A131708 with alternating signs, and the third column is A024493 with alternating signs (both assuming offset 0). (End)

Crossrefs

Programs

Formula

G.f.: x^3/(x+1)/(x^2-x+1). - R. J. Mathar, Nov 14 2007
a(n) = (-A057079(n+1) - (-1)^n)/3. - R. J. Mathar, Jun 13 2011
a(n) = -cos(Pi*(n-1)/3)/3 + sin(Pi*(n-1)/3)/sqrt(3) - (-1)^n/3. - R. J. Mathar, Oct 08 2011
a(n) = ( (-1)^n - (-1)^floor((n+2)/3) )/2. - Bruno Berselli, Jul 09 2013
a(n) + a(n-3) = 0 for n > 3. - Wesley Ivan Hurt, Jun 20 2016

Extensions

Edited by N. J. A. Sloane, Sep 15 2007

A109499 Number of closed walks of length n on the complete graph on 5 nodes from a given node.

Original entry on oeis.org

1, 0, 4, 12, 52, 204, 820, 3276, 13108, 52428, 209716, 838860, 3355444, 13421772, 53687092, 214748364, 858993460, 3435973836, 13743895348, 54975581388, 219902325556, 879609302220, 3518437208884, 14073748835532
Offset: 0

Views

Author

Mitch Harris, Jun 30 2005

Keywords

Crossrefs

Cf. A108020 (bisection), A109502.

Programs

  • GAP
    a:=[1,0];; for n in [3..30] do a[n]:=3*a[n-1]+4*a[n-2]; od; a; # G. C. Greubel, Mar 23 2019
  • Magma
    [(4^n + 4*(-1)^n)/5: n in [0..30]]; // Vincenzo Librandi, Aug 12 2011
    
  • Mathematica
    CoefficientList[Series[(1-3*x)/(1-3*x-4*x^2), {x,0,30}], x] (* or *) LinearRecurrence[{3,4}, {1,0}, 30] (* G. C. Greubel, Dec 30 2017 *)
  • PARI
    a(n)=(4^n+4*(-1)^n)/5 \\ Charles R Greathouse IV, Oct 01 2012
    
  • Sage
    [(4^n+4*(-1)^n)/5 for n in (0..30)] # G. C. Greubel, Mar 23 2019
    

Formula

G.f.: (1 - 3*x)/(1 - 3*x - 4*x^2).
a(n) = (4^n + 4*(-1)^n)/5.
a(n+1) = 4*A015521(n). - Paul Curtz, Nov 01 2009
a(n) = 3*a(n-1) + 4*a(n-1). - G. C. Greubel, Dec 30 2017
a(n) = A108020((n - 1) / 2) = 'ccc...c' (n digits) in base 16, for odd n. - Georg Fischer, Mar 23 2019
E.g.f.: (exp(4*x) + 4*exp(-x))/5. - G. C. Greubel, Mar 23 2019

Extensions

Corrected by Franklin T. Adams-Watters, Sep 18 2006

A115341 a(n) = abs(A154879(n+1)).

Original entry on oeis.org

2, 4, 0, 8, 8, 24, 40, 88, 168, 344, 680, 1368, 2728, 5464, 10920, 21848, 43688, 87384, 174760, 349528, 699048, 1398104, 2796200, 5592408, 11184808, 22369624, 44739240, 89478488, 178956968, 357913944, 715827880, 1431655768, 2863311528
Offset: 0

Views

Author

Roger L. Bagula, Mar 06 2006

Keywords

Comments

General form: a(n)=2^n-a(n-1). - Vladimir Joseph Stephan Orlovsky, Dec 11 2008
For n>=1, a(n) is a(n) is the number of generalized compositions of n+3 when there are i^2-2*i-1 different types of i, (i=1,2,...). - Milan Janjic, Sep 24 2010

Crossrefs

Programs

  • Magma
    [2] cat [(2^(n+1)-8*(-1)^n)/3: n in [1..30]]; // G. C. Greubel, Dec 30 2017
  • Mathematica
    g0[n_] = 2 - Sum[(-1)^(i + 1)/Sqrt[2]^(2*i), {i, 0, n}] f[x_] = ZTransform[g0[n], n, x] g[n_] = InverseZTransform[f[1/x], x, n] a0 = Table[Abs[g[n]], {n, 1, 25}]
    k=0;lst={k};Do[k=2^n-k;AppendTo[lst, k], {n, 3, 5!}];lst (* Vladimir Joseph Stephan Orlovsky, Dec 11 2008 *)
    Table[If[n==0, 2, (2^(n+1)-8*(-1)^n)/3], {n,0,30}] (* G. C. Greubel, Dec 30 2017 *)
  • PARI
    for(n=0,30, print1(if(n==0, 2, (2^(n+1)-8*(-1)^n)/3), ", ")) \\ G. C. Greubel, Dec 30 2017
    

Formula

a(n) = (2^(n+1)-8*(-1)^n)/3, n>0.
a(n) = a(n-1) + 2*a(n-2), n>2.
G.f.: 2+4*x*(1-x)/((1+x)*(1-2*x)).

Extensions

Edited by the Associate Editors of the OEIS, Aug 21 2009

A048573 a(n) = a(n-1) + 2*a(n-2), a(0)=2, a(1)=3.

Original entry on oeis.org

2, 3, 7, 13, 27, 53, 107, 213, 427, 853, 1707, 3413, 6827, 13653, 27307, 54613, 109227, 218453, 436907, 873813, 1747627, 3495253, 6990507, 13981013, 27962027, 55924053, 111848107, 223696213, 447392427, 894784853, 1789569707, 3579139413, 7158278827, 14316557653
Offset: 0

Views

Author

Michael Somos, Jun 17 1999

Keywords

Comments

Number of positive integers requiring exactly n signed bits in the modified non-adjacent form representation. - Ralf Stephan, Aug 02 2003
The n-th entry (n>1) of the sequence is equal to the 1,1-entry of the n-th power of the unnormalized 4 X 4 Haar matrix: [1 1 1 0 / 1 1 -1 0 / 1 1 0 1 / 1 1 0 -1]. - Simone Severini, Oct 27 2004
Pisano period lengths: 1, 1, 6, 2, 2, 6, 6, 2, 18, 2, 10, 6, 12, 6, 6, 2, 8, 18, 18, 2, ... - R. J. Mathar, Aug 10 2012
For n >= 1, a(n) is the number of ways to tile a strip of length n+2 with blue squares and blue and red dominos, with the restriction that the first two tiles must be the same color. - Guanji Chen and Greg Dresden, Jul 15 2024

Examples

			G.f. = 2 + 3*x + 7*x^2 + 13*x^3 + 27*x^4 + 53*x^5 + 107*x^6 + 213*x^7 + 427*x^8 + ...
		

Crossrefs

Programs

  • Magma
    [(5*2^n+(-1)^n)/3: n in [0..35]]; // Vincenzo Librandi, Jul 05 2011
    
  • Mathematica
    LinearRecurrence[{1,2},{2,3},40] (* Harvey P. Dale, Dec 11 2017 *)
  • PARI
    {a(n) = if( n<0, 0, (5*2^n + (-1)^n) / 3)};
    
  • PARI
    {a(n) = if (n<0 ,0, if( n<2, n+2, a(n-1) + 2*a(n-2)))};
    
  • Sage
    [(5*2^n+(-1)^n)/3 for n in range(35)] # G. C. Greubel, Apr 10 2019

Formula

G.f.: (2 + x) / (1 - x - 2*x^2).
a(n) = (5*2^n + (-1)^n) / 3.
a(n) = 2^(n+1) - A001045(n).
a(n) = A084170(n)+1 = abs(A083581(n)-3) = A081254(n+1) - A081254(n) = A084214(n+2)/2.
a(n) = 2*A001045(n+1) + A001045(n) (note that 2 is the limit of A001045(n+1)/A001045(n)). - Paul Barry, Sep 14 2009
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=-3, A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)=-charpoly(A,-1). - Milan Janjic, Jan 27 2010
Equivalently, with different offset, a(n) = b(n+1) with b(0)=1 and b(n) = Sum_{i=0..n-1} (-1)^i (1 + (-1)^i b(i)). - Olivier Gérard, Jul 30 2012
a(n) = A000975(n-2)*10 + 5 + 2*(-1)^(n-2), a(0)=2, a(1)=3. - Yuchun Ji, Mar 18 2019
a(n+1) = Sum_{i=0..n} a(i) + 1 + (1-(-1)^n)/2, a(0)=2. - Yuchun Ji, Apr 10 2019
a(n) = 2^n + J(n+1) = J(n+2) + J(n+1) - J(n), where J is A001045. - Yuchun Ji, Apr 10 2019
a(n) = A001045(n+2) + A078008(n) = A062510(n+1) - A078008(n+1) = (A001045(n+2) + A062510(n+1))/2 = A014551(n) + 2*A001045(n). - Paul Curtz, Jul 14 2021
From Thomas Scheuerle, Jul 14 2021: (Start)
a(n) = A083322(n) + A024493(n).
a(n) = A127978(n) - A102713(n).
a(n) = A130755(n) - A166249(n).
a(n) = A007679(n) + A139763(n).
a(n) = A168642(n) XOR A007283(n).
a(n) = A290604(n) + A083944(n). (End)
From Paul Curtz, Jul 21 2021: (Start)
a(n) = 5*A001045(n) - A280560(n+1) = abs(A140360(n+1)) - A280560(n+1).
a(n) = 2^n + A001045(n+1) = A001045(n+3) - A000079(n).
a(n) = A001045(n+4) - A340627(n). (End)
a(n) = A001045(n+5) - A005010(n).
a(n+1) + a(n) = a(n+2) - a(n) = 5*2^n. - Michael Somos, Feb 22 2023
a(n) = A135318(2*n) + A135318(2*n+1) = A112387(2*n) + A112387(2*n+1). - Paul Curtz, Jun 26 2024
E.g.f.: (cosh(x) + 5*cosh(2*x) - sinh(x) + 5*sinh(2*x))/3. - Stefano Spezia, May 18 2025

Extensions

Formula of Milan Janjic moved here from wrong sequence by Paul D. Hanna, May 29 2010

A087787 a(n) = Sum_{k=0..n} (-1)^(n-k)*A000041(k).

Original entry on oeis.org

1, 0, 2, 1, 4, 3, 8, 7, 15, 15, 27, 29, 48, 53, 82, 94, 137, 160, 225, 265, 362, 430, 572, 683, 892, 1066, 1370, 1640, 2078, 2487, 3117, 3725, 4624, 5519, 6791, 8092, 9885, 11752, 14263, 16922, 20416, 24167, 29007, 34254, 40921, 48213, 57345, 67409
Offset: 0

Views

Author

Vladeta Jovovic, Oct 07 2003

Keywords

Comments

Essentially first differences of A024786 (see the formula). Also, a(n) is the number of 2's in the last section of the set of partitions of n+2 (see A135010). - Omar E. Pol, Sep 10 2008
From Gus Wiseman, May 20 2024: (Start)
Also the number of integer partitions of n containing an even number of ones, ranked by A003159. The a(0) = 1 through a(8) = 15 partitions are:
() . (2) (3) (4) (5) (6) (7) (8)
(11) (22) (32) (33) (43) (44)
(211) (311) (42) (52) (53)
(1111) (222) (322) (62)
(411) (511) (332)
(2211) (3211) (422)
(21111) (31111) (611)
(111111) (2222)
(3311)
(4211)
(22211)
(41111)
(221111)
(2111111)
(11111111)
Also the number of integer partitions of n + 1 containing an odd number of ones, ranked by A036554.
(End)

Crossrefs

The unsigned version is A000070, strict A036469.
For powers of 2 instead number of partitions we have A001045.
The strict or odd version is A025147 or A096765.
The ordered version (compositions instead of partitions) is A078008.
For powers of 2 instead of -1 we have A259401, cf. A259400.
A002865 counts partitions with no ones, column k=0 of A116598.
A072233 counts partitions by sum and length.

Programs

  • Mathematica
    Table[Sum[(-1)^(n-k)*PartitionsP[k], {k,0,n}], {n,0,50}] (* Vaclav Kotesovec, Aug 16 2015 *)
    (* more efficient program *) sig = 1; su = 1; Flatten[{1, Table[sig = -sig; su = su + sig*PartitionsP[n]; Abs[su], {n, 1, 50}]}] (* Vaclav Kotesovec, Nov 06 2016 *)
    Table[Length[Select[IntegerPartitions[n], EvenQ[Count[#,1]]&]],{n,0,30}] (* Gus Wiseman, May 20 2024 *)
  • Python
    from sympy import npartitions
    def A087787(n): return sum(-npartitions(k) if n-k&1 else npartitions(k) for k in range(n+1)) # Chai Wah Wu, Oct 25 2023

Formula

G.f.: 1/(1+x)*1/Product_{k>0} (1-x^k).
a(n) = 1/n*Sum_{k=1..n} (sigma(k)+(-1)^k)*a(n-k).
a(n) = A024786(n+2)-A024786(n+1). - Omar E. Pol, Sep 10 2008
a(n) ~ exp(Pi*sqrt(2*n/3)) / (8*sqrt(3)*n) * (1 + (11*Pi/(24*sqrt(6)) - sqrt(3/2)/Pi)/sqrt(n) - (11/16 + (23*Pi^2)/6912)/n). - Vaclav Kotesovec, Nov 05 2016
a(n) = A000041(n) - a(n-1). - Jon Maiga, Aug 29 2019
Alternating partial sums of A000041. - Gus Wiseman, May 20 2024
Previous Showing 31-40 of 133 results. Next