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 57 results. Next

A003946 Expansion of (1+x)/(1-3*x).

Original entry on oeis.org

1, 4, 12, 36, 108, 324, 972, 2916, 8748, 26244, 78732, 236196, 708588, 2125764, 6377292, 19131876, 57395628, 172186884, 516560652, 1549681956, 4649045868, 13947137604, 41841412812, 125524238436, 376572715308, 1129718145924
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 4.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6. - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001
a(n) is the number of nonreversing random walks of the length of n edges on a two-dimensional square lattice, all beginning at a fixed point P. - Pawel P. Mazur (Pawel.Mazur(AT)pwr.wroc.pl), Apr 06 2005
Binomial transform of {1, 3, 5, 11, 21, 43, ...}, see A001045. Binomial transform is {1, 5, 21, 85, 341, 1365, ...}, see A002450. - Philippe Deléham, Jul 22 2005
For n >= 2, a(n) is equal to the number of functions f:{1,2,...,n+1}->{1,2,3} such that for fixed, different x_1, x_2 in {1,2,...,n} and fixed y_1, y_2 in {1,2,3} we have f(x_1) <> y_1 and f(x_2) <> y_2. - Milan Janjic, Apr 19 2007
Equals row sums of triangle A143865. - Gary W. Adamson, Sep 04 2008
Equals INVERT transform of the odd integers = 1/(1 - x - 3x^2 - 5x^3 - ...). - Gary W. Adamson, Jul 27 2009
a(n) is the number of generalized compositions of n+1 when there are 2 *i-1 different types of the part i, (i=1,2,...). - Milan Janjic, Aug 26 2010
Number of length-n strings of 4 letters with no two adjacent letters identical. The general case (strings of r letters) is the sequence with g.f. (1+x)/(1-(r-1)*x). - Joerg Arndt, Oct 11 2012
The sequence is the INVERTi transform of A015448: (1, 5, 21, 89, 377, ...). - Gary W. Adamson, Aug 06 2016
Let D(m) = {d(m,i)}, i = 1..q, denote the set of the q divisors of a number m, and consider s1(m) and s2(m) the sums of the divisors that are congruent to 1 and 2 (mod 3) respectively. For n > 0, the sequence a(n) lists the numbers m such that s1(m) = 5 and s2(m) = 2. - Michel Lagneau, Feb 09 2017
a(n) is the number of quaternary sequences of length n such that no two consecutive terms have distance 2. - David Nacin, May 31 2017
Also the number of maximal cliques in the n-Sierpinski gasket graph. - Eric W. Weisstein, Dec 01 2017
Number of 3-permutations of n elements avoiding the patterns 231, 321. See Bonichon and Sun. - Michel Marcus, Aug 19 2022

Examples

			G.f. = 1 + 4*x + 12*x^2 + 36*x^3 + 108*x^4 + 324*x^5 + 972*x^6 + 2916*x^7 + ...
		

Crossrefs

Cf. A029653, A143865, column 4 in A265583, A015448.

Programs

Formula

a(n) = floor(4*3^(n-1)). - Michael Somos, Jun 18 2002
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 2. - Philippe Deléham, Jul 10 2005
The Hankel transform of this sequence is [1,-4,0,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
a(n + 1) = (((1 + sqrt(-11))/2)^n + ((1 - sqrt(-11))/2)^n)^2 - (((1 + sqrt(-11))/2)^n - ((1 - sqrt(-11))/2)^n)^2. - Raphie Frank, Dec 07 2015
From Mario C. Enriquez, Apr 01 2017: (Start)
(L(a(n+k)) - 1)/a(n) reduces to the form C/a(n-1), where n > 1, k >= 0, L(a(n)) is the a(n)-th Lucas number and C = (L(a(n+k)) - 1)/3.
(L(a(n+k)) - 1)/3 mod (L(a(n)) - 1)/3 = (L(a(n)) - 1)/3 - 1, where n >= 1, k >= 0 and L(a(n)) is the a(n)-th Lucas number. (End)
E.g.f.: (4*exp(3*x) - 1)/3. - Stefano Spezia, Jan 31 2025

Extensions

Additional comments from Michael Somos, Jun 18 2002
Edited by N. J. A. Sloane, Dec 04 2009

A003948 Expansion of (1+x)/(1-5*x).

Original entry on oeis.org

1, 6, 30, 150, 750, 3750, 18750, 93750, 468750, 2343750, 11718750, 58593750, 292968750, 1464843750, 7324218750, 36621093750, 183105468750, 915527343750, 4577636718750, 22888183593750, 114440917968750, 572204589843750, 2861022949218750, 14305114746093750
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 6.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954, m is 2, 3, 4, 5, 6. - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001
Hamiltonian path in S_4 X P_2n.
For n>=1, a(n+1) is equal to the number of functions f:{1,2,...,n+1}->{1,2,3,4,5,6} such that for fixed, different x_1, x_2,...,x_n in {1,2,...,n+1} and fixed y_1, y_2,...,y_n in {1,2,3,4,5,6} we have f(x_i)<>y_i, (i=1..n). - Milan Janjic, May 10 2007
For n>=1, a(n) equals the numbers of words of length n over the alphabet {0..5} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 30 2017]
a(n) equals the numbers of sequences of length n on {0,...,5} where no two adjacent terms differ by three. - David Nacin, May 30 2017
It appears that these are the only n>1 for which alpha(n)=2n, where alpha(n) is the entry point of n in the Fibonacci sequence, see A001177. - Philippe Schnoebelen, Apr 11 2024

