A002605 -id:A002605 - cp's OEIS frontend

A125145 a(n) = 3a(n-1) + 3a(n-2). a(0) = 1, a(1) = 4.

1, 4, 15, 57, 216, 819, 3105, 11772, 44631, 169209, 641520, 2432187, 9221121, 34959924, 132543135, 502509177, 1905156936, 7222998339, 27384465825, 103822392492, 393620574951, 1492328902329, 5657848431840, 21450532002507

Offset: 0

Tanya Khovanova, Jan 11 2007

Number of aa-avoiding words of length n on the alphabet {a,b,c,d}.
Equals row 3 of the array shown in A180165, the INVERT transform of A028859 and the INVERTi transform of A086347. - Gary W. Adamson, Aug 14 2010
From Tom Copeland, Nov 08 2014: (Start)
This array is one of a family related by compositions of C(x)= [1-sqrt(1-4x)]/2, an o.g.f. for A000108; its inverse Cinv(x) = x(1-x); and the special Mobius transformation P(x,t) = x / (1+t*x) with inverse P(x,-t) in x. Cf. A091867.
O.g.f.: G(x) = P[P[P[-Cinv(-x),-1],-1],-1] = P[-Cinv(-x),-3] = x*(1+x)/[1-3x(1-x)]= x*A125145(x).
Ginv(x) = -C[-P(x,3)] = [-1 + sqrt(1+4x/(1+3x))]/2 = x*A104455(-x).
G(-x) = -x(1-x) * [ 1 - 3*[x*(1+x)] + 3^2*[x*(1+x)]^2 - ...] , and so this array is related to finite differences in the row sums of A030528 * Diag((-3)^1,3^2,(-3)^3,..). (Cf. A146559.)
The inverse of -G(-x) is C[-P(-x,3)]= [1 - sqrt(1-4x/(1-3x))]/2 = x*A104455(x). (End)
Number of 3-compositions of n+1 restricted to parts 1 and 2 (and allowed zeros); see Hopkins & Ouvry reference. - Brian Hopkins, Aug 16 2020

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Jean-Paul Allouche, Jeffrey Shallit, and Manon Stipulanti, Combinatorics on words and generating Dirichlet series of automatic sequences, arXiv:2401.13524 [math.CO], 2025. See p. 14.
Joerg Arndt, Matters Computational (The Fxtbook)
D. Birmajer, J. B. Gil, and M. D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, Example 7.
Martin Burtscher, Igor Szczyrba, and Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
Brian Hopkins and Stéphane Ouvry, Combinatorics of Multicompositions, arXiv:2008.04937 [math.CO], 2020.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (3,3).

Cf. A028859 = a(n+2) = 2 a(n+1) + 2 a(n); A086347 = On a 3 X 3 board, number of n-move routes of chess king ending at a given side cell. a(n) = 4a(n-1) + 4a(n-2).
Cf. A128235.
Cf. A180165, A028859, A086347. - Gary W. Adamson, Aug 14 2010
Cf. A002605, A026150, A030195, A080040, A083337, A106435, A108898.
Cf. A000108, A091867, A125145, A104455, A030528, A146559.

a125145 n = a125145_list !! n
a125145_list =
   1 : 4 : map (* 3) (zipWith (+) a125145_list (tail a125145_list))
-- Reinhard Zumkeller, Oct 15 2011

I:=[1,4]; [n le 2 select I[n] else 3*Self(n-1)+3*Self(n-2): n in [1..40]]; // Vincenzo Librandi, Nov 10 2014

a[0]:=1: a[1]:=4: for n from 2 to 27 do a[n]:=3*a[n-1]+3*a[n-2] od: seq(a[n],n=0..27); # Emeric Deutsch, Feb 27 2007
A125145 := proc(n)
    option remember;
    if n <= 1 then
        op(n+1,[1,4]) ;
    else
        3*(procname(n-1)+procname(n-2)) ;
    end if;
end proc: # R. J. Mathar, Feb 13 2022

nn=23;CoefficientList[Series[(1+x)/(1-3x-3x^2),{x,0,nn}],x] (* Geoffrey Critzer, Feb 09 2014 *)
LinearRecurrence[{3,3},{1,4},30] (* Harvey P. Dale, May 01 2022 *)

G.f.: (1+z)/(1-3z-3z^2). - Emeric Deutsch, Feb 27 2007
a(n) = (5*sqrt(21)/42 + 1/2)*(3/2 + sqrt(21)/2)^n + (-5*sqrt(21)/42 + 1/2)*(3/2 - sqrt(21)/2)^n. - Antonio Alberto Olivares, Mar 20 2008
a(n) = A030195(n)+A030195(n+1). - R. J. Mathar, Feb 13 2022
E.g.f.: exp(3*x/2)*(21*cosh(sqrt(21)*x/2) + 5*sqrt(21)*sinh(sqrt(21)*x/2))/21. - Stefano Spezia, Aug 04 2022
a(n) = (((3 + sqrt(21)) / 2)^(n+2) - ((3 - sqrt(21)) / 2)^(n+2)) / (3 * sqrt(21)). - Werner Schulte, Dec 17 2024

A083100 a(n) = 2a(n-1) + 7a(n-2).

Original entry on oeis.org

1, 9, 25, 113, 401, 1593, 5993, 23137, 88225, 338409, 1294393, 4957649, 18976049, 72655641, 278143625, 1064876737, 4076758849, 15607654857, 59752621657, 228758827313, 875786006225, 3352883803641, 12836269650857, 49142725927201

Offset: 0

Mario Catalani (mario.catalani(AT)unito.it), Apr 23 2003

a(n) = a(n-1) + 8*A015519(n). a(n)/A015519(n+1) converges to sqrt(8).

