A006053 -id:A006053 - cp's OEIS frontend

A053538 Triangle: a(n,m) = ways to place p balls in n slots with m in the rightmost p slots, 0<=p<=n, 0<=m<=n, summed over p, a(n,m)= Sum_{k=0..n} binomial(k,m)*binomial(n-k,k-m), (see program line).

1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 5, 5, 4, 1, 1, 8, 10, 7, 5, 1, 1, 13, 18, 16, 9, 6, 1, 1, 21, 33, 31, 23, 11, 7, 1, 1, 34, 59, 62, 47, 31, 13, 8, 1, 1, 55, 105, 119, 101, 66, 40, 15, 9, 1, 1, 89, 185, 227, 205, 151, 88, 50, 17, 10, 1, 1, 144, 324, 426, 414, 321, 213, 113, 61, 19, 11, 1, 1

Offset: 0

Wouter Meeussen, May 23 2001

Riordan array (1/(1-x-x^2), x(1-x)/(1-x-x^2)). Row sums are A000079. Diagonal sums are A006053(n+2). - Paul Barry, Nov 01 2006
Subtriangle of the triangle given by (0, 1, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 05 2012
Mirror image of triangle in A208342. - Philippe Deléham, Mar 05 2012
A053538 is jointly generated with A076791 as an array of coefficients of polynomials u(n,x):  initially, u(1,x)=v(1,x)=1, for n>1, u(n,x) = x*u(n-1,x) + v(n-1,x) and v(n,x) = u(n-1,x) + v(n-1,x).  See the Mathematica section at A076791. - Clark Kimberling, Mar 08 2012
The matrix inverse starts
    1;
   -1,   1;
   -1,  -1,   1;
    1,  -2,  -1,  1;
    3,   1,  -3, -1,  1;
    1,   6,   1, -4, -1,  1;
   -7,   4,  10,  1, -5, -1,  1;
  -13, -13,   8, 15,  1, -6, -1,  1;
    3, -31, -23, 13, 21,  1, -7, -1, 1; - R. J. Mathar, Mar 15 2013
Also appears to be the number of subsets of {1..n} containing n with k maximal anti-runs of consecutive elements increasing by more than 1. For example, the subset {1,3,6,7,11,12} has maximal anti-runs ((1,3,6),(7,11),(12)) so is counted under a(12,3). For runs instead of anti-runs we get A202064. - Gus Wiseman, Jun 26 2025

			n=4; Table[binomial[k, j]binomial[n-k, k-j], {k, 0, n}, {j, 0, n}] splits {1, 4, 6, 4, 1} into {{1, 0, 0, 0, 0}, {3, 1, 0, 0, 0}, {1, 4, 1, 0, 0}, {0, 0, 3, 1, 0}, {0, 0, 0, 0, 1}} and this gives summed by columns {5, 5, 4, 1, 1}
Triangle begins :
   1;
   1,  1;
   2,  1,  1;
   3,  3,  1, 1;
   5,  5,  4, 1, 1;
   8, 10,  7, 5, 1, 1;
  13, 18, 16, 9, 6, 1, 1;
...
(0, 1, 1, -1, 0, 0, 0, ...) DELTA (1, 0, -1, 1, 0, 0, 0, ...) begins :
  1;
  0,  1;
  0,  1,  1;
  0,  2,  1,  1;
  0,  3,  3,  1, 1;
  0,  5,  5,  4, 1, 1;
  0,  8, 10,  7, 5, 1, 1;
  0, 13, 18, 16, 9, 6, 1, 1;

Alois P. Heinz, Rows n = 0..140, flattened
R. P. Grimaldi, Extraordinary subsets: a generalization, Fib. Quart., 55 (No. 3, 2017), 114-122. See Table 1.

Column k = 1 is A000045.
Row sums are A000079.
Column k = 2 is A010049.
For runs instead of anti-runs we have A202064.
For integer partitions see A268193, strict A384905, runs A116674.
A034839 counts subsets by number of maximal runs.
A384175 counts subsets with all distinct lengths of maximal runs, complement A384176.
A384877 gives lengths of maximal anti-runs in binary indices, firsts A384878.
A384893 counts subsets by number of maximal anti-runs.
Cf. A000071, A001629, A010027, A208342, A384177, A384879, A384889, A384890.

Flat(List([0..12], n-> List([0..n], k-> Sum([0..n], j->  Binomial(j,k)*Binomial(n-j,j-k)) ))); # G. C. Greubel, May 16 2019

[[(&+[Binomial(j,k)*Binomial(n-j,j-k): j in [0..n]]): k in [0..n]]: n in [0..12]]; // G. C. Greubel, May 16 2019

a:= (n, m)-> add(binomial(k, m)*binomial(n-k, k-m), k=0..n):
seq(seq(a(n,m), m=0..n), n=0..12);  # Alois P. Heinz, Sep 19 2013

Table[Sum[Binomial[k, m]*Binomial[n-k, k-m], {k,0,n}], {n,0,12}, {m,0,n}]

{T(n,k) = sum(j=0,n, binomial(j,k)*binomial(n-j,j-k))}; \\ G. C. Greubel, May 16 2019

[[sum(binomial(j,k)*binomial(n-j,j-k) for j in (0..n)) for k in (0..n)] for n in (0..12)] # G. C. Greubel, May 16 2019

From Philippe Deléham, Mar 05 2012: (Start)
T(n,k) = T(n-1,k) + T(n-1,k-1) + T(n-2,k) - T(n-2,k-1), T(0,0) = T(1,0) = T(1,1) = 1 and T(n,k) = 0 if k<0 or if k>n.
G.f.: 1/(1-(1+y)*x-(1-y)*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = A077957(n), A000045(n+1), A000079(n), A001906(n+1), A007070(n), A116415(n), A084326(n+1), A190974(n+1), A190978(n+1), A190984(n+1), A190990(n+1), A190872(n+1) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 respectively. (End)

A120757 Expansion of x^2(2+x)/(1-3x-4*x^2-x^3).

Original entry on oeis.org

0, 2, 7, 29, 117, 474, 1919, 7770, 31460, 127379, 515747, 2088217, 8455018, 34233669, 138609296, 561217582, 2272323599, 9200450421, 37251863241, 150829715006, 610697048403, 2472661868474, 10011603514040, 40536155064419

Offset: 1

Gary W. Adamson & Roger L. Bagula, Jul 01 2006

The (1,1)-entry of the matrix M^n, where M is the 3 X 3 matrix [0,1,1; 1,1,2; 1,2,2].
a(n)/a(n-1) tends to 4.0489173...an eigenvalue of M and a root to the characteristic polynomial x^3 - 3x^2 - 4x - 1.
C(n):=a(n), with a(0):=1 (hence the o.g.f. for C(n) is (1-3*x-2*x^2)/(1-3*x-4*x^2-x^3)), appears in the following formula for the nonnegative powers of rho*sigma, where rho:=2*cos(Pi/7) and sigma:=sin(3*Pi/7)/sin(Pi/7) = rho^2-1 are the ratios of the smaller and larger diagonal length to the side length in a regular 7-gon (heptagon). See the Steinbach reference where the basis <1,rho,sigma> is used in an extension of the rational field. (rho*sigma)^n = C(n) + B(n)*rho + A(n)*sigma,n>=0, with B(n)= |A122600(n-1)|, B(0)=0, and A(n)= A181879(n). For the nonpositive powers see A085810(n)*(-1)^n, A181880(n-2)*(-1)^n and A116423(n+1)*(-1)^(n+1), respectively. See also a comment under A052547.
We have a(n)=cs(3n+1), where the sequence cs(n) and its two conjugate sequences as(n) and bs(n) are defined in the comments to the sequence A214683 (see also A215076, A215100, A006053). We call the sequence a(n) the Ramanujan-type sequence number 5 for the argument 2Pi/7. Since as(3n+1)=bs(3n+1)=0, we obtain the following relation: 49^(1/3)*a(n) = (c(1)/c(4))^(n + 1/3) + (c(4)/c(2))^(n + 1/3) + (c(2)/c(1))^(n + 1/3), where c(j) := Cos(2Pi/7) (for more details and proofs see Witula et al.'s papers). - Roman Witula, Aug 02 2012

			a(7)=1919 because M^7= [1919,3458,4312;3458,6231,7770;4312,7770,9689].

R. Witula, E. Hetmaniok and D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

G. C. Greubel, Table of n, a(n) for n = 1..1000
P. Steinbach, Golden fields: a case for the heptagon, Math. Mag. 70 (1997), no. 1, 22-31, MR 1439165
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,4,1).

Cf. A006053, A214683, A215076, A215100.

a:=[0,2,7]; [ n le 3 select a[n] else 3*Self(n-1) + 4*Self(n-2) + Self(n-3): n in [1..25]]; // Marius A. Burtea, Oct 03 2019

with(linalg): M[1]:=matrix(3,3,[0,1,1,1,1,2,1,2,2]): for n from 2 to 25 do M[n]:=multiply(M[1],M[n-1]) od: seq(M[n][1,1],n=1..25);

LinearRecurrence[{3,4,1},{0,2,7},40] (* Roman Witula, Aug 02 2012 *)

a(n)=([0,1,0; 0,0,1; 1,4,3]^(n-1)*[0;2;7])[1,1] \\ Charles R Greathouse IV, Jun 22 2016

@CachedFunction
def a(n): # a = A120757
    if (n<3): return (0,2,7)[n]
    else: return 3*a(n-1) + 4*a(n-2) + a(n-3)
[a(n) for n in range(40)] # G. C. Greubel, Nov 25 2022

a(n) = 3*a(n-1) + 4*a(n-2) + a(n-3) (follows from the minimal polynomial of the matrix M). See also the o.g.f. given in the name.

Edited by N. J. A. Sloane, Dec 03 2006

New name, old name as comment; o.g.f.; reference.

A214779 a(n) = 3*a(n-2) - a(n-3) with a(0)=-1, a(1)=1, a(2)=-4.

Original entry on oeis.org

-1, 1, -4, 4, -13, 16, -43, 61, -145, 226, -496, 823, -1714, 2965, -5965, 10609, -20860, 37792, -73189, 134236, -257359, 475897, -906313, 1685050, -3194836, 5961463, -11269558, 21079225, -39770137, 74507233, -140389636, 263291836, -495676141, 930265144

Offset: 0

Roman Witula, Jul 28 2012

Ramanujan-type sequence number 2 for argument 2Pi/9 is connected with the sequence A214699 (see also sequences A006053, A214683) - all have "similar" trigonometric description, for example in the case of a(n) the following formula hold true: 9^(1/3)*a(n) = (c(1)/c(4))^(1/3)*c(1)^n + (c(2)/c(1))^(1/3)*c(2)^n + (c(4)/c(2))^(1/3)*c(4)^n = -( (c(2)/c(1))^(1/3)*c(1)^(n+1) + (c(4)/c(2))^(1/3)*c(2)^(n+1) + (c(1)/c(4))^(1/3)*c(4)^(n+1) ), where c(j) := 2*Cos(2Pi*j/9) - for the proof see Witula et al.'s papers.
From a(0),A214699(0),a(2) and c(1)+c(2)+c(4)=0 we deduce
x^3 - 9^(1/3)*x - 1 = (x - (c(1)/c(2))^(1/3))*(x - (c(2)/c(4))^(1/3))*(x - (c(4)/c(1))^(1/3)), and
x^3 - 7*9^(1/3)*x - 1 = (x - (c(1)/c(2))^(1/3)*c(1)^2)*(x - (c(2)/c(4))^(1/3)*c(2)^2)*(x - (c(4)/c(1))^(1/3)*c(4)^2). We note that applying the Newton-Girard formulas to these polynomials two new sequences of real numbers can be discussed: X(n) := (c(1)/c(2))^(n/3) + (c(2)/c(4))^(n/3) + (c(4)/c(1))^(n/3), and Y(n) := ((c(1)/c(2))^(1/3)*c(1)^2)^n + ((c(2)/c(4))^(1/3)*c(2)^2)^n + ((c(4)/c(1))^(1/3)*c(4)^2)^n, where X(n)=9^(1/3)*X(n-2)+X(n-3), X(0)=3, X(1)=0, X(2)=2*9^(1/3), Y(n)=7*9^(1/3)Y(n-2)+Y(n-3), Y(0)=3, Y(1)=0, Y(2)=14*9^(1/3). It could be obtained the following decompositions: X(n) = ax(n) + 9^(1/3)*bx(n) + 81^(1/3)*cx(n), ax(0)=3, bx(0)=cx(0)=ax(1)=bx(1)=cx(1)=ax(2)=bx(2)=0, cx(2)=2, ax(n)=ax(n-3)+9*cx(n-2), bx(n)=bx(n-3)+ax(n-2), cx(n)=cx(n-3)+bx(n-2), and  Y(n) = ay(n) + 9^(1/3)*by(n) + 81^(1/3)*cy(n), ay(0)=3, by(0)=cy(0)=ay(1)=by(1)=cy(1)=ay(2)=cy(2)=0, by(2)=14, ay(n)=ay(n-3)+63*cy(n-2), by(n)=by(n-3)+7*ay(n-2), cy(n)=cy(n-3)+7*by(n-2). All these new sequence of positive integers ax(n),bx(n),...,cy(n) will be presented separately as A214778, A214951, A214954. - Roman Witula, Sep 27 2012
We note that all sums a(n+1) + a(n) are divisible by 3, which easily from recurrence formula for a(n) follows. Then it can be deduced the formula a(n+1) + a(n) = -A214699(n). - Roman Witula, Oct 06 2012

			From a(0)=-1 and A214699(0)=0 we obtain (c(1)/c(4))^(2/3) + (c(2)/c(1))^(2/3) + (c(4)/c(2))^(2/3) = 3*3^(1/3), whereas from a(1)=-1 and A214699(1)=3*3^(1/3) we get (c(1)/c(4))^(2/3)*2c(2) + (c(2)/c(1))^(2/3)*2c(4) + (c(4)/c(2))^(2/3)*2c(1) = 3*3^(1/3).

R. Witula, E. Hetmaniok, D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (0,3,-1).

Cf. A214699, A214778, A214951, A214954, A217053, A217052, A217069.

LinearRecurrence[{0, 3, -1}, {-1, 1, -4}, 40] (* T. D. Noe, Jul 30 2012 *)

G.f.: -(1-x+x^2)/(1-3*x^2+x^3).

A078038 Expansion of (1-x)/(1+x-2*x^2-x^3).

Original entry on oeis.org

1, -2, 4, -7, 13, -23, 42, -75, 136, -244, 441, -793, 1431, -2576, 4645, -8366, 15080, -27167, 48961, -88215, 158970, -286439, 516164, -930072, 1675961, -3019941, 5441791, -9805712, 17669353, -31838986, 57371980, -103380599, 186285573, -335674791, 604865338, -1089929347, 1963985232

Offset: 0

N. J. A. Sloane, Nov 17 2002

From Johannes W. Meijer, May 29 2010: (Start)
The absolute values of the a(n) represent the number of ways White can force checkmate in exactly (n+1) moves, n>=0, ignoring the fifty-move and the triple repetition rules, in the following chess position: White Ka1, Ra8, Bc1, Nb8, pawns a6, a7, b2, c6, d2, f6 and h6; Black Ke8, pawns b3, c7, d3, f7 and h7. (After Noam D. Elkies, see link; diagram 5).
The absolute values of the a(n) represent all paths of length n starting at the third (or fourth) node on the path graph P_6, see the Maple program.
(End)
For n>=1, abs(a(n-1)) is the number of compositions where there is no rise between every second pair of parts, starting with the second and third part; see example. Also, abs(a(n-1)) is the number of compositions of n where there is no fall between every second pair of parts, starting with the second and third part; see example. [Joerg Arndt, May 21 2013]

			From _Joerg Arndt_, May 21 2013: (Start)
There are abs(a(6-1))=23 compositions of 6 where there is no rise between every second pair of parts:
          v   v    <--= no rise over these positions
01:  [ 1 1 1 1 1 1 ]
02:  [ 1 1 1 2 1 ]
03:  [ 1 1 1 3 ]
04:  [ 1 2 1 1 1 ]
05:  [ 1 2 1 2 ]
06:  [ 1 2 2 1 ]
07:  [ 1 3 1 1 ]
08:  [ 1 3 2 ]
09:  [ 1 4 1 ]
10:  [ 1 5 ]
11:  [ 2 1 1 1 1 ]
12:  [ 2 1 1 2 ]
13:  [ 2 2 1 1 ]
14:  [ 2 2 2 ]
15:  [ 2 3 1 ]
16:  [ 2 4 ]
17:  [ 3 1 1 1 ]
18:  [ 3 2 1 ]
19:  [ 3 3 ]
20:  [ 4 1 1 ]
21:  [ 4 2 ]
22:  [ 5 1 ]
23:  [ 6 ]
There are abs(a(6-1))=23 compositions of 6 where there is no fall between every second pair of parts, starting with the second and third part:
          v   v    <--= no fall over these positions
01:  [ 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 2 ]
03:  [ 1 1 1 3 ]
04:  [ 1 1 2 1 1 ]
05:  [ 1 1 2 2 ]
06:  [ 1 1 3 1 ]
07:  [ 1 1 4 ]
08:  [ 1 2 2 1 ]
09:  [ 1 2 3 ]
10:  [ 1 5 ]
11:  [ 2 1 1 1 1 ]
12:  [ 2 1 1 2 ]
13:  [ 2 1 2 1 ]
14:  [ 2 1 3 ]
15:  [ 2 2 2 ]
16:  [ 2 4 ]
17:  [ 3 1 1 1 ]
18:  [ 3 1 2 ]
19:  [ 3 3 ]
20:  [ 4 1 1 ]
21:  [ 4 2 ]
22:  [ 5 1 ]
23:  [ 6 ]
(End)

Noam D. Elkies, New Directions in Enumerative Chess Problems., arXiv:math/0508645 [math.CO]; 2005; The Electronic Journal of Combinatorics, 11 (2), 2004-2005. [From _Johannes W. Meijer_, May 29 2010]
Index entries for linear recurrences with constant coefficients, signature (-1,2,1).

Cf. A028495, A068911, A094790, A287381 (absolute values).

with(GraphTheory): G:= PathGraph(6): A:=AdjacencyMatrix(G): nmax:=36; for n from 0 to nmax do B(n):=A^n; a(n):=add(B(n)[3,k], k=1..6) od: seq(a(n), n=0..nmax); # Johannes W. Meijer, May 29 2010

LinearRecurrence[{-1, 2, 1}, {1, -2, 4}, 40] (* Jean-François Alcover, Jan 08 2019 *)
a[n_]:=Sum[(-(-2)^(n+1)Cos[(Pi r)/7]^n Cot[(Pi r)/14]Sin[(3Pi r)/7])/7,{r,1,5,2}]
Table[a[n],{n,0,40}]//Round (* Herbert Kociemba, Sep 17 2020 *)

Vec((1-x)/(1+x-2*x^2-x^3)+O(x^99)) \\ Charles R Greathouse IV, Sep 25 2012

a(n+3) = -a(n+2)+2*a(n+1)+a(n), a(0)=1, a(1)=-2, a(2)=4. - Wouter Meeussen, Jan 02 2005
a(n) = (-1)^n * (A006053(n+1) + A006053(n+2)). G.f. of |a(n)|: (1+x)/(x^3 - 2*x^2 - x + 1). - Ralf Stephan, Aug 19 2013
a(n) = Sum_{r=1..6} ((-2)^n*(1-(-1)^r)*cos(Pi*r/7)^n*cot(Pi*r/14)*sin(3*Pi*r/7))/7. - Herbert Kociemba, Sep 17 2020

A096976 Number of walks of length n on P_3 plus a loop at the end.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 5, 5, 14, 19, 42, 66, 131, 221, 417, 728, 1341, 2380, 4334, 7753, 14041, 25213, 45542, 81927, 147798, 266110, 479779, 864201, 1557649, 2806272, 5057369, 9112264, 16420730, 29587889, 53317085, 96072133, 173118414, 311945595, 562110290

Offset: 0

Paul Barry, Jul 16 2004

Counts closed walks of length n at the start of P_3 to which a loop has been added at the other extremity. a(n+1) counts walks between the first node and the last. Let A be the adjacency matrix of the graph P_3 with a loop added at the end. A is a 'reverse Jordan matrix' [0,0,1;0,1,1;1,1,0]. a(n) is obtained by taking the (1,1) element of A^n.
Sequence is also related to matrices associated with rhombus substitution tilings showing 7-fold rotational symmetry. Let A_{7,1} be the 3 X 3 unit-primitive matrix (see [Jeffery]) A_{7,1}=[0,1,0; 1,0,1; 0,1,1]; then a(n)=[A_{7,1}^n](1,1). - _L. Edson Jeffery, Jan 05 2012
a(n+2) is the (1,1) element of the n-th power of each of the two 3 X 3 matrices: [0,1,1; 1,0,0; 1,0,1], [0,1,1; 1,1,0; 1,0,0]. - Christopher Hunt Gribble, Apr 03 2014

			G.f. = 1 + x^2 + 2*x^4 + x^5 + 5*x^6 + 5*x^7 + 14*x^8 + 19*x^9 + ... - _Michael Somos_, Dec 12 2023

Michael De Vlieger, Table of n, a(n) for n = 0..3914
Robin Chapman and Nicholas C. Singer, Eigenvalues of a bidiagonal matrix, Amer. Math. Monthly, 111 (2004) 441.
Tomislav Došlić, Mate Puljiz, Stjepan Šebek, and Josip Žubrinić, On a variant of Flory model, arXiv:2210.12411 [math.CO], 2022.
L. E. Jeffery, Unit-primitive matrices
Kai Wang, Fibonacci Numbers And Trigonometric Functions Outline, (2019).
Index entries for linear recurrences with constant coefficients, signature (1,2,-1).

Cf. A006053, A052547, A077998.

LinearRecurrence[{1, 2, -1}, {1, 0, 1}, 60] (* Vladimir Joseph Stephan Orlovsky, Feb 13 2012 *)
a[ n_] := {1, 0, 0} . MatrixPower[{{1, 2, -1}, {1, 0, 0}, {0, 1, 0}}, n] . {1, 1, 3}; (* Michael Somos, Dec 12 2023 *)

{a(n) = [1, 0, 0] * [1, 2, -1; 1, 0, 0; 0, 1, 0]^n * [1, 1, 3]~}; /* Michael Somos, Dec 12 2023 */

G.f. : (1-x-x^2)/(1-x-2x^2+x^3); a(n)=a(n-1)+2a(n-2)-a(n-3).
a(n) = 5a(n-2)-6a(n-4)+a(n-6). - Floor van Lamoen, Nov 02 2005
a(n) = A077998(-n) for all n in Z. - Michael Somos, Dec 12 2023

A215100 a(n) = 3a(n-1) + 4a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=22.

Original entry on oeis.org

2, 5, 22, 88, 357, 1445, 5851, 23690, 95919, 388368, 1572470, 6366801, 25778651, 104375627, 422608286, 1711106017, 6928126822, 28051412820, 113577851765, 459867333397, 1861964820071, 7538941645566, 30524551550379, 123591386053472, 500411306007498, 2026124013786761

Offset: 0

Roman Witula, Aug 03 2012

Ramanujan-type sequence number 4 for the argument 2*Pi/7. We have a(n)=bs(3n+2), where the sequence bs(n) and its two conjugate sequences as(n) and cs(n) are defined in the comments to A214683 (see also A215076, A120757, A006053). Since we also have as(3n+2)=cs(3n+2)=0 from the formula for S(n) (see Comments at A214683) we obtain the relation 7^(1/3)*a(n)= (c(1)/c(4))^(n + 2/3) + (c(4)/c(2))^(n + 2/3) + (c(2)/c(1))^(n + 2/3).

			From 4*a(2) = a(3) = 88 we get 88*7^(1/3) = 4*((c(1)/c(4))^(8/3) + (c(4)/c(2))^(8/3) + (c(2)/c(1))^(8/3))=(c(1)/c(4))^(11/3) + (c(4)/c(2))^(11/3) + (c(2)/c(1))^(11/3).

R. Witula, E. Hetmaniok and D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,4,1).

Cf. A006053, A214683, A215076, A120757.

I:=[2,5,22]; [n le 3 select I[n] else 3*Self(n-1) +4*Self(n-2) +Self(n-3): n in [1..41]]; // G. C. Greubel, Nov 25 2022

LinearRecurrence[{3, 4, 1}, {2, 5, 22}, 40]

Vec((2-x-x^2)/(1-3*x-4*x^2-x^3) + O(x^40)) \\ Michel Marcus, Apr 20 2016

@CachedFunction
def a(n): # a = A215100
    if (n<3): return (2,5,22)[n]
    else: return 3*a(n-1) + 4*a(n-2) + a(n-3)
[a(n) for n in range(41)] # G. C. Greubel, Nov 25 2022

G.f.: (2 - x - x^2)/(1 - 3*x - 4*x^2 - x^3).

More terms from Michel Marcus, Apr 20 2016

A187066 Let i be in {1,2,3} and let r >= 0 be an integer. Let p = {p_1, p_2, p_3} = {-2,0,1}, n=2r+p_i, and define a(-2)=0. Then, a(n)=a(2r+p_i) gives the quantity of H_(7,2,0) tiles in a subdivided H_(7,i,r) tile after linear scaling by the factor x^r, where x=sqrt(2*cos(Pi/7)).

Original entry on oeis.org

1, 0, 0, 1, 2, 1, 1, 3, 5, 4, 5, 9, 14, 14, 19, 28, 42, 47, 66, 89, 131, 155, 221, 286, 417, 507, 728, 924, 1341, 1652, 2380, 2993, 4334, 5373, 7753, 9707, 14041, 17460, 25213, 31501, 45542, 56714, 81927, 102256, 147798, 184183

Offset: 0

L. Edson Jeffery, Mar 09 2011

(Start) See A187067 for supporting theory. Define the matrix
  U_1= (0 1 0)
       (1 0 1)
       (0 1 1).
Let r>=0 and M=(m_(i,j))=(U_1)^r, i,j=1,2,3. Let B_r be the r-th "block" defined by B_r={a(2*r-2),a(2*r),a(2*r+1)} with a(-2)=0. Note that B_r-B_(r-1)-2*B_(r-2)+B_(r-3)={0,0,0}, with B_0={a(-2),a(0),a(1)}={0,1,0}. Let p={p_1,p_2,p_3}=(-2,0,1) and n=2*r+p_i. Then a(n)=a(2*r+p_i)=m_(i,2), where M=(m_(i,j))=(U_1)^r was defined above. Hence the block B_r corresponds component-wise to the second column of M, and a(n)=m_(i,2) gives the quantity of H_(7,2,0) tiles that should appear in a subdivided H_(7,i,r) tile. (End)
Combining blocks A_r, B_r and C_r, from A187065, this sequence and A187067, respectively, as matrix columns [A_r,B_r,C_r] generates the matrix (U_1)^r, and a negative index (-1)*r yields the corresponding inverse [A_(-r),B_(-r),C_(-r)]=(U_1)^(-r) of (U_1)^r. Therefore, the three sequences need not be causal.
Since a(2*r-2)=a(2*(r-1)) for all r, this sequence arises by concatenation of second-column entries m_(2,2) and m_(3,2) from successive matrices M=(U_1)^r.

			Suppose r=3. Then
B_r = B_3 = {a(2*r-2),a(2*r),a(2*r+1)}={a(4),a(6),a(7)} = {2,1,3},
corresponding to the entries in the third column of
M = (U_2)^3 = (0 2 1)
              (2 1 3)
              (1 3 3).
Choose i=2 and set n=2*r+p_i. Then a(n) = a(2*r+p_i) = a(6+0) = a(6) = 1, which equals the entry in row 2 and column 2 of M. Hence a subdivided H_(7,2,3) tile should contain a(6) = m_(2,2) = 1 H_(7,2,0) tiles.

G. C. Greubel, Table of n, a(n) for n = 0..1000
L. Edson Jeffery, Unit-primitive matrices
Roman Witula, Damian Slota and Adam Warzynski, Quasi-Fibonacci Numbers of the Seventh Order, J. Integer Seq., 9 (2006), Article 06.4.3.
Index entries for linear recurrences with constant coefficients, signature (0,1,0,2,0,-1).

Cf. A187065, A187067, A187068, A187069, A187070.

LinearRecurrence[{0,1,0,2,0,-1},{1,0,0,1,2,1},50] (* Harvey P. Dale, Aug 16 2012 *)
CoefficientList[Series[(1 - x^2 + x^3)/(1 - x^2 - 2*x^4 + x^6), {x, 0, 50}], x] (* G. C. Greubel, Oct 20 2017 *)

my(x='x+O('x^50)); Vec((1-x^2+x^3)/(1-x^2-2*x^4+x^6)) \\ G. C. Greubel, Oct 20 2017