Crossrefs

Programs

  • GAP
    Concatenation([1], List([1..30], n-> 6*5^(n-1) )); # G. C. Greubel, Sep 24 2019
  • Magma
    [1] cat [6*5^(n-1): n in [1..30]]; // G. C. Greubel, Sep 24 2019
    
  • Maple
    k := 6; if n = 0 then 1 else k*(k-1)^(n-1); fi;
  • Mathematica
    q = 6; Join[{a = 1}, Table[If[n != 0, a = q*a - a, a = q*a], {n, 0, 25}]] (* and *) Join[{1}, 6*5^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    Join[{1},NestList[5#&,6,30]] (* Harvey P. Dale, Dec 31 2013 *)
    CoefficientList[Series[(1+x)/(1-5x), {x,0,30}], x] (* Michael De Vlieger, Dec 10 2016 *)
  • PARI
    Vec((1+x)/(1-5*x)+O(x^30)) \\ Charles R Greathouse IV, Nov 20 2012
    
  • Sage
    [1]+[6*5^(n-1) for n in (1..30)] # G. C. Greubel, Sep 24 2019
    

Formula

G.f.: (1+x)/(1-5*x).
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 4. - Philippe Deléham, Jul 10 2005
The Hankel transform of this sequence is [1,-6,0,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
a(n) = 6*5^(n-1) for n>0, a(0)=1. - Vincenzo Librandi, Nov 18 2010
G.f.: 2/x - 5 - 8/(x*U(0)) where U(k)= 1 + 2/(3^k - 3^k/(2 + 1 - 12*x*3^k/(6*x*3^k + 1/U(k+1)))) ; (continued fraction, 4-step). - Sergei N. Gladkovskii, Oct 30 2012
E.g.f.: (6*exp(5*x) - 1)/5. - Ilya Gutkovskiy, Dec 10 2016
Sum_{n>=0} 1/a(n) = 29/24. - Bernard Schott, Oct 25 2021

Extensions

Definition corrected by Frans J. Faase, Feb 07 2009
Edited by N. J. A. Sloane, Dec 04 2009

A003950 Expansion of g.f.: (1+x)/(1-7*x).

Original entry on oeis.org

1, 8, 56, 392, 2744, 19208, 134456, 941192, 6588344, 46118408, 322828856, 2259801992, 15818613944, 110730297608, 775112083256, 5425784582792, 37980492079544, 265863444556808, 1861044111897656, 13027308783283592, 91191161482985144, 638338130380896008
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 8.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001.
For n>=1, a(n) equals the number of words of length n on the alphabet {0,1,...,7} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 31 2017]
a(n) is the number of octonary sequences of length n such that no two consecutive terms have distance 4. - David Nacin, May 31 2017

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 6. - Philippe Deléham, Jul 10 2005
From Philippe Deléham, Nov 21 2007: (Start)
a(n) = 8*7^(n-1) for n>=1, a(0)=1 .
G.f.: (1+x)/(1-7x).
The Hankel transform of this sequence is [1,-8,0,0,0,0,0,0,0,0,...]. (End)
a(0)=1, a(1)=8, a(n) = 7*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (8*exp(7*x) - 1)/7. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003952 Expansion of g.f.: (1+x)/(1-9*x).

Original entry on oeis.org

1, 10, 90, 810, 7290, 65610, 590490, 5314410, 47829690, 430467210, 3874204890, 34867844010, 313810596090, 2824295364810, 25418658283290, 228767924549610, 2058911320946490, 18530201888518410, 166771816996665690, 1500946352969991210, 13508517176729920890
Offset: 0

Views

Author

N. J. A. Sloane, Mar 15 1996

Keywords

Comments

Coordination sequence for infinite tree with valency 10.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001
Except 1, all terms are in A033583. - Vincenzo Librandi, May 26 2014
For n>=1, a(n) equals the number of words of length n on alphabet {0,1,...,9} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 31 2017]
a(n) is the number of sequences over the alphabet {0,1,...,9} of length n such that no two consecutive terms have distance 5. - David Nacin, May 31 2017

Crossrefs

Programs

Formula

a(n) = (10*9^n - 0^n)/9. Binomial transform is A000042. - Paul Barry, Jan 29 2004
G.f.: (1+x)/(1-9*x). - Philippe Deléham, Jan 31 2004
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 8. - Philippe Deléham, Jul 10 2005
The Hankel transform of this sequence is: [1,-10,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
E.g.f.: (10*exp(9*x) - 1)/9. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A029653 Numbers in (2,1)-Pascal triangle (by row).

Original entry on oeis.org

1, 2, 1, 2, 3, 1, 2, 5, 4, 1, 2, 7, 9, 5, 1, 2, 9, 16, 14, 6, 1, 2, 11, 25, 30, 20, 7, 1, 2, 13, 36, 55, 50, 27, 8, 1, 2, 15, 49, 91, 105, 77, 35, 9, 1, 2, 17, 64, 140, 196, 182, 112, 44, 10, 1, 2, 19, 81, 204, 336, 378, 294, 156, 54, 11, 1, 2, 21, 100, 285
Offset: 0

Views

Author

Keywords

Comments

Reverse of A029635. Row sums are A003945. Diagonal sums are Fibonacci(n+2) = Sum_{k=0..floor(n/2)} (2n-3k)*C(n-k,n-2k)/(n-k). - Paul Barry, Jan 30 2005
Riordan array ((1+x)/(1-x), x/(1-x)). The signed triangle (-1)^(n-k)T(n,k) or ((1-x)/(1+x), x/(1+x)) is the inverse of A055248. Row sums are A003945. Diagonal sums are F(n+2). - Paul Barry, Feb 03 2005
Row sums = A003945: (1, 3, 6, 12, 24, 48, 96, ...) = (1, 3, 7, 15, 31, 63, 127, ...) - (0, 0, 1, 3, 7, 15, 31, ...); where (1, 3, 7, 15, ...) = A000225. - Gary W. Adamson, Apr 22 2007
Triangle T(n,k), read by rows, given by (2,-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 17 2011
A029653 is jointly generated with A208510 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*v(n-1)x and v(n,x)=u(n-1,x)+x*v(n-1,x)+1. See the Mathematica section. - Clark Kimberling, Feb 28 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 04 2013
The n-th row polynomial is (2 + x)*(1 + x)^(n-1) for n >= 1. More generally, the n-th row polynomial of the Riordan array ( (1-a*x)/(1-b*x), x/(1-b*x) ) is (b - a + x)*(b + x)^(n-1) for n >= 1. - Peter Bala, Feb 25 2018

Examples

			The triangle T(n,k) begins:
n\k 0  1  2   3   4   5   6   7  8  9 10 ...
0:  1
1:  2  1
2:  2  3  1
3:  2  5  4   1
4:  2  7  9   5   1
5:  2  9 16  14   6   1
6:  2 11 25  30  20   7   1
7:  2 13 36  55  50  27   8   1
8:  2 15 49  91 105  77  35   9  1
9:  2 17 64 140 196 182 112  44 10  1
10: 2 19 81 204 336 378 294 156 54 11  1
... Reformatted. - _Wolfdieter Lang_, Jan 09 2015
With the array M(k) as defined in the Formula section, the infinite product M(0)*M(1)*M(2)*... begins
/1        \/1         \/1        \      /1        \
|2 1      ||0 1       ||0 1      |      |2 1      |
|2 1 1    ||0 2 1     ||0 0 1    |... = |2 3 1    |
|2 1 1 1  ||0 2 1 1   ||0 0 2 1  |      |2 5 4 1  |
|2 1 1 1 1||0 2 1 1 1 ||0 0 2 1 1|      |2 7 9 5 1|
|...      ||...       ||...      |      |...      |
- _Peter Bala_, Dec 27 2014
		

References

  • Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.

Crossrefs

(d, 1) Pascal triangles: A007318(d=1), A093560(3), A093561(4), A093562(5), A093563(6), A093564(7), A093565(8), A093644(9), A093645(10).

Programs

  • Haskell
    a029653 n k = a029653_tabl !! n !! k
    a029653_row n = a029653_tabl !! n
    a029653_tabl = [1] : iterate
                   (\xs -> zipWith (+) ([0] ++ xs) (xs ++ [0])) [2, 1]
    -- Reinhard Zumkeller, Dec 16 2013
    
  • Maple
    A029653 :=  proc(n,k)
    if n = 0 then
      1;
    else
      binomial(n-1, k)+binomial(n, k)
    fi
    end proc: # R. J. Mathar, Jun 30 2013
  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + x*v[n - 1, x];
    v[n_, x_] := u[n - 1, x] + x*v[n - 1, x] + 1;
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]  (* A208510 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A029653 *)
    (* Clark Kimberling, Feb 28 2012 *)
  • Python
    from sympy import Poly
    from sympy.abc import x
    def u(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x)
    def v(n, x): return 1 if n==1 else u(n - 1, x) + x*v(n - 1, x) + 1
    def a(n): return Poly(v(n, x), x).all_coeffs()[::-1]
    for n in range(1, 13): print(a(n)) # Indranil Ghosh, May 27 2017
    
  • Python
    from math import comb, isqrt
    def A029653(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),a:=n-comb(r+1,2))*((r<<1)-a)//r if n else 1 # Chai Wah Wu, Nov 12 2024

Formula

T(n, k) = C(n-2, k-1) + C(n-2, k) + C(n-1, k-1) + C(n-1, k) except for n=0.
G.f.: (1 + x + y + xy)/(1 - y - xy). - Ralf Stephan, May 17 2004
T(n, k) = (2n-k)*binomial(n, n-k)/n, n, k > 0. - Paul Barry, Jan 30 2005
Sum_{k=0..n} T(n, k)*x^k gives A003945-A003954 for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10. - Philippe Deléham, Jul 10 2005
T(n, k) = C(n-1, k) + C(n, k). - Philippe Deléham, Jul 10 2005
Equals A097806 * A007318, i.e., the pairwise operator * Pascal's Triangle as infinite lower triangular matrices. - Gary W. Adamson, Apr 22 2007
From Peter Bala, Dec 27 2014: (Start)
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(2 + 5*x + 4*x^2/2! + x^3/3!) = 2 + 7*x + 16*x^2/2! + 30*x^3/3! + 50*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the lower unit triangular array with 1's on the main diagonal and 1's everywhere else below the main diagonal except for the first column which consists of the sequence [1,2,2,2,...]. For k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/ having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. Then the present triangle equals the infinite product M(0)*M(1)*M(2)*... (which is clearly well-defined). See the Example section. (End)

Extensions

More terms from James Sellers

A170732 Expansion of g.f.: (1+x)/(1 - 12*x).

Original entry on oeis.org

1, 13, 156, 1872, 22464, 269568, 3234816, 38817792, 465813504, 5589762048, 67077144576, 804925734912, 9659108818944, 115909305827328, 1390911669927936, 16690940039135232, 200291280469622784, 2403495365635473408, 28841944387625680896, 346103332651508170752
Offset: 0

Views

Author

N. J. A. Sloane, Dec 05 2009

Keywords

Comments

For n >= 1, a(n) equals the number of words of length n-1 on the alphabet {0,1,...,12} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015

Crossrefs

Programs

  • GAP
    k:=13;; Concatenation([1], List([1..25], n-> k*(k-1)^(n-1) )); # G. C. Greubel, Sep 24 2019
  • Magma
    k:=13; [1] cat [k*(k-1)^(n-1): n in [1..25]]; // G. C. Greubel, Sep 24 2019
    
  • Maple
    k:=13; seq(`if`(n=0, 1, k*(k-1)^(n-1)), n = 0..25); # G. C. Greubel, Sep 24 2019
  • Mathematica
    Join[{1}, 13*12^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    Join[{1},NestList[12#&,13,20]] (* Harvey P. Dale, Nov 24 2024 *)
  • PARI
    a(n)=if(n,13*12^(n-1),1) \\ Charles R Greathouse IV, Jul 01 2016
    
  • Python
    for i in range(1001):print(i,13*12**(i-1) if i>0 else 1) # Kenny Lau, Aug 01 2017
    
  • Sage
    k=13; [1]+[k*(k-1)^(n-1) for n in (1..25)] # G. C. Greubel, Sep 24 2019
    

Formula

a(0)=1; for n > 0, a(n) = 13*12^(n-1). - Vincenzo Librandi, Dec 05 2009
E.g.f.: (13*exp(12*x) - 1)/12. - G. C. Greubel, Sep 24 2019

A170733 Expansion of g.f.: (1+x)/(1-13*x).

Original entry on oeis.org

1, 14, 182, 2366, 30758, 399854, 5198102, 67575326, 878479238, 11420230094, 148462991222, 1930018885886, 25090245516518, 326173191714734, 4240251492291542, 55123269399790046, 716602502197270598, 9315832528564517774, 121105822871338731062, 1574375697327403503806
Offset: 0

Views

Author

N. J. A. Sloane, Dec 04 2009

Keywords

Comments

For n>=1, a(n) equals the numbers of words of length n-1 on alphabet {0,1,...,13} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015

Crossrefs

Programs

  • GAP
    k:=14;; Concatenation([1], List([1..25], n-> k*(k-1)^(n-1) )); # G. C. Greubel, Sep 24 2019
  • Magma
    k:=14; [1] cat [k*(k-1)^(n-1): n in [1..25]]; // G. C. Greubel, Sep 24 2019
    
  • Maple
    k:=14; seq(`if`(n=0, 1, k*(k-1)^(n-1)), n = 0..25); # G. C. Greubel, Sep 24 2019
  • Mathematica
    Join[{1}, 14*13^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    CoefficientList[Series[(1+x)/(1-13x), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 10 2012 *)
    Join[{1},NestList[13#&,14,20]] (* Harvey P. Dale, Oct 09 2017 *)
  • PARI
    vector(26, n, k=14; if(n==1, 1, k*(k-1)^(n-2))) \\ G. C. Greubel, Sep 24 2019
    
  • Sage
    k=14; [1]+[k*(k-1)^(n-1) for n in (1..25)] # G. C. Greubel, Sep 24 2019
    

Formula

a(n) = Sum_{k=0..n} A097805(n,k)*(-1)^(n-k)*14^k. - Philippe Deléham, Dec 04 2009
a(0) = 1; for n>0, a(n) = 14*13^(n-1). - Vincenzo Librandi, Dec 05 2009
a(0)=1, a(1)=14, a(n) = 13*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (14*exp(13*x) - 1)/13. - G. C. Greubel, Sep 24 2019

A170734 Expansion of g.f.: (1+x)/(1-14*x).

Original entry on oeis.org

1, 15, 210, 2940, 41160, 576240, 8067360, 112943040, 1581202560, 22136835840, 309915701760, 4338819824640, 60743477544960, 850408685629440, 11905721598812160, 166680102383370240, 2333521433367183360, 32669300067140567040, 457370200939967938560, 6403182813159551139840
Offset: 0

Views

Author

N. J. A. Sloane, Dec 04 2009

Keywords

Comments

For n>=1, a(n) equals the numbers of words of length n-1 on alphabet {0,1,...,14} with no two adjacent letters identical. -Milan Janjic, Jan 31 2015

Crossrefs

Programs

  • GAP
    k:=15;; Concatenation([1], List([1..25], n-> k*(k-1)^(n-1) )); # G. C. Greubel, Sep 24 2019
  • Magma
    k:=15; [1] cat [k*(k-1)^(n-1): n in [1..25]]; // G. C. Greubel, Sep 24 2019
    
  • Maple
    k:=15; seq(`if`(n=0, 1, k*(k-1)^(n-1)), n = 0..25); # G. C. Greubel, Sep 24 2019
  • Mathematica
    Join[{1}, 15*14^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    CoefficientList[Series[(1+x)/(1-14x), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 10 2012 *)
  • PARI
    vector(26, n, k=15; if(n==1, 1, k*(k-1)^(n-2))) \\ G. C. Greubel, Sep 24 2019
    
  • Sage
    k=15; [1]+[k*(k-1)^(n-1) for n in (1..25)] # G. C. Greubel, Sep 24 2019
    

Formula

a(n) = Sum_{k=0..n} A097805(n,k)*(-1)^(n-k)*15^k. - Philippe Deléham, Dec 04 2009
a(0) = 1; for n>0, a(n) = 15*14^(n-1). - Vincenzo Librandi, Dec 05 2009
a(0) = 1, a(1) = 15, a(n) = 14*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (15*exp(14*x) -1)/14. - G. C. Greubel, Sep 24 2019

A167916 Number of reduced words of length n in Coxeter group on 12 generators S_i with relations (S_i)^2 = (S_i S_j)^16 = I.

Original entry on oeis.org

1, 12, 132, 1452, 15972, 175692, 1932612, 21258732, 233846052, 2572306572, 28295372292, 311249095212, 3423740047332, 37661140520652, 414272545727172, 4556998002998892, 50126978032987746, 551396758362864480
Offset: 0

Views

Author

John Cannon and N. J. A. Sloane, Dec 03 2009

Keywords

Comments

The initial terms coincide with those of A003954, although the two sequences are eventually different.
Computed with MAGMA using commands similar to those used to compute A154638.

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Integers(), 30); Coefficients(R!( (1+x)*(1-x^16)/(1-11*x+65*x^16-55*x^17) )); // G. C. Greubel, Nov 10 2023
    
  • Mathematica
    CoefficientList[Series[(1+t)*(1-t^16)/(1-11*t+65*t^16-55*t^17), {t,0,50}], t] (* G. C. Greubel, Jul 01 2016; Nov 10 2023 *)
    coxG[{16,55,-10}] (* The coxG program is at A169452 *) (* G. C. Greubel, Nov 10 2023 *)
  • SageMath
    def A167916_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( (1+x)*(1-x^16)/(1-11*x+65*x^16-55*x^17) ).list()
    A167916_list(30) # G. C. Greubel, Nov 10 2023

Formula

G.f.: (t^16 + 2*t^15 + 2*t^14 + 2*t^13 + 2*t^12 + 2*t^11 + 2*t^10 + 2*t^9 + 2*t^8 + 2*t^7 + 2*t^6 + 2*t^5 + 2*t^4 + 2*t^3 + 2*t^2 + 2*t + 1)/( 55*t^16 - 10*t^15 - 10*t^14 - 10*t^13 - 10*t^12 - 10*t^11 - 10*t^10 - 10*t^9 - 10*t^8 - 10*t^7 - 10*t^6 - 10*t^5 - 10*t^4 - 10*t^3 - 10*t^2 - 10*t + 1).
From G. C. Greubel, Nov 10 2023: (Start)
a(n) = 10*Sum_{j=1..15} a(n-j) - 55*a(n-16).
G.f.: (1+x)*(1-x^16)/(1 - 11*x + 65*x^16 - 55*x^17). (End)

A158497 Triangle T(n,k) formed by the coordination sequences and the number of leaves for trees.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 1, 3, 6, 12, 1, 4, 12, 36, 108, 1, 5, 20, 80, 320, 1280, 1, 6, 30, 150, 750, 3750, 18750, 1, 7, 42, 252, 1512, 9072, 54432, 326592, 1, 8, 56, 392, 2744, 19208, 134456, 941192, 6588344, 1, 9, 72, 576, 4608, 36864, 294912, 2359296, 18874368, 150994944, 1, 10, 90, 810, 7290, 65610, 590490, 5314410, 47829690, 430467210, 3874204890
Offset: 0

Views

Author

Thomas Wieder, Mar 20 2009

Keywords

Comments

Consider the k-fold Cartesian products CP(n,k) of the vector A(n) = [1, 2, 3, ..., n].
An element of CP(n,k) is a n-tuple T_t of the form T_t = [i_1, i_2, i_3, ..., i_k] with t=1, .., n^k.
We count members T of CP(n,k) which satisfy some condition delta(T_t), so delta(.) is an indicator function which attains values of 1 or 0 depending on whether T_t is to be counted or not; the summation sum_{CP(n,k)} delta(T_t) over all elements T_t of CP produces the count.
For the triangle here we have delta(T_t) = 0 if for any two i_j, i_(j+1) in T_t one has i_j = i_(j+1): T(n,k) = Sum_{CP(n,k)} delta(T_t) = Sum_{CP(n,k)} delta(i_j = i_(j+1)).
The test on i_j > i_(j+1) generates A158498. One gets the Pascal triangle A007318 if the indicator function tests whether for any two i_j, i_(j+1) in T_t one has i_j >= i_(j+1).
Use of other indicator functions can also calculate the Bell numbers A000110, A000045 or A000108.

Examples

			Array, A(n, k) = n*(n-1)^(k-1) for n > 1, A(n, k) = 1 otherwise, begins as:
  1,  1,   1,    1,     1,      1,       1,        1,        1, ... A000012;
  1,  1,   1,    1,     1,      1,       1,        1,        1, ... A000012;
  1,  2,   2,    2,     2,      2,       2,        2,        2, ... A040000;
  1,  3,   6,   12,    24,     48,      96,      192,      384, ... A003945;
  1,  4,  12,   36,   108,    324,     972,     2916,     8748, ... A003946;
  1,  5,  20,   80,   320,   1280,    5120,    20480,    81920, ... A003947;
  1,  6,  30,  150,   750,   3750,   18750,    93750,   468750, ... A003948;
  1,  7,  42,  252,  1512,   9072,   54432,   326592,  1959552, ... A003949;
  1,  8,  56,  392,  2744,  19208,  134456,   941192,  6588344, ... A003950;
  1,  9,  72,  576,  4608,  36864,  294912,  2359296, 18874368, ... A003951;
  1, 10,  90,  810,  7290,  65610,  590490,  5314410, 47829690, ... A003952;
  1, 11, 110, 1100, 11000, 110000, 1100000, 11000000, ............. A003953;
  1, 12, 132, 1452, 15972, 175692, 1932612, 21258732, ............. A003954;
  1, 13, 156, 1872, 22464, 269568, 3234816, 38817792, ............. A170732;
  ... ;
The triangle begins as:
  1
  1, 1;
  1, 2,  2;
  1, 3,  6,  12;
  1, 4, 12,  36,  108;
  1, 5, 20,  80,  320,  1280;
  1, 6, 30, 150,  750,  3750,  18750;
  1, 7, 42, 252, 1512,  9072,  54432, 326592;
  1, 8, 56, 392, 2744, 19208, 134456, 941192, 6588344;
  ...;
T(3,3) = 12 counts the triples (1,2,1), (1,2,3), (1,3,1), (1,3,2), (2,1,2), (2,1,3), (2,3,1), (2,3,2), (3,1,2), (3,1,3), (3,2,1), (3,2,3) out of a total of 3^3 = 27 triples in the CP(3,3).
		

Crossrefs

Array rows n: A170733 (n=14), ..., A170769 (n=50).
Columns k: A000012(n) (k=0), A000027(n) (k=1), A002378(n-1) (k=2), A011379(n-1) (k=3), A179824(n) (k=4), A101362(n-1) (k=5), 2*A168351(n-1) (k=6), 2*A168526(n-1) (k=7), 2*A168635(n-1) (k=8), 2*A168675(n-1) (k=9), 2*A170783(n-1) (k=10), 2*A170793(n-1) (k=11).
Diagonals k: A055897 (k=n), A055541 (k=n-1), A373395 (k=n-2), A379612 (k=n-3).
Sums: (-1)^n*A065440(n) (signed row).

Programs

  • Magma
    A158497:= func< n,k | k le 1 select n^k else n*(n-1)^(k-1) >;
    [A158497(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Mar 18 2025
    
  • Mathematica
    A158497[n_, k_]:= If[n<2 || k==0, 1, n*(n-1)^(k-1)];
    Table[A158497[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Mar 18 2025 *)
  • SageMath
    def A158497(n,k): return n^k if k<2 else n*(n-1)^(k-1)
    print(flatten([[A158497(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Mar 18 2025

Formula

T(n, k) = (n-1)^(k-1) + (n-1)^k = n*A079901(n-1,k-1), k > 0.
Sum_{k=0..n} T(n,k) = (n*(n-1)^n - 2)/(n-2), n > 2.

Extensions

Edited by R. J. Mathar, Mar 31 2009
More terms added by G. C. Greubel, Mar 18 2025
Showing 1-10 of 57 results. Next