a(n-1) is the number of compositions of n when there is 1 type of 1 and 8 types of other natural numbers. - Milan Janjic, Aug 13 2010

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (2, 7).

Essentially a duplicate of A084058.

The following sequences (and others) belong to the same family: A001333, A000129, A026150, A002605, A046717, A015518, A084057, A063727, A002533, A002532, A083098, A083099, A083100, A015519.

I:=[1,1]; [n le 2 select I[n] else 2*Self(n-1) + 7*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 08 2018

CoefficientList[Series[(1 + 7 x)/(1 - 2 x - 7 x^2), {x, 0, 25}], x] (* Or *) a[n_] := Simplify[((1 + Sqrt[8])^n + (1 - Sqrt[8])^n)/2]; Array[a, 25, 0] (* Or *) LinearRecurrence[{2, 7}, {1, 1}, 28] (* Or *) Table[ MatrixPower[{{1, 2}, {4, 1}}, n][[1, 1]], {n, 0, 25}] (* Robert G. Wilson v, Sep 18 2013 *)

a(n)=([0,1; 7,2]^n*[1;9])[1,1] \\ Charles R Greathouse IV, Apr 06 2016

x='x+O('x^30); Vec((1+7*x)/(1-2*x-7*x^2)) \\ G. C. Greubel, Jan 08 2018