Recurrence: a(n) = a(n-2)+2*a(n-4)-a(n-6).
a(2*n) = A052547(n), a(2n+1) = A006053(n+1).
G.f.: (1-x^2+x^3)/(1-x^2-2*x^4+x^6).
Closed-form: a(n) = -(1/14)*[[X_1+Y_1*(-1)^(n-1)]*[(w_2)^2-(w_3)^2]*(w_1)^(n-1)+[X_2+Y_2*(-1)^(n-1)]*[(w_3)^2-(w_1)^2]*(w_2)^(n-1)+[X_3+Y_3*(-1)^(n-1)]*[(w_1)^2-(w_2)^2]*(w_3)^(n-1)], where w_k = sqrt[2*(-1)^(k-1)*cos(k*Pi/7)], X_k = (w_k)^5-(w_k)^3+(w_k)^2 and Y_k = -(w_k)^5+(w_k)^3+(w_k)^2, k=1,2,3.

A215560 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=a(1)=3, a(2)=101.

Original entry on oeis.org

3, 3, 101, 444, 5981, 38468, 390974, 2948431, 26868565, 216624495, 1888775906, 15657923053, 134074085330, 1124375492334, 9556192325235, 80523923708399, 682280993578341, 5760499663646612, 48746948619251921, 411906111379078256, 3483838470286469746, 29447943482916260935

