cp's OEIS Frontend

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

Previous Showing 31-40 of 146 results. Next

A073133 Table by antidiagonals of T(n,k) = n*T(n,k-1) + T(n,k-2) starting with T(n,1) = 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 1, 3, 5, 3, 1, 4, 10, 12, 5, 1, 5, 17, 33, 29, 8, 1, 6, 26, 72, 109, 70, 13, 1, 7, 37, 135, 305, 360, 169, 21, 1, 8, 50, 228, 701, 1292, 1189, 408, 34, 1, 9, 65, 357, 1405, 3640, 5473, 3927, 985, 55, 1, 10, 82, 528, 2549, 8658, 18901, 23184, 12970, 2378, 89
Offset: 1

Views

Author

Henry Bottomley, Jul 16 2002

Keywords

Comments

Columns of the array are generated from Fibonacci polynomials f(x). They are: (1), (x), (x^2 + 1), (x^3 + 2x), (x^4 + 3x^2 + 1), (x^5 + 4x^3 + 3x), (x^6 + 5x^4 + 6x^2 +1), ... If column headings start 0, 1, 2, ... then the terms in the n-th column are generated from the n-th degree Fibonacci polynomial. For example, column 5 (8, 70, 360, ...) is generated from f(x), x = 1,2,3,...; fifth-degree polynomial x^5 + 4x^3 + 3x; e.g., f(2) = 70 = 2^5 + 4*8 + 3*2. - Gary W. Adamson, Apr 02 2006
The ratio of two consecutive entries of the sequence in the n-th row approaches (n + sqrt(n^2 + 4))/2. Example: The sequence beginning (1, 3, 10, 33, ...) tends to 3.302775... = (3 + sqrt(13))/2. - Gary W. Adamson, Aug 12 2013
As to the array sequences, (n+1)-th sequence is the INVERT transform of the n-th sequence. - Gary W. Adamson, Aug 20 2013
The array can be extended infinitely above the Fibonacci row by taking successive INVERTi transforms, resulting in:
...
1, -2, 5, -12, 29, -70, ...
1, -1, 2, -3, 5, -8, ...
l, 0, 1, 0, 1, 0, ...
1, 1, 2, 3, 5, 8, ...
1, 2, 5, 12, 29, 70, ...
...
This results in an infinite array in which sequences above the (1, 0, 1, 0, ...) are reflections of the sequences below, except for the alternate signs. Any sequence in the (+ sign) row starting (1, n, ...) is the (2*n-th) INVERT transform of the same sequence but with alternate signs. Example: (1, 2, 5, 12, ...) is the (2*2) = fourth INVERT transform of (1, -2, 5, -12, ...) by inspection. Conjecture: This "reflection" principle will result from taking successive INVERT transforms of any aerated sequence starting 1, ... and with positive signs. Likewise, the rows above the aerated sequence are successive INVERTi transforms of the aerated sequence. - Gary W. Adamson, Jul 14 2019
From Michael A. Allen, Feb 21 2023: (Start)
Row n is the n-metallonacci sequence.
T(n,k) is the number of tilings of a (k-1)-board (a board with dimensions (k-1) X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are n kinds of squares available. (End)

Examples

			Table begins:
  1, 1,  2,  3,   5,    8,   13, ...
  1, 2,  5, 12,  29,   70,  169, ...
  1, 3, 10, 33, 109,  360, 1189, ...
  1, 4, 17, 72, 305, 1292, 5473, ... etc.
		

Crossrefs

Programs

  • GAP
    T:= function(n,k)
        if k<0 then return 0;
        elif k=1 then return 1;
        else return n*T(n,k-1) + T(n,k-2);
        fi;
      end;
    Flat(List([1..15], n-> List([1..n], k-> T(n-k+1,k) ))); # G. C. Greubel, Aug 12 2019
  • Maple
    A073133 := proc(n,k)
        option remember;
        if k <= 1 then
            k;
        else
            n*procname(n,k-1)+procname(n,k-2) ;
        end if;
    end proc:
    seq(seq( A073133(d-k,k),k=1..d-1),d=2..13) ; # R. J. Mathar, Aug 16 2019
  • Mathematica
    T[n_, 1]:= 1; T[n_, k_]:= T[n, k] = If[k<0, 0, n*T[n, k-1] + T[n, k-2]]; Table[T[n-k+1, k], {n, 15}, {k, n}]//Flatten (* G. C. Greubel, Aug 12 2019 *)
  • PARI
    T(n,k) = if(k==1, 1, k<0, 0, n*T(n,k-1)+T(n,k-2));
    for(n=1,15, for(k=1,n, print1(T(n-k+1,k), ", "))) \\ G. C. Greubel, Aug 12 2019
    
  • Sage
    def T(n, k):
        if (k<0): return 0
        elif (k==1): return 1
        else: return n*T(n, k-1) + T(n, k-2)
    [[T(n-k+1, k) for k in (1..n)] for n in (1..15)] # G. C. Greubel, Aug 12 2019
    

Formula

T(n, k) = A073134(n, k) + 2*A073135(n, k-2) = Sum_{j=0..k-1} abs(A049310(k-1, j)*n^j).
T(n,k) = [[0,1; 1,n]^{k+1}]{1,1}, n,k in {1,2,...}. - _L. Edson Jeffery, Sep 23 2012
G.f. for row n: x/(1-n*x-x^2). - L. Edson Jeffery, Aug 28 2013

A037027 Skew Fibonacci-Pascal triangle read by rows.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 3, 5, 3, 1, 5, 10, 9, 4, 1, 8, 20, 22, 14, 5, 1, 13, 38, 51, 40, 20, 6, 1, 21, 71, 111, 105, 65, 27, 7, 1, 34, 130, 233, 256, 190, 98, 35, 8, 1, 55, 235, 474, 594, 511, 315, 140, 44, 9, 1, 89, 420, 942, 1324, 1295, 924, 490, 192, 54, 10, 1, 144, 744, 1836
Offset: 0

Views

Author

Floor van Lamoen, Jan 01 1999

Keywords

Comments

T(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (0,1), (1,0), (2,0). - Joerg Arndt, Jun 30 2011
T(n,k) is the number of lattice paths of length n, starting from the origin and ending at (n,k), using horizontal steps H=(1,0), up steps U=(1,1) and down steps D=(1,-1), never containing UUU, DD, HD. For instance, for n=4 and k=2, we have the paths; HHUU, HUHU, HUUH, UHHU, UHUH, UUHH, UUDU, UDUU, UUUD. - Emanuele Munarini, Mar 15 2011
Row sums form Pell numbers A000129, T(n,0) forms Fibonacci numbers A000045, T(n,1) forms A001629. T(n+k,n-k) is polynomial sequence of degree k.
T(n,k) gives a convolved Fibonacci sequence (A001629, A001872, etc.).
As a Riordan array, this is (1/(1-x-x^2),x/(1-x-x^2)). An interesting factorization is (1/(1-x^2),x/(1-x^2))*(1/(1-x),x/(1-x)) [abs(A049310) times A007318]. Diagonal sums are the Jacobsthal numbers A001045(n+1). - Paul Barry, Jul 28 2005
T(n,k) = T'(n+1,k+1), T' given by [0, 1, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 19 2005
Equals A049310 * A007318 as infinite lower triangular matrices. - Gary W. Adamson, Oct 28 2007
This triangle may also be obtained from the coefficients of the Morgan-Voyce polynomials defined by: Mv(x, n) = (x + 1)*Mv(x, n - 1) + Mv(x, n - 2). - Roger L. Bagula, Apr 09 2008
Row sums are A000129. - Roger L. Bagula, Apr 09 2008
Absolute value of coefficients of the characteristic polynomial of tridiagonal matrices with 1's along the main diagonal, and i's along the superdiagonal and the subdiagonal (where i=sqrt(-1), see Mathematica program). - John M. Campbell, Aug 23 2011
A037027 is jointly generated with A122075 as an array of coefficients of polynomials v(n,x): initially, u(1,x)=v(1,x)=1; for n>1, u(n,x)=u(n-1,x)+(x+1)*v(n-1)x and v(n,x)=u(n-1,x)+x*v(n-1,x). See the Mathematica section at A122075. - Clark Kimberling, Mar 05 2012
For a closed-form formula for arbitrary left and right borders of Pascal like triangle see A228196. - Boris Putievskiy, Aug 18 2013
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 09 2013
Row n, for n>=0, shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(n)*x^n which is the denominator of the n-th convergent of the continued fraction [x+1, x+1, x+1, ...]; see A230000. - Clark Kimberling, Nov 13 2013
T(n,k) is the number of ternary words of length n having k letters 2 and avoiding a runs of odd length for the letter 0. - Milan Janjic, Jan 14 2017
Let T(m, n, k) be an m-bonacci Pascal's triangle, where T(m, n, 0) gives the values of F(m, n), the n-th m-bonacci number, and T(m, n, k) gives the values for the k-th convolution of F(m, n). Then the classic Pascal triangle is T(1, n, k) and this sequence is T(2, n, k). T(m, n, k) is the number of compositions of n using only the positive integers 1, 1' and 2 through m, with the part 1' used exactly k times. G.f. for k-th column of T(m, n, k): x/(1 - x - x^2 - ... - x^m)^k. The row sum for T(m, n, k) is the number of compositions of n using only the positive integers 1, 1' and 2 through m. G.f. for row sum of T(m, n, k): 1/(1 - 2x - x^2 - ... - x^m). - Gregory L. Simay, Jul 24 2021

Examples

			Ratio of row polynomials R(3)/R(2) = (3 + 5*x + 3*x^2 + x^3)/(2 + 2*x + x^2) = [1+x; 1+x, 1+x].
Triangle begins:
                                 1;
                              1,    1;
                           2,    2,    1;
                        3,    5,    3,    1;
                     5,   10,    9,    4,    1;
                  8,   20,   22,   14,    5,    1;
              13,   38,   51,   40,   20,    6,    1;
           21,   71,  111,  105,   65,   27,    7,    1;
        34,  130,  233,  256,  190,   98,   35,    8,    1;
     55,  235,  474,  594,  511,  315,  140,   44,    9,    1;
  89,  420,  942, 1324, 1295,  924,  490,  192,   54,   10,    1;
		

Crossrefs

A038112(n) = T(2n, n). A038137 is reflected version. Maximal row entries: A038149.
Diagonal differences are in A055830. Vertical sums are in A091186.
Some other Fibonacci-Pascal triangles: A027926, A036355, A074829, A105809, A109906, A111006, A114197, A162741, A228074.

Programs

  • Haskell
    a037027 n k = a037027_tabl !! n !! k
    a037027_row n = a037027_tabl !! n
    a037027_tabl = [1] : [1,1] : f [1] [1,1] where
       f xs ys = ys' : f ys ys' where
         ys' = zipWith3 (\u v w -> u + v + w) (ys ++ [0]) (xs ++ [0,0]) ([0] ++ ys)
    -- Reinhard Zumkeller, Jul 07 2012
  • Maple
    T := (n,k) -> `if`(n=0,1,binomial(n,k)*hypergeom([(k-n)/2, (k-n+1)/2], [-n], -4)): seq(seq(simplify(T(n,k)), k=0..n), n=0..10); # Peter Luschny, Apr 25 2016
    # Uses function PMatrix from A357368. Adds a row above and a column to the left.
    PMatrix(10, n -> combinat:-fibonacci(n)); # Peter Luschny, Oct 07 2022
  • Mathematica
    Mv[x, -1] = 0; Mv[x, 0] = 1; Mv[x, 1] = 1 + x; Mv[x_, n_] := Mv[x, n] = ExpandAll[(x + 1)*Mv[x, n - 1] + Mv[x, n - 2]]; Table[ CoefficientList[ Mv[x, n], x], {n, 0, 10}] // Flatten (* Roger L. Bagula, Apr 09 2008 *)
    Abs[Flatten[Table[CoefficientList[CharacteristicPolynomial[Array[KroneckerDelta[#1,#2]+KroneckerDelta[#1,#2+1]*I+KroneckerDelta[#1,#2-1]*I&,{n,n}],x],x],{n,1,20}]]] (* John M. Campbell, Aug 23 2011 *)
    T[n_, k_] := Binomial[n, k] Hypergeometric2F1[(k-n)/2, (k-n+1)/2, -n, -4];
    Table[T[n, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 16 2019, after Peter Luschny *)
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, if( n==0 && k==0, 1, T(n-1, k) + T(n-1, k-1) + T(n-2, k)))}; /* Michael Somos, Sep 29 2003 */
    
  • PARI
    T(n,k)=if(nPaul D. Hanna, Feb 27 2004
    

Formula

T(n, m) = T'(n-1, m) + T'(n-2, m) + T'(n-1, m-1), where T'(n, m) = T(n, m) for n >= 0 and 0< = m <= n and T'(n, m) = 0 otherwise.
G.f.: 1/(1 - y - y*z - y^2).
G.f. for k-th column: x/(1-x-x^2)^k.
T(n, m) = Sum_{k=0..n-m} binomial(m+k, m)*binomial(k, n-k-m), n >= m >= 0, otherwise 0. - Wolfdieter Lang, Jun 17 2002
T(n, m) = ((n-m+1)*T(n, m-1) + 2*(n+m)*T(n-1, m-1))/(5*m), n >= m >= 1; T(n, 0)= A000045(n+1); T(n, m)= 0 if n < m. - Wolfdieter Lang, Apr 12 2000
Chebyshev coefficient triangle (abs(A049310)) times Pascal's triangle (A007318) as product of lower triangular matrices. T(n, k) = Sum_{j=0..n} binomial((n+j)/2, j)*(1+(-1)^(n+j))*binomial(j, k)/2. - Paul Barry, Dec 22 2004
Let R(n) = n-th row polynomial in x, with R(0)=1, then R(n+1)/R(n) equals the continued fraction [1+x;1+x, ...(1+x) occurring (n+1) times ..., 1+x] for n >= 0. - Paul D. Hanna, Feb 27 2004
T(n,k) = Sum_{j=0..n} binomial(n-j,j)*binomial(n-2*j,k); in Egorychev notation, T(n,k) = res_w(1-w-w^2)^(-k-1)*w^(-n+k+1). - Paul Barry, Sep 13 2006
Sum_{k=0..n} T(n,k)*x^k = A000045(n+1), A000129(n+1), A006190(n+1), A001076(n+1), A052918(n), A005668(n+1), A054413(n), A041025(n), A099371(n+1), A041041(n), A049666(n+1), A041061(n), A140455(n+1), A041085(n), A154597(n+1), A041113(n) for x = 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 respectively. - Philippe Deléham, Nov 29 2009
T((m+1)*n+r-1, m*n+r-1)*r/(m*n+r) = Sum_{k=1..n} k/n*T((m+1)*n-k-1, m*n-1)*(r+k,r), n >= m > 1.
T(n-1,m-1) = (m/n)*Sum_{k=1..n-m+1} k*A000045(k)*T(n-k-1,m-2), n >= m > 1. - Vladimir Kruchinin, Mar 17 2011
T(n,k) = binomial(n,k)*hypergeom([(k-n)/2, (k-n+1)/2], [-n], -4) for n >= 1. - Peter Luschny, Apr 25 2016

Extensions

Examples from Paul D. Hanna, Feb 27 2004

A172236 Array A(n,k) = n*A(n,k-1) + A(n,k-2) read by upward antidiagonals, starting A(n,0) = 0, A(n,1) = 1.

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 5, 3, 0, 1, 4, 10, 12, 5, 0, 1, 5, 17, 33, 29, 8, 0, 1, 6, 26, 72, 109, 70, 13, 0, 1, 7, 37, 135, 305, 360, 169, 21, 0, 1, 8, 50, 228, 701, 1292, 1189, 408, 34, 0, 1, 9, 65, 357, 1405, 3640, 5473, 3927, 985, 55, 0, 1, 10, 82, 528, 2549, 8658, 18901, 23184, 12970, 2378, 89, 0, 1, 11, 101, 747, 4289, 18200, 53353, 98145, 98209, 42837, 5741, 144
Offset: 1

Views

Author

Roger L. Bagula, Jan 29 2010

Keywords

Comments

Equals A073133 with an additional column A(.,0).
If the first column and top row are deleted, antidiagonal reading yields A118243.
Adding a top row of 1's and antidiagonal reading downwards yields A157103.
Antidiagonal sums are 0, 1, 2, 5, 12, 32, 93, 297, 1035, 3911, 15917, 69350, ....
From Jianing Song, Jul 14 2018: (Start)
All rows have strong divisibility, that is, gcd(A(n,k_1), A(n,k_2)) = A(n,gcd(k_1,k_2)) holds for all k_1, k_2 >= 0.
Let E(n,m) be the smallest number l such that m divides A(n,l), we have: for odd primes p that are not divisible by n^2 + 4, E(n,p) divides p - ((n^2+4)/p) if p == 3 (mod 4) and (p - ((n^2+4)/p))/2 if p == 1 (mod 4). E(n,p) = p for odd primes p that are divisible by n^2 + 4. E(n,2) = 2 for even n and 3 for odd n. Here ((n^2+4)/p) is the Legendre symbol. A prime p such that p^2 divides T(n,E(n,p)) is called an n-Wall-Sun-Sun prime.
E(n,p^e) <= p^(e-1)*E(n,p) for all primes p. If p^2 does not divide A(n, E(n,p)), then E(n,p^e) = p^(e-1)*E(n,p) for all exponent e. Specially, for primes p >= 5 that are divisible by n^2 + 4, p^2 is never divisible by A(n,p), so E(n,p^e) = p^e as described above. E(n,m_1*m_2) = lcm(E(n,m_1),E(n,m_2)) if gcd(m_1,m_2) = 1.
Let pi(n,m) be the Pisano period of A(n, k) modulo m, i.e, the smallest number l such that A(n, k+1) == T(n,k) (mod m) holds for all k >= 0, we have: for odd primes p that are not divisible by n^2 - 4, pi(n,p) divides p - 1 if ((n^2+4)/p) = 1 and 2(p+1) if ((n^2+4)/p) = -1. pi(n,p) = 4p for odd primes p that are divisible by n^2 + 4. pi(n,2) = 2 even n and 3 for odd n.
pi(n,p^e) <= p^(e-1)*pi(n,p) for all primes p. If p^2 does not divide A(n, E(n,p)), then pi(n,p^e) = p^(e-1)*pi(n,p) for all exponent e. Specially, for primes p >= 5 that are divisible by n^2 + 4, p^2 is never divisible by A(n, p), so pi(n,p^e) = 4p^e as described above. pi(n,m_1*m_2) = lcm(pi(n,m_1), pi(n,m_2)) if gcd(m_1,m_2) = 1.
If n != 2, the largest possible value of pi(n,m)/m is 4 for even n and 6 for odd n. For even n, pi(n,p^e) = 4p^e; for odd n, pi(n,2p^e) = 12p^e, where p is any odd prime factor of n^2 + 4. For n = 2 it is 8/3, obtained by m = 3^e.
Let z(n,m) be the number of zeros in a period of A(n, k) modulo m, i.e., z(n,m) = pi(n,m)/E(n,m), then we have: z(n,p) = 4 for odd primes p that are divisible by n^2 + 4. For other odd primes p, z(n,p) = 4 if E(n,p) is odd; 1 if E(n,p) is even but not divisible by 4; 2 if E(n,p) is divisible by 4; see the table below. z(n,2) = z(n,4) = 1.
Among all values of z(n,p) when p runs through all odd primes that are not divisible by n^2 + 4, we have:
((n^2+4)/p)...p mod 8....proportion of 1.....proportion of 2.....proportion of 4
......1..........1......1/6 (conjectured)...2/3 (conjectured)...1/6 (conjectured)*
......1..........5......1/2 (conjectured)...........0...........1/2 (conjectured)*
......1.........3,7.............1...................0...................0
.....-1.........1,5.............0...................0...................1
.....-1.........3,7.............0...................1...................0
* The result is that among all odd primes that are not divisible by n^2 + 4, 7/24 of them are with z(n,p) = 1, 5/12 are with z(n,p) = 2 and 7/24 are with z(n,p) = 4 if n^2 + 4 is a twice a square; 1/3 of them are with z(n,p) = 1, 1/3 are with z(n,p) = 2 and 1/3 are with z(n,p) = 4 otherwise. [Corrected by Jianing Song, Jul 06 2019]
z(n,p^e) = z(n,p) for all odd primes p; z(n,2^e) = 1 for even n and 2 for odd n, e >= 3.
(End)
From Michael A. Allen, Mar 06 2023: (Start)
Removing the first (n=0) row of A352361 gives this sequence.
Row n is the n-metallonacci sequence.
A(n,k) is (for k>0) the number of tilings of a (k-1)-board (a board with dimensions (k-1) X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are n kinds of squares available. (End)

Examples

			The array, A(n, k), starts in row n = 1 with columns k >= 0 as
  0      1      1      2      3      5      8
  0      1      2      5     12     29     70
  0      1      3     10     33    109    360
  0      1      4     17     72    305   1292
  0      1      5     26    135    701   3640
  0      1      6     37    228   1405   8658
  0      1      7     50    357   2549  18200
  0      1      8     65    528   4289  34840
  0      1      9     82    747   6805  61992
  0      1     10    101   1020  10301 104030
  0      1     11    122   1353  15005 166408
Antidiagonal triangle, T(n, k), begins as:
  0;
  0, 1;
  0, 1, 1;
  0, 1, 2,  2;
  0, 1, 3,  5,   3;
  0, 1, 4, 10,  12,   5;
  0, 1, 5, 17,  33,  29,    8;
  0, 1, 6, 26,  72, 109,   70,   13;
  0, 1, 7, 37, 135, 305,  360,  169,  21;
  0, 1, 8, 50, 228, 701, 1292, 1189, 408, 34;
		

Crossrefs

Rows n include: A000045 (n=1), A000129 (n=2), A006190 (n=3), A001076 (n=4), A052918 (n=5), A005668 (n=6), A054413 (n=7), A041025 (n=8), A099371 (n=9), A041041 (n=10), A049666 (n=11), A041061 (n=12), A140455 (n=13), A041085 (n=14), A154597 (n=15), A041113 (n=16), A178765 (n=17), A041145 (n=18), A243399 (n=19), A041181 (n=20). (Note that there are offset shifts for rows n = 5, 7, 8, 10, 12, 14, 16..20.)
Columns k include: A000004 (k=0), A000012 (k=1), A000027 (k=2), A002522 (k=3), A054602 (k=4), A057721 (k=5), A124152 (k=6).
Entry points for A(n,k) modulo m: A001177 (n=1), A214028 (n=2), A322907 (n=3).
Pisano period for A(n,k) modulo m: A001175 (n=1), A175181 (n=2), A175182 (n=3), A175183 (n=4), A175184 (n=5), A175185 (n=6).
Number of zeros in a period for A(n,k) modulo m: A001176 (n=1), A214027 (n=2), A322906 (n=3).
Sums include: A304357, A304359.
Similar to: A073133.

Programs

  • Magma
    A172236:= func< n,k | k le 1 select k else Evaluate(DicksonSecond(k-1,-1), n-k) >;
    [A172236(n,k): k in [0..n-1], n in [1..13]]; // G. C. Greubel, Sep 29 2024
    
  • Mathematica
    A172236[n_,k_]:=Fibonacci[k, n-k];
    Table[A172236[n, k], {n,15}, {k,0,n-1}]//Flatten
  • PARI
    A(n, k) = if (k==0, 0, if (k==1, 1, n*A(n, k-1) + A(n, k-2)));
    tabl(nn) = for(n=1, nn, for (k=0, nn, print1(A(n, k), ", ")); print); \\ Jianing Song, Jul 14 2018 (program from Michel Marcus; see also A316269)
    
  • PARI
    A(n, k) = ([n, 1; 1, 0]^k)[2, 1] \\ Jianing Song, Nov 10 2018
    
  • SageMath
    def A172236(n,k): return sum(binomial(k-j-1,j)*(n-k)^(k-2*j-1) for j in range(1+(k-1)//2))
    flatten([[A172236(n,k) for k in range(n)] for n in range(1,14)]) # G. C. Greubel, Sep 29 2024

Formula

A(n,k) = (((n + sqrt(n^2 + 4))/2)^k - ((n-sqrt(n^2 + 4))/2)^k)/sqrt(n^2 + 4), n >= 1, k >= 0. - Jianing Song, Jun 27 2018
For n >= 1, Sum_{i=1..k} 1/A(n,2^i) = ((u^(2^k-1) + v^(2^k-1))/(u + v)) * (1/A(n,2^k)), where u = (n + sqrt(n^2 + 4))/2, v = (n - sqrt(n^2 + 4))/2 are the two roots of the polynomial x^2 - n*x - 1. As a result, Sum_{i>=1} 1/A(n,2^i) = (n^2 + 4 - n*sqrt(n^2 + 4))/(2*n). - Jianing Song, Apr 21 2019
From G. C. Greubel, Sep 29 2024: (Start)
A(n, k) = F_{k}(n) (Fibonacci polynomials F_{n}(x)) (array).
T(n, k) = F_{k}(n-k) (antidiagonal triangle).
Sum_{k=0..n-1} T(n, k) = A304357(n) - (1-(-1)^n)/2.
Sum_{k=0..n-1} (-1)^k*T(n, k) = (-1)*A304359(n) + (1-(-1)^n)/2.
T(2*n, n) = A084844(n).
T(2*n+1, n+1) = A084845(n). (End)

Extensions

More terms from Jianing Song, Jul 14 2018

A291219 p-INVERT of (0,1,0,1,0,1,...), where p(S) = 1 - S - S^3.

Original entry on oeis.org

1, 1, 3, 5, 11, 21, 42, 83, 163, 323, 635, 1255, 2473, 4880, 9625, 18985, 37451, 73869, 145715, 287421, 566954, 1118331, 2205947, 4351307, 8583091, 16930447, 33395857, 65874464, 129939569, 256310161, 505580371, 997274197, 1967156763, 3880282533, 7653987242
Offset: 0

Views

Author

Clark Kimberling, Aug 24 2017

Keywords

Comments

Suppose s = (c(0), c(1), c(2), ...) is a sequence and p(S) is a polynomial. Let S(x) = c(0)*x + c(1)*x^2 + c(2)*x^3 + ... and T(x) = (-p(0) + 1/p(S(x)))/x. The p-INVERT of s is the sequence t(s) of coefficients in the Maclaurin series for T(x). Taking p(S) = 1 - S gives the "INVERT" transform of s, so that p-INVERT is a generalization of the "INVERT" transform (e.g., A033453).
In the following guide to p-INVERT sequences using s = (1,0,1,0,1,...) = A000035, in some cases t(1,0,1,0,1,...) is a shifted version of the indicated sequence.
p(S) t(1,0,1,0,1,...)
1 - S A000045 (Fibonacci numbers)
1 - S^2 A147600
1 - S^3 A291217
1 - S^5 A291218
1 - S - S^2 A289846
1 - S - S^3 A291219
1 - S - S^4 A291220
1 - S^3- S^6 A291221
1 - S^2- S^3 A291222
1 - S^3- S^4 A291223
1 - 2S A052542
1 - 3S A006190
(1 - S)^2 A239342
(1 - S)^3 A276129
(1 - S)^4 A291224
(1 - S)^5 A291225
(1 - S)^6 A291226
1 - S - 2 S^2 A291227
1 - 2 S - 2 S^2 A291228
1 - 3 S - 2 S^2 A060801
(1 - S)(1 - 2 S) A291229
(1 - S)(1 - 2 S)(1 - 3 S) A291230
(1 - S)(1 - 2 S)(1 - 3 S)( 1 - 4 S) A291231
(1 - 2 S)^2 A291264
(1 - 3 S)^2 A291232
1 - S - S^2 - S^3 A291233
1 - S - S^2 - S^3 - S^4 A291234
1 - S - S^2 - S^3 - S^4 - S^5 A291235
(1 - S)(1 - 3 S) A291236
(1 - S)(1 - 2S)( 1 - 4S) A291237
(1 - S)^2 (1 - 2S) A291238
(1 - S^2) (1 - 2S) A291239
(1 - S^3)^2 A291240
1 - S - S^2 + S^3 A291241
1 - 2 S - S^2 + S^3 A291242
1 - 3 S + S^2 A291243
1 - 4 S + S^2 A291244
1 - 5 S + S^2 A291245
1 - 6 S + S^2 A291246
1 - S - S^2 - S^3 + S^4 A291247
1 - S - S^2 - S^3 - S^4 + S^5 A291248
1 - S - S^2 - S^3 + S^4 + S^5 A291249
1 - S - 2 S^2 + 2 S^3 A291250
1 - 3 S^2 + 2 S^3 A291251 (includes negative terms)
(1 - S^3)^3 A291252
(1 - S - S^2)^2 A291253
(1 - 2 S - S^2)^2 A291254
(1 - S - 2 S^2)^2 A291255

Crossrefs

Programs

  • Magma
    I:=[1,1,3,5,11,21]; [n le 6 select I[n] else Self(n-1)+3*Self(n-2)-Self(n-3)-3*Self(n-4)+Self(n-5)+Self(n-6): n in [1..45]]; // Vincenzo Librandi, Aug 25 2017
  • Mathematica
    z = 60; s = x/(1 - x^2); p = 1 - s - s^3;
    Drop[CoefficientList[Series[s, {x, 0, z}], x], 1]  (* A000035 *)
    Drop[CoefficientList[Series[1/p, {x, 0, z}], x], 1]  (* A291219 *)
    LinearRecurrence[{1, 3, -1, -3, 1, 1}, {1, 1, 3, 5, 11, 21}, 50] (* Vincenzo Librandi, Aug 25 2017 *)

Formula

G.f.: -(1 - x^2 + x^4)/(-1 + x + 3*x^2 - x^3 - 3*x^4 + x^5 + x^6).
a(n) = a(n-1) + 3*a(n-2) - a(n-3) - 3*a(n-4) + a(n-5) + a(n-6) for n >= 7.

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

Original entry on oeis.org

2, 3, 11, 36, 119, 393, 1298, 4287, 14159, 46764, 154451, 510117, 1684802, 5564523, 18378371, 60699636, 200477279, 662131473, 2186871698, 7222746567, 23855111399, 78788080764, 260219353691, 859446141837, 2838557779202
Offset: 0

Views

Author

Keywords

Comments

For more information about this type of recurrence follow the Khovanova link and see A086902 and A054413. - Johannes W. Meijer, Jun 12 2010

References

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

Crossrefs

Programs

  • Haskell
    a006497 n = a006497_list !! n
    a006497_list = 2 : 3 : zipWith (+) (map (* 3) $ tail a006497_list) a006497_list
    -- Reinhard Zumkeller, Feb 19 2011
    
  • Magma
    [ n eq 1 select 2 else n eq 2 select 3 else 3*Self(n-1)+Self(n-2): n in [1..30] ]; // Vincenzo Librandi, Aug 20 2011
    
  • Maple
    a:= n-> (<<0|1>, <1|3>>^n. <<2, 3>>)[1, 1]:
    seq(a(n), n=0..30);  # Alois P. Heinz, Jan 26 2018
  • Mathematica
    Table[LucasL[n, 3], {n, 0, 30}] (* Zerinvary Lajos, Jul 09 2009 *)
    LucasL[Range[0, 30], 3] (* Eric W. Weisstein, Apr 17 2018 *)
    LinearRecurrence[{3,1},{2,3},30] (* Harvey P. Dale, Feb 17 2020 *)
  • PARI
    my(x='x+O('x^30)); Vec((2-3*x)/(1-3*x-x^2)) \\ G. C. Greubel, Jul 05 2017
    
  • PARI
    apply( {A006497(n)=[2,3]*([0,1;1,3]^n)[,1]}, [0..30]) \\ M. F. Hasler, Mar 06 2020
  • Sage
    [lucas_number2(n,3,-1) for n in range(0, 30)] # Zerinvary Lajos, Apr 30 2009
    

Formula

G.f.: (2-3*x)/(1-3*x-x^2). - Simon Plouffe in his 1992 dissertation
From Gary W. Adamson, Jun 15 2003: (Start)
a(n) = ((3 + sqrt(13))/2)^n + ((3 - sqrt(13))/2)^n. See bronze mean (A098316).
A006190(n-2) + A006190(n) = a(n-1).
a(n)^2 - 13*A006190(n)^2 = 4(-1)^n. (End)
From Paul Barry, Nov 15 2003: (Start)
E.g.f.: 2*exp(3*x/2)*cosh(sqrt(13)*x/2).
a(n) = 2^(1-n)*Sum_{k=0..floor(n/2)} C(n, 2*k)* (13)^k * 3^(n-2*k).
a(n) = 2*T(n, 3i/2)*(-i)^n with T(n, x) Chebyshev's polynomials of the first kind (see A053120) and i^2=-1. (End)
From Hieronymus Fischer, Jan 02 2009: (Start)
fract(((3+sqrt(13))/2)^n) = (1/2)*(1+(-1)^n) - (-1)^n*((3+sqrt(13))/2)^(-n) = (1/2)*(1+(-1)^n) - ((3-sqrt(13))/2)^n.
See A001622 for a general formula concerning the fractional parts of powers of numbers x>1, which satisfy x-x^(-1)=floor(x).
a(n) = round(((3+sqrt(13))/2)^n) for n > 0. (End)
From Johannes W. Meijer, Jun 12 2010: (Start)
a(2n+1) = 3*A097783(n), a(2n) = A057076(n).
a(3n+1) = A041018(5n), a(3n+2) = A041018(5n+3) and a(3n+3) = 2*A041018(5n+4).
Limit_{k -> infinity} a(n+k)/a(k) = (a(n) + A006190(n)*sqrt(13))/2.
Limit_{n -> infinity} a(n)/A006190(n) = sqrt(13).
(End)
a(n) = sqrt(13*(A006190(n))^2 + 4*(-1)^n). - Vladimir Shevelev, Mar 13 2013
G.f.: G(0), where G(k) = 1 + 1/(1 - (x*(13*k-9))/((x*(13*k+4)) - 6/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 15 2013
a(n) = [x^n] ( (1 + 3*x + sqrt(1 + 6*x + 13*x^2))/2 )^n for n >= 1. - Peter Bala, Jun 23 2015
a(n) = Lucas(n,3), Lucas polynomials, L(n,x), evaluated at x=3. - G. C. Greubel, Jun 06 2019
a(n) = 2 * Sum_{k=0..n-2} A168561(n-2,k)*3^k + 3 * Sum_{k=0..n-1} A168561(n-1,k)*3^k, n>0. - R. J. Mathar, Feb 14 2024
a(n) = 2*A006190(n+1) - 3*A006190(n). - R. J. Mathar, Feb 14 2024
a(2*n+1) = 3 + 3*Sum_{k=1..n} a(2*k). - Greg Dresden and Canran Wang, Jul 11 2024
From Peter Bala, Jul 14 2025: (Start)
The following series telescope (Cf. A000032):
For k >= 1, Sum_{n >= 1} (-1)^((k+1)*(n+1)) * a(2*n*k)/(a((2*n-1)*k)*a((2*n+1)*k)) = 1/a(k)^2.
For positive even k, Sum_{n >= 1} 1/(a(k*n) - (a(k) + 2)/a(k*n)) = 1/(a(k) - 2) and
Sum_{n >= 1} (-1)^(n+1)/(a(k*n) + (a(k) - 2)/a(k*n)) = 1/(a(k) + 2).
For positive odd k, Sum_{n >= 1} 1/(a(k*n) - (-1)^n*(a(2*k) + 2)/a(k*n)) = (a(k) + 2)/(2*(a(2*k) - 2)) and
Sum_{n >= 1} (-1)^(n+1)/(a(k*n) - (-1)^n*(a(2*k) + 2)/a(k*n)) = (a(k) - 2)/(2*(a(2*k) - 2)). (End)

Extensions

Definition completed by M. F. Hasler, Mar 06 2020

A352361 Array read by ascending antidiagonals. A(n, k) = Fibonacci(k, n), where Fibonacci(n, x) are the Fibonacci polynomials.

Original entry on oeis.org

0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 1, 2, 2, 0, 0, 1, 3, 5, 3, 1, 0, 1, 4, 10, 12, 5, 0, 0, 1, 5, 17, 33, 29, 8, 1, 0, 1, 6, 26, 72, 109, 70, 13, 0, 0, 1, 7, 37, 135, 305, 360, 169, 21, 1, 0, 1, 8, 50, 228, 701, 1292, 1189, 408, 34, 0, 0, 1, 9, 65, 357, 1405, 3640, 5473, 3927, 985, 55, 1
Offset: 0

Views

Author

Peter Luschny, Mar 18 2022

Keywords

Comments

From Michael A. Allen, Mar 26 2023: (Start)
Row n is the n-metallonacci sequence for n>0.
A(n,k), for n > 0 and k > 0, is the number of tilings of a (k-1)-board (a board with dimensions (k-1) X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are n kinds of squares available. (End)

Examples

			Array, A(n,k), starts:
  n\k 0, 1, 2,  3,   4,    5,     6,      7,       8,        9, ...
  -------------------------------------------------------------------------
  [0] 0, 1, 0,  1,   0,    1,     0,      1,       0,        1, ... A000035;
  [1] 0, 1, 1,  2,   3,    5,     8,     13,      21,       34, ... A000045;
  [2] 0, 1, 2,  5,  12,   29,    70,    169,     408,      985, ... A000129;
  [3] 0, 1, 3, 10,  33,  109,   360,   1189,    3927,    12970, ... A006190;
  [4] 0, 1, 4, 17,  72,  305,  1292,   5473,   23184,    98209, ... A001076;
  [5] 0, 1, 5, 26, 135,  701,  3640,  18901,   98145,   509626, ... A052918;
  [6] 0, 1, 6, 37, 228, 1405,  8658,  53353,  328776,  2026009, ... A005668;
  [7] 0, 1, 7, 50, 357, 2549, 18200, 129949,  927843,  6624850, ... A054413;
  [8] 0, 1, 8, 65, 528, 4289, 34840, 283009, 2298912, 18674305, ... A041025;
  [9] 0, 1, 9, 82, 747, 6805, 61992, 564733, 5144589, 46866034, ... A099371;
      |  |  |  | A054602 |   A124152;
      |  |  |  A002522   A057721;
      |  |  A001477;
      |  A000012;
      A000004;
Antidiagonals, T(n, k), begin as:
  0;
  0, 1;
  0, 1, 0;
  0, 1, 1,  1;
  0, 1, 2,  2,   0;
  0, 1, 3,  5,   3,   1;
  0, 1, 4, 10,  12,   5,   0;
  0, 1, 5, 17,  33,  29,   8,   1;
  0, 1, 6, 26,  72, 109,  70,  13,  0;
  0, 1, 7, 37, 135, 305, 360, 169, 21, 1;
		

Crossrefs

Other versions of this array are A073133, A157103, A172236.
Rows n: A000035 (n=0), A000045 (n=1), A000129 (n=2), A006190 (n=3), A001076 (n=4), A052918 (n=5), A005668 (n=6), A054413 (n=7), A041025 (n=8), A099371 (n=9).
Columns k: A000004 (k=0), A000012 (k=1), A001477 (k=2), A002522 (k=3), A054602 (k=4), A057721 (k=5), A124152 (k=6).
Cf. A084844 (main diagonal), A352362 (Lucas polynomials), A350470 (Jacobsthal polynomials).
Sums include: A304357 (row sums), A304359.
Cf. A084845.

Programs

  • Magma
    A352361:= func< n, k | k le 1 select k else Evaluate(DicksonSecond(k-1, -1), n-k) >;
    [A352361(n, k): k in [0..n], n in [0..13]]; // G. C. Greubel, Sep 29 2024
    
  • Maple
    seq(seq(combinat:-fibonacci(k, n - k), k = 0..n), n = 0..11);
  • Mathematica
    Table[Fibonacci[k, n - k], {n, 0, 9}, {k, 0, n}] // Flatten
    (* or *)
    A[n_, k_] := With[{s = Sqrt[n^2 + 4]}, ((n + s)^k - (n - s)^k) / (2^k*s)];
    Table[Simplify[A[n, k]], {n, 0, 9}, {k, 0, 9}] // TableForm
  • PARI
    A(n, k) = ([1, k; 1, k-1]^n)[2, 1] ;
    export(A)
    for(k = 0, 9, print(parvector(10, n, A(n - 1, k))))
    
  • SageMath
    def A352361(n, k): return lucas_number1(k,n-k,-1)
    flatten([[A352361(n, k) for k in range(n+1)] for n in range(14)]) # G. C. Greubel, Sep 29 2024

Formula

A(n, k) = Sum_{j=0..floor((k-1)/2)} binomial(k-j-1, j)*n^(k-2*j-1).
A(n, k) = ((n + s)^k - (n - s)^k) / (2^k*s) where s = sqrt(n^2 + 4).
A(n, k) = [x^k] (x / (1 - n*x - x^2)).
A(n, k) = n^(k-1)*hypergeom([1 - k/2, 1/2 - k/2], [1 - k], -4/n^2) for n,k >= 1.
A(n, n) = T(2*n, n) = A084844(n).
From G. C. Greubel, Sep 29 2024: (Start)
T(n, k) = A(n-k, k) (antidiagonal triangle).
T(2*n+1, n+1) = A084845(n).
Sum_{k=0..n} T(n, k) = A304357(n) (row sums).
Sum_{k=0..n} (-1)^k*T(n, k) = (-1)*A304359(n). (End)

A054413 a(n) = 7*a(n-1) + a(n-2), with a(0)=1 and a(1)=7.

Original entry on oeis.org

1, 7, 50, 357, 2549, 18200, 129949, 927843, 6624850, 47301793, 337737401, 2411463600, 17217982601, 122937341807, 877779375250, 6267392968557, 44749530155149, 319514104054600, 2281348258537349, 16288951913816043, 116304011655249650, 830417033500563593
Offset: 0

Views

Author

Henry Bottomley, May 10 2000

Keywords

Comments

In general, sequences with recurrence a(n) = k*a(n-1) + a(n-2) and a(0)=1 (and a(-1)=0) have the generating function 1/(1-k*x-x^2). If k is odd (k>=3) they satisfy a(3n) = b(5n), a(3n+1) = b(5n+3), a(3n+2) = 2*b(5n+4) where b(n) is the sequence of denominators of continued fraction convergents to sqrt(k^2+4). [If k is even then a(n) is the sequence of denominators of continued fraction convergents to sqrt(k^2/4+1).]
a(p-1) == 53^((p-1)/2) (mod p), for odd primes p. - Gary W. Adamson, Feb 22 2009 [See A087475 for more info about this congruence. - Jason Yuen, Apr 05 2025]
From Johannes W. Meijer, Jun 12 2010: (Start)
For the sequence given above k=7 which implies that it is associated with A041091.
For a similar statement about sequences with recurrence a(n) = k*a(n-1) + a(n-2) but with a(0) = 2, and a(-1) = 0, see A086902; a sequence that is associated with A041090.
For more information follow the Khovanova link and see A087130, A140455 and A178765.
(End)
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 7's along the main diagonal and 1's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 08 2011
a(n) equals the number of words of length n on alphabet {0,1,...,7} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
From Michael A. Allen, Feb 21 2023: (Start)
Also called the 7-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence.
a(n) is the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are 7 kinds of squares available. (End)

Crossrefs

Row n=7 of A073133, A172236 and A352361.
Cf. A099367 (squares).

Programs

Formula

a(3n) = A041091(5n), a(3n+1) = A041091(5n+3), a(3n+2) = 2*A041091(5n+4).
G.f.: 1/(1 - 7x - x^2).
a(n) = U(n, 7*i/2)*(-i)^n with i^2=-1 and Chebyshev's U(n, x/2) = S(n, x) polynomials. See A049310.
a(n) = F(n, 7), the n-th Fibonacci polynomial evaluated at x=7. - T. D. Noe, Jan 19 2006
From Sergio Falcon, Sep 24 2007: (Start)
a(n) = (sigma^n - (-sigma)^(-n))/(sqrt(53)) with sigma = (7+sqrt(53))/2;
a(n) = Sum_{i=0..floor((n-1)/2)} binomial(n-1-i,i)*7^(n-1-2i). (End)
a(n) = ((7 + sqrt(53))^n - (7 - sqrt(53))^n)/(2^n*sqrt(53)). Offset 1. a(3)=50. - Al Hakanson (hawkuu(AT)gmail.com), Jan 17 2009
From Johannes W. Meijer, Jun 12 2010: (Start)
a(2n+1) = 7*A097836(n), a(2n) = A097838(n).
Lim_{k->oo} a(n+k)/a(k) = (A086902(n) + A054413(n-1)*sqrt(53))/2.
Lim_{n->oo} A086902(n)/A054413(n-1) = sqrt(53).
(End)
Sum_{n>=0} (-1)^n/(a(n)*a(n+1)) = (sqrt(53)-7)/2. - Vladimir Shevelev, Feb 23 2013
From Kai Wang, Feb 24 2020: (Start)
Sum_{m>=0} 1/(a(m)*a(m+2)) = 1/49.
Sum_{m>=0} 1/(a(2*m)*a(2*m+2)) = (sqrt(53)-7)/14.
In general, for sequences with recurrence f(n)= k*f(n-1)+f(n-2) and f(0)=1,
Sum_{m>=0} 1/(f(m)*f(m+2)) = 1/(k^2).
Sum_{m>=0} 1/(f(2*m)*f(2*m+2)) = (sqrt(k^2+4) - k)/(2*k). (End)
E.g.f.: (1/53)*exp(7*x/2)*(53*cosh(sqrt(53)*x/2) + 7*sqrt(53)*sinh(sqrt(53)*x/2)). - Stefano Spezia, Feb 26 2020
G.f.: x/(1 - 7*x - x^2) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (m*k + 7 - m + x)/(1 + m*k*x) ) for arbitrary m (a telescoping series). - Peter Bala, May 08 2024

Extensions

Formula corrected by Johannes W. Meijer, May 30 2010, Jun 02 2010
Extended by T. D. Noe, May 23 2011

A005668 Denominators of continued fraction convergents to sqrt(10).

Original entry on oeis.org

0, 1, 6, 37, 228, 1405, 8658, 53353, 328776, 2026009, 12484830, 76934989, 474094764, 2921503573, 18003116202, 110940200785, 683644320912, 4212806126257, 25960481078454, 159975692596981, 985814636660340, 6074863512559021, 37434995712014466, 230684837784645817
Offset: 0

Views

Author

Keywords

Comments

a(2*n+1) with b(2*n+1) := A005667(2*n+1), n >= 0, give all (positive integer) solutions to Pell equation b^2 - 10*a^2 = -1, a(2*n) with b(2*n) := A005667(2*n), n>=1, give all (positive integer) solutions to Pell equation b^2 - 10*a^2 = +1 (cf. Emerson reference).
Bisection: a(2*n)= 6*S(n-1,2*19) = 6*A078987(n-1), n >= 0 and a(2*n+1) = A097315(n), n >= 0, with S(n,x) Chebyshev's polynomials of the second kind. S(-1,x)=0. See A049310.
Sqrt(10) = 6/2 + 6/37 + 6/(37*1405) + 6/(1405*53353) + ... - Gary W. Adamson, Dec 21 2007
a(p) == 40^((p-1)/2) mod p, for odd primes p. - Gary W. Adamson, Feb 22 2009
For n>=2, a(n) equals the permanent of the (n-1)X(n-1) tridiagonal matrix with 6's along the main diagonal and 1's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 08 2011
For n>=1, a(n) equals the number of words of length n-1 on alphabet {0,1,...,6} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
From Michael A. Allen, Feb 15 2023: (Start)
Also called the 6-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are 6 kinds of squares available. (End)

Examples

			G.f. = x + 6*x^2 + 37*x^3 + 228*x^4 + 1405*x^5 + 8658*x^6 + 53353*x^7 + ...
		

References

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

Crossrefs

Row n=6 of A073133, A172236, A352361.

Programs

  • Magma
    [n le 2 select n-1 else 6*Self(n-1)+Self(n-2): n in [1..25]]; // Vincenzo Librandi, Feb 23 2013
    
  • Maple
    evalf(sqrt(10),200); convert(%,confrac,fractionlist); fractionlist;
    A005668:=-z/(-1+6*z+z**2); - Simon Plouffe in his 1992 dissertation.
    a := n -> `if`(n<2,n,6^(n-1)*hypergeom([1-n/2,(1-n)/2], [1-n], -1/9)):
    seq(simplify(a(n)), n=0..23); # Peter Luschny, Jun 28 2017
  • Mathematica
    LinearRecurrence[{6,1}, {0,1}, 30] (* Vincenzo Librandi, Feb 23 2013 *)
    a[ n_] := (-I)^(n - 1) ChebyshevU[ n - 1, 3 I]; (* Michael Somos, May 28 2014 *)
    a[ n_] := MatrixPower[ {{0, 1}, {1, 6}}, n + 1][[1, 1]]; (* Michael Somos, May 28 2014 *)
    Fibonacci[Range[0,30],6] (* G. C. Greubel, Jun 06 2019 *)
    Join[{0},Convergents[Sqrt[10],30]//Denominator] (* Harvey P. Dale, Dec 28 2022 *)
  • PARI
    {a(n) = ([0, 1; 1, 6]^(n+1)) [1, 1]}; /* Michael Somos, May 28 2014 */
    
  • PARI
    {a(n) = (-I)^(n-1) * polchebyshev( n-1, 2, 3*I)}; /* Michael Somos, May 28 2014 */
  • Sage
    from sage.combinat.sloane_functions import recur_gen3; it = recur_gen3(0,1,6,6,1,0); [next(it) for i in range(1,22)] # Zerinvary Lajos, Jul 09 2008
    
  • Sage
    [lucas_number1(n,6,-1) for n in range(0, 21)]# Zerinvary Lajos, Apr 24 2009
    

Formula

G.f.: x / (1 - 6*x - x^2).
a(n) = 6*a(n-1) + a(n-2).
a(n) = ((-i)^(n-1))*S(n-1, 3*i) with S(n, x) Chebyshev's polynomials of the second kind (see A049310) and i^2=-1.
a(n) = F(n, 6), the n-th Fibonacci polynomial evaluated at x=6. - T. D. Noe, Jan 19 2006
From Sergio Falcon, Sep 24 2007: (Start)
a(n) = ((3+sqrt(10))^n - (3-sqrt(10))^n)/(2*sqrt(10)).
a(n) = Sum_{i=0..floor((n-1)/2)} binomial(n-1-i,i)*6^(n-1-2*i). (End)
Sum_{n>=1}(-1)^(n-1)/(a(n)*a(n+1)) = sqrt(10) - 3. - Vladimir Shevelev, Feb 23 2013
a(-n) = -(-1)^n * a(n). - Michael Somos, May 28 2014
a(n) = 6^(n-1)*hypergeom([1-n/2, (1-n)/2], [1-n], -1/9) for n >= 2. - Peter Luschny, Jun 28 2017
G.f.: x/(1 - 6*x - x^2) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (m*k + 6 - m + x)/(1 + m*k*x) ) for arbitrary m (a telescoping series). - Peter Bala, May 08 2024

Extensions

Chebyshev comments from Wolfdieter Lang, Jan 21 2003

A027465 Cube of lower triangular normalized binomial matrix.

Original entry on oeis.org

1, 3, 1, 9, 6, 1, 27, 27, 9, 1, 81, 108, 54, 12, 1, 243, 405, 270, 90, 15, 1, 729, 1458, 1215, 540, 135, 18, 1, 2187, 5103, 5103, 2835, 945, 189, 21, 1, 6561, 17496, 20412, 13608, 5670, 1512, 252, 24, 1, 19683, 59049, 78732, 61236, 30618, 10206, 2268
Offset: 0

Views

Author

Keywords

Comments

Rows of A013610 reversed. - Michael Somos, Feb 14 2002
Row sums are powers of 4 (A000302), antidiagonal sums are A006190 (a(n) = 3*a(n-1) + a(n-2)). - Gerald McGarvey, May 17 2005
Triangle of coefficients in expansion of (3+x)^n.
Also: Pure Galton board of scheme (3,1). Also: Multiplicity (number) of pairs of n-dimensional binary vectors with dot product (overlap) k. There are 2^n = A000079(n) binary vectors of length n and 2^(2n) = 4^n = A000302(n) different pairs to form dot products k = Sum_{i=1..n} v[i]*u[i] between these, 0 <= k <= n. (Since dot products are symmetric, there are only 2^n*(2^n-1)/2 different non-ordered pairs, actually.) - R. J. Mathar, Mar 17 2006
Mirror image of A013610. - Zerinvary Lajos, Nov 25 2007
T(i,j) is the number of i-permutations of 4 objects a,b,c,d, with repetition allowed, containing j a's. - Zerinvary Lajos, Dec 21 2007
The antidiagonals of the sequence formatted as a square array (see Examples section) and summed with alternating signs gives a bisection of Fibonacci sequence, A001906. Example: 81-(27-1)=55. Similar rule applied to rows gives A000079. - Mark Dols, Sep 01 2009
Triangle T(n,k), read by rows, given by (3,0,0,0,0,0,0,0,...)DELTA (1,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 09 2011
T(n,k) = binomial(n,k)*3^(n-k), the number of subsets of [2n] with exactly k symmetric pairs, where elements i and j of [2n] form a symmetric pair if i+j=2n+1. Equivalently, if n couples attend a (ticketed) event that offers door prizes, then the number of possible prize distributions that have exactly k couples as dual winners is T(n,k). - Dennis P. Walsh, Feb 02 2012
T(n,k) is the number of ordered pairs (A,B) of subsets of {1,2,...,n} such that the intersection of A and B contains exactly k elements. For example, T(2,1) = 6 because we have ({1},{1}); ({1},{1,2}); ({2},{2}); ({2},{1,2}); ({1,2},{1}); ({1,2},{2}). Sum_{k=0..n} T(n,k)*k = A002697(n) (see comment there by Ross La Haye). - Geoffrey Critzer, Sep 04 2013
Also the convolution triangle of A000244. - Peter Luschny, Oct 09 2022

Examples

			Example: n = 3 offers 2^3 = 8 different binary vectors (0,0,0), (0,0,1), ..., (1,1,0), (1,1,1). a(3,2) = 9 of the 2^4 = 64 pairs have overlap k = 2: (0,1,1)*(0,1,1) = (1,0,1)*(1,0,1) = (1,1,0)*(1,1,0) = (1,1,1)*(1,1,0) = (1,1,1)*(1,0,1) = (1,1,1)*(0,1,1) = (0,1,1)*(1,1,1) = (1,0,1)*(1,1,1) = (1,1,0)*(1,1,1) = 2.
For example, T(2,1)=6 since there are 6 subsets of {1,2,3,4} that have exactly 1 symmetric pair, namely, {1,4}, {2,3}, {1,2,3}, {1,2,4}, {1,3,4}, and {2,3,4}.
The present sequence formatted as a triangular array:
     1
     3     1
     9     6     1
    27    27     9     1
    81   108    54    12    1
   243   405   270    90   15    1
   729  1458  1215   540  135   18   1
  2187  5103  5103  2835  945  189  21  1
  6561 17496 20412 13608 5670 1512 252 24 1
  ...
A013610 formatted as a triangular array:
  1
  1  3
  1  6   9
  1  9  27   27
  1 12  54  108   81
  1 15  90  270  405   243
  1 18 135  540 1215  1458   729
  1 21 189  945 2835  5103  5103  2187
  1 24 252 1512 5670 13608 20412 17496 6561
   ...
A099097 formatted as a square array:
      1     0     0    0   0 0 0 0 0 0 0 ...
      3     1     0    0   0 0 0 0 0 0 ...
      9     6     1    0   0 0 0 0 0 ...
     27    27     9    1   0 0 0 0 ...
     81   108    54   12   1 0 0 ...
    243   405   270   90  15 1 ...
    729  1458  1215  540 135 ...
   2187  5103  5103 2835 ...
   6561 17496 20412 ...
  19683 59049 ...
  59049 ...
		

Crossrefs

Programs

  • Haskell
    a027465 n k = a027465_tabl !! n !! k
    a027465_row n = a027465_tabl !! n
    a027465_tabl = iterate (\row ->
       zipWith (+) (map (* 3) (row ++ [0])) (map (* 1) ([0] ++ row))) [1]
    -- Reinhard Zumkeller, May 26 2013
  • Maple
    for i from 0 to 12 do seq(binomial(i, j)*3^(i-j), j = 0 .. i) od; # Zerinvary Lajos, Nov 25 2007
    # Uses function PMatrix from A357368. Adds column 1, 0, 0, ... to the left.
    PMatrix(10, n -> 3^(n-1)); # Peter Luschny, Oct 09 2022
  • Mathematica
    t[n_, k_] := Binomial[n, k]*3^(n-k); Table[t[n, n-k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Sep 19 2012 *)
  • PARI
    {T(n, k) = polcoeff( (3 + x)^n, k)}; /* Michael Somos, Feb 14 2002 */
    

Formula

Numerators of lower triangle of (b^2)[ i, j ] where b[ i, j ] = binomial(i-1, j-1)/2^(i-1) if j <= i, 0 if j > i.
Triangle whose (i, j)-th entry is binomial(i, j)*3^(i-j).
a(n, m) = 4^(n-1)*Sum_{j=m..n} b(n, j)*b(j, m) = 3^(n-m)*binomial(n-1, m-1), n >= m >= 1; a(n, m) := 0, n < m. G.f. for m-th column: (x/(1-3*x))^m (m-fold convolution of A000244, powers of 3). - Wolfdieter Lang, Feb 2006
G.f.: 1 / (1 - x(3+y)).
a(n,k) = 3*a(n-1,k) + a(n-1,k-1) - R. J. Mathar, Mar 17 2006
From the formalism of A133314, the e.g.f. for the row polynomials of A027465 is exp(x*t)*exp(3x). The e.g.f. for the row polynomials of the inverse matrix is exp(x*t)*exp(-3x). p iterates of the matrix give the matrix with e.g.f. exp(x*t)*exp(p*3x). The results generalize for 3 replaced by any number. - Tom Copeland, Aug 18 2008
T(n,k) = A164942(n,k)*(-1)^k. - Philippe Deléham, Oct 09 2011
Let P and P^T be the Pascal matrix and its transpose and H = P^3 = A027465. Then from the formalism of A132440 and A218272,
exp[x*z/(1-3z)]/(1-3z) = exp(3z D_z z) e^(x*z)= exp(3D_x x D_x) e^(z*x)
= (1 z z^2 z^3 ...) H (1 x x^2/2! x^3/3! ...)^T
= (1 x x^2/2! x^3/3! ...) H^T (1 z z^2 z^3 ...)^T = Sum_{n>=0} (3z)^n L_n(-x/3), where D is the derivative operator and L_n(x) are the regular (not normalized) Laguerre polynomials. - Tom Copeland, Oct 26 2012
E.g.f. for column k: x^k/k! * exp(3x). - Geoffrey Critzer, Sep 04 2013

A157103 Array A(n, k) = Fibonacci(n+1, k), with A(n, 0) = A(n, n) = 1, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 3, 5, 3, 1, 1, 5, 12, 10, 4, 1, 1, 8, 29, 33, 17, 5, 1, 1, 13, 70, 109, 72, 26, 6, 1, 1, 21, 169, 360, 305, 135, 37, 7, 1, 1, 34, 408, 1189, 1292, 701, 228, 50, 8, 1, 1, 55, 985, 3927, 5473, 3640, 1405, 357, 65, 9, 1
Offset: 0

Views

Author

Gary W. Adamson, Feb 22 2009

Keywords

Comments

From Michael A. Allen, Mar 30 2023: (Start)
Column k is the k-metallonacci sequence for k > 0.
T(n,k) is, for n > 0 and k > 0, the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are k kinds of squares available. (End)

Examples

			Array begins:
  1,  1,   1,    1,     1,     1,      1,      1, ... (A000012);
  1,  1,   2,    3,     4,     5,      6,      7, ... (A000027);
  1,  2,   5,   10,    17,    26,     37,     50, ... (A002522);
  1,  3,  12,   33,    72,   135,    228,    357, ...;
  1,  5,  29,  109,   305,   701,   1405,   2549, ...;
  1,  8,  70,  360,  1292,  3640,   8658,  18200, ...;
  1, 13, 169, 1189,  5473, 18901,  53353, 129949, ...;
  1, 21, 408, 3927, 23184, 98145, 328776, 927843, ...;
  ...
First few rows of the triangle:
  1;
  1,   1;
  1,   1,    1;
  1,   2,    2,     1;
  1,   3,    5,     3,     1;
  1,   5,   12,    10,     4,     1;
  1,   8,   29,    33,    17,     5,     1;
  1,  13,   70,   109,    72,    26,     6,     1;
  1,  21,  169,   360,   305,   135,    37,     7,    1;
  1,  34,  408,  1189,  1292,   701,   228,    50,    8,   1;
  1,  55,  985,  3927,  5473,  3640,  1405,   357,   65,   9,   1;
  1,  89, 2378, 12970, 23184, 18901,  8658,  2549,  528,  82,  10,  1;
  1, 144, 5741, 42837, 98209, 98145, 53353, 18200, 4289, 747, 101, 11, 1;
  ...
Example: Column 3 = (1, 3, 10, 33, 109, 360, ...) = A006190.
		

Crossrefs

Essentially the transpose of A073133, A172236, A352361.

Programs

  • Magma
    A157103:= func< n,k | k eq 0 or k eq n select 1 else Evaluate(DicksonSecond(n, -1), k) >;
    [A157103(n-k, k): k in [0..n], n in [0..15]]; // G. C. Greubel, Jan 11 2022
    
  • Maple
    A157103 := proc(n,k)
        if k = 0 then
            1;
        else
            mul(k-2*I*cos(l*Pi/(n+1)),l=1..n) ;
            combine(%,trig) ;
            round(%) ;
        end if;
    end proc:
    seq( seq(A157103(d-k,k),k=0..d),d=0..12) ; # R. J. Mathar, Feb 27 2023
  • Mathematica
    (* First program *)
    T[, 0]=1; T[n, n_]=1; T[, ]=0;
    T[n_, k_] /; 0 <= k <= n := k T[n-1, k] + T[n-2, k];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}]//Flatten (* Jean-François Alcover, Aug 07 2018 *)
    (* Second program *)
    T[n_, k_]:= If[k==0 || k==n, 1, Fibonacci[n-k+1, k]];
    Table[T[n, k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Jan 11 2022 *)
  • Sage
    def A157103(n,k): return 1 if (k==0 or k==n) else lucas_number1(n+1, k, -1)
    flatten([[A157103(n-k, k) for k in (0..n)] for n in (0..10)]) # G. C. Greubel, Jan 11 2022

Formula

A(n, k) = Fibonacci(n+1, k), with A(n, 0) = A(n, n) = 1 (array).
A(n, 1) = A000045(n+1).
T(n, k) = k*T(n-1, k) + T(n-2, k) with T(n, 0) = T(n, n) = 1 (triangle).
From G. C. Greubel, Jan 11 2022: (Start)
T(n, k) = Fibonacci(n-k+1, k), with T(n, 0) = T(n, n) = 1.
T(2*n, n) = A084845(n) for n >= 1, with T(0, 0) = 1.
T(2*n+1, n+1) = A084844(n). (End)

Extensions

Edited by G. C. Greubel, Jan 11 2022
Previous Showing 31-40 of 146 results. Next