G.f.: (1+7*x)/(1-2*x-7*x^2).
a(n) = binomial transform of 1,8,8,64,64,512. - Al Hakanson (hawkuu(AT)gmail.com), Aug 17 2009
If p[1]=1, and p[i]=8,(i>1), and if A is Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)=det A. - Milan Janjic, Apr 29 2010
G.f.: G(0)/(2*x) - 1/x, where G(k)= 1 + 1/(1 - x*(8*k-1)/(x*(8*k+7) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 28 2013

A080599 Expansion of e.g.f.: 2/(2-2*x-x^2).

Original entry on oeis.org

1, 1, 3, 12, 66, 450, 3690, 35280, 385560, 4740120, 64751400, 972972000, 15949256400, 283232149200, 5416632421200, 110988861984000, 2425817682288000, 56333385828720000, 1385151050307024000, 35950878932544576000, 982196278209226080000, 28175806418228108640000

Offset: 0

Detlef Pauly (dettodet(AT)yahoo.de), Feb 24 2003

Number of ordered partitions of {1,..,n} with at most 2 elements per block. - Bob Proctor, Apr 18 2005
In other words, number of preferential arrangements of n things (see A000670) in which each clump has size 1 or 2. - N. J. A. Sloane, Apr 13 2014
Recurrences (of the hypergeometric type of the Jovovic formula) mean: multiplying the sequence vector from the left with the associated matrix of the recurrence coefficients (here: an infinite lower triangular matrix with the natural numbers in the main diagonal and the triangular series in the subdiagonal) recovers the sequence up to an index shift. In that sense, this sequence here and many other sequences of the OEIS are eigensequences. - Gary W. Adamson, Feb 14 2011
Number of intervals in the weak (Bruhat) order of S_n that are Boolean algebras. - Richard Stanley, May 09 2011
a(n) = D^n(1/(1-x)) evaluated at x = 0, where D is the operator sqrt(1+2*x)*d/dx. Cf. A000085, A005442 and A052585. - Peter Bala, Dec 07 2011
From Gus Wiseman, Jul 04 2020: (Start)
Also the number of (1,1,1)-avoiding or cubefree sequences of length n covering an initial interval of positive integers. For example, the a(0) = 1 through a(3) = 12 sequences are:
  ()  (1)  (11)  (112)
           (12)  (121)
           (21)  (122)
                 (123)
                 (132)
                 (211)
                 (212)
                 (213)
                 (221)
                 (231)
                 (312)
                 (321)
(End)

			From _Gus Wiseman_, Jul 04 2020: (Start)
The a(0) = 1 through a(3) = 12 ordered set partitions with block sizes <= 2 are:
  {}  {{1}}  {{1,2}}    {{1},{2,3}}
             {{1},{2}}  {{1,2},{3}}
             {{2},{1}}  {{1,3},{2}}
                        {{2},{1,3}}
                        {{2,3},{1}}
                        {{3},{1,2}}
                        {{1},{2},{3}}
                        {{1},{3},{2}}
                        {{2},{1},{3}}
                        {{2},{3},{1}}
                        {{3},{1},{2}}
                        {{3},{2},{1}}
(End)

Alois P. Heinz, Table of n, a(n) for n = 0..200
S. Alex Bradt, Jennifer Elder, Pamela E. Harris, Gordon Rojas Kirby, Eva Reutercrona, Yuxuan (Susan) Wang, and Juliet Whidden, Unit interval parking functions and the r-Fubini numbers, arXiv:2401.06937 [math.CO], 2024. See page 10.
Jennifer Elder, Pamela E. Harris, Jan Kretschmann, and J. Carlos Martínez Mori, Boolean intervals in the weak order of S_n, arXiv:2306.14734 [math.CO], 2023.
Laura Gellert and Raman Sanyal, On degree sequences of undirected, directed, and bidirected graphs, arXiv preprint arXiv:1512.08448 [math.CO], 2015.
Hannah Golab, Pattern avoidance in Cayley permutations, Master's Thesis, Northern Arizona Univ. (2024). See p. 36.
Harri Hakula, Helmut Harbrecht, Vesa Kaarnioja, Frances Y. Kuo, and Ian H. Sloan, Uncertainty quantification for random domains using periodic random variables, arXiv:2210.17329 [math.NA], 2022.
Dixy Msapato, Counting the number of tau-exceptional sequences over Nakayama algebras, arXiv:2002.12194 [math.RT], 2020.
Robert A. Proctor, Let's Expand Rota's Twelvefold Way For Counting Partitions!, arXiv:math/0606404 [math.CO], 2006-2007.
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.
Index entries for related partition-counting sequences

Cf. A000085, A000670, A002530, A002605, A005442, A052585, A052611.
Cf. A090932, A102726, A106356, A232464, A333755, A335455, A335456.
Column k=2 of A276921.
Cubefree numbers are A004709.
(1,1)-avoiding patterns are A000142.
(1,1,1)-avoiding compositions are A232432.
(1,1,1)-matching patterns are A335508.
(1,1,1)-avoiding permutations of prime indices are A335511.
(1,1,1)-avoiding compositions are ranked by A335513.
(1,1,1,1)-avoiding patterns are A189886.

[n le 2 select 1 else (n-1)*Self(n-1) + Binomial(n-1,2)*Self(n-2): n in [1..31]]; // G. C. Greubel, Jan 31 2023

a:= n-> n! *(Matrix([[1,1], [1/2,0]])^n)[1,1]:
seq(a(n), n=0..40);  # Alois P. Heinz, Jun 01 2009
a:= gfun:-rectoproc({a(n) = n*a(n-1)+(n*(n-1)/2)*a(n-2),a(0)=1,a(1)=1},a(n),remember):
seq(a(n), n=0..40); # Robert Israel, Nov 01 2015

Table[n!*SeriesCoefficient[-2/(-2+2*x+x^2),{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 13 2012 *)
Round@Table[n! ((1+Sqrt[3])^(n+1) - (1-Sqrt[3])^(n+1))/(2^(n+1) Sqrt[3]), {n, 0, 20}] (* Vladimir Reshetnikov, Oct 31 2015 *)

Vec(serlaplace((2/(2-2*x-x^2) + O(x^30)))) \\ Michel Marcus, Nov 02 2015

A002605=BinaryRecurrenceSequence(2,2,0,1)
def A080599(n): return factorial(n)*A002605(n+1)/2^n
[A080599(n) for n in range(41)] # G. C. Greubel, Jan 31 2023

a(n) = n*a(n-1) + (n*(n-1)/2)*a(n-2). - Vladeta Jovovic, Aug 22 2003
E.g.f.: 1/(1-x-x^2/2). - Richard Stanley, May 09 2011
a(n) ~ n!*((1+sqrt(3))/2)^(n+1)/sqrt(3). - Vaclav Kotesovec, Oct 13 2012
a(n) = n!*((1+sqrt(3))^(n+1) - (1-sqrt(3))^(n+1))/(2^(n+1)*sqrt(3)). - Vladimir Reshetnikov, Oct 31 2015
a(n) = A090932(n) * A002530(n+1). - Robert Israel, Nov 01 2015

A214992 Power ceiling-floor sequence of (golden ratio)^4.

Original entry on oeis.org

7, 47, 323, 2213, 15169, 103969, 712615, 4884335, 33477731, 229459781, 1572740737, 10779725377, 73885336903, 506417632943, 3471038093699, 23790849022949, 163064905066945, 1117663486445665, 7660579500052711

Offset: 0

Clark Kimberling, Nov 08 2012, Jan 24 2013

Let f = floor and c = ceiling.  For x > 1, define four sequences as functions of x, as follows:
p1(0) = f(x), p1(n) = f(x*p1(n-1));
p2(0) = f(x), p2(n) = c(x*p2(n-1)) if n is odd and p2(n) = f(x*p1(n-1)) if n is even;
p3(0) = c(x), p3(n) = f(x*p3(n-1)) if n is odd and p3(n) = c(x*p3(n-1)) if n is even;
p4(0) = c(x), p4(n) = c(x*p4(n-1)).
The present sequence is given by a(n) = p3(n).
Following the terminology at A214986, call the four sequences power floor, power floor-ceiling, power ceiling-floor, and power ceiling sequences.  In the table below, a sequence is identified with an A-numbered sequence if they appear to agree except possibly for initial terms.  Notation: S(t)=sqrt(t), r = (1+S(5))/2 = golden ratio, and Limit = limit of p3(n)/p2(n).
x ......p1..... p2..... p3..... p4.......Limit
r.......A000045 A000045 A000045 A000045..r
r^2.....A001519 A001654 A061646 A001906..-1+S(5)
r^3.....A024551 A001076 A015448 A049652..-1+S(5)
r^4.....A049685 A157335 A214992 A004187..-19+9*S(5)
r^5.....A214993 A049666 A015457 A214994...(-9+5*S(5))/2
r^6.....A007805 A156085 A214995 A049660..-151+68*S(5)
1+S(2)..A024537 A000129 A001333 A048739..S(2)
2+S(2)..A007052 A214996 A214997 A007070..(1+S(2))/2
1+S(3)..A057960 A002605 A028859 A077846..(1+S(3))/2
2+S(3)..A001835 A109437 A214998 A001353..-4+3*S(3)
S(5)....A214999 A215091 A218982 A218983..1.26879683...
2+S(5)..A024551 A001076 A015448 A049652..-1+S(5)
2+S(6)..A218984 A090017 A123347 A218985..S(3/2)
2+S(7)..A218986 A015530 A126473 A218987..(1+S(7))/3
2+S(8)..A218988 A057087 A086347 A218989..(1+S(2))/2
3+S(8)..A001653 A084158 A218990 A001109..-13+10*S(2)
3+S(10).A218991 A005668 A015451 A218992..-2+S(10)
...
Properties of p1, p2, p3, p4:
(1) If x > 2, the terms of p2 and p3 interlace: p2(0) < p3(0) < p2(1) < p3(1) < p2(2) < p3(2)... Also, p1(n) <= p2(n) <= p3(n) <= p4(n) <= p1(n+1) for all x>0 and n>=0.
(2) If x > 2, the limits L(x) = limit(p/x^n) exist for the four functions p(x), and L1(x) <= L2(x) <= L3(x) <= L4 (x). See the Mathematica programs for plots of the four functions; one of them also occurs in the Odlyzko and Wilf article, along with a discussion of the special case x = 3/2.
(3) Suppose that x = u + sqrt(v) where v is a nonsquare positive integer.  If u = f(x) or u = c(x), then p1, p2, p3, p4 are linear recurrence sequences.  Is this true for sequences p1, p2, p3, p4 obtained from x = (u + sqrt(v))^q for every positive integer q?
(4) Suppose that x is a Pisot-Vijayaraghavan number. Must p1, p2, p3, p4 then be linearly recurrent?  If x is also a quadratic irrational b + c*sqrt(d), must the four limits L(x) be in the field Q(sqrt(d))?
(5) The Odlyzko and Wilf article (page 239) raises three interesting questions about the power ceiling function; it appears that they remain open.

			a(0) = ceiling(r) = 7, where r = ((1+sqrt(5))/2)^4 = 6.8...; a(1) = floor(7*r) = 47; a(2) = ceiling(47) = 323.

Clark Kimberling, Table of n, a(n) for n = 0..250
A. M. Odlyzko and H. S. Wilf, Functional iteration and the Josephus problem, Glasgow Math. J. 33, 235-240, 1991.
Index entries for linear recurrences with constant coefficients, signature (6,6,-1).
Index entries for sequences related to the Josephus Problem.

Cf. A001622, A214986, A049685, A157335, A004187.

(* Program 1.  A214992 and related sequences *)
x = GoldenRatio^4; z = 30; (* z = # terms in sequences *)
z1 = 100; (* z1 = # digits in approximations *)
f[x_] := Floor[x]; c[x_] := Ceiling[x];
p1[0] = f[x]; p2[0] = f[x]; p3[0] = c[x]; p4[0] = c[x];
p1[n_] := f[x*p1[n - 1]]
p2[n_] := If[Mod[n, 2] == 1, c[x*p2[n - 1]], f[x*p2[n - 1]]]
p3[n_] := If[Mod[n, 2] == 1, f[x*p3[n - 1]], c[x*p3[n - 1]]]
p4[n_] := c[x*p4[n - 1]]
Table[p1[n], {n, 0, z}]  (* A049685 *)
Table[p2[n], {n, 0, z}]  (* A157335 *)
Table[p3[n], {n, 0, z}]  (* A214992 *)
Table[p4[n], {n, 0, z}]  (* A004187 *)
Table[p4[n] - p1[n], {n, 0, z}]  (* A004187 *)
Table[p3[n] - p2[n], {n, 0, z}]  (* A098305 *)
(* Program 2.  Plot of power floor and power ceiling functions, p1(x) and p4(x) *)
f[x_] := f[x] = Floor[x]; c[x_] := c[x] = Ceiling[x];
p1[x_, 0] := f[x]; p1[x_, n_] := f[x*p1[x, n - 1]];
p4[x_, 0] := c[x]; p4[x_, n_] := c[x*p4[x, n - 1]];
Plot[Evaluate[{p1[x, 10]/x^10, p4[x, 10]/x^10}], {x, 2, 3}, PlotRange -> {0, 4}]
(* Program 3. Plot of power floor-ceiling and power ceiling-floor functions, p2(x) and p3(x) *)
f[x_] := f[x] = Floor[x]; c[x_] := c[x] = Ceiling[x];
p2[x_, 0] := f[x]; p3[x_, 0] := c[x];
p2[x_, n_] := If[Mod[n, 2] == 1, c[x*p2[x, n - 1]], f[x*p2[x, n - 1]]]
p3[x_, n_] := If[Mod[n, 2] == 1, f[x*p3[x, n - 1]], c[x*p3[x, n - 1]]]
Plot[Evaluate[{p2[x, 10]/x^10, p3[x, 10]/x^10}], {x, 2, 3}, PlotRange -> {0, 4}]

a(n) = floor(r*a(n-1)) if n is odd and a(n) = ceiling(r*a(n-1)) if n is even, where a(0) = ceiling(r), r = (golden ratio)^4 = (7 + sqrt(5))/2.
a(n) = 6*a(n-1) + 6*a(n-2) - a(n-3).
G.f.: (7 + 5*x - x^2)/((1 + x)*(1 - 7*x + x^2)).
a(n) = (10*(-2)^n+(10+3*sqrt(5))*(7-3*sqrt(5))^(n+2)+(10-3*sqrt(5))*(7+3*sqrt(5))^(n+2))/(90*2^n). - Bruno Berselli, Nov 14 2012
a(n) = 7*A157335(n) + 5*A157335(n-1) - A157335(n-2). - R. J. Mathar, Feb 05 2020
E.g.f.: exp(-x)*(5 + 2*exp(9*x/2)*(155*cosh(3*sqrt(5)*x/2) + 69*sqrt(5)*sinh(3*sqrt(5)*x/2)))/45. - Stefano Spezia, Oct 28 2024

A063967 Triangle read by rows, T(n,k) = T(n-1,k) + T(n-2,k) + T(n-1,k-1) + T(n-2,k-1) and T(0,0) = 1.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 3, 7, 5, 1, 5, 15, 16, 7, 1, 8, 30, 43, 29, 9, 1, 13, 58, 104, 95, 46, 11, 1, 21, 109, 235, 271, 179, 67, 13, 1, 34, 201, 506, 705, 591, 303, 92, 15, 1, 55, 365, 1051, 1717, 1746, 1140, 475, 121, 17, 1, 89, 655, 2123, 3979, 4759, 3780, 2010, 703, 154, 19, 1

Offset: 0

Henry Bottomley, Sep 05 2001

			T(3,1) = T(2,1) + T(1,1) + T(2,0) + T(1,0) = 3 + 1 + 2 + 1 = 7.
Triangle begins:
   1,
   1,   1,
   2,   3,   1,
   3,   7,   5,   1,
   5,  15,  16,   7,   1,
   8,  30,  43,  29,   9,   1,
  13,  58, 104,  95,  46,  11,  1,
  21, 109, 235, 271, 179,  67, 13,  1,
  34, 201, 506, 705, 591, 303, 92, 15, 1

Reinhard Zumkeller, Rows n = 0..120 of triangle, flattened
E. Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, Advances in Mathematics, 34 (2005) pp. 101-122.
Emanuele Munarini, A generalization of André-Jeannin's symmetric identity, Pure Mathematics and Applications (2018) Vol. 27, No. 1, 98-118.

Row sums are A002605.
Columns include: A000045(n+1), A023610(n-1).
Main diagonal: A000012, a(n, n-1) = A005408(n-1).
Matrix inverse: A091698, matrix square: A091700.
Cf. A321620.
Sum_{k=0..n} x^k*T(n,k) is (-1)^n*A057086(n) (x=-11), (-1)^n*A057085(n+1) (x=-10), (-1)^n*A057084(n) (x=-9), (-1)^n*A030240(n) (x=-8), (-1)^n*A030192(n) (x=-7), (-1)^n*A030191(n) (x=-6), (-1)^n*A001787(n+1) (x=-5), A000748(n) (x=-4), A108520(n) (x=-3), A049347(n) (x=-2), A000007(n) (x=-1), A000045(n) (x=0), A002605(n) (x=1), A030195(n+1) (x=2), A057087(n) (x=3), A057088(n) (x=4), A057089(n) (x=5), A057090(n) (x=6), A057091(n) (x=7), A057092(n) (x=8), A057093(n) (x=9). - Philippe Deléham, Nov 03 2006

a063967_tabl = [1] : [1,1] : f [1] [1,1] where
   f us vs = ws : f vs ws where
     ws = zipWith (+) ([0] ++ us ++ [0]) $
          zipWith (+) (us ++ [0,0]) $ zipWith (+) ([0] ++ vs) (vs ++ [0])
-- Reinhard Zumkeller, Apr 17 2013

T[n_, k_] := Sum[Binomial[j, n - j]*Binomial[j, k], {j, 0, n}]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Oct 11 2017, after Paul Barry *)
(* Function RiordanSquare defined in A321620. *)
RiordanSquare[1/(1 - x - x^2), 11] // Flatten (* Peter Luschny, Nov 27 2018 *)

G.f.: 1/(1-x*(1+x)*(1+y)). - Vladeta Jovovic, Oct 11 2003
Riordan array (1/(1-x-x^2), x(1+x)/(1-x-x^2)). The inverse of the signed version (1/(1+x-x^2),x(1-x)/(1+x-x^2)) is abs(A091698). - Paul Barry, Jun 10 2005
T(n, k) = Sum_{j=0..n} C(j, n-j)C(j, k). - Paul Barry, Nov 09 2005
Diagonal sums are A002478. - Paul Barry, Nov 09 2005
A026729*A007318 as infinite lower triangular matrices. - Philippe Deléham, Dec 11 2008
Central coefficients T(2*n,n) are A137644. - Paul Barry, Apr 15 2010
Product of Riordan arrays (1, x(1+x))*(1/(1-x), x/(1-x)), that is, A026729*A007318. - Paul Barry, Mar 14 2011
Triangle T(n,k), read by rows, given by (1,1,-1,0,0,0,0,0,0,0,...) DELTA (1,0,0,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 12 2011

A073370 Convolution triangle of A001045(n+1) (generalized (1,2)-Fibonacci), n>=0.

Original entry on oeis.org

1, 1, 1, 3, 2, 1, 5, 7, 3, 1, 11, 16, 12, 4, 1, 21, 41, 34, 18, 5, 1, 43, 94, 99, 60, 25, 6, 1, 85, 219, 261, 195, 95, 33, 7, 1, 171, 492, 678, 576, 340, 140, 42, 8, 1, 341, 1101, 1692, 1644, 1106, 546, 196, 52, 9, 1

Offset: 0

Wolfdieter Lang, Aug 02 2002

The g.f. for the row polynomials P(n,x) = Sum_{m=0..n} T(n,m)*x^m is 1/(1-(1+x+2*z)*z). See Shapiro et al. reference and comment under A053121 for such convolution triangles.
Riordan array (1/(1-x-2*x^2), x/(1-x-2*x^2)). - Paul Barry, Mar 15 2005
Subtriangle (obtained by dropping the first column) of the triangle given by (0, 1, 2, -2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Feb 19 2013
The number of ternary words of length n having k letters equal 2 and 0,1 avoid runs of odd lengths. - Milan Janjic, Jan 14 2017

			Triangle begins as:
    1;
    1,   1;
    3,   2,   1;
    5,   7,   3,   1;
   11,  16,  12,   4,   1;
   21,  41,  34,  18,   5,   1;
   43,  94,  99,  60,  25,   6,   1;
   85, 219, 261, 195,  95,  33,   7,   1;
  171, 492, 678, 576, 340, 140,  42,   8,   1;
The triangle (0, 1, 2, -2, 0, 0, ...) DELTA (1, 0, 0, 0, 0, ...) begins:
  1;
  0,  1;
  0,  1,  1;
  0,  3,  2,  1;
  0,  5,  7,  3,  1;
  0, 11, 16, 12,  4,  1;
  0, 21, 41, 34, 18,  5,  1; - _Philippe Deléham_, Feb 19 2013

G. C. Greubel, Rows n = 0..50 of the triangle, flattened
W. Lang, First 10 rows.
Milan Janjić, Words and Linear Recurrences, J. Int. Seq. 21 (2018), #18.1.4.

Columns: A001045 (k=0), A073371 (k=1), A073372 (k=2), A073373 (k=3), A073374 (k=4), A073375 (k=5), A073376 (k=6), A073377 (k=7), A073378 (k=8), A073379 (k=9).
Cf. A002605 (row sums), A006130 (diagonal sums), A073399, A073400.
Cf. A000045, A077957.

A073370:= func< n,k | (&+[Binomial(n-j,k)*Binomial(n-k-j,j)*2^j: j in [0..Floor((n-k)/2)]]) >;
[A073370(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 01 2022

# Uses function PMatrix from A357368. Adds a row above and a column to the left.
PMatrix(10, n -> (2^n - (-1)^n) / 3); # Peter Luschny, Oct 07 2022

T[n_, k_]:= T[n, k]= Sum[Binomial[n-j,k]*Binomial[n-k-j,j]*2^j, {j,0,Floor[(n- k)/2]}];
Table[T[n, k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Oct 01 2022 *)

def A073370(n,k): return binomial(n,k)*sum( 2^j * binomial(2*j,j) * binomial(n-k,2*j)/binomial(n,j) for j in range(1+(n-k)//2))
flatten([[A073370(n,k) for k in range(n+1)] for n in range(12)]) # G. C. Greubel, Oct 01 2022

T(n, m) = Sum_{k=0..floor((n-m)/2)} binomial(n-k, m)*binomial(n-m-k, k)*2^k, if n > m, else 0.
Sum_{k=0..n} T(n, k) = A002605(n+1).
T(n, m) = (1*(n-m+1)*T(n, m-1) + 2*2*(n+m)*T(n-1, m-1))/((1^2 + 4*2)*m), n >= m >= 1, T(n, 0) = A001045(n+1), n >= 0, else 0.
T(n, m) = (p(m-1, n-m)*1*(n-m+1)*T(n-m+1) + q(m-1, n-m)*2*(n-m+2)*T(n-m))/(m!*9^m), n >= m >= 1, with T(n) = T(n, m=0) = A001045(n+1), else 0; p(k, n) = Sum_{j=0..k} (A(k, j)*n^(k-j) and q(k, n) = Sum_{j=0..k} B(k, j)*n^(k-j), with the number triangles A(k, m) = A073399(k, m) and B(k, m) = A073400(k, m).
G.f.: 1/(1-(1+2*x)*x)^(m+1) = 1/((1+x)*(1-2*x))^(m+1), m >= 0, for column m (without leading zeros).
T(n, 0) = A001045(n), T(1, 1) = 1, T(n, k) = 0 if k>n, T(n, k) = T(n-1, k-1) + 2*T(n-2, k) + T(n-1, k) otherwise. - Paul Barry, Mar 15 2005
G.f.: (1+x)*(1-2*x)/(1-x-2*x^2-x*y) for the triangle including the 1, 0, 0, 0, 0, ... column. - R. J. Mathar, Aug 11 2015
From Peter Bala, Oct 07 2019: (Start)
Recurrence for row polynomials: R(n,x) = (1 + x)*R(n-1,x) + 2*R(n-2,x) with R(0,x) = 1 and R(1,x) = 1 + x.
The row reverse polynomial x^n*R(n,1/x) is equal to the numerator polynomial of the finite continued fraction 1 + x/(1 - 2*x/(1 + ... + x/(1 - 2*x/(1)))) (with 2*n partial numerators). Cf. A110441. (End)
From G. C. Greubel, Oct 01 2022: (Start)
T(n, k) = binomial(n,k)*Sum_{j=0..floor((n-k)/2)} 2^j*binomial(2*j, j)*binomial(n-k, 2*j)/binomial(n, j).
T(n, k) = binomial(n, k)*Hypergeometric2F1([(k-n)/2, (k-n+1)/2], [-2*n], -8).
Sum_{k=0..n} (-1)^k * T(n, k) = A077957(n).
Sum_{k=0..floor(n/2)} T(n-k, k) = A006130(n).
Sum_{k=0..floor(n/2)} (-1)^k * T(n-k, k) = A000045(n+1). (End)

A038208 Triangle whose (i,j)-th entry is binomial(i,j)*2^i.

Original entry on oeis.org

1, 2, 2, 4, 8, 4, 8, 24, 24, 8, 16, 64, 96, 64, 16, 32, 160, 320, 320, 160, 32, 64, 384, 960, 1280, 960, 384, 64, 128, 896, 2688, 4480, 4480, 2688, 896, 128, 256, 2048, 7168, 14336, 17920, 14336, 7168, 2048, 256, 512, 4608, 18432, 43008, 64512, 64512, 43008, 18432, 4608, 512

Offset: 0

N. J. A. Sloane

Triangle obtained from expansion of (2 + 2*x)^n.

			    1;
    2,    2;
    4,    8,    4;
    8,   24,   24,     8;
   16,   64,   96,    64,    16;
   32,  160,  320,   320,   160,   32;
   64,  384,  960,  1280,   960,   384,   64;
  128,  896, 2688,  4480,  4480,  2688,  896,  128;
  256, 2048, 7168, 14336, 17920, 14336, 7168, 2048, 256;

Harvey P. Dale, Table of n, a(n) for n = 0..1000
B. N. Cyvin et al., Isomer enumeration of unbranched catacondensed polygonal systems with pentagons and heptagons, Match, No. 34 (Oct 1996), pp. 109-121.
S. Elizalde, E. Deutsch, The degree of asymmetry of a sequence, Enum. Combinat. Applic. 2 (2022) no 1 #S2R7, corollary 2.2 at m=2, then dropping each 2nd row
Franck Ramaharo, Statistics on some classes of knot shadows, arXiv:1802.07701 [math.CO], 2018.
Franck Ramaharo, A generating polynomial for the pretzel knot, arXiv:1805.10680 [math.CO], 2018.
Franck Ramaharo, A bracket polynomial for 2-tangle shadows, arXiv:2002.06672 [math.CO], 2020.

Cf. A000079, A000302 (row sums), A002605 (diagonal sums), A027306.

[Binomial(n,k)*2^n: k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 17 2018

nn=8;Map[Select[#,#>0&]&,Transpose[Table[Range[0,nn]!CoefficientList[Series[2^k x^k/k! Exp[2x],{x,0,nn}],x],{k,0,nn}]]]//Grid (* Geoffrey Critzer, Feb 13 2014 *)
Flatten[Table[Binomial[i,j]2^i,{i,0,10},{j,0,i}]] (* Harvey P. Dale, May 28 2015 *)

for(n=0,10, for(k=0,n, print1(binomial(n,k)*2^n, ", "))) \\ G. C. Greubel, Oct 17 2018

flatten([[binomial(n,k)*2^n for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Mar 21 2022

E.g.f. for column k: 2^k*x^k/k!*exp(2*x). - Geoffrey Critzer, Feb 13 2014
From G. C. Greubel, Mar 21 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n).
Sum_{k=0..n} T(n, k) = A000302(n).
Sum_{k=0..floor(n/2)} T(n-k, k) = A002605(n+1).
Sum_{k=0..floor(n/2)} T(n, k) = 2^n*A027306(n). (End)

A028860 a(n+2) = 2a(n+1) + 2a(n); a(0) = -1, a(1) = 1.

Original entry on oeis.org

-1, 1, 0, 2, 4, 12, 32, 88, 240, 656, 1792, 4896, 13376, 36544, 99840, 272768, 745216, 2035968, 5562368, 15196672, 41518080, 113429504, 309895168, 846649344, 2313089024, 6319476736, 17265131520, 47169216512, 128868696064, 352075825152, 961889042432, 2627929735168

Offset: 0

N. J. A. Sloane

a(n+1) is the top left entry of the n-th power of the 3 X 3 matrix [0, 1, 1; 1, 1, 1; 1, 1, 1]. - R. J. Mathar, Feb 04 2014

(A002605, a(.+1)) is the canonical basis of the space of linear recurrent sequences with signature (2, 2), i.e., any sequence s(n) = 2(s(n-1) + s(n-2)) is given by s = s(0)*A002605 + s(1)*a(.+1). - M. F. Hasler, Aug 06 2018

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Martin Burtscher, Igor Szczyrba, and Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 924
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (2,2).

Cf. A002605, A026150, A030195, A080040, A083337, A106435, A108898, A125145.

a:=[-1,1];; for n in [3..30] do a[n]:=2*a[n-1]+2*a[n-2]; od; a; # Muniru A Asiru, Aug 07 2018

a028860 n = a028860_list !! n
a028860_list =
   -1 : 1 : map (* 2) (zipWith (+) a028860_list (tail a028860_list))
-- Reinhard Zumkeller, Oct 15 2011

I:=[-1,1]; [n le 2 select I[n] else 2*Self(n-1)+2*Self(n-2): n in [1..40]]; // Vincenzo Librandi, Aug 13 2018

seq(coeff(series((3*x-1)/(1-2*x-2*x^2), x,n+1),x,n),n=0..30); # Muniru A Asiru, Aug 07 2018

(With a different offset) M = {{0, 2}, {1, 2}} v[1] = {0, 1} v[n_] := v[n] = M.v[n - 1] a = Table[Abs[v[n][[1]]], {n, 1, 50}] (* Roger L. Bagula, May 29 2005 *)
LinearRecurrence[{2,2},{-1,1},40] (* Harvey P. Dale, Dec 13 2012 *)
CoefficientList[Series[(-3 x + 1)/(2 x^2 + 2 x - 1), {x, 0, 27}], x] (* Robert G. Wilson v, Aug 07 2018 *)

apply( A028860(n)=([2,2;1,0]^n)[2,]*[1,-1]~, [0..30]) \\ 15% faster than (A^n*[1,-1]~)[2]. - M. F. Hasler, Aug 06 2018

A028860 = BinaryRecurrenceSequence(2,2,-1,1)
[A028860(n) for n in range(51)] # G. C. Greubel, Dec 08 2022

a(n) = 4*A028859(n-4), for n > 3.
From R. J. Mathar, Nov 27 2008: (Start)
G.f.: -(1 - 3*x)/(1 - 2*x - 2*x^2).
a(n) = 3*A002605(n-1) - A002605(n). (End)
a(n) = det A, where A is the Hessenberg matrix of order n+1 defined by: A[i,j] = p(j - i + 1) (i <= j), A[i,j] = -1 (i = j + 1), A[i,j] = 0 otherwise, with p(i) = fibonacci(2i - 4). - Milan Janjic, May 08 2010, edited by M. F. Hasler, Aug 06 2018
a(n) = (2*sqrt(3) - 3)/6*(1 + sqrt(3))^n - (2*sqrt(3) + 3)/6*(1 - sqrt(3))^n. - Sergei N. Gladkovskii, Jul 18 2012
a(n) = 2*A002605(n-2) for n >= 2. - M. F. Hasler, Aug 06 2018
E.g.f.: exp(x)*(2*sqrt(3)*sinh(sqrt(3)*x) - 3*cosh(sqrt(3)*x))/3. - Franck Maminirina Ramaharo, Nov 11 2018

Edited by N. J. A. Sloane, Apr 11 2009

Edited and initial values added in definition by M. F. Hasler, Aug 06 2018

A367298 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 2 + 4x, p(n,x) = up(n-1,x) + vp(n-2,x) for n >= 3, where u = p(2,x), v = 1 - 2x - x^2.

Original entry on oeis.org

1, 2, 4, 5, 14, 15, 12, 48, 76, 56, 29, 148, 326, 372, 209, 70, 436, 1212, 1904, 1718, 780, 169, 1242, 4169, 8228, 10191, 7642, 2911, 408, 3456, 13576, 32176, 49992, 51488, 33112, 10864, 985, 9448, 42492, 117304, 218254, 281976, 249612, 140712, 40545

Offset: 1

Clark Kimberling, Nov 26 2023

Because (p(n,x)) is a strong divisibility sequence, for each integer k, the sequence (p(n,k)) is a strong divisibility sequence of integers.

			First eight rows:
    1
    2     4
    5    14     15
   12    48     76     56
   29   148    326    372    209
   70   436   1212   1904   1718   780
  169  1242   4169   8228  10191  7642    2911
  408  3456  13576  32176  49992  51488  33112  10864
Row 4 represents the polynomial p(4,x) = 12 + 48*x + 76*x^2 + 56*x^3, so (T(4,k)) = (12,48,76,56), k=0..3.

Rigoberto Flórez, Robinson Higuita, and Antara Mukherjee, Characterization of the strong divisibility property for generalized Fibonacci polynomials, Integers, 18 (2018), Paper No. A14.

Cf. A000129 (column 1), A001353 (p(n,n-1)), A154244 (row sums, p(n,1)), A002605 (alternating row sums, p(n,-1)), A190989 (p(n,2)), A005668 (p(n,-2)), A190869 (p(n,-3)), A094440, A367208, A367209, A367210, A367211, A367297, A367299, A367300, A367301.

p[1, x_] := 1; p[2, x_] := 2 + 4 x; u[x_] := p[2, x]; v[x_] := 1 - 2 x - x^2;
p[n_, x_] := Expand[u[x]*p[n - 1, x] + v[x]*p[n - 2, x]]
Grid[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]
Flatten[Table[CoefficientList[p[n, x], x], {n, 1, 10}]]

p(n,x) = u*p(n-1,x) + v*p(n-2,x) for n >= 3, where p(1,x) = 1, p(2,x) = 2 + 4*x, u = p(2,x), and v = 1 - 2*x - x^2.

p(n,x) = k*(b^n - c^n), where k = -(1/sqrt(8 + 8*x + 12*x^2)), b = (1/2)*(4*x + 2 + 1/k), c = (1/2)*(4*x + 2 - 1/k).

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

Original entry on oeis.org

0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 0, 0, 0, -1, 0, -1, 0

Offset: 0

Paul Barry, Jul 14 2005

Transform of A002605 by the Riordan array A102587. Denominator is the 12th cyclotomic polynomial.

Antti Karttunen, Table of n, a(n) for n = 0..65537
Index entries for linear recurrences with constant coefficients, signature (0,1,0,-1).

Cf. A002605, A010892, A014021, A102587, A322796.

A110161:= func< n | KroneckerSymbol(12, n) >;
[A110161(n): n in [0..120]]; // G. C. Greubel, Oct 23 2024

a[ n_] := JacobiSymbol[ 12, n]; (* Michael Somos, Jan 29 2015 *)
LinearRecurrence[{0,1,0,-1},{0,1,0,0},110] (* Harvey P. Dale, Jul 11 2015 *)

{a(n) = kronecker( 12, n)}; /* Michael Somos, Jun 11 2007 */

def A110161(n): return kronecker(12, n)
[A110161(n) for n in range(121)] # G. C. Greubel, Oct 23 2024

Periodic of length 12: 0, 1, 0, 0, 0, -1, 0, -1, 0, 0, 0, 1. - T. D. Noe, Dec 12 2006
From Michael Somos, Jun 11 2007: (Start)
Euler transform of length 12 sequence [0, 0, 0, -1, 0, -1, 0, 0, 0, 0, 0, 1].
a(n) is multiplicative with a(2^e) = a(3^e) = 0^e, a(p^e) = 1 if p == 1, 11 (mod 12), a(p^e) = (-1)^e if p == 5, 7 (mod 12).
a(n) = a(-n) = -a(n + 6) for all n in Z.
G.f.: x * (1 - x^4) * (1 - x^6) / (1 - x^12). (End)
a(2*n - 1) = A010892(n). - Michael Somos, Jan 29 2015
a(n) = A014021(n+1). - R. J. Mathar, Nov 13 2023

Corrected by T. D. Noe, Dec 12 2006

cp's OEIS Frontend

A125145 a(n) = 3a(n-1) + 3a(n-2). a(0) = 1, a(1) = 4.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Formula

A083100 a(n) = 2*a(n-1) + 7*a(n-2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

A080599 Expansion of e.g.f.: 2/(2-2*x-x^2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A214992 Power ceiling-floor sequence of (golden ratio)^4.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A063967 Triangle read by rows, T(n,k) = T(n-1,k) + T(n-2,k) + T(n-1,k-1) + T(n-2,k-1) and T(0,0) = 1.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A073370 Convolution triangle of A001045(n+1) (generalized (1,2)-Fibonacci), n>=0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

SageMath

Formula

A038208 Triangle whose (i,j)-th entry is binomial(i,j)*2^i.

Views

A083100 a(n) = 2a(n-1) + 7a(n-2).

A028860 a(n+2) = 2a(n+1) + 2a(n); a(0) = -1, a(1) = 1.

A367298 Triangular array T(n,k), read by rows: coefficients of strong divisibility sequence of polynomials p(1,x) = 1, p(2,x) = 2 + 4x, p(n,x) = up(n-1,x) + vp(n-2,x) for n >= 3, where u = p(2,x), v = 1 - 2x - x^2.