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 41-50 of 196 results. Next

A352523 Number of integer compositions of n with exactly k nonfixed points (parts not on the diagonal).

Original entry on oeis.org

1, 1, 0, 0, 2, 0, 1, 1, 2, 0, 0, 4, 2, 2, 0, 0, 5, 5, 4, 2, 0, 1, 3, 12, 8, 6, 2, 0, 0, 7, 14, 19, 14, 8, 2, 0, 0, 8, 21, 33, 32, 22, 10, 2, 0, 0, 9, 30, 54, 63, 54, 32, 12, 2, 0, 1, 6, 47, 80, 116, 116, 86, 44, 14, 2, 0, 0, 11, 53, 129, 194, 229, 202, 130, 58, 16, 2, 0
Offset: 0

Views

Author

Gus Wiseman, Mar 26 2022

Keywords

Comments

A nonfixed point in a composition c is an index i such that c_i != i.

Examples

			Triangle begins:
   1
   1   0
   0   2   0
   1   1   2   0
   0   4   2   2   0
   0   5   5   4   2   0
   1   3  12   8   6   2   0
   0   7  14  19  14   8   2   0
   0   8  21  33  32  22  10   2   0
   0   9  30  54  63  54  32  12   2   0
   1   6  47  80 116 116  86  44  14   2   0
   ...
For example, row n = 6 counts the following compositions (empty column indicated by dot):
  (123)  (6)   (24)    (231)    (2112)   (21111)    .
         (15)  (33)    (312)    (2121)   (111111)
         (42)  (51)    (411)    (3111)
               (114)   (1113)   (11112)
               (132)   (1122)   (11121)
               (141)   (1311)   (11211)
               (213)   (2211)
               (222)   (12111)
               (321)
               (1131)
               (1212)
               (1221)
		

Crossrefs

Column k = 0 is A010054.
Row sums are A011782.
The version for permutations is A098825.
The corresponding rank statistic is A352513.
Column k = 1 is A352520.
A238349 and A238350 count comps by fixed points, first col A238351, rank stat A352512.
A352486 gives the nonfixed points of A122111, counted by A330644.
A352521 counts comps by strong nonexcedances, first A219282, stat A352514.
A352522 counts comps by weak nonexcedances, first col A238874, stat A352515.
A352524 counts comps by strong excedances, first col A008930, stat A352516.
A352525 counts comps by weak excedances, first col A177510, stat A352517.