Offset: 0

Roman Witula, Aug 16 2012

The Ramanujan-type sequence number 6 for the argument 2Pi/7 (see also A214683, A215112, A006053, A006054, A215076, A215100, A120757 for the numbers: 1, 1a, 2, 2a, 3, 4 and 5 respectively).
The sequence a(n) is one of the three special sequences (the remaining two are A215569 and A215572) connected with the following recurrence relation: T(n):=49^(1/3)*T(n-2)+T(n-3), with T(0)=3, T(1)=0, and T(2)=2*49^(1/3) - see the comments to A214683.
It can be proved that
T(n) = (c(1)^4/c(2))^(n/3) + (c(2)^4/c(4))^(n/3) + (c(4)^4/c(1))^(n/3), where c(j):=2*cos(2*Pi*j/7), and the following decomposition hold true:
T(n) = at(n) + bt(n)*7^(1/3) + ct(n)*49^(1/3), where sequences at(n), bt(n), and ct(n) satisfy the following system of recurrence equations: at(n)=7*bt(n-2)+at(n-3),
bt(n)=7*ct(n-2)+bt(n-3), ct(n)=at(n-2)+ct(n-3), with at(0)=3, at(1)=at(2)=bt(0)=bt(1)=bt(2)=ct(0)=ct(1)=0, ct(2)=2 - for details see the first Witula reference.
It follows that a(n)=at(3*n), bt(3*n)=ct(3*n)=0.
Every difference of the form a(n)-a(n-2)-a(n-3) is divisible by 3. Because the difference a(n+1)-a(n) is congruent to the difference a(n-4)-a(n-2) modulo 3 we easily deduce that a(6)-a(5) and a(7)-a(6)-2 are both divisible by 3.

