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 21-30 of 215 results. Next

A002697 a(n) = n*4^(n-1).

Original entry on oeis.org

0, 1, 8, 48, 256, 1280, 6144, 28672, 131072, 589824, 2621440, 11534336, 50331648, 218103808, 939524096, 4026531840, 17179869184, 73014444032, 309237645312, 1305670057984, 5497558138880, 23089744183296
Offset: 0

Views

Author

Keywords

Comments

Coefficient of x^(2n-2) in Chebyshev polynomial T(2n) is -a(n).
Let M_n be the n X n matrix m_(i,j) = 1 + 2*abs(i-j); then det(M_n) = (-1)^(n-1)*a(n-1). - Benoit Cloitre, May 28 2002
Number of subsequences 00 in all words of length n+1 on the alphabet {0,1,2,3}. Example: a(2)=8 because we have 000,001,002,003,100,200,300 (the other 57=A125145(3) words of length 3 have no subsequences 00). a(n) = Sum_{k=0..n} k*A128235(n+1, k). - Emeric Deutsch, Feb 27 2007
Let P(A) be the power set of an n-element set A. Then a(n) = the sum of the size of the symmetric difference of x and y for every subset {x,y} of P(A). - Ross La Haye, Dec 30 2007 (See the comment from Bernard Schott below.)
Let P(A) be the power set of an n-element set A and B be the Cartesian product of P(A) with itself. Then remove (y,x) from B when (x,y) is in B and x != y and call this R35. Then a(n) = the sum of the size of the symmetric difference of x and y for every (x,y) of R35. [proposed edit of comment just above; by Ross La Haye]
The numbers in this sequence are the Wiener indices of the graphs of n-cubes (Boolean hypercubes). For example, the 3-cube is the graph of the standard cube whose Wiener index is 48. - K.V.Iyer, Feb 26 2009
From Gary W. Adamson, Sep 06 2009: (Start)
Starting (1, 8, 48, ...) = 4th binomial transform of [1, 4, 0, 0, 0, ...].
Equals the sum of terms in 2^n X 2^n semi-magic square arrays in which each row and column is composed of a binomial frequency of terms in the set (1, 3, 5, 7, ...).
The first few such arrays = [1] [1,3; 3,1]; /Q.
[1, 3, 5, 3;
3, 1, 3, 5;
5, 3, 1, 3;
3, 5, 3, 1]
(sum of terms = 48, with a binomial frequency of (1, 2, 1) as to (1, 3, 5) in each row and column)
[1, 3, 5, 3, 5, 7, 5, 3;
3, 1, 3, 5, 7, 5, 3, 5;
5, 3, 1, 3, 5, 3, 5, 7;
3, 5, 3, 1, 3, 5, 7, 5;
5, 7, 5, 3, 1, 3, 5, 3;
7, 5, 3, 5, 3, 1, 3, 5;
5, 3, 5, 7, 5, 3, 1, 3;
3, 5, 7, 5, 3, 5, 3, 1]
(sum of terms = 256, with each row and column composed of one 1, three 3's, three 5's, and one 7)
... (End)
Let P(A) be the power set of an n-element set A and B be the Cartesian product of P(A) with itself. Then a(n) = the sum of the size of the intersection of x and y for every (x,y) of B. - Ross La Haye, Jan 05 2013
Following the last comment of Ross, A002699 is the similar sequence when "intersection" is replaced by "symmetric difference" and A212698 is the similar sequence when "intersection" is replaced by "union". - Bernard Schott, Jan 04 2013
Also, following the first comment of Ross, A082134 is the similar sequence when "symmetric difference" is replaced by "intersection" and A133224 is the similar sequence when "symmetric difference" is replaced by "union". - Bernard Schott, Jan 15 2013
Let [n] denote the set {1,2,3,...,n} and denote an n-permutation of the elements of [n] by p = p(1)p(2)p(3)...p(n), where p(i) is the i-th entry in the linear order given by p. Then (p(i),p(j)) is an inversion of p if i < j but p(i) > p(j). Denote the number of inversions of p by inv(p) and call a 2n-permutation p = p(1)p(2)...p(2n) 2-ordered if p(1) < p(3) < ... < p(2n-1) and p(2) < p(4) < ... < p(2n). Then Sum(inv(p)) = n*4^(n-1), where the sum is taken over all 2-ordered 2n-permutations of p. See Bona reference below. - Ross La Haye, Jan 21 2014
Sum over all peaks of Dyck paths of semilength n of the product of the x and y coordinates. - Alois P. Heinz, May 29 2015
Sum of the number of all edges over all j-dimensional subcubes of the boolean hypercube graph of dimension n, Q_n, for all j, so a(n) = Sum_{j=1..n} binomial(n,j)*2^(n-j) * j*2^(j-1). - Constantinos Kourouzides, Mar 24 2024

Examples

			From _Bernard Schott_, Jan 04 2013: (Start)
See the comment about intersection of X and Y.
If A={b,c}, then in P(A) we have:
{b}Inter{b}={b},
{b}Inter{b,c}={b},
{c}Inter{c}={c},
{c}Inter{b,c}={c},
{b,c}Inter{b}={b},
{b,c}Inter{c}={c},
{b,c}Inter{b,c}={b,c}
and : #{b}+ #{b}+ #{c}+ #{c}+ #{b}+ #{c}+ #{b,c} = 8 = 2*4^(2-1) = a(2).
The other intersections are empty.
(End)
		

References

  • Miklos Bona, Combinatorics of Permutations, Chapman and Hall/CRC, 2004, pp. 1, 43, 64.
  • C. Lanczos, Applied Analysis. Prentice-Hall, Englewood Cliffs, NJ, 1956, p. 516.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

Formula