Programs

  • Maple
    b:= proc(n, i) option remember; expand(`if`(n=0, 1,
          add(`if`(i=j, 1, x)*b(n-j, i+1), j=1..n)))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n, 1)):
    seq(T(n), n=0..12);  # Alois P. Heinz, Mar 19 2025
  • Mathematica
    pnq[y_]:=Length[Select[Range[Length[y]],#!=y[[#]]&]];
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],pnq[#]==k&]],{n,0,9},{k,0,n}]
  • PARI
    T_xy(max_row) = {my(N=max_row+1, x='x+O('x^N), h= sum(i=0, N, prod(j=1, i, y*(x/(1-x)-x^j)+x^j))); vector(N, n, my(r=Vecrev(polcoeff(h, n-1))); if(n<2, r, concat(r,[0])))}
    T_xy(10) \\ John Tyler Rascoe, Mar 21 2025

Formula

G.f.: Sum_{i>=0} Product_{j=1..i} y*(x/(1-x) - x^j) + x^j. - John Tyler Rascoe, Mar 19 2025

A069723 a(n) = 2^(n-1)*binomial(2*n-3, n-1).

Original entry on oeis.org

1, 2, 12, 80, 560, 4032, 29568, 219648, 1647360, 12446720, 94595072, 722362368, 5538111488, 42600857600, 328635187200, 2541445447680, 19696202219520, 152935217233920, 1189496134041600, 9265548833587200, 72271280901980160, 564404288948797440, 4412615349963325440
Offset: 1

Views

Author

Valery A. Liskovets, Apr 07 2002

Keywords

Comments

Number of rooted unicursal planar maps with n edges and two vertices of valency 1 (unicursal means that exactly two vertices are of odd valency; there is an Eulerian path).

Crossrefs

Main diagonal of array A082137.

Programs

  • Maple
    Z:=(1-sqrt(1-z))*8^n/sqrt(1-z)/2: Zser:=series(Z, z=0, 33): seq(coeff(Zser, z, n), n=0..20); # Zerinvary Lajos, Jan 16 2007
  • Mathematica
    Table[2^(n - 1) * Binomial[2*n - 3, n - 1], {n, 1,50}] (* G. C. Greubel, Jan 15 2017 *)
  • Sage
    # Assuming offset 0:
    A069723  = lambda n: (rising_factorial(n, n)/factorial(n)) << n
    [A069723(n) for n in (0..20)] # Peter Luschny, Nov 30 2014

Formula

a(n) = A069722(n)/2, n>1.
G.f.: 4*x/(sqrt(1-8*x) * (1-sqrt(1-8*x))). - Paul Barry, Sep 06 2004
With offset 0: a(n) = (0^n + 2^n*binomial(2n, n))/2. - Paul Barry, Sep 24 2004
D-finite with recurrence (-n+1)*a(n) + 4*(2*n-3)*a(n-1) = 0. - R. J. Mathar, Dec 03 2012
With offset 0: a(n) = 2^n*rf(n,n)/n! = 2^n*A088218(n), where rf denotes the rising factorial. - Peter Luschny, Nov 30 2014
a(n) = Sum_{k=0..n} binomial(n+k-1,k)*binomial(2*n-1, n-k). - Vladimir Kruchinin, Nov 11 2016
a(n) ~ 2^(3*n-4)/sqrt(Pi*n). - Ilya Gutkovskiy, Nov 11 2016
From Amiram Eldar, Jan 16 2024: (Start)
Sum_{n>=1} 1/a(n) = 9/7 + 16*arcsin(1/(2*sqrt(2)))/(7*sqrt(7)).
Sum_{n>=1} (-1)^(n+1)/a(n) = 7/9 - 8*log(2)/27. (End)

A180281 Triangle read by rows: T(n,k) = number of arrangements of n indistinguishable balls in n boxes with the maximum number of balls in any box equal to k.

Original entry on oeis.org

1, 1, 2, 1, 6, 3, 1, 18, 12, 4, 1, 50, 50, 20, 5, 1, 140, 195, 90, 30, 6, 1, 392, 735, 392, 147, 42, 7, 1, 1106, 2716, 1652, 672, 224, 56, 8, 1, 3138, 9912, 6804, 2970, 1080, 324, 72, 9, 1, 8952, 35850, 27600, 12825, 4950, 1650, 450, 90, 10, 1, 25652, 128865, 110715, 54450, 22022, 7865, 2420, 605, 110, 11
Offset: 1

Views

Author

R. H. Hardin, Aug 24 2010

Keywords

Comments

To clarify a slight ambiguity in the definition, the heaviest box in such an arrangement should contain exactly k balls. - Gus Wiseman, Sep 22 2016

Examples

			The T(4,2)=18 arrangements are {0022, 0112, 0121, 0202, 0211, 0220, 1012, 1021, 1102, 1120, 1201, 1210, 2002, 2011, 2020, 2101, 2110, 2200}.
Triangle starts
  1
  1   2
  1   6   3
  1  18  12  4
  1  50  50 20  5
  1 140 195 90 30 6
  ...
		

Crossrefs

Row sums give A088218.
T(n,ceiling(n/2)) gives A318160.
T(2n,n) gives A318161.
T(2n-1,n) gives A318161.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, 1,
          `if`(i=0, 0, add(b(n-j, i-1, k), j=0..min(n, k))))
        end:
    T:= (n, k)-> b(n$2, k)-b(n$2, k-1):
    seq(seq(T(n,k), k=1..n), n=1..12);  # Alois P. Heinz, Aug 16 2018
    # second Maple program:
    T:= (n, k)-> coeff(series(((x^(k+1)-1)/(x-1))^n
                 -((x^k-1)/(x-1))^n, x, n+1), x, n):
    seq(seq(T(n, k), k=1..n), n=1..12);  # Alois P. Heinz, Aug 17 2018
  • Mathematica
    T[n_,k_]:=Select[Tuples[Range[0,k],n],And[Max[#]===k,Total[#]===n]&]; (* Gus Wiseman, Sep 22 2016 *)
    SequenceForm@@@T[4,2] (* example *)
    Join@@Table[Length[T[n,k]],{n,1,6},{k,1,n}] (* sequence *)
    (* Second program: *)
    b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i == 0, 0, Sum[b[n-j, i-1, k], {j, 0, Min[n, k]}]]];
    T[n_, k_] := b[n, n, k] - b[n, n, k-1];
    Table[Table[T[n, k], {k, 1, n}], {n, 1, 12}] // Flatten (* Jean-François Alcover, Aug 28 2022, after Alois P. Heinz *)

Formula

Empirical: right half of table, T(n,k) = n*binomial(2*n-k-2,n-2) for 2*k > n; also, T(n,2) = Sum_{j=1..n} binomial(n,j)*binomial(n-j,j) = 2*A097861(n). - Robert Gerbicz in the Sequence Fans Mailing List
From Alois P. Heinz, Aug 17 2018: (Start)
T(n,k) = [x^n] ((x^(k+1)-1)/(x-1))^n - ((x^k-1)/(x-1))^n.
T(n,k) = A305161(n,k) - A305161(n,k-1). (End)

A345914 Numbers k such that the k-th composition in standard order (row k of A066099) has reverse-alternating sum >= 0.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 19, 20, 21, 22, 24, 26, 27, 28, 30, 31, 32, 35, 36, 37, 38, 40, 41, 42, 43, 44, 46, 47, 48, 50, 51, 52, 53, 54, 55, 56, 58, 59, 60, 61, 62, 63, 64, 67, 69, 70, 72, 73, 74, 76, 79, 80, 82, 83, 84, 86, 87, 88
Offset: 1

Views

Author

Gus Wiseman, Jul 04 2021

Keywords

Comments

The reverse-alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(k-i) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The sequence of terms together with the corresponding compositions begins:
     0: ()           19: (3,1,1)        40: (2,4)
     1: (1)          20: (2,3)          41: (2,3,1)
     2: (2)          21: (2,2,1)        42: (2,2,2)
     3: (1,1)        22: (2,1,2)        43: (2,2,1,1)
     4: (3)          24: (1,4)          44: (2,1,3)
     6: (1,2)        26: (1,2,2)        46: (2,1,1,2)
     7: (1,1,1)      27: (1,2,1,1)      47: (2,1,1,1,1)
     8: (4)          28: (1,1,3)        48: (1,5)
    10: (2,2)        30: (1,1,1,2)      50: (1,3,2)
    11: (2,1,1)      31: (1,1,1,1,1)    51: (1,3,1,1)
    12: (1,3)        32: (6)            52: (1,2,3)
    13: (1,2,1)      35: (4,1,1)        53: (1,2,2,1)
    14: (1,1,2)      36: (3,3)          54: (1,2,1,2)
    15: (1,1,1,1)    37: (3,2,1)        55: (1,2,1,1,1)
    16: (5)          38: (3,1,2)        56: (1,1,4)
		

Crossrefs

The version for prime indices is A000027, counted by A000041.
These compositions are counted by A116406.
The case of non-Heinz numbers of partitions is A119899, counted by A344608.
The version for Heinz numbers of partitions is A344609, counted by A344607.
These are the positions of terms >= 0 in A344618.
The version for unreversed alternating sum is A345913.
The opposite (k <= 0) version is A345916.
The strict (k > 0) case is A345918.
The complement is A345920, counted by A294175.
A011782 counts compositions.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A236913 counts partitions of 2n with reverse-alternating sum <= 0.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A344610 counts partitions by sum and positive reverse-alternating sum.
A344611 counts partitions of 2n with reverse-alternating sum >= 0.
A345197 counts compositions by sum, length, and alternating sum.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
Compositions of n, 2n, or 2n+1 with alternating/reverse-alternating sum k:
- k = 0: counted by A088218, ranked by A344619/A344619.
- k = 1: counted by A000984, ranked by A345909/A345911.
- k = -1: counted by A001791, ranked by A345910/A345912.
- k = 2: counted by A088218, ranked by A345925/A345922.
- k = -2: counted by A002054, ranked by A345924/A345923.
- k >= 0: counted by A116406, ranked by A345913/A345914.
- k <= 0: counted by A058622(n-1), ranked by A345915/A345916.
- k > 0: counted by A027306, ranked by A345917/A345918.
- k < 0: counted by A294175, ranked by A345919/A345920.
- k != 0: counted by A058622, ranked by A345921/A345921.
- k even: counted by A081294, ranked by A053754/A053754.
- k odd: counted by A000302, ranked by A053738/A053738.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[ Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    sats[y_]:=Sum[(-1)^(i-Length[y])*y[[i]],{i,Length[y]}];
    Select[Range[0,100],sats[stc[#]]>=0&]

A345915 Numbers k such that the k-th composition in standard order (row k of A066099) has alternating sum <= 0.

Original entry on oeis.org

0, 3, 6, 10, 12, 13, 15, 20, 24, 25, 27, 30, 36, 40, 41, 43, 46, 48, 49, 50, 51, 53, 54, 55, 58, 60, 61, 63, 72, 80, 81, 83, 86, 92, 96, 97, 98, 99, 101, 102, 103, 106, 108, 109, 111, 116, 120, 121, 123, 126, 136, 144, 145, 147, 150, 156, 160, 161, 162, 163
Offset: 1

Views

Author

Gus Wiseman, Jul 08 2021

Keywords

Comments

The alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The sequence of terms together with the corresponding compositions begins:
     0: ()
     3: (1,1)
     6: (1,2)
    10: (2,2)
    12: (1,3)
    13: (1,2,1)
    15: (1,1,1,1)
    20: (2,3)
    24: (1,4)
    25: (1,3,1)
    27: (1,2,1,1)
    30: (1,1,1,2)
    36: (3,3)
    40: (2,4)
    41: (2,3,1)
		

Crossrefs

The version for Heinz numbers of partitions is A028260 (counted by A027187).
These compositions are counted by A058622.
These are the positions of terms <= 0 in A124754.
The reverse-alternating version is A345916.
The opposite (k >= 0) version is A345917.
The strictly negative (k < 0) version is A345919.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A011782 counts compositions.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A236913 counts partitions of 2n with reverse-alternating sum <= 0.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A345197 counts compositions by sum, length, and alternating sum.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
Compositions of n, 2n, or 2n+1 with alternating/reverse-alternating sum k:
- k = 0: counted by A088218, ranked by A344619/A344619.
- k = 1: counted by A000984, ranked by A345909/A345911.
- k = -1: counted by A001791, ranked by A345910/A345912.
- k = 2: counted by A088218, ranked by A345925/A345922.
- k = -2: counted by A002054, ranked by A345924/A345923.
- k >= 0: counted by A116406, ranked by A345913/A345914.
- k <= 0: counted by A058622(n-1), ranked by A345915/A345916.
- k > 0: counted by A027306, ranked by A345917/A345918.
- k < 0: counted by A294175, ranked by A345919/A345920.
- k != 0: counted by A058622, ranked by A345921/A345921.
- k even: counted by A081294, ranked by A053754/A053754.
- k odd: counted by A000302, ranked by A053738/A053738.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
    Select[Range[0,100],ats[stc[#]]<=0&]

A345916 Numbers k such that the k-th composition in standard order (row k of A066099) has reverse-alternating sum <= 0.

Original entry on oeis.org

0, 3, 5, 9, 10, 13, 15, 17, 18, 23, 25, 29, 33, 34, 36, 39, 41, 43, 45, 46, 49, 50, 53, 55, 57, 58, 61, 63, 65, 66, 68, 71, 75, 77, 78, 81, 85, 89, 90, 95, 97, 98, 103, 105, 109, 113, 114, 119, 121, 125, 129, 130, 132, 135, 136, 139, 141, 142, 145, 147, 149
Offset: 1

Views

Author

Gus Wiseman, Jul 08 2021

Keywords

Comments

The reverse-alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(k-i) y_i.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The sequence of terms together with the corresponding compositions begins:
     0: ()
     3: (1,1)
     5: (2,1)
     9: (3,1)
    10: (2,2)
    13: (1,2,1)
    15: (1,1,1,1)
    17: (4,1)
    18: (3,2)
    23: (2,1,1,1)
    25: (1,3,1)
    29: (1,1,2,1)
    33: (5,1)
    34: (4,2)
    36: (3,3)
		

Crossrefs

The version for Heinz numbers of partitions is A000290.
These compositions are counted by A058622.
These are the positions of terms <= 0 in A344618.
The opposite (k >= 0) version is A345914.
The version for unreversed alternating sum is A345915.
The strictly negative (k < 0) version is A345920.
A011782 counts compositions.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A236913 counts partitions of 2n with reverse-alternating sum <= 0.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A344611 counts partitions of 2n with reverse-alternating sum >= 0.
A345197 counts compositions by sum, length, and alternating sum.
Standard compositions: A000120, A066099, A070939, A228351, A124754, A344618.
Compositions of n, 2n, or 2n+1 with alternating/reverse-alternating sum k:
- k = 0: counted by A088218, ranked by A344619/A344619.
- k = 1: counted by A000984, ranked by A345909/A345911.
- k = -1: counted by A001791, ranked by A345910/A345912.
- k = 2: counted by A088218, ranked by A345925/A345922.
- k = -2: counted by A002054, ranked by A345924/A345923.
- k >= 0: counted by A116406, ranked by A345913/A345914.
- k <= 0: counted by A058622(n-1), ranked by A345915/A345916.
- k > 0: counted by A027306, ranked by A345917/A345918.
- k < 0: counted by A294175, ranked by A345919/A345920.
- k != 0: counted by A058622, ranked by A345921/A345921.
- k even: counted by A081294, ranked by A053754/A053754.
- k odd: counted by A000302, ranked by A053738/A053738.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    sats[y_]:=Sum[(-1)^(i-Length[y])*y[[i]],{i,Length[y]}];
    Select[Range[0,100],sats[stc[#]]<=0&]

A177510 Number of compositions (p0, p1, p2, ...) of n with pi - p0 <= i and pi >= p0.

Original entry on oeis.org

1, 1, 2, 3, 5, 8, 14, 25, 46, 87, 167, 324, 634, 1248, 2466, 4887, 9706, 19308, 38455, 76659, 152925, 305232, 609488, 1217429, 2432399, 4860881, 9715511, 19421029, 38826059, 77626471, 155211785, 310357462, 620608652, 1241046343, 2481817484, 4963191718, 9925669171, 19850186856, 39698516655, 79394037319
Offset: 0

Views

Author

Mats Granvik, Dec 11 2010

Keywords

Comments

a(0)=1, otherwise row sums of A179748.
For n>=1 cumulative sums of A008930.
a(n) is proportional to A048651*A000079. The error (a(n)-A048651*A000079) divided by sequence A186425 tends to the golden ratio A001622. This can be seen when using about 1000 decimals of the constant A048651 = 0.2887880950866024212... - [Mats Granvik, Jan 01 2015]
From Gus Wiseman, Mar 31 2022: (Start)
Also the number of integer compositions of n with exactly one part on or above the diagonal. For example, the a(1) = 1 through a(5) = 8 compositions are:
(1) (2) (3) (4) (5)
(11) (21) (31) (41)
(111) (112) (212)
(211) (311)
(1111) (1112)
(1121)
(2111)
(11111)
(End)

Examples

			From _Joerg Arndt_, Mar 24 2014: (Start)
The a(7) = 25 such compositions are:
01:  [ 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 2 ]
03:  [ 1 1 1 1 2 1 ]
04:  [ 1 1 1 1 3 ]
05:  [ 1 1 1 2 1 1 ]
06:  [ 1 1 1 2 2 ]
07:  [ 1 1 1 3 1 ]
08:  [ 1 1 1 4 ]
09:  [ 1 1 2 1 1 1 ]
10:  [ 1 1 2 1 2 ]
11:  [ 1 1 2 2 1 ]
12:  [ 1 1 2 3 ]
13:  [ 1 1 3 1 1 ]
14:  [ 1 1 3 2 ]
15:  [ 1 2 1 1 1 1 ]
16:  [ 1 2 1 1 2 ]
17:  [ 1 2 1 2 1 ]
18:  [ 1 2 1 3 ]
19:  [ 1 2 2 1 1 ]
20:  [ 1 2 2 2 ]
21:  [ 1 2 3 1 ]
22:  [ 2 2 3 ]
23:  [ 2 3 2 ]
24:  [ 3 4 ]
25:  [ 7 ]
(End)
		

Crossrefs

Cf. A238859 (compositions with subdiagonal growth), A238876 (partitions with subdiagonal growth), A001227 (partitions into distinct parts with subdiagonal growth).
Cf. A238860 (partitions with superdiagonal growth), A238861 (compositions with superdiagonal growth), A000009 (partitions into distinct parts have superdiagonal growth by definition).
The version for partitions is A001477, strong A002620.
The version for permutations is A057427, strong A000295.
The opposite version is A238874, first column of A352522.
The version for fixed points is A240736, nonfixed A352520.
The strong version is A351983, column k=1 of A352524.
This is column k = 1 of A352525.
A238349 counts compositions by fixed points, first col A238351.
A352517 counts weak excedances of standard compositions.

Programs

  • Maple
    A179748 := proc(n,k) option remember; if k= 1 then 1; elif k> n then 0 ; else add( procname(n-i,k-1),i=1..k-1) ; end if; end proc:
    A177510 := proc(n) add(A179748(n,k),k=1..n) ;end proc:
    seq(A177510(n),n=1..20) ; # R. J. Mathar, Dec 14 2010
  • Mathematica
    Clear[t, nn]; nn = 39; t[n_, 1] = 1; t[n_, k_] := t[n, k] = If[n >= k, Sum[t[n - i, k - 1], {i, 1, k - 1}], 0]; Table[Sum[t[n, k], {k, 1, n}], {n, 1, nn}] (* Mats Granvik, Jan 01 2015 *)
    pdw[y_]:=Length[Select[Range[Length[y]],#<=y[[#]]&]]; Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],pdw[#]==1&]],{n,0,10}] (* Gus Wiseman, Mar 31 2022 *)
  • PARI
    N=66; q='q+O('q^N); Vec( 1 + q/(1-q) * sum(n=0, N, q^n * prod(k=1, n, (1-q^k)/(1-q) ) ) ) \\ Joerg Arndt, Mar 24 2014
  • Sage
    @CachedFunction
    def T(n, k): # A179748
        if n == 0:  return int(k==0);
        if k == 1:  return int(n>=1);
        return sum( T(n-i, k-1) for i in [1..k-1] );
    # to display triangle A179748 including column zero = [1,0,0,0,...]:
    #for n in [0..10]: print([ T(n,k) for k in [0..n] ])
    def a(n): return sum( T(n,k) for k in [0..n] )
    print([a(n) for n in [0..66]])
    # Joerg Arndt, Mar 24 2014
    

Formula

G.f.: 1 + q/(1-q) * sum(n>=0, q^n * prod(k=1..n, (1-q^k)/(1-q) ) ). [Joerg Arndt, Mar 24 2014]

Extensions

New name and a(0) = 1 prepended, Joerg Arndt, Mar 24 2014

A039004 Numbers whose base-4 representation has the same number of 1's and 2's.

Original entry on oeis.org

0, 3, 6, 9, 12, 15, 18, 24, 27, 30, 33, 36, 39, 45, 48, 51, 54, 57, 60, 63, 66, 72, 75, 78, 90, 96, 99, 102, 105, 108, 111, 114, 120, 123, 126, 129, 132, 135, 141, 144, 147, 150, 153, 156, 159, 165, 177, 180, 183, 189, 192, 195, 198, 201, 204, 207, 210, 216, 219
Offset: 1

Views

Author

Keywords

Comments

Numbers such that sum (-1)^k*b(k) = 0 where b(k)=k-th binary digit of n (see A065359). - Benoit Cloitre, Nov 18 2003
Conjecture: a(C(2n,n)-1) = 4^n - 1. (A000984 is C(2n,n)). - Gerald McGarvey, Nov 18 2007
From Russell Jay Hendel, Jun 23 2015: (Start)
We prove the McGarvey conjecture (A) a(e(n,n)-1) = 4^n-1, with e(n,m) = A034870(n,m) = binomial(2n,m), the even rows of Pascal's triangle. By the comment from Hendel in A034870, we have the function s(n,k) = #{n-digit, base-4 numbers with n-k more 1-digits than 2-digits}. As shown in A034870, (B) #s(n,k)= e(n,k) with # indicating cardinality, that is, e(n,k) = binomial(2n,k) gives the number of n-digit, base-4 numbers with n-k more 1-digits than 2-digits.
We now show that (B) implies (A). By definition, s(n,n) contains the e(n,n) = binomial(2n,n) numbers with an equal number of 1-digits and 2-digits. The biggest n-digit, base-4 number is 333...3 (n copies of 3). Since 333...33 has zero 1-digits and zero 2-digits it follows that 333...333 is a member of s(n,n) and hence it is the biggest member of s(n,n). But 333...333 (n copies of 3) in base 4 has value 4^n-1. Since A039004 starts with index 0 (that is, 0 is the 0th member of A039004), it immediately follows that 4^n-1 is the (e(n,n)-1)st member of A039004, proving the McGarvey conjecture. (End)
Also numbers whose alternating sum of binary expansion is 0, i.e., positions of zeros in A345927. These are numbers whose binary expansion has the same number of 1's at even positions as at odd positions. - Gus Wiseman, Jul 28 2021

Crossrefs

A subset of A001969 (evil numbers).
A base-2 version is A031443 (digitally balanced numbers).
Positions of 0's in A065359 and A345927.
Positions of first appearances are A086893.
The version for standard compositions is A344619.
A000120 and A080791 count binary digits, with difference A145037.
A003714 lists numbers with no successive binary indices.
A011782 counts compositions.
A030190 gives the binary expansion of each nonnegative integer.
A070939 gives the length of an integer's binary expansion.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A101211 lists run-lengths in binary expansion:
- row-lengths: A069010
- reverse: A227736
- ones only: A245563
A138364 counts compositions with alternating sum 0:
- bisection: A001700/A088218
- complement: A058622
A328594 lists numbers whose binary expansion is aperiodic.
A345197 counts compositions by length and alternating sum.

Programs

  • Fortran
    c See link in A139351.
  • Maple
    N:= 1000: # to get all terms up to N, which should be divisible by 4
    B:= Array(0..N-1):
    d:= ceil(log[4](N));
    S:= Array(0..N-1,[seq(op([0,1,-1,0]),i=1..N/4)]):
    for i from 1 to d do
      B:= B + S;
      S:= Array(0..N-1,i-> S[floor(i/4)]);
    od:
    select(t -> B[t]=0, [$0..N-1]); # Robert Israel, Jun 24 2015
  • Mathematica
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
    Select[Range[0,100],ats[IntegerDigits[#,2]]==0&] (* Gus Wiseman, Jul 28 2021 *)
  • PARI
    for(n=0,219,if(sum(i=1,length(binary(n)),(-1)^i*component(binary(n),i))==0,print1(n,",")))
    

Formula

Conjecture: there is a constant c around 5 such that a(n) is asymptotic to c*n. - Benoit Cloitre, Nov 24 2002
That conjecture is false. The number of members of the sequence from 0 to 4^d-1 is binomial(2d,d) which by Stirling's formula is asymptotic to 4^d/sqrt(Pi*d). If Cloitre's conjecture were true we would have 4^d-1 asymptotic to c*4^d/sqrt(Pi*d), a contradiction. - Robert Israel, Jun 24 2015

A060150 a(0) = 1; for n > 0, binomial(2n-1, n-1)^2.

Original entry on oeis.org

1, 1, 9, 100, 1225, 15876, 213444, 2944656, 41409225, 590976100, 8533694884, 124408576656, 1828114918084, 27043120090000, 402335398890000, 6015361252737600, 90324408810638025, 1361429497505672100, 20589520178326522500, 312321918272897610000
Offset: 0

Views

Author

N. J. A. Sloane, Apr 10 2001

Keywords

Comments

Number of square lattice walks that start at (0,0) and end at (1,0) after 2n-1 steps, free to pass through (1,0) at intermediate steps. - Steven Finch, Dec 20 2001
Number of paths of length n connecting two neighboring nodes in optimal chordal graph of degree 4, G(2*d(G)^2+2*d(G)+1,2d(G)+1), of diameter d(G). - B. Dubalski (dubalski(AT)atr.bydgoszcz.pl), Feb 05 2002
a(n) is the number of ways to place n red balls and n blue balls into n distinguishable boxes with no restrictions on the number of balls put in a box. - Geoffrey Critzer, Jul 08 2013
The number of square lattice walks of n steps that start at the origin and end at (k,0) is zero if n-k is odd and [binomial(n,(n-k)/2)]^2 if n-k is even. - R. J. Mathar, Sep 28 2020

References

  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, 1994 Addison-Wesley company, Inc.
  • A. P. Prudnikov, Yu. A. Brychkov and O.I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", New York, Gordon and Breach Science Publishers, 1986-1992, Eq. (5.1.29.2)
  • K. A. Ross and C. R. B. Wright, Discrete Mathematics, 1992 Prentice Hall Inc.

Crossrefs

Programs

  • Maple
    seq(coeff(series(EllipticK(4*sqrt(x))/(2*Pi) + 3/4, x=0, n+1), x, n), n=0..30);  # Mark van Hoeij, Apr 30 2013
  • Mathematica
    Table[Binomial[2n-1,n]^2,{n,0,19}] (* Geoffrey Critzer, Jul 08 2013 *)
  • PARI
    a(n)=if(n<2, 1, binomial(2*n-1,n-1)^2)
    
  • PARI
    for (n=0, 200, if (n==0, a=1, a=binomial(2*n - 1, n - 1)^2); write("b060150.txt", n, " ", a)) \\ Harry J. Smith, Jul 02 2009

Formula

a(n) = A088218(n)^2.
a(n) = A002894(n)/4 for n>0.
G.f.: 1 + (1/AGM(1, sqrt(1-16*x))-1)/4. - Michael Somos, Dec 12 2002
G.f. = 1 + (K(16x)-1)/4 = 1 + Sum_{k>0} q^k/(1+q^(2k)) where K(16x) is the complete Elliptic integral of the first kind at 16x=k^2 and q is the nome. - Michael Somos, May 09 2005
G.f.: 1 + x*3F2((1, 3/2, 3/2); (2, 2))(16*x). - Olivier Gérard, Feb 16 2011
E.g.f.: Sum_{n>0} a(n)*x^(2n-1)/(2n-1)! = BesselI(0, 2x)*BesselI(1, 2x) . - Michael Somos, Jun 22 2005
D-finite with recurrence n^2*a(n) -4*(2*n-1)^2*a(n-1)=0. - R. J. Mathar, Jul 26 2014
From Seiichi Manyama, Oct 19 2016: (Start)
Let the number of multisets of length k on n symbols be denoted by ((n, k)) = binomial(n+k-1, k).
a(n) = (Sum_{0 <= k <= n} binomial(n, k)^2 * ((2*n, n - k)))/3 for n > 0. (End)
a(n) ~ 4^(2*n-1)/(Pi*n). - Ilya Gutkovskiy, Oct 19 2016
For n >= 1, a(n) = 1/n * Sum_{k = 0..n-1} (n + 2*k)*binomial(n+k-1, k)^2 = ( 1/(4*n) * Sum_{k = 0..n} (n + 2*k)*binomial(-n+k-1, k)^2 )^2. - Peter Bala, Nov 02 2024

A100100 Triangle T(n,k) = binomial(2*n-k-1, n-k) read by rows.

Original entry on oeis.org

1, 1, 1, 3, 2, 1, 10, 6, 3, 1, 35, 20, 10, 4, 1, 126, 70, 35, 15, 5, 1, 462, 252, 126, 56, 21, 6, 1, 1716, 924, 462, 210, 84, 28, 7, 1, 6435, 3432, 1716, 792, 330, 120, 36, 8, 1, 24310, 12870, 6435, 3003, 1287, 495, 165, 45, 9, 1, 92378, 48620, 24310, 11440, 5005, 2002
Offset: 0

Views

Author

Paul Barry, Nov 08 2004

Keywords

Comments

Sometimes called a Catalan triangle, although there are many other triangles that carry that name - see A009766, A008315, A028364, A033184, A053121, A059365, A062103.
Number of nodes of outdegree k in all ordered trees with n edges. Equivalently, number of ascents of length k in all Dyck paths of semilength n. Example: T(3,2) = 3 because the Dyck paths of semilength 3 are UDUDUD, UD(UU)DD, (UU)DDUD, (UU)DUDD and UUUDDD, where U = (1,1), D = (1,-1), the ascents of length 2 being shown between parentheses. - Emeric Deutsch, Nov 19 2006
Riordan array (f(x), x*g(x)) where f(x) is the g.f. of A088218 and g(x) is the g.f. of A000108. - Philippe Deléham, Jan 23 2010
T(n,k) is the number of nonnegative paths of upsteps U = (1,1) and downsteps D = (1,-1) of length 2*n with k returns to ground level, the horizontal line through the initial vertex. Example: T(2,1) = 2 counts UDUU, UUDD. Also, T(n,k) = number of these paths whose last descent has length k, that is, k downsteps follow the last upstep. Example: T(2,1) = 2 counts UUUD, UDUD. - David Callan, Nov 21 2011
Belongs to the hitting-time subgroup of the Riordan group. Multiplying this triangle by the square Pascal matrix gives A092392 read as a square array. See the example below. - Peter Bala, Nov 03 2015

Examples

			From _Paul Barry_, Mar 15 2010: (Start)
Triangle begins in row n=0 with columns 0<=k<=n as:
    1;
    1,   1;
    3,   2,   1;
   10,   6,   3,  1;
   35,  20,  10,  4,  1;
  126,  70,  35, 15,  5, 1;
  462, 252, 126, 56, 21, 6, 1;
Production matrix begins
  1, 1;
  2, 1, 1;
  3, 1, 1, 1;
  4, 1, 1, 1, 1;
  5, 1, 1, 1, 1, 1;
  6, 1, 1, 1, 1, 1, 1;
  7, 1, 1, 1, 1, 1, 1, 1;
(End)
A092392 as a square array = A100100 * square Pascal matrix:
/1   1  1  1 ...\   / 1          \/1 1  1  1 ...\
|2   3  4  5 ...|   | 1 1        ||1 2  3  4 ...|
|6  10 15 21 ...| = | 3 2 1      ||1 3  6 10 ...|
|20 35 56 84 ...|   |10 6 3 1    ||1 4 10 20 ...|
|70 ...         |   |35 ...      ||1 ...        |
- _Peter Bala_, Nov 03 2015
		

Crossrefs

Row sums are A000984. Equivalent to A092392, to which A088218 has been added as a first column. Columns include A088218, A000984, A001700, A001791, A002054, A002694. Diagonal sums are A100217. Matrix inverse is A100218.
Cf. A059481 (mirrored). Cf. A033184, A094527, A113955.

Programs

  • Haskell
    a100100 n k = a100100_tabl !! n !! n
    a100100_row n = a100100_tabl !! n
    a100100_tabl = [1] : f a092392_tabl where
       f (us : wss'@(vs : wss)) = (vs !! 1 : us) : f wss'
    -- Reinhard Zumkeller, Jan 15 2014
    
  • Magma
    /* As triangle */ [[Binomial(2*n - k - 1, n - k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Nov 21 2018
  • Maple
    A100100 := proc(n,k)
        binomial(2*n-k-1,n-1) ;
    end proc:
    seq(seq(A100100(n,k),k=0..n),n=0..10) ; # R. J. Mathar, Feb 06 2015
  • Mathematica
    Flatten[Table[Binomial[2 n - k - 1, n - k], {n, 0, 11}, {k, 0, n}]] (* Vincenzo Librandi, Nov 21 2018 *)
  • PARI
    T(n,k)=binomial(2*n-k-1,n-k) \\ Charles R Greathouse IV, Jan 16 2012
    

Formula

From Peter Bala, Sep 06 2015: (Start)
Matrix product A094527 * P^(-1) = A113955 * P^(-2), where P denotes Pascal's triangle A007318.
Essentially, the logarithmic derivative of A033184. (End)
Let column(k) = [T(n, k), n >= k], then the generating function for column(k) is (2/(sqrt(1-4*x)+1))^(k-1)/sqrt(1-4*x). - Peter Luschny, Mar 19 2021
O.g.f. row polynomials R(n, x) := Sum_{k=0..n} T(n, k)*x^k, i.e. o.g.f. of the triangle, G(z,x) = 1/((2 - c(z))*(1 - x*z*c(z))), with c the o.g.f. of A000108 (Catalan). See the Riordan coomment by Philippe Deléham above. - Wolfdieter Lang, Apr 06 2021
Previous Showing 41-50 of 196 results. Next