R. Witula, E. Hetmaniok, D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012

Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,46,1).

Cf. A214683, A215112, A006053, A006054, A215076, A215100, A120757, A214699.

LinearRecurrence[{3, 46, 1}, {3, 3, 101}, 50]

Vec((3-6*x-46*x^2)/(1-3*x-46*x^2-x^3) + O(x^40)) \\ Michel Marcus, Apr 20 2016

a(n) = (c(1)^4/c(2))^n + (c(2)^4/c(4))^n + (c(4)^4/c(1))^n, where c(j) = 2*cos(2*Pi*j/7).

G.f.: (3-6*x-46*x^2)/(1-3*x-46*x^2-x^3).

More terms from Michel Marcus, Apr 20 2016

A215572 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=106.

Original entry on oeis.org

2, 5, 106, 550, 6531, 44999, 435973, 3384404, 30252969, 246877464, 2135653370, 17793576423, 151867661753, 1276243154087, 10832435479322, 91356359187721, 773637352766062, 6534137016412674, 55281085635664595, 467187197014742851, 3951025667301212597, 33398969150217473532

Offset: 0

Roman Witula, Aug 16 2012

The Ramanujan-type sequence number 8 for the argument 2Pi/7 (see also A214683, A215112, A006053, A006054, A215076, A215100, A120757, A215560, A215569 for the numbers: 1, 1a, 2, 2a, 3-7 respectively). The sequence a(n) is one of the three special sequences (the remaining two are A215560 and A215569) connected with the following recurrence relation:
(c(1)^4/c(2))^(n/3) + (c(2)^4/c(4))^(n/3) + (c(4)^4/c(1))^(n/3) = at(n) + bt(n)*7^(1/3) + ct(n)*49^(1/3), where c(j):=2*cos(2*Pi*j/7), and the sequences at(n), bt(n), and ct(n) are defined in comments to A215560 (see also A215569). It follows that a(n)=ct(3*n+2), at(3*n+2)=bt(3*n+2)=0, which implies the first formula below.
We note that if a(n), a(n+1) and a(n+2) are all odd for some n in N then a(n+3) is even, a(n+4) is odd, a(n+5) and a(n+6) are both even, and the numbers a(n+7), a(n+8), a(n+9) are all odd again. In consequence, this situation hold for every n of the form 7*k+4, k=0,1,..., in the other words cyclical through all sequence a(n), n=4,5,... (from n=1 whenever we start from odd-even-even sequence).

			From 4*a(1)+5*a(2)=a(3) we obtain 4*((c(1)^4/c(2))^(5/3) + (c(2)^4/c(4))^(5/3) + (c(4)^4/c(1))^(5/3)) + 5*((c(1)^4/c(2))^(8/3) + (c(2)^4/c(4))^(8/3) + (c(4)^4/c(1))^(8/3)) = (4 + 5*c(1)^4/c(2))*((c(1)^4/c(2))^(5/3) + (4 + 5*c(2)^4/c(4))*((c(2)^4/c(4))^(5/3) + (4 + 5*c(4)^4/c(1))*((c(4)^4/c(1))^(5/3) = (c(1)^4/c(2))^(11/3) + (c(2)^4/c(4))^(11/3) + (c(4)^4/c(1))^(11/3) = 550*49^(1/3).

R. Witula, E. Hetmaniok, D. Slota, Sums of the powers of any order roots taken from the roots of a given polynomial, Proceedings of the Fifteenth International Conference on Fibonacci Numbers and Their Applications, Eger, Hungary, 2012.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Roman Witula, Ramanujan Type Trigonometric Formulas: The General Form for the Argument 2*Pi/7, Journal of Integer Sequences, Vol. 12 (2009), Article 09.8.5.
Roman Witula, Full Description of Ramanujan Cubic Polynomials, Journal of Integer Sequences, Vol. 13 (2010), Article 10.5.7.
Roman Witula, Ramanujan Cubic Polynomials of the Second Kind, Journal of Integer Sequences, Vol. 13 (2010), Article 10.7.5.
Roman Witula, Ramanujan Type Trigonometric Formulae, Demonstratio Math. 45 (2012) 779-796.
Index entries for linear recurrences with constant coefficients, signature (3,46,1).

Cf. A214683, A215112, A006053, A006054, A215076, A215100, A120757, A215560, A215569.

LinearRecurrence[{3,46,1}, {2,5,106}, 50]
CoefficientList[Series[(2 - x - x^2)/(1 - 3*x - 46*x^2 - x^3), {x,0,50}], x] (* G. C. Greubel, Apr 16 2017 *)

Vec((2-x-x^2)/(1-3*x-46*x^2-x^3) + O(x^40)) \\ Michel Marcus, Apr 20 2016

49^(1/3)*a(n) = (c(1)^4/c(2))^(n+2/3) + (c(2)^4/c(4))^(n+2/3) + (c(4)^4/c(1))^(n+2/3) = (c(1)*(c(1)/c(2))^(1/3))^(3*n+2) + (c(2)*(c(2)/c(4))^(1/3))^(3*n+2) + (c(4)*(c(4)/c(1))^(1/3))^(3*n+2).

G.f.: (2-x-x^2)/(1-3*x-46*x^2-x^3).

More terms from Michel Marcus, Apr 20 2016

A096975 Trace sequence of a path graph plus loop.

Original entry on oeis.org

3, 1, 5, 4, 13, 16, 38, 57, 117, 193, 370, 639, 1186, 2094, 3827, 6829, 12389, 22220, 40169, 72220, 130338, 234609, 423065, 761945, 1373466, 2474291, 4459278, 8034394, 14478659, 26088169, 47011093, 84708772, 152642789, 275049240

Offset: 0

Paul Barry, Jul 16 2004

Let A be the adjacency matrix of the graph P_3 with a loop added at the end. Then a(n) = trace(A^n). A is a 'reverse Jordan matrix' [0,0,1;0,1,1;1,1,0]. a(n) = abs(A094648(n)).
From L. Edson Jeffery, Mar 22 2011: (Start)
Let A be the unit-primitive matrix (see [Jeffery])
A = A_(7,1) =
(0 1 0)
(1 0 1)
(0 1 1).
Then a(n) = Trace(A^n). (End)

Michael De Vlieger, Table of n, a(n) for n = 0..3910
A. Akbary, Q. Wang, A generalized Lucas sequence and permutations binomials, Proc. Am. Math. Soc. 134 (2006) 15-22, sequence a(n) with l=7.
Robin Chapman and Nicholas C. Singer, Eigenvalues of a bidiagonal matrix, Amer. Math. Monthly, 111 (2004), p. 441.
Tomislav Došlić, Mate Puljiz, Stjepan Šebek, and Josip Žubrinić, On a variant of Flory model, arXiv:2210.12411 [math.CO], 2022.
L. E. Jeffery, Unit-primitive matrix
Genki Shibukawa, New identities for some symmetric polynomials and their applications, arXiv:1907.00334 [math.CA], 2019.
Q. Wang, On generalized Lucas sequences, Contemp. Math. 531 (2010) 127-141, Table 1 (k=3).
Index entries for linear recurrences with constant coefficients, signature (1,2,-1).

Cf. A006053, A052547, A096976.

A033304(n) = a(-1-n). - Michael Somos, Aug 03 2006

CoefficientList[Series[(3 - 2 x - 2 x^2)/(1 - x - 2 x^2 + x^3), {x, 0, 33}], x] (* Michael De Vlieger, Aug 21 2019 *)

{a(n)=if(n>=0, n+=1; polsym(x^3-x^2-2*x+1,n-1)[n], n=1-n; polsym(1-x-2*x^2+x^3,n-1)[n])} /* Michael Somos, Aug 03 2006 */

a(n)=trace([0,1,0;1,0,1;0,1,1]^n); /* Joerg Arndt, Apr 30 2011 */

G.f.: (3-2*x-2*x^2)/(1-x-2*x^2+x^3);
a(n) = a(n-1) + 2*a(n-2) - a(n-3);
a(n) = (2*sqrt(7)*sin(atan(sqrt(3)/9)/3)/3+1/3)^n + (1/3-2*sqrt(7)*sin(atan(sqrt(3)/9)/3+Pi/3)/3)^n + (2*sqrt(7)*cos(acot(-sqrt(3)/9)/3)/3+1/3)^n.
a(n) = 2^n*((cos(Pi/7))^n+(cos(3*Pi/7))^n+(cos(5*Pi/7))^n). - Vladimir Shevelev, Aug 25 2010
a(n) = (-1)^n*A094648(n). - R. J. Mathar, Nov 05 2024

cp's OEIS Frontend

A053538 Triangle: a(n,m) = ways to place p balls in n slots with m in the rightmost p slots, 0<=p<=n, 0<=m<=n, summed over p, a(n,m)= Sum_{k=0..n} binomial(k,m)*binomial(n-k,k-m), (see program line).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A120757 Expansion of x^2*(2+x)/(1-3*x-4*x^2-x^3).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

A214779 a(n) = 3*a(n-2) - a(n-3) with a(0)=-1, a(1)=1, a(2)=-4.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

A078038 Expansion of (1-x)/(1+x-2*x^2-x^3).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A096976 Number of walks of length n on P_3 plus a loop at the end.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A215100 a(n) = 3*a(n-1) + 4*a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=22.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

A120757 Expansion of x^2(2+x)/(1-3x-4*x^2-x^3).

A215100 a(n) = 3a(n-1) + 4a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=22.

A187066 Let i be in {1,2,3} and let r >= 0 be an integer. Let p = {p_1, p_2, p_3} = {-2,0,1}, n=2r+p_i, and define a(-2)=0. Then, a(n)=a(2r+p_i) gives the quantity of H_(7,2,0) tiles in a subdivided H_(7,i,r) tile after linear scaling by the factor x^r, where x=sqrt(2*cos(Pi/7)).

A215560 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=a(1)=3, a(2)=101.

A215572 a(n) = 3a(n-1) + 46a(n-2) + a(n-3) with a(0)=2, a(1)=5, a(2)=106.