a(n) = n*4^(n-1).
G.f.: x/(1-4x)^2. a(n+1) is the convolution of powers of 4 (A000302). - Wolfdieter Lang, May 16 2003
Third binomial transform of n. E.g.f.: x*exp(4x). - Paul Barry, Jul 22 2003
a(n) = Sum_{k=0..n} k*binomial(2*n, 2*k). - Benoit Cloitre, Jul 30 2003
For n>=0, a(n+1) = Sum_{i+j+k+l=n} binomial(2i, i)*binomial(2j, j)*binomial(2k, k)*binomial(2l, l). - Philippe Deléham, Jan 22 2004
a(n) = Sum_{k=0..n} 4^(n-k)*binomial(n-k+1, k)*binomial(1, (k+1)/2)*(1-(-1)^k)/2. - Paul Barry, Oct 15 2004
Sum_{n>0} 1/a(n) = 8*log(2) - 4*log(3). - Jaume Oliver Lafont, Sep 11 2009
a(0) = 0, a(n) = 4*a(n-1) + 4^(n-1). - Vincenzo Librandi, Dec 31 2010
a(n+1) is the convolution of A000984 with A002457. - Rui Duarte, Oct 08 2011
a(0) = 0, a(1) = 1, a(n) = 8*a(n-1) - 16*a(n-2). - Harvey P. Dale, Jan 18 2012
a(n) = A002699(n)/2 = A212698(n)/3. - Bernard Schott, Jan 04 2013
G.f.: W(0)*x/2 , where W(k) = 1 + 1/( 1 - 4*x*(k+2)/( 4*x*(k+2) + (k+1)/W(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 19 2013
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(5/4). - Amiram Eldar, Oct 28 2020
a(n) = (1/2)*Sum_{k=0..n} k*binomial(2*n, k). Compare this with the formula of Benoit Cloitre above. - Wolfdieter Lang, Nov 12 2021
a(n) = (-1)^(n-1)*det(M(n)) for n > 0, where M(n) is the n X n symmetric Toeplitz matrix whose first row consists of 1, 3, ..., 2*n-1. - Stefano Spezia, Aug 04 2022

A124734 Table with all compositions sorted first by total, then by length and finally lexicographically.

Original entry on oeis.org

1, 2, 1, 1, 3, 1, 2, 2, 1, 1, 1, 1, 4, 1, 3, 2, 2, 3, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 5, 1, 4, 2, 3, 3, 2, 4, 1, 1, 1, 3, 1, 2, 2, 1, 3, 1, 2, 1, 2, 2, 2, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 6, 1, 5, 2, 4, 3, 3, 4, 2, 5, 1, 1, 1, 4, 1, 2, 3, 1, 3, 2, 1, 4, 1, 2, 1
Offset: 1

Views

Author

Keywords

Comments

This is similar to the Abramowitz and Stegun ordering for partitions (see A036036). The standard ordering for compositions is A066099, which is more similar to the Mathematica partition ordering (A080577).
This can be regarded as a table in two ways: with each composition as a row, or with the compositions of each integer as a row. The first way has A124736 as row lengths and A070939 as row sums; the second has A001792 as row lengths and A001788 as row sums.
This sequence includes every finite sequence of positive integers.

Examples

			The table starts:
1
2; 1 1
3; 1 2; 2 1; 1 1 1
4; 1 3; 2 2; 3 1; 1 1 2; 1 2 1; 2 1 1; 1 1 1 1;
		

Crossrefs

Programs

  • Mathematica
    Table[Sort@Flatten[Permutations /@ IntegerPartitions@n, 1], {n, 8}] // Flatten (* Robert Price, Jun 13 2020 *)

A008949 Triangle read by rows of partial sums of binomial coefficients: T(n,k) = Sum_{i=0..k} binomial(n,i) (0 <= k <= n); also dimensions of Reed-Muller codes.

Original entry on oeis.org

1, 1, 2, 1, 3, 4, 1, 4, 7, 8, 1, 5, 11, 15, 16, 1, 6, 16, 26, 31, 32, 1, 7, 22, 42, 57, 63, 64, 1, 8, 29, 64, 99, 120, 127, 128, 1, 9, 37, 93, 163, 219, 247, 255, 256, 1, 10, 46, 130, 256, 382, 466, 502, 511, 512, 1, 11, 56, 176, 386, 638, 848, 968, 1013, 1023, 1024, 1, 12, 67, 232, 562, 1024, 1486, 1816, 1981, 2036, 2047, 2048
Offset: 0

Views

Author

Keywords

Comments

The second-left-from-middle column is A000346: T(2n+2, n) = A000346(n). - Ed Catmur (ed(AT)catmur.co.uk), Dec 09 2006
T(n,k) is the maximal number of regions into which n hyperplanes of co-dimension 1 divide R^k (the Cake-Without-Icing numbers). - Rob Johnson, Jul 27 2008
T(n,k) gives the number of vertices within distance k (measured along the edges) of an n-dimensional unit cube, (i.e., the number of vertices on the hypercube graph Q_n whose distance from a reference vertex is <= k). - Robert Munafo, Oct 26 2010
A triangle formed like Pascal's triangle, but with 2^n for n >= 0 on the right border instead of 1. - Boris Putievskiy, Aug 18 2013
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 04 2013
Consider each "1" as an apex of two sequences: the first is the set of terms in the same row as the "1", but the rightmost term in the row repeats infinitely. Example: the row (1, 4, 7, 8) becomes (1, 4, 7, 8, 8, 8, ...). The second sequence begins with the same "1" but is the diagonal going down and to the right, thus: (1, 5, 16, 42, 99, 219, 466, ...). It appears that for all such sequence pairs, the binomial transform of the first, (1, 4, 7, 8, 8, 8, ...) in this case; is equal to the second: (1, 5, 16, 42, 99, ...). - Gary W. Adamson, Aug 19 2015
Let T* be the infinite tree with root 0 generated by these rules: if p is in T*, then p+1 is in T* and x*p is in T*. Let q(n) be the sum of polynomials in the n-th generation of T*. For n >= 0, row n of A008949 gives the coefficients of q(n+1); e.g., (row 3) = (1, 4, 7, 8) matches x^3 + 4*x^2 + 7*x + 9, which is the sum of the 8 polynomials in the 4th generation of T*. - Clark Kimberling, Jun 16 2016
T(n,k) is the number of subsets of [n]={1,...,n} of at most size k. Equivalently, T(n,k) is the number of subsets of [n] of at least size n-k. Counting the subsets of at least size (n-k) by conditioning on the largest element m of the smallest (n-k) elements of such a subset provides the formula T(n,k) = Sum_{m=n-k..n} C(m-1,n-k-1)*2^(n-m), and, by letting j=m-n+k, we obtain T(n,k) = Sum_{j=0..k} C(n+j-k-1,j)*2^(k-j). - Dennis P. Walsh, Sep 25 2017
If the interval of integers 1..n is shifted up or down by k, making the new interval 1+k..n+k or 1-k..n-k, then T(n-1,n-1-k) (= 2^(n-1)-T(n-1,k-1)) is the number of subsets of the new interval that contain their own cardinal number as an element. - David Pasino, Nov 01 2018

Examples

			Triangle begins:
  1;
  1,  2;
  1,  3,  4;
  1,  4,  7,   8;
  1,  5, 11,  15,  16;
  1,  6, 16,  26,  31,  32;
  1,  7, 22,  42,  57,  63,  64;
  1,  8, 29,  64,  99, 120, 127, 128;
  1,  9, 37,  93, 163, 219, 247, 255,  256;
  1, 10, 46, 130, 256, 382, 466, 502,  511,  512;
  1, 11, 56, 176, 386, 638, 848, 968, 1013, 1023, 1024;
  ...
		

References

  • F. J. MacWilliams and N. J. A. Sloane, The Theory of Error-Correcting Codes, Elsevier-North Holland, 1978, p. 376.

Crossrefs

Row sums sequence is A001792.
T(n, m)= A055248(n, n-m).

Programs

  • GAP
    T:=Flat(List([0..11],n->List([0..n],k->Sum([0..k],j->Binomial(n+j-k-1,j)*2^(k-j))))); # Muniru A Asiru, Nov 25 2018
    
  • Haskell
    a008949 n k = a008949_tabl !! n !! k
    a008949_row n = a008949_tabl !! n
    a008949_tabl = map (scanl1 (+)) a007318_tabl
    -- Reinhard Zumkeller, Nov 23 2012
    
  • Magma
    [[(&+[Binomial(n,j): j in [0..k]]): k in [0..n]]: n in [0..12]]; // G. C. Greubel, Nov 25 2018
    
  • Maple
    A008949 := proc(n,k) local i; add(binomial(n,i),i=0..k) end; # Typo corrected by R. J. Mathar, Oct 26 2010
  • Mathematica
    Table[Length[Select[Subsets[n], (Length[ # ] <= k) &]], {n, 0, 12}, {k, 0, n}] // Grid (* Geoffrey Critzer, May 13 2009 *)
    Flatten[Accumulate/@Table[Binomial[n,i],{n,0,20},{i,0,n}]] (* Harvey P. Dale, Aug 08 2015 *)
    T[ n_, k_] := If[ n < 0 || k > n, 0, Binomial[n, k] Hypergeometric2F1[1, -k, n + 1 - k, -1]]; (* Michael Somos, Aug 05 2017 *)
  • PARI
    A008949(n)=T8949(t=sqrtint(2*n-sqrtint(2*n)),n-t*(t+1)/2)
    T8949(r,c)={ 2*c > r || return(sum(i=0,c,binomial(r,i))); 1<M. F. Hasler, May 30 2010
    
  • PARI
    {T(n, k) = if(k>n, 0, sum(i=0, k, binomial(n, i)))}; /* Michael Somos, Aug 05 2017 */
    
  • PARI
    row(n) = my(v=vector(n+1, k, binomial(n,k-1))); vector(#v, k, sum(i=1, k, v[i])); \\ Michel Marcus, Apr 13 2025
    
  • Sage
    [[sum(binomial(n,j) for j in range(k+1)) for k in range(n+1)] for n in range(12)] # G. C. Greubel, Nov 25 2018

Formula

From partial sums across rows of Pascal triangle A007318.
T(n, 0) = 1, T(n, n) = 2^n, T(n, k) = T(n-1, k-1) + T(n-1, k), 0 < k < n.
G.f.: (1 - x*y)/((1 - y - x*y)*(1 - 2*x*y)). - Antonio Gonzalez (gonfer00(AT)gmail.com), Sep 08 2009
T(2n,n) = A032443(n). - Philippe Deléham, Sep 16 2009
T(n,k) = 2 T(n-1,k-1) + binomial(n-1,k) = 2 T(n-1,k) - binomial(n-1,k). - M. F. Hasler, May 30 2010
T(n,k) = binomial(n,n-k)* 2F1(1, -k; n+1-k; -1). - Olivier Gérard, Aug 02 2012
For a closed-form formula for arbitrary left and right borders of Pascal like triangle see A228196. - Boris Putievskiy, Aug 18 2013
T(n,floor(n/2)) = A027306(n). - Reinhard Zumkeller, Nov 14 2014
T(n,n) = 2^n, otherwise for 0 <= k <= n-1, T(n,k) = 2^n - T(n,n-k-1). - Bob Selcoe, Mar 30 2017
For fixed j >= 0, lim_{n -> oo} T(n+1,n-j+1)/T(n,n-j) = 2. - Bob Selcoe, Apr 03 2017
T(n,k) = Sum_{j=0..k} C(n+j-k-1,j)*2^(k-j). - Dennis P. Walsh, Sep 25 2017

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Mar 23 2000

A083710 Number of integer partitions of n with a part dividing all the other parts.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 11, 12, 20, 25, 37, 43, 70, 78, 114, 143, 196, 232, 330, 386, 530, 641, 836, 1003, 1340, 1581, 2037, 2461, 3127, 3719, 4746, 5605, 7038, 8394, 10376, 12327, 15272, 17978, 22024, 26095, 31730, 37339, 45333, 53175, 64100, 75340, 90138
Offset: 0

Views

Author

N. J. A. Sloane, Jun 16 2003

Keywords

Comments

Since the summand (part) which divides all the other summands is necessarily the smallest, an equivalent definition is: "Number of partitions of n such that smallest part divides every part." - Joerg Arndt, Jun 08 2009
The first few partitions that fail the criterion are 5=3+2, 7=5+2=4+3=3+2+2. So a(5) = A000041(5) - 1 = 6, a(7) = A000041(7) - 3 = 12. - Vladeta Jovovic, Jun 17 2003
Starting with offset 1 = inverse Mobius transform (A051731) of the partition numbers, A000041. - Gary W. Adamson, Jun 08 2009

Examples

			From _Gus Wiseman_, Apr 18 2021: (Start)
The a(1) = 1 through a(7) = 12 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (21)   (22)    (41)     (33)      (61)
             (111)  (31)    (221)    (42)      (331)
                    (211)   (311)    (51)      (421)
                    (1111)  (2111)   (222)     (511)
                            (11111)  (321)     (2221)
                                     (411)     (3211)
                                     (2211)    (4111)
                                     (3111)    (22111)
                                     (21111)   (31111)
                                     (111111)  (211111)
                                               (1111111)
(End)
		

References

  • L. M. Chawla, M. O. Levan and J. E. Maxfield, On a restricted partition function and its tables, J. Natur. Sci. and Math., 12 (1972), 95-101.

Crossrefs

Cf. A000041, A051731. - Gary W. Adamson, Jun 08 2009
The case with no 1's is A083711.
The strict case is A097986.
The version for "divisible by" instead of "dividing" is A130689.
The case where there is also a part divisible by all the others is A130714.
The complement of these partitions is counted by A338470.
The Heinz numbers of these partitions are dense, complement of A342193.
The case where there is also no part divisible by all the others is A343345.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A018818 counts partitions into divisors (strict: A033630).
A167865 counts strict chains of divisors > 1 summing to n.

Programs

  • Maple
    with(combinat): with(numtheory): a := proc(n) c := 0: l := sort(convert(divisors(n), list)): for i from 1 to nops(l)-0 do c := c+numbpart(l[i]-1) od: RETURN(c): end: for j from 0 to 60 do printf(`%d, `, a(j)) od: # Zerinvary Lajos, Apr 14 2007
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And@@IntegerQ/@(#/Min@@#)&]],{n,0,30}] (* Gus Wiseman, Apr 18 2021 *)

Formula

Equals left border of triangle A137587 starting (1, 2, 3, 5, 6, 11, ...). - Gary W. Adamson, Jan 27 2008
G.f.: 1 + Sum_{n>=1} x^n/eta(x^n). The g.f. for partitions into parts that are a multiple of n is x^n/eta(x^n), now sum over n. - Joerg Arndt, Jun 08 2009
Gary W. Adamson's comment is equivalent to the formula a(n) = Sum_{d|n} p(d-1) where p(i) = number of partitions of i (A000041(i)). Hence A083710 has g.f. Sum_{d>=1} p(d-1)*x^d/(1-x^d), - N. J. A. Sloane, Jun 08 2009

Extensions

More terms from Vladeta Jovovic, Jun 17 2003
Name shortened by Gus Wiseman, Apr 18 2021

A228369 Triangle read by rows in which row n lists the compositions (ordered partitions) of n in lexicographic order.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 1, 2, 2, 1, 3, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 3, 2, 1, 1, 2, 2, 3, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 1, 1, 2, 2, 1, 3, 1, 1, 4, 2, 1, 1, 1, 2, 1, 2, 2, 2, 1, 2, 3, 3, 1, 1, 3, 2, 4, 1, 5, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Omar E. Pol, Aug 28 2013

Keywords

Comments

The representation of the compositions (for fixed n) is as lists of parts, the order between individual compositions (for the same n) is lexicographic. - Joerg Arndt, Sep 02 2013
The equivalent sequence for partitions is A026791.
Row n has length A001792(n-1).
Row sums give A001787, n >= 1.
The m-th composition has length A008687(m+1), m >= 1. - Andrey Zabolotskiy, Jul 19 2017

Examples

			Illustration of initial terms:
-----------------------------------
n  j       Diagram   Composition j
-----------------------------------
.               _
1  1           |_|   1;
.             _ _
2  1         | |_|   1, 1,
2  2         |_ _|   2;
.           _ _ _
3  1       | | |_|   1, 1, 1,
3  2       | |_ _|   1, 2,
3  3       |   |_|   2, 1,
3  4       |_ _ _|   3;
.         _ _ _ _
4  1     | | | |_|   1, 1, 1, 1,
4  2     | | |_ _|   1, 1, 2,
4  3     | |   |_|   1, 2, 1,
4  4     | |_ _ _|   1, 3,
4  5     |   | |_|   2, 1, 1,
4  6     |   |_ _|   2, 2,
4  7     |     |_|   3, 1,
4  8     |_ _ _ _|   4;
.
Triangle begins:
[1];
[1,1],[2];
[1,1,1],[1,2],[2,1],[3];
[1,1,1,1],[1,1,2],[1,2,1],[1,3],[2,1,1],[2,2],[3,1],[4];
[1,1,1,1,1],[1,1,1,2],[1,1,2,1],[1,1,3],[1,2,1,1],[1,2,2],[1,3,1],[1,4],[2,1,1,1],[2,1,2],[2,2,1],[2,3],[3,1,1],[3,2],[4,1],[5];
...
		

Crossrefs

Programs

  • Haskell
    a228369 n = a228369_list !! (n - 1)
    a228369_list = concatMap a228369_row [1..]
    a228369_row 0 = []
    a228369_row n
      | 2^k == 2 * n + 2 = [k - 1]
      | otherwise        = a228369_row (n `div` 2^k) ++ [k] where
        k = a007814 (n + 1) + 1
    -- Peter Kagey, Jun 27 2016
    
  • Mathematica
    Table[Sort[Join@@Permutations/@IntegerPartitions[n],OrderedQ[PadRight[{#1,#2}]]&],{n,5}] (* Gus Wiseman, Dec 14 2017 *)
  • PARI
    gen_comp(n)=
    {  /* Generate compositions of n as lists of parts (order is lex): */
        my(ct = 0);
        my(m, z, pt);
        \\ init:
        my( a = vector(n, j, 1) );
        m = n;
        while ( 1,
            ct += 1;
            pt = vector(m, j, a[j]);
            /* for A228369  print composition: */
            for (j=1, m, print1(pt[j],", ") );
    \\        /* for A228525 print reversed (order is colex): */
    \\        forstep (j=m, 1, -1, print1(pt[j],", ") );
            if ( m<=1,  return(ct) );  \\ current is last
            a[m-1] += 1;
            z = a[m] - 2;
            a[m] = 1;
            m += z;
        );
        return(ct);
    }
    for(n=1, 12, gen_comp(n) );
    \\ Joerg Arndt, Sep 02 2013
    
  • Python
    a = [[[]], [[1]]]
    for s in range(2, 9):
        a.append([])
        for k in range(1, s+1):
            for ss in a[s-k]:
                a[-1].append([k]+ss)
    print(a)
    # Andrey Zabolotskiy, Jul 19 2017

A006234 a(n) = n*3^(n-4).

Original entry on oeis.org

1, 4, 15, 54, 189, 648, 2187, 7290, 24057, 78732, 255879, 826686, 2657205, 8503056, 27103491, 86093442, 272629233, 860934420, 2711943423, 8523250758, 26732013741, 83682825624, 261508830075, 815907549834, 2541865828329
Offset: 3

Views

Author

Keywords

Comments

For n >= 1 a(n) is also the determinant of the n-3 X n-3 matrix with 4's on the diagonal and 1's elsewhere. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 06 2001
a(n+3) = det(M(n)) where M(n) is the n X n matrix with m(i,i) = 4, m(i,j) = i/j for i != j. - Benoit Cloitre, Feb 01 2003
Main diagonal of array defined by m(1,j) = j; m(i,1) = i and m(i,j) = m(i-1,j) + 2*m(i-1,j-1). - Benoit Cloitre, Jun 13 2003
a(n+3) is the number of words of length n on {A, B, C, D} with no D appearing anywhere to the right of an A. - Rob Pratt, Aug 04 2004
Number of spanning trees in the book graph of order n-2, i.e., S_{n-2} X P_2 (S_k = the star graph on k nodes) (conjectured). This conjecture is true - see Doslic (2013). - N. J. A. Sloane, Dec 28 2013
Conjecture: a(n+2) is the total number of parts used in the compositions of n if the parts can be runs of any length from 1 to n, and contain any integers from 1 to n. (The number of such compositions is given by A000244(n-1).) - Gregory L. Simay, May 27 2017
a(n+3) is the number of words of length n defined on 4 letters where one of the letters is used at most once. - Enrique Navarrete, Mar 14 2024

Examples

			For n=3, the total number of parts is (3+2)3^(3+2-4)=(5)(3)=15 (each part indicated by "[]"): [3]; [2,1]; [1,2]; [2],[1]; [1],[2]; [1,1,1]; [1,1],[1]; [1],[1,1]; [1],[1],[1]. Note that these 15 parts are arranged into 9 = A000244(3-1)compositions. - _Gregory L. Simay_, May 27 2017
		

References

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

Crossrefs

Binomial transform of A001792.

Programs

Formula

G.f.: (1-2*x)/(1-3*x)^2. - Simon Plouffe in his 1992 dissertation.
a(n+3) = Sum_{k=0..n} A112626(n, k). - Ross La Haye, Jan 11 2006
G.f.: Hypergeometric2F1([1,4],[3],3*x). - R. J. Mathar, Aug 09 2015
From Amiram Eldar, Jan 18 2021: (Start)
Sum_{n>=1} 1/a(n) = 81*log(3/2).
Sum_{n>=1} (-1)^(n+1)/a(n) = 81*log(4/3). (End)
E.g.f.: x*(exp(3*x) - 3*x - 1)/27. - Stefano Spezia, Mar 04 2023
E.g.f. (with offset 0): exp(3*x)*(1+x). - Enrique Navarrete, Mar 14 2024

A228525 Triangle read by rows in which row n lists the compositions (ordered partitions) of n in colexicographic order.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 2, 2, 2, 1, 3, 4, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 2, 1, 2, 2, 1, 1, 3, 1, 4, 1, 1, 1, 1, 2, 2, 1, 2, 1, 2, 2, 3, 2, 1, 1, 3, 2, 3, 1, 4, 5, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Omar E. Pol, Aug 24 2013

Keywords

Comments

The representation of the compositions (for fixed n) is as lists of parts, the order between individual compositions (for the same n) is co-lexicographic. [Joerg Arndt, Sep 02 2013]
The equivalent sequence for partitions is A211992.
Row n has length A001792(n-1).
Row sums give A001787, n >= 1.

Examples

			Illustration of initial terms:
---------------------------------
n  j     Diagram     Composition
---------------------------------
.         _
1  1     |_|         1;
.         _ _
2  1     |_| |       1, 1,
2  2     |_ _|       2;
.         _ _ _
3  1     |_| | |     1, 1, 1,
3  2     |_ _| |     2, 1,
3  3     |_|   |     1, 2,
3  4     |_ _ _|     3;
.         _ _ _ _
4  1     |_| | | |   1, 1, 1, 1,
4  2     |_ _| | |   2, 1, 1,
4  3     |_|   | |   1, 2, 1,
4  4     |_ _ _| |   3, 1,
4  5     |_| |   |   1, 1, 2,
4  6     |_ _|   |   2, 2,
4  7     |_|     |   1, 3,
4  8     |_ _ _ _|   4;
.
Triangle begins:
[1];
[1,1],[2];
[1,1,1],[2,1],[1,2],[3];
[1,1,1,1],[2,1,1],[1,2,1],[3,1],[1,1,2],[2,2],[1,3],[4];
[1,1,1,1,1],[2,1,1,1],[1,2,1,1],[3,1,1],[1,1,2,1],[2,2,1],[1,3,1],[4,1],[1,1,1,2],[2,1,2],[1,2,2],[3,2],[1,1,3],[2,3],[1,4],[5];
		

Crossrefs

Programs

  • PARI
    gen_comp(n)=
    {  /* Generate compositions of n as lists of parts (order is lex): */
        my(ct = 0);
        my(m, z, pt);
        \\ init:
        my( a = vector(n, j, 1) );
        m = n;
        while ( 1,
            ct += 1;
            pt = vector(m, j, a[j]);
    \\        /* for A228369  print composition: */
    \\        for (j=1, m, print1(pt[j],", ") );
            /* for A228525 print reversed (order is colex): */
            forstep (j=m, 1, -1, print1(pt[j],", ") );
            if ( m<=1,  return(ct) );  \\ current is last
            a[m-1] += 1;
            z = a[m] - 2;
            a[m] = 1;
            m += z;
        );
        return(ct);
    }
    for(n=1, 12, gen_comp(n) );
    \\ Joerg Arndt, Sep 02 2013

A045891 First differences of A045623.

Original entry on oeis.org

1, 1, 3, 7, 16, 36, 80, 176, 384, 832, 1792, 3840, 8192, 17408, 36864, 77824, 163840, 344064, 720896, 1507328, 3145728, 6553600, 13631488, 28311552, 58720256, 121634816, 251658240, 520093696, 1073741824, 2214592512, 4563402752
Offset: 0

Views

Author

Keywords

Comments

Let M_n be the n X n matrix m_(i,j) = 3 + abs(i-j), then det(M_n) =(-1)^(n+1)*a(n+1). - Benoit Cloitre, May 28 2002
If X_1, X_2, ..., X_n are 2-blocks of a (2n+3)-set X then, for n>=1, a(n+2) is the number of (n+1)-subsets of X intersecting each X_i, (i=1..n). - Milan Janjic, Nov 18 2007
Equals row sums of triangle A152194. - Gary W. Adamson, Nov 28 2008
An elephant sequence, see A175655. For the central square 16 A[5] vectors, with decimal values between 19 and 400, lead to this sequence (without the first leading 1). For the corner squares these vectors lead to the companion sequence A045623. - Johannes W. Meijer, Aug 15 2010
a(n) is the total number of runs of 1 in the compositions of n+1. For example, a(3) = A045623(3) - A045623(2) = 12 - 5 = 7 runs of only 1 in the compositions of 4, enumerated "()" as follows: 3,(1); (1),3; 2,(1,1);(1),2,(1); (1,1),2; (1,1,1,1). More generally, the total number of runs of only part k in the compositions of n+k is A045623(n) - A045623(n-k). - Gregory L. Simay, May 02 2017
This is essentially the p-INVERT of (1,1,1,1,1,...) for p(S) = 1 - S - S^2 + S^3; see A291000. - Clark Kimberling, Aug 24 2017

Examples

			G.f. = 1 + x + 3*x^2 + 7*x^3 + 16*x^4 + 36*x^5 + 80*x^6 + ... - _Michael Somos_, Mar 26 2022
		

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n-2} (k+3)*binomial(n-2,k) for n >= 2. - N. J. A. Sloane, Jan 30 2008
a(n) = (n+4)*2^(n-3), n >= 2, with a(0) = a(1) = 1.
G.f.: (1-x)^3/(1-2*x)^2.
Equals binomial transform of A027656.
Starting 1, 3, 7, 16, ... this is ((n+5)*2^n - 0^n)/4, the binomial transform of (1, 2, 2, 3, 3, ...). - Paul Barry, May 20 2003
From Paul Barry, Nov 29 2004: (Start)
a(n) = ((n+4)*2^(n-1) + 3*C(0, n) - C(1, n))/4;
a(n) = Sum_{k=0..floor(n/2)} C(n, 2*k)*(k+1). (End)
a(n) = A045623(n-1) + 2^(n-2) = A034007(n+1) - 2^(n-2) for n>=2. - Philippe Deléham, Apr 20 2009
G.f.: 1 + Q(0)*x/(1-x)^2, where Q(k)= 1 + (k+1)*x/(1 - x - x*(1-x)/(x + (k+1)*(1-x)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 25 2013
a(n) = Sum_{k=0..n} (k+1)*C(n-2,n-k). Peter Luschny, Apr 20 2015
From Amiram Eldar, Jan 13 2021: (Start)
Sum_{n>=0} 1/a(n) = 128*log(2) - 1292/15.
Sum_{n>=0} (-1)^n/a(n) = 782/15 - 128*log(3/2). (End)
E.g.f.: (2 - x + exp(2*x)*(2 + x))/4. - Stefano Spezia, Mar 26 2022

A055248 Triangle of partial row sums of triangle A007318(n,m) (Pascal's triangle). Triangle A008949 read backwards. Riordan (1/(1-2x), x/(1-x)).

Original entry on oeis.org

1, 2, 1, 4, 3, 1, 8, 7, 4, 1, 16, 15, 11, 5, 1, 32, 31, 26, 16, 6, 1, 64, 63, 57, 42, 22, 7, 1, 128, 127, 120, 99, 64, 29, 8, 1, 256, 255, 247, 219, 163, 93, 37, 9, 1, 512, 511, 502, 466, 382, 256, 130, 46, 10, 1, 1024, 1023, 1013, 968, 848, 638, 386, 176, 56, 11, 1
Offset: 0

Views

Author

Wolfdieter Lang, May 26 2000

Keywords

Comments

In the language of the Shapiro et al. reference (also given in A053121) such a lower triangular (ordinary) convolution array, considered as matrix, belongs to the Riordan-group. The g.f. for the row polynomials p(n,x) (increasing powers of x) is 1/((1-2*z)*(1-x*z/(1-z))).
Binomial transform of the all 1's triangle: as a Riordan array, it factors to give (1/(1-x),x/(1-x))(1/(1-x),x). Viewed as a number square read by antidiagonals, it has T(n,k) = Sum_{j=0..n} binomial(n+k,n-j) and is then the binomial transform of the Whitney square A004070. - Paul Barry, Feb 03 2005
Riordan array (1/(1-2x), x/(1-x)). Antidiagonal sums are A027934(n+1), n >= 0. - Paul Barry, Jan 30 2005; edited by Wolfdieter Lang, Jan 09 2015
Eigensequence of the triangle = A005493: (1, 3, 10, 37, 151, 674, ...); row sums of triangles A011971 and A159573. - Gary W. Adamson, Apr 16 2009
Read as a square array, this is the generalized Riordan array ( 1/(1 - 2*x), 1/(1 - x) ) as defined in the Bala link (p. 5), which factorizes as ( 1/(1 - x), x/(1 - x) )*( 1/(1 - x), x )*( 1, 1 + x ) = P*U*transpose(P), where P denotes Pascal's triangle, A007318, and U is the lower unit triangular array with 1's on or below the main diagonal. - Peter Bala, Jan 13 2016

Examples

			The triangle a(n,m) begins:
n\m    0    1    2   3   4   5   6   7  8  9 10 ...
0:     1
1:     2    1
2:     4    3    1
3:     8    7    4   1
4:    16   15   11   5   1
5:    32   31   26  16   6   1
6:    64   63   57  42  22   7   1
7:   128  127  120  99  64  29   8   1
8:   256  255  247 219 163  93  37   9  1
9:   512  511  502 466 382 256 130  46 10  1
10: 1024 1023 1013 968 848 638 386 176 56 11  1
... Reformatted. - _Wolfdieter Lang_, Jan 09 2015
Fourth row polynomial (n=3): p(3,x)= 8 + 7*x + 4*x^2 + x^3.
The matrix inverse starts
   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;
   ...
which may be related to A029653. - _R. J. Mathar_, Mar 29 2013
From _Peter Bala_, Dec 23 2014: (Start)
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     |
|4 3 1   ||0 2 1     ||0 0 1    |... = |4  5 1   |
|8 7 4 1 ||0 4 3 1   ||0 0 2 1  |      |8 19 9 1 |
|...     ||0 8 7 4 1 ||0 0 4 3 1|      |...      |
|...     ||...       ||...      |      |         |
= A143494. (End)
Matrix factorization of square array as P*U*transpose(P):
/1      \ /1        \ /1 1 1 1 ...\    /1  1  1  1 ...\
|1 1     ||1 1       ||0 1 2 3 ... |   |2  3  4  5 ... |
|1 2 1   ||1 1 1     ||0 0 1 3 ... | = |4  7 11 16 ... |
|1 3 3 1 ||1 1 1 1   ||0 0 0 1 ... |   |8 15 26 42 ... |
|...     ||...       ||...         |   |...            |
- _Peter Bala_, Jan 13 2016
		

Crossrefs

Column sequences: A000079 (powers of 2, m=0), A000225 (m=1), A000295 (m=2), A002662 (m=3), A002663 (m=4), A002664 (m=5), A035038 (m=6), A035039 (m=7), A035040 (m=8), A035041 (m=9), A035042 (m=10).
Row sums: A001792(n) = A055249(n, 0).
Alternating row sums: A011782.
Cf. A011971, A159573. - Gary W. Adamson, Apr 16 2009

Programs

  • Haskell
    a055248 n k = a055248_tabl !! n !! k
    a055248_row n = a055248_tabl !! n
    a055248_tabl = map reverse a008949_tabl
    -- Reinhard Zumkeller, Jun 20 2015
  • Maple
    T := (n,k) -> 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n + 1], [n-k + 2], 1/2).
    seq(seq(simplify(T(n,k)), k=0..n),n=0..10); # Peter Luschny, Oct 10 2019
  • Mathematica
    a[n_, m_] := Sum[ Binomial[n, m + j], {j, 0, n}]; Table[a[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013, after Paul Barry *)
    T[n_, k_] := Binomial[n, k] * Hypergeometric2F1[1, k - n, k + 1, -1];
    Flatten[Table[T[n, k], {n, 0, 7}, {k, 0, n}]]  (* Peter Luschny, Oct 06 2023 *)

Formula

a(n, m) = A008949(n, n-m), if n > m >= 0.
a(n, m) = Sum_{k=m..n} A007318(n, k) (partial row sums in columns m).
Column m recursion: a(n, m) = Sum_{j=m..n-1} a(j, m) + A007318(n, m) if n >= m >= 0, a(n, m) := 0 if n
G.f. for column m: (1/(1-2*x))*(x/(1-x))^m, m >= 0.
a(n, m) = Sum_{j=0..n} binomial(n, m+j). - Paul Barry, Feb 03 2005
Inverse binomial transform (by columns) of A112626. - Ross La Haye, Dec 31 2006
T(2n,n) = A032443(n). - Philippe Deléham, Sep 16 2009
From Peter Bala, Dec 23 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)*(8 + 7*x + 4*x^2/2! + x^3/3!) = 8 + 15*x + 26*x^2/2! + 42*x^3/3! + 64*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the present triangle. 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. The infinite product M(0)*M(1)*M(2)*..., which is clearly well-defined, is equal to A143494 (but with a different offset). See the Example section. Cf. A106516. (End)
a(n,m) = Sum_{p=m..n} 2^(n-p)*binomial(p-1,m-1), n >= m >= 0, else 0. - Wolfdieter Lang, Jan 09 2015
T(n, k) = 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n+1], [n-k+2], 1/2). - Peter Luschny, Oct 10 2019
T(n, k) = binomial(n, k)*hypergeom([1, k - n], [k + 1], -1). - Peter Luschny, Oct 06 2023
n-th row polynomial R(n, x) = (2^n - x*(1 + x)^n)/(1 - x). These polynomials can be used to find series acceleration formulas for the constants log(2) and Pi. - Peter Bala, Mar 03 2025

A059570 Number of fixed points in all 231-avoiding involutions in S_n.

Original entry on oeis.org

1, 2, 6, 14, 34, 78, 178, 398, 882, 1934, 4210, 9102, 19570, 41870, 89202, 189326, 400498, 844686, 1776754, 3728270, 7806066, 16311182, 34020466, 70837134, 147266674, 305718158, 633805938, 1312351118, 2714180722, 5607318414, 11572550770, 23860929422
Offset: 1

Author

Emeric Deutsch, Feb 16 2001

Keywords

Comments

Number of odd parts in all compositions (ordered partitions) of n: a(3)=6 because in 3=2+1=1+2=1+1+1 we have 6 odd parts. Number of even parts in all compositions (ordered partitions) of n+1: a(3)=6 because in 4=3+1=1+3=2+2=2+1+1=1+2+1=1+1+2=1+1+1+1 we have 6 even parts.
Convolved with (1, 2, 2, 2, ...) = A001787: (1, 4, 12, 32, 80, ...). - Gary W. Adamson, May 23 2009
An elephant sequence, see A175654. For the corner squares 36 A[5] vectors, with decimal values between 15 and 480, lead to this sequence. For the central square these vectors lead to the companion sequence 4*A172481, for n>=-1. - Johannes W. Meijer, Aug 15 2010
a(n) is the total number of runs of equal parts in the compositions of n. a(5) = 34 because there are 34 runs of equal parts in the compositions of 5, with parentheses enclosing each run: (5), (4)(1), (1)(4), (3)(2), (2)(3), (3)(1,1), (1)(3)(1), (1,1)(3), (2,2)(1), (2)(1)(2), (1)(2,2), (2)(1,1,1), (1)(2)(1,1), (1,1)(2)(1), (1,1,1)(2), (1,1,1,1,1). - Gregory L. Simay, Apr 28 2017
a(n) - a(n-2) is the number of 1's in all compositions of n and more generally, the number of k's in all compositions of n+k-1. - Gregory L. Simay, May 01 2017

Examples

			a(3) = 6 because in the 231-avoiding involutions of {1,2,3}, i.e., in 123, 132, 213, 321, we have altogether 6 fixed points (3+1+1+1).
		

Programs

  • Magma
    [(3*n+4)*2^n/18-2*(-1)^n/9: n in [1..40]]; // Vincenzo Librandi, May 01 2017
  • Mathematica
    LinearRecurrence[{3,0,-4},{1,2,6},30] (* Harvey P. Dale, Dec 29 2013 *)
    Table[(3 n + 4) 2^n/18 - 2 (-1)^n/9, {n, 30}] (* Vincenzo Librandi, May 01 2017 *)

Formula

a(n) = (3*n+4)*2^n/18 - 2*(-1)^n/9.
G.f.: z*(1-z)/((1+z)*(1-2*z)^2).
a(n) = Sum_{j=0..n} Sum_{k=0..n} binomial(n-k, k+j)*2^k. - Paul Barry, Aug 29 2004
a(n) = Sum_{k=0..n+1} (-1)^(k+1)*binomial(n+1, k+j)*A001045(k). - Paul Barry, Jan 30 2005
Convolution of "Expansion of (1-x)/(1-x-2*x^2)" (A078008) with "Powers of 2" (A000079), treating the result as if offset=1. - Graeme McRae, Jul 12 2006
Convolution of "Difference sequence of A045623" (A045891) with "Positive integers repeated" (A008619), treating the result as if offset=1. - Graeme McRae, Jul 12 2006
a(n) = 3*a(n-1)-4*a(n-3); a(1)=1,a(2)=2,a(3)=6. - Philippe Deléham, Aug 30 2006
Equals row sums of A128255. (1, 2, 6, 14, 34, ...) - (0, 0, 1, 2, 6, 14, 34, ...) = A045623: (1, 2, 5, 12, 28, 64, ...). - Gary W. Adamson, Feb 20 2007
Equals triangle A059260 * [1, 2, 3, ...] as a vector. - Gary W. Adamson, Mar 06 2012
a(n) + a(n-1) = A001792(n-1). - Gregory L. Simay, Apr 30 2017
a(n) - a(n-2) = A045623(n-1). - Gregory L. Simay, May 01 2017
a(n) = A045623(n-1) + A045623(n-3) + A045623(n-5) + ... - Gregory L. Simay, Feb 19 2018
a(n) = A225084(2n,n). - Alois P. Heinz, Aug 30 2018

Extensions

More terms from Eugene McDonnell (eemcd(AT)mac.com), Jan 13 2005
Previous Showing 21-30 of 215 results. Next