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-10 of 14 results. Next

A030528 Triangle read by rows: a(n,k) = binomial(k,n-k).

Original entry on oeis.org

1, 1, 1, 0, 2, 1, 0, 1, 3, 1, 0, 0, 3, 4, 1, 0, 0, 1, 6, 5, 1, 0, 0, 0, 4, 10, 6, 1, 0, 0, 0, 1, 10, 15, 7, 1, 0, 0, 0, 0, 5, 20, 21, 8, 1, 0, 0, 0, 0, 1, 15, 35, 28, 9, 1, 0, 0, 0, 0, 0, 6, 35, 56, 36, 10, 1, 0, 0, 0, 0, 0, 1, 21, 70, 84, 45, 11, 1, 0, 0, 0, 0, 0, 0, 7, 56, 126, 120, 55, 12, 1
Offset: 1

Views

Author

Keywords

Comments

A convolution triangle of numbers obtained from A019590.
a(n,m) := s1(-1; n,m), a member of a sequence of triangles including s1(0; n,m)= A023531(n,m) (unit matrix) and s1(2; n,m)= A007318(n-1,m-1) (Pascal's triangle).
The signed triangular matrix a(n,m)*(-1)^(n-m) is the inverse matrix of the triangular Catalan convolution matrix A033184(n+1,m+1), n >= m >= 0, with A033184(n,m) := 0 if n
Riordan array (1+x, x(1+x)). The signed triangle is the Riordan array (1-x,x(1-x)), inverse to (c(x),xc(x)) with c(x) g.f. for A000108. - Paul Barry, Feb 02 2005 [with offset 0]
Also, a(n,k)=number of compositions of n into k parts of 1's and 2's. Example: a(6,4)=6 because we have 2211, 2121, 2112, 1221, 1212 and 1122. - Emeric Deutsch, Apr 05 2005 [see MacMahon and Riordan. - Wolfdieter Lang, Jul 27 2023]
Subtriangle of A026729. - Philippe Deléham, Aug 31 2006
a(n,k) is the number of length n-1 binary sequences having no two consecutive 0's with exactly k-1 1's. Example: a(6,4)=6 because we have 01011, 01101, 01110, 10101, 10110, 11010. - Geoffrey Critzer, Jul 22 2013
Mirrored, shifted Fibonacci polynomials of A011973. The polynomials (illustrated below) of this entry have the property that p(n,t) = t * [p(n-1,t) + p(n-2,t)]. The additive properties of Pascal's triangle (A007318) are reflected in those of these polynomials, as can be seen in the Example Section below and also when the o.g.f. G(x,t) below is expanded as the series x*(1+x) + t * [x*(1+x)]^2 + t^2 * [x*(1+x)]^3 + ... . See also A053122 for a relation to Cartan matrices. - Tom Copeland, Nov 04 2014
Rows of this entry appear as columns of an array for an infinitesimal generator presented in the Copeland link. - Tom Copeland, Dec 23 2015
For n >= 2, the n-th row is also the coefficients of the vertex cover polynomial of the (n-1)-path graph P_{n-1}. - Eric W. Weisstein, Apr 10 2017
With an additional initial matrix element a_(0,0) = 1 and column of zeros a_(n,0) = 0 for n > 0, these are antidiagonals read from bottom to top of the numerical coefficients of the Maurer-Cartan form matrix of the Leibniz group L^(n)(1,1) presented on p. 9 of the Olver paper, which is generated as exp[c. * M] with (c.)^n = c_n and M the Lie infinitesimal generator A218272. Cf. A011973. And A169803. - Tom Copeland, Jul 02 2018

Examples

			Triangle starts:
  [ 1]  1
  [ 2]  1   1
  [ 3]  0   2   1
  [ 4]  0   1   3   1
  [ 5]  0   0   3   4   1
  [ 6]  0   0   1   6   5   1
  [ 7]  0   0   0   4  10   6   1
  [ 8]  0   0   0   1  10  15   7   1
  [ 9]  0   0   0   0   5  20  21   8   1
  [10]  0   0   0   0   1  15  35  28   9   1
  [11]  0   0   0   0   0   6  35  56  36  10   1
  [12]  0   0   0   0   0   1  21  70  84  45  11   1
  [13]  0   0   0   0   0   0   7  56 126 120  55  12   1
  ...
From _Tom Copeland_, Nov 04 2014: (Start)
For quick comparison to other polynomials:
  p(1,t) = 1
  p(2,t) = 1 + 1 t
  p(3,t) = 0 + 2 t + 1 t^2
  p(4,t) = 0 + 1 t + 3 t^2 + 1 t^3
  p(5,t) = 0 + 0   + 3 t^2 + 4 t^3 +  1 t^4
  p(6,t) = 0 + 0   + 1 t^2 + 6 t^3 +  5 t^4 +  1 t^5
  p(7,t) = 0 + 0   + 0     + 4 t^3 + 10 t^4 +  6 t^5 + 1 t^6
  p(8,t) = 0 + 0   + 0     + 1 t^3 + 10 t^4 + 15 t^5 + 7 t^6 + 1 t^7
  ...
Reading along columns gives rows for Pascal's triangle. (End)
		

References

  • P. A. MacMahon, Combinatory Analysis, Two volumes (bound as one), Chelsea Publishing Company, New York, 1960, Vol. I, nr. 124, p. 151.
  • John Riordan, An Introduction to Combinatorial Analysis, John Wiley & Sons, London, 1958. eq. (35), p.124, 11. p. 154.

Crossrefs

Row sums A000045(n+1) (Fibonacci). a(n, 1)= A019590(n) (Fermat's last theorem). Cf. A049403.

Programs

  • Magma
    /* As triangle */ [[Binomial(k, n-k): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Nov 05 2014
  • Maple
    for n from 1 to 12 do seq(binomial(k,n-k),k=1..n) od; # yields sequence in triangular form - Emeric Deutsch, Apr 05 2005
  • Mathematica
    nn=10;CoefficientList[Series[(1+x)/(1-y x - y x^2),{x,0,nn}],{x,y}]//Grid (* Geoffrey Critzer, Jul 22 2013 *)
    Table[Binomial[k, n - k], {n, 13}, {k, n}] // Flatten (* Michael De Vlieger, Dec 23 2015 *)
    CoefficientList[Table[x^(n/2 - 1) Fibonacci[n + 1, Sqrt[x]], {n, 10}],
       x] // Flatten (* Eric W. Weisstein, Apr 10 2017 *)

Formula

a(n, m) = 2*(2*m-n+1)*a(n-1, m)/n + m*a(n-1, m-1)/n, n >= m >= 1; a(n, m) := 0, n
G.f. for m-th column: (x*(1+x))^m.
As a number triangle with offset 0, this is T(n, k) = Sum_{i=0..n} (-1)^(n+i)*binomial(n, i)*binomial(i+k+1, 2k+1). The antidiagonal sums give the Padovan sequence A000931(n+5). Inverse binomial transform of A078812 (product of lower triangular matrices). - Paul Barry, Jun 21 2004
G.f.: (1 + x)/(1 - y*x - y*x^2). - Geoffrey Critzer, Jul 22 2013 [offset 0] [with offset 1: g.f. of row polynomials in y: x*(1+x)*y/(1 - x*(1+x)*y). - Wolfdieter Lang, Jul 27 2023]
From Tom Copeland, Nov 04 2014: (Start)
O.g.f: G(x,t) = x*(1+x) / [1 - t*x*(1+x)] = -P[Cinv(-x),t], where P(x,t)= x / (1 + t*x) and Cinv(x)= x*(1-x) are the compositional inverses in x of Pinv(x,t) = -P(-x,t) = x / (1 - t*x) and C(x) = [1-sqrt(1-4*x)]/2, an o.g.f. for the shifted Catalan numbers A000108.
Therefore, Ginv(x,t) = -C[Pinv(-x,t)] = {-1 + sqrt[1 + 4*x/(1+t*x)]}/2, which is -A124644(-x,t).
This places this array in a family of arrays related by composition of P and C and their inverses and interpolation by t, such as A091867 and A104597, and associated to the Catalan, Motzkin, Fine, and Fibonacci numbers. Cf. A104597 (polynomials shifted in t) A125145, A146559, A057078, A000045, A155020, A125145, A039717, A001792, A057862, A011973, A115139. (End)

Extensions

More terms from Emeric Deutsch, Apr 05 2005

A028859 a(n+2) = 2*a(n+1) + 2*a(n); a(0) = 1, a(1) = 3.

Original entry on oeis.org

1, 3, 8, 22, 60, 164, 448, 1224, 3344, 9136, 24960, 68192, 186304, 508992, 1390592, 3799168, 10379520, 28357376, 77473792, 211662336, 578272256, 1579869184, 4316282880, 11792304128, 32217174016, 88018956288, 240472260608, 656982433792, 1794909388800, 4903783645184, 13397386067968
Offset: 0

Keywords

Comments

Number of words of length n without adjacent 0's from the alphabet {0,1,2}. For example, a(2) counts 01,02,10,11,12,20,21,22. - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Jun 12 2001
Individually, both this sequence and A002605 are convergents to 1+sqrt(3). Mutually, both sequences are convergents to 2+sqrt(3) and 1+sqrt(3)/2. - Klaus E. Kastberg (kastberg(AT)hotkey.net.au), Nov 04 2001 [Can someone clarify what is meant by the obscure second phrase, "Mutually..."? - M. F. Hasler, Aug 06 2018]
Add a loop at two vertices of the graph C_3=K_3. a(n) counts walks of length n+1 between these vertices. - Paul Barry, Oct 15 2004
Prefaced with a 1 as (1 + x + 3x^2 + 8x^3 + 22x^4 + ...) = 1 / (1 - x - 2x^2 - 3x^3 - 5x^4 - 8x^5 - 13x^6 - 21x^7 - ...). - Gary W. Adamson, Jul 28 2009
Equals row 2 of the array in A180165, and the INVERTi transform of A125145. - Gary W. Adamson, Aug 14 2010
Pisano period lengths: 1, 1, 3, 1, 24, 3, 48, 1, 9, 24, 10, 3, 12, 48, 24, 1, 144, 9, 180, 24, .... - R. J. Mathar, Aug 10 2012
Also the number of independent vertex sets and vertex covers in the n-centipede graph. - Eric W. Weisstein, Sep 21 2017
From Gus Wiseman, May 19 2020: (Start)
Conjecture: Also the number of length n + 1 sequences that cover an initial interval of positive integers and whose non-adjacent parts are weakly decreasing. For example, (3,2,3,1,2) has non-adjacent pairs (3,3), (3,1), (3,2), (2,1), (2,2), (3,2), all of which are weakly decreasing, so is counted under a(11). The a(1) = 1 through a(3) = 8 sequences are:
(1) (11) (111)
(12) (121)
(21) (211)
(212)
(221)
(231)
(312)
(321)
The case of compositions is A333148, or A333150 for strict compositions, or A333193 for strictly decreasing parts. A version for ordered set partitions is A332872. Standard composition numbers of these compositions are A334966. Unimodal normal sequences are A227038. See also: A001045, A001523, A032020, A100471, A100881, A115981, A329398, A332836, A332872.
(End)
Number of 2-compositions of n+1 restricted to parts 1 and 2 (and allowed zeros); see Hopkins & Ouvry reference. - Brian Hopkins, Aug 16 2020
The number of ternary strings of length n not containing 00. Complement of A186244. - R. J. Mathar, Feb 13 2022

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see p. 73).

Crossrefs

Cf. A155020 (same sequence with term 1 prepended).
Cf. A002605.

Programs

  • Haskell
    a028859 n = a028859_list !! n
    a028859_list =
       1 : 3 : map (* 2) (zipWith (+) a028859_list (tail a028859_list))
    -- Reinhard Zumkeller, Oct 15 2011
    
  • Maple
    a[0]:=1:a[1]:=3:for n from 2 to 24 do a[n]:=2*a[n-1]+2*a[n-2] od: seq(a[n],n=0..24); # Emeric Deutsch
  • Mathematica
    a[n_]:=(MatrixPower[{{1,3},{1,1}},n].{{2},{1}})[[2,1]]; Table[a[n],{n,0,40}] (* Vladimir Joseph Stephan Orlovsky, Feb 20 2010 *)
    Table[2^(n - 1) Hypergeometric2F1[(1 - n)/2, -n/2, -n, -2], {n, 20}] (* Eric W. Weisstein, Jun 14 2017 *)
    LinearRecurrence[{2, 2}, {1, 3}, 20] (* Eric W. Weisstein, Jun 14 2017 *)
  • PARI
    a(n)=([1,3;1,1]^n*[2;1])[2,1] \\ Charles R Greathouse IV, Mar 27 2012
    
  • PARI
    A028859(n)=([1,1]*[2,2;1,0]^n)[1] \\ M. F. Hasler, Aug 06 2018

Formula

a(n) = a(n-1) + A052945(n) = A002605(n) + A002605(n-1).
G.f.: -(x+1)/(2*x^2+2*x-1).
a(n) = [(1+sqrt(3))^(n+2)-(1-sqrt(3))^(n+2)]/(4*sqrt(3)). - Emeric Deutsch, Feb 01 2005
If p[i]=fibonacci(i+1) and if A is the 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, May 08 2010
a(n) = 3^n - A186244(n). - Toby Gottfried, Mar 07 2013
E.g.f.: exp(x)*(cosh(sqrt(3)*x) + 2*sinh(sqrt(3)*x)/sqrt(3)). - Stefano Spezia, Mar 02 2024

Extensions

Definition completed by M. F. Hasler, Aug 06 2018

A064613 Second binomial transform of the Catalan numbers.

Original entry on oeis.org

1, 3, 10, 37, 150, 654, 3012, 14445, 71398, 361114, 1859628, 9716194, 51373180, 274352316, 1477635912, 8016865533, 43773564294, 240356635170, 1326359740956, 7351846397334, 40913414754324, 228508350629892
Offset: 0

Author

Karol A. Penson, Sep 24 2001

Keywords

Comments

Exponential convolution of Catalan numbers and powers of 2. - Vladeta Jovovic, Dec 03 2004
Hankel transform of this sequence gives A000012 = [1,1,1,1,1,...]. - Philippe Deléham, Oct 24 2007
a(n) is the number of Motzkin paths of length n in which the (1,0)-steps at level 0 come in 3 colors and those at a higher level come in 4 colors. Example: a(3)=37 because, denoting U=(1,1), H=(1,0), and D=(1,-1), we have 3^3 = 27 paths of shape HHH, 3 paths of shape HUD, 3 paths of shape UDH, and 4 paths of shape UHD. - Emeric Deutsch, May 02 2011
a(n) is the number of Schroeder paths of semilength n in which the (2,0)-steps come in 2 colors and having no (2,0)-steps at levels 1,3,5,... - José Luis Ramírez Ramírez, Mar 30 2013
From Tom Copeland, Nov 08 2014: (Start)
This array is one of a family of Catalan arrays related by compositions of the special fractional linear (Möbius) transformations P(x,t)=x/(1-t*x); its inverse Pinv(x,t) = P(x,-t); and an o.g.f. of the Catalan numbers A000108, C(x) = [1-sqrt(1-4x)]/2; and its inverse Cinv(x) = x*(1-x). (Cf A126930.)
O.g.f.: G(x) = C[P[P(x,-1),-1]] = C[P(x,-2)] = (1-sqrt(1-4*x/(1-2*x)))/2 = x*A064613(x).
Ginv(x) = Pinv[Cinv(x),-2] = P[Cinv(x),2] = x(1-x)/[1+2x(1-x)] = (x-x^2)/[1+2(x-x^2)] = x - 3 x^2 + 8 x^3 - ... is -A155020(-x) ignoring first term there. (Cf. A146559, A125145.)(End)

Crossrefs

Programs

  • Magma
    I:=[3,10]; [1] cat [n le 2 select I[n] else ((8*n-2)*Self(n-1)-(12*n-12)*Self(n-2))div (n+1): n in [1..30]]; // Vincenzo Librandi, Jan 23 2017
  • Mathematica
    CoefficientList[Series[(1-Sqrt[(1-6*x)/(1-2*x)])/2/x, {x, 0, 20}], x] (* Vaclav Kotesovec, Jun 29 2013 *)
    a[n_] := 2^n Hypergeometric2F1[1/2, -n, 2, -2];
    Array[a, 22, 0] (* Peter Luschny, Jan 27 2020 *)
  • PARI
    x='x+O('x^66); Vec((1-sqrt((1-6*x)/(1-2*x)))/(2*x)) /* Joerg Arndt, Mar 31 2013 */
    

Formula

a(n) = Sum_{k=0..n} binomial(n, k)*binomial(2*k, k)*2^(n-k)/(k+1).
a(n) = 2^n*hypergeom([1/2, -n], [2], -2).
G.f.: (1-sqrt((1-6*x)/(1-2*x)))/(2*x). - Vladeta Jovovic, May 03 2003
With offset 1: a(1) = 1, a(n) = 2^(n-1) + Sum_{i=1..n-1} a(i)*a(n-i). - Benoit Cloitre, Mar 16 2004
D-finite with recurrence (n+1)*a(n) = (8*n-2)*a(n-1) - (12*n-12)*a(n-2). - Vladeta Jovovic, Jul 16 2004
E.g.f.: exp(4*x)*(BesselI(0, 2*x) - BesselI(1, 2*x)). - Vladeta Jovovic, Dec 03 2004
Inverse binomial transform of A104455. - Philippe Deléham, Nov 30 2007
G.f.: 1/(1-3*x-x^2/(1-4*x-x^2/(1-4*x-x^2/(1-4*x-x^2/(1-... (continued fraction). - Paul Barry, Jul 02 2009
a(n) = Sum_{0<=k<=n} A052179(n,k)*(-1)^k. - Philippe Deléham, Nov 28 2009
From Gary W. Adamson, Jul 21 2011: (Start)
a(n) = the upper left term in M^n, M = an infinite square production matrix as follows:
3, 1, 0, 0, ...
1, 3, 1, 0, ...
1, 1, 3, 1, ...
1, 1, 1, 3, ...
... (End)
a(n) ~ 2^(n-3/2)*3^(n+3/2)/(n^(3/2)*sqrt(Pi)). - Vaclav Kotesovec, Jun 29 2013
G.f. A(x) satisfies: A(x) = 1/(1 - 2*x) + x * A(x)^2. - Ilya Gutkovskiy, Jun 30 2020

Extensions

Name clarified using a comment of Vladeta Jovovic by Peter Bala, Jan 27 2020

A057960 Number of base-5 (n+1)-digit numbers starting with a zero and with adjacent digits differing by one or less.

Original entry on oeis.org

1, 2, 5, 13, 35, 95, 259, 707, 1931, 5275, 14411, 39371, 107563, 293867, 802859, 2193451, 5992619, 16372139, 44729515, 122203307, 333865643, 912137899, 2492007083, 6808289963, 18600594091, 50817768107, 138836724395, 379308985003, 1036291418795, 2831200807595
Offset: 0

Author

Henry Bottomley, May 18 2001

Keywords

Comments

Or, number of three-choice paths along a corridor of width 5 and length n, starting from one side.
If b(n) is the number of three-choice paths along a corridor of width 5 and length n, starting from any of the five positions at the beginning of the corridor, then b(n) = a(n+2) for n >= 0. - Pontus von Brömssen, Sep 06 2021

Examples

			a(6) = 259 since a(5) = 21 + 30 + 25 + 14 + 5 so a(6) = (21+30) + (21 + 30 + 25) + (30+25+14) + (25+14+5) + (14+5) = 51 + 76 + 69 + 44 + 19.
		

Crossrefs

The "three-choice" comes in the recurrence b(n+1, i) = b(n, i-1) + b(n, i) + b(n, i+1) if 1 <= i <= 5. Narrower corridors produce A000012, A000079, A000129, A001519. An infinitely wide corridor (i.e., just one wall) would produce A005773. Two-choice corridors are A000124, A000125, A000127.
Cf. A038754, A052948, A155020 (first differences), A188866.

Programs

  • Maple
    with(combstruct): ZL0:=S=Prod(Sequence(Prod(a, Sequence(b))), b): ZL1:=Prod(begin_blockP, Z, end_blockP): ZL2:=Prod(begin_blockLR, Z, Sequence(Prod(mu_length, Z), card>=1), end_blockLR): ZL3:=Prod(begin_blockRL, Sequence(Prod(mu_length, Z), card>=1), Z, end_blockRL):Q:=subs([a=Union(ZL1, ZL2, ZL3), b=ZL3], ZL0), begin_blockP=Epsilon, end_blockP=Epsilon, begin_blockLR=Epsilon, end_blockLR=Epsilon, begin_blockRL=Epsilon, end_blockRL=Epsilon, mu_length=Epsilon:temp15:=draw([S, {Q}, unlabelled], size=15):seq(count([S, {Q}, unlabelled], size=n+2), n=0..28); # Zerinvary Lajos, Mar 08 2008
  • Mathematica
    Join[{a=1,b=2},Table[c=(a+b)*2-1;a=b;b=c,{n,0,50}]] (* Vladimir Joseph Stephan Orlovsky, Nov 22 2010 *)
    CoefficientList[Series[(1-x-x^2)/((1-x)*(1-2*x-2*x^2)),{x,0,100}],x] (* Vincenzo Librandi, Aug 13 2012 *)
  • Python
    from functools import cache
    @cache
    def B(n, j):
        if not 0 <= j < 5:
            return 0
        if n == 0:
            return j == 0
        return B(n - 1, j - 1) + B(n - 1, j) + B(n - 1, j + 1)
    def A057960(n):
        return sum(B(n, j) for j in range(5))
    print([A057960(n) for n in range(30)]) # Pontus von Brömssen, Sep 06 2021

Formula

a(n) = Sum_{0 <= i <= 6} b(n, i) where b(n, 0) = b(n, 6) = 0, b(0, 1) = 1, b(0, n) = 0 if n <> 1 and b(n+1, i) = b(n, i-1) + b(n, i) + b(n, i+1) if 1 <= i <= 5.
a(n) = 3*a(n-1) - 2*a(n-3) = 2*A052948(n) - A052948(n-2).
a(n) = ceiling((1+sqrt(3))^(n+2)/12). - Mitch Harris, Apr 26 2006
a(n) = floor(a(n-1)*(a(n-1) + 1/2)/a(n-2)). - Franklin T. Adams-Watters and Max Alekseyev, Apr 25 2006
a(n) = floor(a(n-1)*(1+sqrt(3))). - Philippe Deléham, Jul 25 2003
From Paul Barry, Sep 16 2003: (Start)
G.f.: (1-x-x^2)/((1-x)*(1-2*x-2*x^2));
a(n) = 1/3 + (2+sqrt(3))*(1+sqrt(3))^n/6 + (2-sqrt(3))*(1-sqrt(3))^n/6.
Binomial transform of A038754 (with extra leading 1). (End)
More generally, it appears that a(base,n) = a(base-1,n) + 3^(n-1) for base >= n; a(base,n) = a(base-1,n) + 3^(n-1)-2 when base = n-1. - R. H. Hardin, Dec 26 2006
a(n) = A188866(4,n-1) for n >= 2. - Pontus von Brömssen, Sep 06 2021
a(n) = 2*a(n-1) + 2*a(n-2) - 1 for n >= 2, a(0) = 1, a(1) = 2. - Philippe Deléham, Mar 01 2024
E.g.f.: exp(x)*(1 + 2*cosh(sqrt(3)*x) + sqrt(3)*sinh(sqrt(3)*x))/3. - Stefano Spezia, Mar 02 2024

Extensions

This is the result of merging two identical entries submitted by Henry Bottomley and R. H. Hardin. - N. J. A. Sloane, Aug 14 2012
Name clarified by Pontus von Brömssen, Sep 06 2021

A155116 a(n) = 3*a(n-1) + 3*a(n-2), n>2, a(0)=1, a(1)=2, a(2)=8.

Original entry on oeis.org

1, 2, 8, 30, 114, 432, 1638, 6210, 23544, 89262, 338418, 1283040, 4864374, 18442242, 69919848, 265086270, 1005018354, 3810313872, 14445996678, 54768931650, 207644784984, 787241149902, 2984657804658, 11315696863680, 42901064005014
Offset: 0

Author

Philippe Deléham, Jan 20 2009

Keywords

Comments

From Johannes W. Meijer, Aug 14 2010: (Start)
A berserker sequence, see A180140 and A180147. For the central square 16 A[5] vectors with decimal values between 3 and 384 lead to this sequence. These vectors lead for the corner squares to A123620 and for the side squares to A180142.
This sequence belongs to a family of sequences with GF(x)=(1-(2*k-1)*x-k*x^2)/(1-3*x+(k-4)*x^2). Berserker sequences that are members of this family are A000007 (k=2), A155116 (k=1; this sequence), A000302 (k=0), 6*A179606 (k=-1; with leading 1 added) and 2*A180141 (k=-2; n>=1 and a(0)=1). Some other members of this family are (-2)*A003688 (k=3; with leading 1 added), (-4)*A003946 (k=4; with leading 1 added), (-6)*A002878 (k=5; with leading 1 added) and (-8)*A033484 (k=6; with leading 1 added).
Inverse binomial transform of A101368 (without the first leading 1).
(End)

Crossrefs

Sequences of the form a(n) = m*(a(n-1) + a(n-2)) with a(0)=1, a(1) = m-1, a(2) = m^2 -1: A155020 (m=2), this sequence (m=3), A155117 (m=4), A155119 (m=5), A155127 (m=6), A155130 (m=7), A155132 (m=8), A155144 (m=9), A155157 (m=10).

Programs

  • Magma
    m:=3; [1] cat [n le 2 select (m-1)*(m*n-(m-1)) else m*(Self(n-1) + Self(n-2)): n in [1..30]]; // G. C. Greubel, Mar 25 2021
    
  • Mathematica
    With[{m=3}, LinearRecurrence[{m, m}, {1, m-1, m^2-1}, 30]] (* G. C. Greubel, Mar 25 2021 *)
  • PARI
    Vec((1-x-x^2)/(1-3*x-3*x^2)+O(x^99)) \\ Charles R Greathouse IV, Jan 12 2012
    
  • Sage
    m=3; [1]+[-(m-1)*(sqrt(m)*i)^(n-2)*chebyshev_U(n, -sqrt(m)*i/2) for n in (1..30)] # G. C. Greubel, Mar 25 2021

Formula

G.f.: (1-x-x^2)/(1-3*x-3*x^2).
a(n) = 2*A125145(n-1), n>=1 .
a(n) = ( (2+4*A)*A^(-n-1) + (2+4*B)*B^(-n-1) )/21 with A=(-3+sqrt(21))/6 and B=(-3-sqrt(21))/6 for n>=1 with a(0)=1. [corrected by Johannes W. Meijer, Aug 12 2010]
Contribution from Johannes W. Meijer, Aug 14 2010: (Start)
a(n) = A123620(n)/2 for n>=1.
(End)
a(n) = (1/3)*[n=0] - 2*(sqrt(3)*i)^(n-2)*ChebyshevU(n, -sqrt(3)*i/2). - G. C. Greubel, Mar 25 2021

A155117 a(n) = 4*a(n-1) + 4*a(n-2), n>2, a(0)=1, a(1)=3, a(2)=15.

Original entry on oeis.org

1, 3, 15, 72, 348, 1680, 8112, 39168, 189120, 913152, 4409088, 21288960, 102792192, 496324608, 2396467200, 11571167232, 55870537728, 269766819840, 1302549430272, 6289265000448, 30367257722880, 146626090893312
Offset: 0

Author

Philippe Deléham, Jan 20 2009

Keywords

Crossrefs

Sequences of the form a(n) = m*(a(n-1) + a(n-2)) with a(0)=1, a(1) = m-1, a(2) = m^2 -1: A155020 (m=2), A155116 (m=3), this sequence (m=4), A155119 (m=5), A155127 (m=6), A155130 (m=7), A155132 (m=8), A155144 (m=9), A155157 (m=10).
Cf. A000129.

Programs

  • Magma
    m:=4; [1] cat [n le 2 select (m-1)*(m*n-(m-1)) else m*(Self(n-1) + Self(n-2)): n in [1..30]]; // G. C. Greubel, Mar 25 2021
    
  • Maple
    1,seq(simplify(-3*(2*I)^(n-2)*ChebyshevU(n, -I)), n = 1..30); # G. C. Greubel, Mar 25 2021
  • Mathematica
    With[{m=4}, LinearRecurrence[{m, m}, {1, m-1, m^2-1}, 30]] (* G. C. Greubel, Mar 25 2021 *)
  • Sage
    m=4; [1]+[-(m-1)*(sqrt(m)*i)^(n-2)*chebyshev_U(n, -sqrt(m)*i/2) for n in (1..30)] # G. C. Greubel, Mar 25 2021

Formula

G.f.: (1-x-x^2)/(1-4*x-4*x^2) . a(n)=3*A086347(n), n>=1 .
From G. C. Greubel, Mar 25 2021: (Start)
a(n) = (1/4)*[n=0] - 3*(2*i)^(n-2)*ChebyshevU(n, -i).
a(n) = (1/4)*[n=0] + 3*2^(n-2)*P_{n+1}, where P_{n} = A000129(n) (Pell numbers). (End)

A155119 a(n) = 5*a(n-1) + 5*a(n-2), n > 2, a(0)=1, a(1)=4, a(2)=24.

Original entry on oeis.org

1, 4, 24, 140, 820, 4800, 28100, 164500, 963000, 5637500, 33002500, 193200000, 1131012500, 6621062500, 38760375000, 226907187500, 1328337812500, 7776225000000, 45522814062500, 266495195312500, 1560090046875000, 9132926210937500
Offset: 0

Author

Philippe Deléham, Jan 20 2009

Keywords

Crossrefs

Sequences of the form a(n) = m*(a(n-1) + a(n-2)) with a(0)=1, a(1) = m-1, a(2) = m^2 -1: A155020 (m=2), A155116 (m=3), A155117 (m=4), this sequence (m=5), A155127 (m=6), A155130 (m=7), A155132 (m=8), A155144 (m=9), A155157 (m=10).

Programs

  • Magma
    m:=5; [1] cat [n le 2 select (m-1)*(m*n-(m-1)) else m*(Self(n-1) + Self(n-2)): n in [1..30]]; // G. C. Greubel, Mar 25 2021
    
  • Mathematica
    With[{m=5}, LinearRecurrence[{m, m}, {1, m-1, m^2-1}, 30]] (* G. C. Greubel, Mar 25 2021 *)
  • Sage
    m=5; [1]+[-(m-1)*(sqrt(m)*i)^(n-2)*chebyshev_U(n, -sqrt(m)*i/2) for n in (1..30)] # G. C. Greubel, Mar 25 2021

Formula

G.f.: (1 - x - x^2) / (1 - 5*x - 5*x^2).
a(n) = (1/5)*[n=0] - 4*(sqrt(5)*i)^(n-2)*ChebyshevU(n, -sqrt(5)*i/2). - G. C. Greubel, Mar 25 2021
E.g.f.: (3 + 4*exp(5*x/2)*(3*cosh(3*sqrt(5)*x/2) + sqrt(5)*sinh(3*sqrt(5)*x/2)))/15. - Stefano Spezia, May 31 2023

Extensions

a(20) corrected and a(21) from Sean A. Irvine, May 19 2019

A155127 a(n) = 6*a(n-1) + 6*a(n-2), n>2, a(0)=1, a(1)=5, a(2)=35.

Original entry on oeis.org

1, 5, 35, 240, 1650, 11340, 77940, 535680, 3681720, 25304400, 173916720, 1195326720, 8215460640, 56464724160, 388081108800, 2667274997760, 18332136639360, 125996469822720, 865971638772480, 5951808651571200
Offset: 0

Author

Philippe Deléham, Jan 20 2009

Keywords

Crossrefs

Sequences of the form a(n) = m*(a(n-1) + a(n-2)) with a(0)=1, a(1) = m-1, a(2) = m^2 -1: A155020 (m=2), A155116 (m=3), A155117 (m=4), A155119 (m=5), this sequence (m=6), A155130 (m=7), A155132 (m=8), A155144 (m=9), A155157 (m=10).

Programs

  • Magma
    m:=6; [1] cat [n le 2 select (m-1)*(m*n-(m-1)) else m*(Self(n-1) + Self(n-2)): n in [1..30]]; // G. C. Greubel, Mar 25 2021
    
  • Maple
    m:=6; 1,seq(simplify((1-m)*(sqrt(m)*I)^(n-2)*ChebyshevU(n, -I*sqrt(m)/2)), n = 1..30); # G. C. Greubel, Mar 25 2021
  • Mathematica
    LinearRecurrence[{6,6},{1,5,35},20] (* Harvey P. Dale, Apr 14 2015 *)
  • Sage
    m=6; [1]+[-(m-1)*(sqrt(m)*i)^(n-2)*chebyshev_U(n, -sqrt(m)*i/2) for n in (1..30)] # G. C. Greubel, Mar 25 2021

Formula

G.f.: (1-x-x^2)/(1-6*x-6*x^2) .
a(n) = (1/6)*[n=0] - 5*(sqrt(6)*i)^(n-2)*ChebyshevU(n, -sqrt(6)*i/2). - G. C. Greubel, Mar 25 2021

A155130 a(n) = 7*a(n-1) + 7*a(n-2), n>2, a(0)=1, a(1)=6, a(2)=48.

Original entry on oeis.org

1, 6, 48, 378, 2982, 23520, 185514, 1463238, 11541264, 91031514, 718009446, 5663286720, 44669073162, 352326519174, 2778969146352, 21919069658682, 172886271635238, 1363637389057440, 10755665624848746, 84835121097343302
Offset: 0

Author

Philippe Deléham, Jan 20 2009

Keywords

Crossrefs

Sequences of the form a(n) = m*(a(n-1) + a(n-2)) with a(0)=1, a(1) = m-1, a(2) = m^2 -1: A155020 (m=2), A155116 (m=3), A155117 (m=4), A155119 (m=5), A155127 (m=6), this sequence (m=7), A155132 (m=8), A155144 (m=9), A155157 (m=10).

Programs

  • Magma
    m:=7; [1] cat [n le 2 select (m-1)*(m*n-(m-1)) else m*(Self(n-1) + Self(n-2)): n in [1..30]]; // G. C. Greubel, Mar 25 2021
    
  • Maple
    m:= 7; 1,seq(simplify((1-m)*(sqrt(m)*I)^(n-2)*ChebyshevU(n, -I*sqrt(m)/2)), n = 1..30); # G. C. Greubel, Mar 25 2021
  • Mathematica
    LinearRecurrence[{7,7},{1,6,48},30] (* Harvey P. Dale, Mar 11 2018 *)
  • Sage
    m=7; [1]+[-(m-1)*(sqrt(m)*i)^(n-2)*chebyshev_U(n, -sqrt(m)*i/2) for n in (1..30)] # G. C. Greubel, Mar 25 2021

Formula

G.f.: (1-x-x^2)/(1-7*x-7*x^2) .
a(n) = (1/7)*[n=0] - 6*(sqrt(7)*i)^(n-2)*ChebyshevU(n, -sqrt(7)*i/2). - G. C. Greubel, Mar 25 2021

A155132 a(n) = 8*a(n-1) + 8*a(n-2), n > 2, a(0)=1, a(1)=7, a(2)=63.

Original entry on oeis.org

1, 7, 63, 560, 4984, 44352, 394688, 3512320, 31256064, 278147072, 2475225088, 22026977280, 196017618944, 1744356769792, 15522995109888, 138138815037440, 1229294481178624, 10939466369728512, 97350086807257088
Offset: 0

Author

Philippe Deléham, Jan 20 2009

Keywords

Crossrefs

Sequences of the form a(n) = m*(a(n-1) + a(n-2)) with a(0)=1, a(1) = m-1, a(2) = m^2 -1: A155020 (m=2), A155116 (m=3), A155117 (m=4), A155119 (m=5), A155127 (m=6), A155130 (m=7), this sequence (m=8), A155144 (m=9), A155157 (m=10).

Programs

  • Magma
    [1] cat [n le 2 select 7*(8*n-7) else 8*(Self(n-1) + Self(n-2)): n in [1..30]]; // G. C. Greubel, Mar 24 2021
    
  • Mathematica
    With[{m=8}, LinearRecurrence[{m, m}, {1, m-1, m^2-1}, 30]] (* G. C. Greubel, Mar 24 2021 *)
  • Sage
    [1]+[-7*(2*sqrt(2)*i)^(n-2)*chebyshev_U(n, -sqrt(2)*i) for n in (1..30)] # G. C. Greubel, Mar 24 2021

Formula

G.f.: (1-x-x^2)/(1-8*x-8*x^2) .
a(n) = (1/8)*[n=0] - 7*(2*sqrt(2)*i)^(n-2)*ChebyshevU(n, -sqrt(2)*I). - G. C. Greubel, Mar 24 2021
Showing 1-10 of 14 results. Next