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

A126120 Catalan numbers (A000108) interpolated with 0's.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 5, 0, 14, 0, 42, 0, 132, 0, 429, 0, 1430, 0, 4862, 0, 16796, 0, 58786, 0, 208012, 0, 742900, 0, 2674440, 0, 9694845, 0, 35357670, 0, 129644790, 0, 477638700, 0, 1767263190, 0, 6564120420, 0, 24466267020, 0, 91482563640, 0, 343059613650, 0
Offset: 0

Views

Author

Philippe Deléham, Mar 06 2007

Keywords

Comments

Inverse binomial transform of A001006.
The Hankel transform of this sequence gives A000012 = [1,1,1,1,1,...].
Counts returning walks (excursions) of length n on a 1-d integer lattice with step set {+1,-1} which stay in the chamber x >= 0. - Andrew V. Sutherland, Feb 29 2008
Moment sequence of the trace of a random matrix in G=USp(2)=SU(2). If X=tr(A) is a random variable (A distributed according to the Haar measure on G) then a(n) = E[X^n]. - Andrew V. Sutherland, Feb 29 2008
Essentially the same as A097331. - R. J. Mathar, Jun 15 2008
Number of distinct proper binary trees with n nodes. - Chris R. Sims (chris.r.sims(AT)gmail.com), Jun 30 2010
-a(n-1), with a(-1):=0, n>=0, is the Z-sequence for the Riordan array A049310 (Chebyshev S). For the definition see that triangle. - Wolfdieter Lang, Nov 04 2011
See A180874 (also A238390 and A097610) and A263916 for relations to the general Bell A036040, cycle index A036039, and cumulant expansion polynomials A127671 through the Faber polynomials. - Tom Copeland, Jan 26 2016
A signed version is generated by evaluating polynomials in A126216 that are essentially the face polynomials of the associahedra. This entry's sequence is related to an inversion relation on p. 34 of Mizera, related to Feynman diagrams. - Tom Copeland, Dec 09 2019

Examples

			G.f. = 1 + x^2 + 2*x^4 + 5*x^6 + 14*x^8 + 42*x^10 + 132*x^12 + 429*x^14 + ...
From _Gus Wiseman_, Nov 14 2022: (Start)
The a(0) = 1 through a(8) = 14 ordered binary rooted trees with n + 1 nodes (ranked by A358375):
  o  .  (oo)  .  ((oo)o)  .  (((oo)o)o)  .  ((((oo)o)o)o)
                 (o(oo))     ((o(oo))o)     (((o(oo))o)o)
                             ((oo)(oo))     (((oo)(oo))o)
                             (o((oo)o))     (((oo)o)(oo))
                             (o(o(oo)))     ((o((oo)o))o)
                                            ((o(o(oo)))o)
                                            ((o(oo))(oo))
                                            ((oo)((oo)o))
                                            ((oo)(o(oo)))
                                            (o(((oo)o)o))
                                            (o((o(oo))o))
                                            (o((oo)(oo)))
                                            (o(o((oo)o)))
                                            (o(o(o(oo))))
(End)
		

References

  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Ch. 49, Hemisphere Publishing Corp., 1987.

Crossrefs

Cf. A126216.
The unordered version is A001190, ranked by A111299.
These trees (ordered binary rooted) are ranked by A358375.

Programs

  • Magma
    &cat [[Catalan(n), 0]: n in [0..30]]; // Vincenzo Librandi, Jul 28 2016
    
  • Maple
    with(combstruct): grammar := { BB = Sequence(Prod(a,BB,b)), a = Atom, b = Atom }: seq(count([BB,grammar], size=n),n=0..47); # Zerinvary Lajos, Apr 25 2007
    BB := {E=Prod(Z,Z), S=Union(Epsilon,Prod(S,S,E))}: ZL:=[S,BB,unlabeled]: seq(count(ZL, size=n), n=0..45); # Zerinvary Lajos, Apr 22 2007
    BB := [T,{T=Prod(Z,Z,Z,F,F), F=Sequence(B), B=Prod(F,Z,Z)}, unlabeled]: seq(count(BB, size=n+1), n=0..45); # valid for n> 0. # Zerinvary Lajos, Apr 22 2007
    seq(n!*coeff(series(hypergeom([],[2],x^2),x,n+2),x,n),n=0..45); # Peter Luschny, Jan 31 2015
    # Using function CompInv from A357588.
    CompInv(48, n -> ifelse(irem(n, 2) = 0, 0, (-1)^iquo(n-1, 2))); # Peter Luschny, Oct 07 2022
  • Mathematica
    a[n_?EvenQ] := CatalanNumber[n/2]; a[n_] = 0; Table[a[n], {n, 0, 45}] (* Jean-François Alcover, Sep 10 2012 *)
    a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ BesselI[ 1, 2 x] / x, {x, 0, n}]]; (* Michael Somos, Mar 19 2014 *)
    bot[n_]:=If[n==1,{{}},Join@@Table[Tuples[bot/@c],{c,Table[{k,n-k-1},{k,n-1}]}]];
    Table[Length[bot[n]],{n,10}] (* Gus Wiseman, Nov 14 2022 *)
    Riffle[CatalanNumber[Range[0,50]],0,{2,-1,2}] (* Harvey P. Dale, May 28 2024 *)
  • Python
    from math import comb
    def A126120(n): return 0 if n&1 else comb(n,m:=n>>1)//(m+1) # Chai Wah Wu, Apr 22 2024
  • Sage
    def A126120_list(n) :
        D = [0]*(n+2); D[1] = 1
        b = True; h = 2; R = []
        for i in range(2*n-1) :
            if b :
                for k in range(h,0,-1) : D[k] -= D[k-1]
                h += 1; R.append(abs(D[1]))
            else :
                for k in range(1,h, 1) : D[k] += D[k+1]
            b = not b
        return R
    A126120_list(46) # Peter Luschny, Jun 03 2012
    

Formula

a(2*n) = A000108(n), a(2*n+1) = 0.
a(n) = A053121(n,0).
(1/Pi) Integral_{0 .. Pi} (2*cos(x))^n *2*sin^2(x) dx. - Andrew V. Sutherland, Feb 29 2008
G.f.: (1 - sqrt(1 - 4*x^2)) / (2*x^2) = 1/(1-x^2/(1-x^2/(1-x^2/(1-x^2/(1-... (continued fraction). - Philippe Deléham, Nov 24 2009
G.f. A(x) satisfies A(x) = 1 + x^2*A(x)^2. - Vladimir Kruchinin, Feb 18 2011
E.g.f.: I_1(2x)/x Where I_n(x) is the modified Bessel function. - Benjamin Phillabaum, Mar 07 2011
Apart from the first term the e.g.f. is given by x*HyperGeom([1/2],[3/2,2], x^2). - Benjamin Phillabaum, Mar 07 2011
a(n) = Integral_{x=-2..2} x^n*sqrt((2-x)*(2+x))/(2*Pi) dx. - Peter Luschny, Sep 11 2011
E.g.f.: E(0)/(1-x) where E(k) = 1-x/(1-x/(x-(k+1)*(k+2)/E(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 05 2013
G.f.: 3/2- sqrt(1-4*x^2)/2 = 1/x^2 + R(0)/x^2, where R(k) = 2*k-1 - x^2*(2*k-1)*(2*k+1)/R(k+1); (continued fraction). - Sergei N. Gladkovskii, Oct 28 2013 (warning: this is not the g.f. of this sequence, R. J. Mathar, Sep 23 2021)
G.f.: 1/Q(0), where Q(k) = 2*k+1 + x^2*(1-4*(k+1)^2)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Jan 09 2014
a(n) = n!*[x^n]hypergeom([],[2],x^2). - Peter Luschny, Jan 31 2015
a(n) = 2^n*hypergeom([3/2,-n],[3],2). - Peter Luschny, Feb 03 2015
a(n) = ((-1)^n+1)*2^(2*floor(n/2)-1)*Gamma(floor(n/2)+1/2)/(sqrt(Pi)* Gamma(floor(n/2)+2)). - Ilya Gutkovskiy, Jul 23 2016
D-finite with recurrence (n+2)*a(n) +4*(-n+1)*a(n-2)=0. - R. J. Mathar, Mar 21 2021
From Peter Bala, Feb 03 2024: (Start)
a(n) = 2^n * Sum_{k = 0..n} (-2)^(-k)*binomial(n, k)*Catalan(k+1).
G.f.: 1/(1 + 2*x) * c(x/(1 + 2*x))^2 = 1/(1 - 2*x) * c(-x/(1 - 2*x))^2 = c(x^2), where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108. (End)

Extensions

An erroneous comment removed by Tom Copeland, Jul 23 2016

A036043 Irregular triangle read by rows: row n (n >= 0) gives number of parts in all partitions of n (in Abramowitz and Stegun order).

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 3, 1, 2, 2, 3, 4, 1, 2, 2, 3, 3, 4, 5, 1, 2, 2, 2, 3, 3, 3, 4, 4, 5, 6, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 5, 5, 6, 7, 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 6, 6, 7, 8, 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 6, 6, 6, 7, 7, 8, 9
Offset: 0

Views

Author

Keywords

Comments

The sequence of row lengths of this array is p(n) = A000041(n) (partition numbers).
The sequence of row sums is A006128(n).
The number of times k appears in row n is A008284(n,k). - Franklin T. Adams-Watters, Jan 12 2006
The next level (row) gets created from each node by adding one or two more nodes. If a single node is added, its value is one more than the value of its parent. If two nodes are added, the first is equal in value to the parent and the value of the second is one more than the value of the parent. See A128628. - Alford Arnold, Mar 27 2007
The 1's in the (flattened) sequence mark the start of a new row, the value that precedes the 1 equals the row number minus one. (I.e., the 1 preceded by a 0 is the start of row 1, the 1 preceded by a 6 is the start of row 7, etc.) - M. F. Hasler, Jun 06 2018
Also the maximum part in the n-th partition in graded lexicographic order (sum/lex, A193073). - Gus Wiseman, May 24 2020

Examples

			0;
1;
1, 2;
1, 2, 3;
1, 2, 2, 3, 4;
1, 2, 2, 3, 3, 4, 5;
1, 2, 2, 2, 3, 3, 3, 4, 4, 5, 6;
1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 5, 5, 6, 7;
		

References

  • Abramowitz and Stegun, Handbook, p. 831, column labeled "m".

Crossrefs

Row lengths are A000041.
Partition lengths of A036036 and A334301.
The version not sorted by length is A049085.
The generalization to compositions is A124736.
The Heinz number of the same partition is A334433.
The number of distinct elements in the same partition is A334440.
The maximum part of the same partition is A334441.
Lexicographically ordered reversed partitions are A026791.
Lexicographically ordered partitions are A193073.

Programs

  • Maple
    with(combinat): nmax:=9: for n from 1 to nmax do y(n):=numbpart(n): P(n):=sort(partition(n)): for k from 1 to y(n) do B(k) := P(n)[k] od: for k from 1 to y(n) do s:=0: j:=0: while sJohannes W. Meijer, Jun 21 2010, revised Nov 29 2012
    # alternative implementation based on A119441 by R. J. Mathar, Jul 12 2013
    A036043 := proc(n,k)
        local pi;
        pi := ASPrts(n)[k] ;
        nops(pi) ;
    end proc:
    for n from 1 to 10 do
        for k from 1 to A000041(n) do
            printf("%d,",A036043(n,k)) ;
        end do:
        printf("\n") ;
    end do:
  • Mathematica
    Table[Length/@Sort[IntegerPartitions[n]],{n,0,30}] (* Gus Wiseman, May 22 2020 *)
  • PARI
    A036043(n,k)=#partitions(n)[k] \\ M. F. Hasler, Jun 06 2018
    
  • SageMath
    def A036043_row(n):
        return [len(p) for k in (0..n) for p in Partitions(n, length=k)]
    for n in (0..10): print(A036043_row(n)) # Peter Luschny, Nov 02 2019

Formula

a(n) = A001222(A334433(n)). - Gus Wiseman, May 22 2020

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Jun 17 2001
a(0) inserted by Franklin T. Adams-Watters, Jun 24 2014
Incorrect formula deleted by M. F. Hasler, Jun 06 2018

A048996 Irregular triangle read by rows. Preferred multisets: numbers refining A007318 using format described in A036038.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 3, 1, 1, 2, 2, 3, 3, 4, 1, 1, 2, 2, 1, 3, 6, 1, 4, 6, 5, 1, 1, 2, 2, 2, 3, 6, 3, 3, 4, 12, 4, 5, 10, 6, 1, 1, 2, 2, 2, 1, 3, 6, 6, 3, 3, 4, 12, 6, 12, 1, 5, 20, 10, 6, 15, 7, 1, 1, 2, 2, 2, 2, 3, 6, 6, 3, 3, 6, 1, 4, 12, 12, 12, 12, 4, 5, 20, 10, 30, 5, 6, 30, 20, 7, 21, 8, 1
Offset: 0

Views

Author

Keywords

Comments

This array gives in row n>=1 the multinomial numbers (call them M_0 numbers) m!/product((a_j)!,j=1..n) with the exponents of the partitions of n with number of parts m:=sum(a_j,j=1..n), given in the Abramowitz-Stegun order. See p. 831 of the given reference. See also the arrays for the M_1, M_2 and M_3 multinomial numbers A036038, A036039 and A036040 (or A080575).
For a signed version see A111786.
These M_0 multinomial numbers give the number of compositions of n >= 1 with parts corresponding to the partitions of n (in A-St order). See an n = 5 example below. The triangle with the summed entries of like number of parts m is A007318(n-1, m-1) (Pascal). - Wolfdieter Lang, Jan 29 2021

Examples

			Table starts:
[1]
[1]
[1, 1]
[1, 2, 1]
[1, 2, 1, 3, 1]
[1, 2, 2, 3, 3, 4, 1]
[1, 2, 2, 1, 3, 6, 1, 4, 6,  5, 1]
[1, 2, 2, 2, 3, 6, 3, 3, 4, 12, 4, 5, 10, 6, 1]
.
T(5,6) = 4 because there are four multisets using the first four digits {0,1,2,3}: 32100, 32110, 32210 and 33210
T(5,6) = 4 because there are 4 compositions of 5 that can be formed from the partition 2+1+1+1. - _Geoffrey Critzer_, May 19 2013
These 4 compositions 2+1+1+1, 1+2+1+1, 1+1+2+1 and 1+1+1+2 of 5 correspond to the 4 set partitions of [5] :={1,2,3,4,5}, with 4 blocks of consecutive numbers, namely {1,2},{3},{4},{5} and {1},{2,3},{4},{5} and {1},{2},{3,4},{5} and {1},{2},{3},{4,5}. - _Wolfdieter Lang_, May 30 2018
		

Crossrefs

Cf. A000670, A007318, A036035, A036038, A019538, A115621, A309004, A000079 (row sums), A000041 (row lengths).

Programs

  • Maple
    nmax:=9: with(combinat): for n from 1 to nmax do P(n):=sort(partition(n)): for r from 1 to numbpart(n) do B(r):=P(n)[r] od: for m from 1 to numbpart(n) do s:=0: j:=0: while sA036040(n, m) := (add(q(t), t=1..n))!/(mul(q(t)!, t=1..n)); od: od: seq(seq(A036040(n, m), m=1..numbpart(n)), n=1..nmax); # Johannes W. Meijer, Jul 14 2016
  • PARI
    C(sig)={my(S=Set(sig)); (#sig)!/prod(k=1, #S, (#select(t->t==S[k], sig))!)}
    Row(n)={apply(C, [Vecrev(p) | p<-partitions(n)])}
    { for(n=0, 7, print(Row(n))) } \\ Andrew Howroyd, Oct 18 2020
  • SageMath
    from collections import Counter
    def ASPartitions(n, k):
        Q = [p.to_list() for p in Partitions(n, length=k)]
        for q in Q: q.reverse()
        return sorted(Q)
    def A048996_row(n):
        h = lambda p: product(map(factorial, Counter(p).values()))
        return [factorial(len(p))//h(p) for k in (0..n) for p in ASPartitions(n, k)]
    for n in (1..10): print(A048996_row(n)) # Peter Luschny, Nov 02 2019 [corrected on notice from Sean A. Irvine, Apr 30 2022]
    

Formula

T(n,k) = A036040(n,k) * Factorial(A036043(n,k)) / A036038(n,k) = A049019(n,k) / A036038(n,k).
If the n-th partition is P, a(n) is the multinomial coefficient of the signature of P. - Franklin T. Adams-Watters, May 30 2006
T(n,k) = A309004(A036035(n,k)). - Andrew Howroyd, Oct 19 2020

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Jun 17 2001
a(0)=1 prepended by Andrew Howroyd, Oct 19 2020

A080510 Triangle read by rows: T(n,k) gives the number of set partitions of {1,...,n} with maximum block length k.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 9, 4, 1, 1, 25, 20, 5, 1, 1, 75, 90, 30, 6, 1, 1, 231, 420, 175, 42, 7, 1, 1, 763, 2016, 1015, 280, 56, 8, 1, 1, 2619, 10024, 6111, 1890, 420, 72, 9, 1, 1, 9495, 51640, 38010, 12978, 3150, 600, 90, 10, 1, 1, 35695, 276980, 244035, 91938, 24024, 4950, 825, 110, 11, 1
Offset: 1

Views

Author

Wouter Meeussen, Mar 22 2003

Keywords

Comments

Row sums are A000110 (Bell numbers). Second column is A001189 (Degree n permutations of order exactly 2).
From Peter Luschny, Mar 09 2009: (Start)
Partition product of Product_{j=0..n-1} ((k + 1)*j - 1) and n! at k = -1, summed over parts with equal biggest part (see the Luschny link).
Underlying partition triangle is A036040.
Same partition product with length statistic is A008277.
Diagonal a(A000217) = A000012.
Row sum is A000110. (End)
From Gary W. Adamson, Feb 24 2011: (Start)
Construct an array in which the n-th row is the partition function G(n,k), where G(n,1),...,G(n,6) = A000012, A000085, A001680, A001681, A110038, A148092, with the first few rows
1, 1, 1, 1, 1, 1, 1, ... = A000012
1, 2, 4, 10, 26, 76, 232, ... = A000085
1, 2, 5, 14, 46, 166, 652, ... = A001680
1, 2, 5, 15, 51, 196, 827, ... = A001681
1, 2 5 15 52 202 869, ... = A110038
1, 2, 5 15 52 203 876, ... = A148092
...
Rows tend to A000110, the Bell numbers. Taking finite differences from the top, then reorienting, we obtain triangle A080510.
The n-th row of the array is the eigensequence of an infinite lower triangular matrix with n diagonals of Pascal's triangle starting from the right and the rest zeros. (End)

Examples

			T(4,3) = 4 since there are 4 set partitions with longest block of length 3: {{1},{2,3,4}}, {{1,3,4},{2}}, {{1,2,3},{4}} and {{1,2,4},{3}}.
Triangle begins:
  1;
  1,    1;
  1,    3,     1;
  1,    9,     4,    1;
  1,   25,    20,    5,    1;
  1,   75,    90,   30,    6,   1;
  1,  231,   420,  175,   42,   7,  1;
  1,  763,  2016, 1015,  280,  56,  8,  1;
  1, 2619, 10024, 6111, 1890, 420, 72,  9,  1;
  ...
		

Crossrefs

Columns k=1..10 give: A000012 (for n>0), A001189, A229245, A229246, A229247, A229248, A229249, A229250, A229251, A229252. - Alois P. Heinz, Sep 17 2013
T(2n,n) gives A276961.
Take differences along rows of A229223. - N. J. A. Sloane, Jan 10 2018

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
           add(b(n-i*j, i-1) *n!/i!^j/(n-i*j)!/j!, j=0..n/i)))
        end:
    T:= (n, k)-> b(n, k) -b(n, k-1):
    seq(seq(T(n, k), k=1..n), n=1..12);  # Alois P. Heinz, Apr 20 2012
  • Mathematica
    << DiscreteMath`NewCombinatorica`; Table[Length/@Split[Sort[Max[Length/@# ]&/@SetPartitions[n]]], {n, 12}]
    (* Second program: *)
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i<1, 0, Sum[b[n-i*j, i-1]*n!/i!^j/(n-i*j)!/j!, {j, 0, n/i}]]]; T[n_, k_] := b[n, k]-b[n, k-1]; Table[Table[T[n, k], {k, 1, n}], {n, 1, 12}] // Flatten (* Jean-François Alcover, Feb 25 2014, after Alois P. Heinz *)

Formula

E.g.f. for k-th column: exp(exp(x)*GAMMA(k, x)/(k-1)!-1)*(exp(x^k/k!)-1). - Vladeta Jovovic, Feb 04 2005
From Peter Luschny, Mar 09 2009: (Start)
T(n,0) = [n = 0] (Iverson notation) and for n > 0 and 1 <= m <= n.
T(n,m) = Sum_{a} M(a)|f^a| where a = a_1,...,a_n such that
1*a_1 + 2*a_2 + ... + n*a_n = n and max{a_i} = m, M(a) = n!/(a_1!*...*a_n!),
f^a = (f_1/1!)^a_1*...*(f_n/n!)^a_n and f_n = Product_{j=0..n-1} (-1) = (-1)^n. (End)
From Ludovic Schwob, Jan 15 2022: (Start)
T(2n,n) = C(2n,n)*(A000110(n)-1/2) for n>0.
T(n,m) = C(n,m)*A000110(n-m) for 2m > n > 0. (End)

A106400 Thue-Morse sequence: let A_k denote the first 2^k terms; then A_0 = 1 and for k >= 0, A_{k+1} = A_k B_k, where B_k is obtained from A_k by interchanging 1's and -1's.

Original entry on oeis.org

1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, -1, 1, -1, -1, 1, -1, 1, 1, -1, -1, 1, 1, -1, 1, -1, -1, 1, 1, -1, -1, 1, -1, 1
Offset: 0

Views

Author

Michael Somos, May 02 2005

Keywords

Comments

See A010060, the main entry for the Thue-Morse sequence, for additional information. - N. J. A. Sloane, Aug 13 2014
a(A000069(n)) = -1; a(A001969(n)) = +1. - Reinhard Zumkeller, Apr 29 2012
Partial sums of every third terms give A005599. - Reinhard Zumkeller, May 26 2013
Fixed point of the morphism 1 --> 1,-1 and -1 --> -1,1. - Robert G. Wilson v, Apr 07 2014
Fibbinary numbers (A003714) gives the numbers n for which a(n) = A132971(n). - Antti Karttunen, May 30 2017

Examples

			G.f. = 1 - x - x^2 + x^3 - x^4 + x^5 + x^6 - x^7 - x^8 + x^9 + x^10 + ...
The first 2^2 = 4 terms are 1, -1, -1, 1. Exchanging 1 and -1 gives -1, 1, 1, -1, which are a(4) through a(7). - _Michael B. Porter_, Jul 29 2016
		

Crossrefs

Convolution inverse of A018819.
Partial sums of A292118.

Programs

  • Haskell
    import Data.List (transpose)
    a106400 n = a106400_list !! n
    a106400_list =  1 : concat
       (transpose [map negate a106400_list, tail a106400_list])
    -- Reinhard Zumkeller, Apr 29 2012
    
  • Magma
    [1-2*(&+Intseq(n,2) mod(2)): n in [0..100]]; // Vincenzo Librandi, Sep 01 2015
    
  • Maple
    A106400 := proc(n)
            1-2*A010060(n) ;
    end proc: # R. J. Mathar, Jul 22 2012
    subs("0"=1,"1"=-1, StringTools:-Explode(StringTools:-ThueMorse(1000))); # Robert Israel, Sep 01 2015
    # third Maple program:
    a:= n-> (-1)^add(i, i=Bits[Split](n)):
    seq(a(n), n=0..120);  # Alois P. Heinz, Apr 13 2020
  • Mathematica
    tm[0] = 0; tm[n_?EvenQ] := tm[n/2]; tm[n_] := 1 - tm[(n-1)/2]; Table[(-1)^tm[n], {n, 0, 101}] (* Jean-François Alcover, Oct 24 2013 *)
    Nest[ Flatten[# /. {1 -> {1, -1}, -1 -> {-1, 1}}] &, {1}, 7] (* Robert G. Wilson v, Apr 07 2014 *)
    Table[Coefficient[Product[1 - x^(2^k), {k, 0, Log2[n + 1]}], x, n], {n, 0, 20}] (* Vladimir Reshetnikov, Nov 11 2016 *)
    (-1)^ThueMorse[Range[0,100]] (* Paolo Xausa, Dec 18 2023 *)
  • PARI
    {a(n) = if( n<1, n>=0, a(n\2) * (-1)^(n%2))};
    
  • PARI
    {a(n) = my(A, m); if( n<1, n==0, m=1; A = 1 + O(x); while( m<=n, m*=2; A = subst(A, x, x^2) * (1-x)); polcoeff(A, n))};
    
  • PARI
    a(n) = { 1 - 2 * (hammingweight(n) % 2) };  \\ Gheorghe Coserea, Aug 30 2015
    
  • PARI
    apply( {A106400(n)=(-1)^hammingweight(n)}, [0..99]) \\ M. F. Hasler, Feb 07 2020
    
  • Python
    def aupto(nn):
        A = [1]
        while len(A) < nn+1: A += [-i for i in A]
        return A[:nn+1]
    print(aupto(101)) # Michael S. Branicky, Jun 26 2022
    
  • Python
    def A106400(n): return -1 if n.bit_count()&1 else 1 # Chai Wah Wu, Mar 01 2023

Formula

a(n) = (-1)^A010060(n).
a(n) = (-1)^wt(n), where wt(n) is the binary weight of n, A000120(n).
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = v^3 - 2*u*v*w + u^2*w.
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u6*u1^3 - 3*u6*u2*u1^2 + 3*u6*u2^2*u1 - u3*u2^3.
Euler transform of sequence b(n) where b(2^k) = -1 and zero otherwise.
G.f.: Product_{k>=0} (1 - x^(2^k)) = A(x) = (1-x) * A(x^2).
a(n) = B_n(-A038712(1)*0!, ..., -A038712(n)*(n-1)!)/n!, where B_n(x_1, ..., x_n) is the n-th complete Bell polynomial. See the Wikipedia link for complete Bell polynomials , and A036040 for the coefficients of these partition polynomials. - Gevorg Hmayakyan, Jul 10 2016 (edited by - Wolfdieter Lang, Aug 31 2016)
a(n) = A008836(A005940(1+n)). [Analogous to Liouville's lambda] - Antti Karttunen, May 30 2017
a(n) = (-1)^A309303(n), see the closed form (5) in the MathWorld link. - Vladimir Reshetnikov, Jul 23 2019

A134278 A certain partition array in Abramowitz-Stegun order (A-St order), called M_3(6).

Original entry on oeis.org

1, 6, 1, 66, 18, 1, 1056, 264, 108, 36, 1, 22176, 5280, 3960, 660, 540, 60, 1, 576576, 133056, 95040, 43560, 15840, 23760, 3240, 1320, 1620, 90, 1, 17873856, 4036032, 2794176, 2439360, 465696, 665280, 304920, 249480, 36960, 83160, 22680, 2310
Offset: 1

Views

Author

Wolfdieter Lang, Nov 13 2007

Keywords

Comments

For the A-St order of partitions see the Abramowitz-Stegun reference given in A117506.
Partition number array M_3(6), the k=6 member in the family of a generalization of the multinomial number arrays M_3 = M_3(1) = A036040.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42, ...].
The S2(6,n,m):=A049385(n,m) numbers (generalized Stirling2 numbers) are obtained by summing in row n all numbers with the same part number m. In the same manner the S2(n,m) (Stirling2) numbers A008277 are obtained from the partition array M_3= A036040.
a(n,k) enumerates unordered forests of increasing 6-ary trees related to the k-th partition of n in the A-St order. The m-forest is composed of m such trees, with m the number of parts of the partition.

Examples

			[1]; [6,1]; [66,18,1]; [1056,264,108,36,1]; [22176,5280,3960,660,540,60,1]; ...
There are a(4,3) = 108 = 3*6^2 unordered 2-forests with 4 vertices, composed of two 6-ary increasing trees, each with two vertices: there are 3 increasing labelings (1,2)(3,4); (1,3)(2,4); (1,4)(2,3) and each tree comes in six versions from the 6-ary structure.
		

Crossrefs

Cf. A049412 (row sums, also of triangle A049385).
Cf. A134273 (M_3(5) partition array).

Formula

a(n,k) = n!*Product_{j=1..n} (S2(6,j,1)/j!)^e(n,k,j)/e(n,k,j)! with S2(6,n,1) = A049385(n,1) = A008548(n) = (5*n-4)(!^5) (quintuple- or 5-factorials) and the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n. Exponents 0 can be omitted due to 0!=1.

A060540 Square array read by antidiagonals downwards: T(n,k) = (n*k)!/(k!^n*n!), (n>=1, k>=1), the number of ways of dividing nk labeled items into n unlabeled boxes with k items in each box.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 10, 15, 1, 1, 35, 280, 105, 1, 1, 126, 5775, 15400, 945, 1, 1, 462, 126126, 2627625, 1401400, 10395, 1, 1, 1716, 2858856, 488864376, 2546168625, 190590400, 135135, 1, 1, 6435, 66512160, 96197645544, 5194672859376, 4509264634875, 36212176000, 2027025, 1
Offset: 1

Views

Author

Henry Bottomley, Apr 02 2001

Keywords

Comments

The Copeland link gives the associations of this entry with the operator calculus of Appell Sheffer polynomials, the combinatorics of simple set partitions encoded in the Faa di Bruno formula for composition of analytic functions (formal Taylor series), the Pascal matrix, and the geometry of the n-dimensional simplices (hypertriangles, or hypertetrahedra). These, in turn, are related to simple instances of the application of the exponential formula / principle / schema giving the number of not-necessarily-connected objects composed from an ensemble of connected objects. - Tom Copeland, Jun 09 2021

Examples

			Array begins:
  1,   1,       1,          1,             1,                 1, ...
  1,   3,      10,         35,           126,               462, ...
  1,  15,     280,       5775,        126126,           2858856, ...
  1, 105,   15400,    2627625,     488864376,       96197645544, ...
  1, 945, 1401400, 2546168625, 5194672859376, 11423951396577720, ...
  ...
		

Crossrefs

Main diagonal is A057599.
Related to A057599, see also A096126 and A246048.
Cf. A060358, A361948 (includes row/col 0).
Cf. A000217, A000292, A000332, A000389, A000579, A000580, A007318, A036040, A099174, A133314, A132440, A135278 (associations in Copeland link).

Programs

  • Mathematica
    T[n_, k_] := (n*k)!/(k!^n*n!);
    Table[T[n-k+1, k], {n, 1, 10}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Jun 29 2018 *)
  • PARI
    { i=0; for (m=1, 20, for (n=1, m, k=m - n + 1; write("b060540.txt", i++, " ", (n*k)!/(k!^n*n!))); ) } \\ Harry J. Smith, Jul 06 2009

Formula

T(n,k) = (n*k)!/(k!^n*n!) = T(n-1,k)*A060543(n,k) = A060538(n,k)/k!.
T(n,k) = Product_{j=2..n} binomial(j*k-1,k-1). - M. F. Hasler, Aug 22 2014

Extensions

Definition reworded by M. F. Hasler, Aug 23 2014

A117506 Irregular triangle read by rows: dimensions of the irreducible representations of the symmetric group S_n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 3, 2, 3, 1, 1, 4, 5, 6, 5, 4, 1, 1, 5, 9, 5, 10, 16, 5, 10, 9, 5, 1, 1, 6, 14, 14, 15, 35, 21, 21, 20, 35, 14, 15, 14, 6, 1, 1, 7, 20, 28, 14, 21, 64, 70, 56, 42, 35, 90, 56, 70, 14, 35, 64, 28, 21, 20, 7, 1
Offset: 0

Views

Author

Wolfdieter Lang, Apr 13 2006

Keywords

Comments

The n-th row has partition(n) = A000041(n) entries.
Also the numbers of standard Young tableaux for Young diagrams (or partitions).
Also "generalized" Catalan numbers. For a partition of n, n=(n_1+...+n_d), this is the number of integral lattice paths from (0,...,0) to (n_1,...,n_d) such that for any point p=(p_1,...p_d) on such a path p_i is never less than p_j whenever iGraham H. Hawkes, Jul 05 2013
The irreducible representations of S_n correspond to Young diagrams or partitions.
Partitions of n are ordered according to Abramowitz-Stegun (A-St) (see the reference, pp. 831-2). In contrast to A-St, a partition has nondecreasing parts (reverse notation of A-St).
The dimension of a representation of S_n corresponding to a Young diagram or partition is a(n,k) for the k-th partition of n in this A-St order.
One could call these numbers a(n,k) M_4 (similar to M_0, M_1, M_2, M_3 given in A111786, A036038, A036039, A036040, respectively).
From Wolfdieter Lang, Oct 09 2015: (Start)
The first formula given below appears in A. Young, Q.S.A. III, PLMS 28 (1928) 255-292 (third paper on "On Quantitative Substitutional Analysis"), Theorem II on p. 260, and he calls it f; see the collected papers (CP) reference, p. 357. Note the shorthand notation for the products; see Q.S.A. II, PLMS 34 (1902) 361-397, p. 366, CP, p. 97, for the explicit one.
This formula also can be found in the Glass-Ng link, Theorem 1, p. 702, using the Vandermonde determinant in the numerator and re-indexing the denominator.
The product of the hook length numbers, called H(n, k) in this formula below, is found in A263003(n, k).
The squared row entries sum to n!. See A. Young, Q.S.A. II (see above), pp. 367-368, CP pp. 98-99. Also Q.S.A. III, p. 265, CP p. 362.
(End)

Examples

			[1];
[1];
[1, 1];
[1, 2, 1];
[1, 3, 2, 3, 1];
[1, 4, 5, 6, 5, 4, 1];
[1, 5, 9, 5, 10, 16, 5, 10, 9, 5, 1];...
a(4,4)=3 because the 4th partition of n=4 in A-St order is [2,1,1],
and H(4,4)=(4!*2!*1!)/Vandermonde([4,2,1]) = (4!*2)/6 =4*2, hence
4!/H(4,4) = 3.
a(4,4)=3 because the hook lengths of the Young diagram of [2,1,1] are [4, 1; 2; 1], hence 4!/(4*1*2*1) = 3.
The sum of the squared entries of each row gives n!: n = 5: 2*(1^1 + 4^2 + 5^2) + 6^2 = 120 = 5!. - _Wolfdieter Lang_, Oct 09 2015
		

References

  • G. de B. Robinson (ed.), The Collected Papers of Alfred Young 1873-1940, University of Toronto Press, 1977.
  • G. B. Wybourne, Symmetry principles and atomic spectroscopy, Wiley, New York, 1970, p. 9.

Crossrefs

Cf. A000041, A000085 (row sums), A060240 (rows sorted), A263003.
Cf. A067924.

Programs

  • Maple
    h:= l-> (n-> mul(mul(1+l[i]-j+add(`if`(l[k]>=j, 1, 0),
                 k=i+1..n), j=1..l[i]), i=1..n))(nops(l)):
    g:= (n, i, l)-> `if`(n=0 or i=1, [h([l[], 1$n])],
        [g(n, i-1, l)[], g(n-i, min(n-i, i), [l[], i])[]]):
    T:= n-> map(x-> n!/x, g(n$2, []))[]:
    seq(T(n), n=0..10);  # Alois P. Heinz, Nov 05 2015
  • Mathematica
    h[l_List] := Function[n, Product[Product[1 + l[[i]] - j + Sum[If[l[[k]] >= j, 1, 0], {k, i+1, n}], {j, 1, l[[i]]}], {i, 1, n}]][Length[l]]; g[n_, i_, l_List] := If[n==0 || i==1, Join[{h[Join[l, Array[1&, n]]]}], If[i<1, {}, Join[{g[n, i-1, l]}, If[i>n, {}, g[n-i, i, Join[l, {i}]]]]]] // Flatten; T[n_] := n!/ g[n, n, {}]; Table[T[n], {n, 0, 10}] // Flatten (* Jean-François Alcover, Dec 19 2015, after Alois P. Heinz *)

Formula

a(n,k) = n!/H(n,k) with H(n,k):= Product_{i=1..m(n,k)} (x_i)!/Det(x_i^(m(n,k)-j)) with the Vandermonde determinant for the variables x_i:=lambda(n,k)_i + m(n,k)-i, i,j=1..m(n,k) if m(n,k) is the number of parts of the k-th partition of n, called lambda(n,k), in the A-St order (see above). Lambda(n,k)_i denotes the i-th part of the partition lambda(n,k), sorted in decreasing order (this is the reverse of the A-St notation).
a(n,k) = n!/Product_{j=1..n}(h(n,k,j) with the hook numbers h(n,k,j) of the Young diagram of the partition lambda(n,k) in the A-St order. See the link for 'hook length formula'.

Extensions

Row n=0 prepended by Alois P. Heinz, Nov 05 2015

A263916 Coefficients of the Faber partition polynomials.

Original entry on oeis.org

-1, -2, 1, -3, 3, -1, -4, 4, 2, -4, 1, -5, 5, 5, -5, -5, 5, -1, -6, 6, 6, -6, 3, -12, 6, -2, 9, -6, 1, -7, 7, 7, -7, 7, -14, 7, -7, -7, 21, -7, 7, -14, 7, -1, -8, 8, 8, -8, 8, -16, 8, 4, -16, -8, 24, -8, -8, 12, 24, -32, 8, 2, -16, 20, -8, 1
Offset: 1

Views

Author

Tom Copeland, Oct 29 2015

Keywords

Comments

The coefficients of the Faber polynomials F(n,b(1),b(2),...,b(n)) (Bouali, p. 52) in the order of the partitions of Abramowitz and Stegun. Compare with A115131 and A210258.
These polynomials occur in discussions of the Virasoro algebra, univalent function spaces and the Schwarzian derivative, symmetric functions, and free probability theory. They are intimately related to symmetric functions, free probability, and Appell sequences through the raising operator R = x - d log(H(D))/dD for the Appell sequence inverse pair associated to the e.g.f.s H(t)e^(xt) (cf. A094587) and (1/H(t))e^(xt) with H(0)=1.
Instances of the Faber polynomials occur in discussions of modular invariants and modular functions in the papers by Asai, Kaneko, and Ninomiya, by Ono and Rolen, and by Zagier. - Tom Copeland, Aug 13 2019
The Faber polynomials, denoted by s_n(a(t)) where a(t) is a formal power series defined by a product formula, are implicitly defined by equation 13.4 on p. 62 of Hazewinkel so as to extract the power sums of the reciprocals of the zeros of a(t). This is the Newton identity expressing the power sum symmetric polynomials in terms of the elementary symmetric polynomials/functions. - Tom Copeland, Jun 06 2020
From Tom Copeland, Oct 15 2020: (Start)
With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an
A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!
B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n
C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.
Expansions of log(f(x)) are given in
I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials
II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.
Expansions of exp(f(x)-1) are given in
III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind
IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials
V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.
Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)

Examples

			F(1,b1) = - b1
F(2,b1,b2) = -2 b2 + b1^2
F(3,b1,b2,b3) = -3 b3 + 3 b1 b2 - b1^3
F(4,b1,...) = -4 b4 + 4 b1 b3 + 2 b2^2  - 4 b1^2 b2 + b1^4
F(5,...) = -5 b5 + 5 b1 b4 + 5 b2 b3 - 5 b1^2 b3 - 5 b1 b2^2 + 5 b1^3 b2 - b1^5
------------------------------
IF(1,b1) = -b1
IF(2,b1,,b2) = -b2 + b1^2
IF(3,b1,b2,b3) = -2 b3 + 3 b1 b2 - b1^3
IF(4,b1,...) = -6 b4 + 8 b1 b3 + 3 b2^2  - 6 b1^2 b2 + b1^4
IF(5,...) = -24 b5 + 30 b1 b4 + 20 b2 b3 - 20 b1^2 b3 - 15 b1 b2^2 + 10 b1^3 b2 - b1^5
------------------------------
For 1/(1+x)^2 = 1- 2x + 3x^2 - 4x^3 + 5x^4 - ..., F(n,-2,3,-4,...) = (-1)^(n+1) 2.
------------------------------
F(n,x,2x,...,nx), F(n,-x,2x,-3x,...,(-1)^n n*x), and F(n,(2-x),1,0,0,...) are related to the Chebyshev polynomials through A127677 and A111125. See also A110162, A156308, A208513, A217476, and A220668.
------------------------------
For b1 = p, b2 = q, and all other indeterminates 0, see A113279 and A034807.
For b1 = -y, b2 = 1 and all other indeterminates 0, see A127672.
		

References

  • H. Airault, "Symmetric sums associated to the factorization of Grunsky coefficients," in Groups and Symmetries: From Neolithic Scots to John McKay, CRM Proceedings and Lecture Notes: Vol. 47, edited by J. Harnad and P. Winternitz, American Mathematical Society, 2009.
  • D. Bleeker and B. Booss, Index Theory with Applications to Mathematics and Physics, International Press, 2013, (see section 16.7 Characteristic Classes and Curvature).
  • M. Hazewinkel, Formal Groups and Applications, Academic Press, New York San Francisco London, 1978, p. 120.
  • F. Hirzebruch, Topological methods in algebraic geometry. Second, corrected printing of the third edition. Die Grundlehren der Mathematischen Wissenschaften, Band 131 Springer-Verlag, Berlin Heidelberg New York, 1978, p. 11 and 92.
  • D. Knutson, λ-Rings and the Representation Theory of the Symmetric Group, Lect. Notes in Math. 308, Springer-Verlag, 1973, p. 35.
  • D. Yau, Lambda-Rings, World Scientific Publishing Co., Singapore, 2010, p. 45.

Crossrefs

Programs

  • Mathematica
    F[0] = 1; F[1] = -b[1]; F[2] = b[1]^2 - 2 b[2]; F[n_] := F[n] = -b[1] F[n - 1] - Sum[b[n - k] F[k], {k, 1, n - 2}] - n b[n] // Expand;
    row[n_] := (List @@ F[n]) /. b[_] -> 1 // Reverse;
    Table[row[n], {n, 1, 8}] // Flatten // Rest (* Jean-François Alcover, Jun 12 2017 *)

Formula

-log(1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n.
-d(1 + b(1) x + b(2) x^2 + ...)/dx / (1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) x^(n-1).
F(n,b(1),...,b(n)) = -n*b(n) - Sum_{k=1..n-1} b(n-k)*F(k,b(1),...,b(k)).
Umbrally, with B(x) = 1 + b(1) x + b(2) x^2 + ..., B(x) = exp[log(1-F.x)] and 1/B(x) = exp[-log(1-F.x)], establishing a connection to the e.g.f. of A036039 and the symmetric polynomials.
The Stirling partition polynomials of the first kind St1(n,b1,b2,...,bn;-1) = IF(n,b1,b2,...,bn) (cf. the Copeland link Lagrange a la Lah, signed A036039, and p. 184 of Airault and Bouali), i.e., the cyclic partition polynomials for the symmetric groups, and the Faber polynomials form an inverse pair for isolating the indeterminates in their definition, for example, F(3,IF(1,b1),IF(2,b1,b2)/2!,IF(3,b1,b2,b3)/3!)= b3, with bk = b(k), and IF(3,F(1,b1),F(2,b1,b2),F(3,b1,b2,b3))/3!= b3.
The polynomials specialize to F(n,t,t,...) = (1-t)^n - 1.
See Newton Identities on Wikipedia on relation between the power sum symmetric polynomials and the complete homogeneous and elementary symmetric polynomials for an expression in multinomials for the coefficients of the Faber polynomials.
(n-1)! F(n,x[1],x[2]/2!,...,x[n]/n!) = - p_n(x[1],...,x[n]), where p_n are the cumulants of A127671 expressed in terms of the moments x[n]. - Tom Copeland, Nov 17 2015
-(n-1)! F(n,B(1,x[1]),B(2,x[1],x[2])/2!,...,B(n,x[1],...,x[n])/n!) = x[n] provides an extraction of the indeterminates of the complete Bell partition polynomials B(n,x[1],...,x[n]) of A036040. Conversely, IF(n,-x[1],-x[2],-x[3]/2!,...,-x[n]/(n-1)!) = B(n,x[1],...,x[n]). - Tom Copeland, Nov 29 2015
For a square matrix M, determinant(I - x M) = exp[-Sum_{k>0} (trace(M^k) x^k / k)] = Sum_{n>0} [ P_n(-trace(M),-trace(M^2),...,-trace(M^n)) x^n/n! ] = 1 + Sum_{n>0} (d[n] x^n), where P_n(x[1],...,x[n]) are the cycle index partition polynomials of A036039 and d[n] = P_n(-trace(M),-trace(M^2),...,-trace(M^n)) / n!. Umbrally, det(I - x M)= exp[log(1 - b. x)] = exp[P.(-b_1,..,-b_n)x] = 1 / (1-d.x), where b_k = tr(M^k). Then F(n,d[1],...,d[n]) = tr[M^n]. - Tom Copeland, Dec 04 2015
Given f(x) = -log(g(x)) = -log(1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n, action on u_n = F(n,b(1),...,b(n)) with A133932 gives the compositional inverse finv(x) of f(x), with F(1,b(1)) not equal to zero, and f(g(finv(x))) = f(e^(-x)). Note also that exp(f(x)) = 1 / g(x) = exp[Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n] implies relations among A036040, A133314, A036039, and the Faber polynomials. - Tom Copeland, Dec 16 2015
The Dress and Siebeneicher paper gives combinatorial interpretations and various relations that the Faber polynomials must satisfy for integral values of its arguments. E.g., Eqn. (1.2) p. 2 implies [2 * F(1,-1) + F(2,-1,b2) + F(4,-1,b2,b3,b4)] mod(4) = 0. This equation implies that [F(n,b1,b2,...,bn)-(-b1)^n] mod(n) = 0 for n prime. - Tom Copeland, Feb 01 2016
With the elementary Schur polynomials S(n,a_1,a_2,...,a_n) = Lah(n,a_1,a_2,...,a_n) / n!, where Lah(n,...) are the refined Lah polynomials of A130561, F(n,S(1,a_1),S(2,a_1,a_2),...,S(n,a_1,...,a_n)) = -n * a_n since sum_{n > 0} a_n x^n = log[sum{n >= 0} S(n,a_1,...,a_n) x^n]. Conversely, S(n,-F(1,a_1),-F(2,a_1,a_2)/2,...,-F(n,a_1,...,a_n)/n) = a_n. - Tom Copeland, Sep 07 2016
See Corollary 3.1.3 on p. 38 of Ardila and Copeland's two MathOverflow links to relate the Faber polynomials, with arguments being the signed elementary symmetric polynomials, to the logarithm of determinants, traces of powers of an adjacency matrix, and number of walks on graphs. - Tom Copeland, Jan 02 2017
The umbral inverse polynomials IF appear on p. 19 of Konopelchenko as partial differential operators. - Tom Copeland, Nov 19 2018

Extensions

More terms from Jean-François Alcover, Jun 12 2017

A055302 Triangle of number of labeled rooted trees with n nodes and k leaves, n >= 1, 1 <= k <= n.

Original entry on oeis.org

1, 2, 0, 6, 3, 0, 24, 36, 4, 0, 120, 360, 140, 5, 0, 720, 3600, 3000, 450, 6, 0, 5040, 37800, 54600, 18900, 1302, 7, 0, 40320, 423360, 940800, 588000, 101136, 3528, 8, 0, 362880, 5080320, 16087680, 15876000, 5143824, 486864, 9144, 9, 0, 3628800
Offset: 1

Views

Author

Christian G. Bower, May 11 2000

Keywords

Comments

Beginning with the second row, dividing each row by n gives the mirror of row n-1 of A141618. Under the exponential transform, the mirror of A141618 is generated, relating the number of connected graphs here to the number of disconnected graphs associated with A141618 (cf. A127671 and A036040). - Tom Copeland, Oct 25 2014

Examples

			Triangle begins
     1,
     2,     0;
     6,     3,     0;
    24,    36,     4,     0;
   120,   360,   140,     5,    0;
   720,  3600,  3000,   450,    6, 0;
  5040, 37800, 54600, 18900, 1302, 7, 0;
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 313.

Crossrefs

Row sums give A000169. Columns 1 through 12: A000142, A055303-A055313. Cf. A055314.
Cf. A248120 for a natural refinement.

Programs

  • Maple
    T:= (n, k)-> (n!/k!)*Stirling2(n-1, n-k):
    seq(seq(T(n, k), k=1..n), n=1..10);  # Alois P. Heinz, Nov 13 2013
  • Mathematica
    Table[Table[n!/k! StirlingS2[n-1,n-k], {k,1,n}], {n,0,10}]//Grid  (* Geoffrey Critzer, Dec 01 2012 *)
  • PARI
    A055302(n,k)=n!/k!*stirling(n-1, n-k,2);
    for(n=1,10,for(k=1,n,print1(A055302(n,k),", "));print());
    \\ Joerg Arndt, Oct 27 2014

Formula

E.g.f. (relative to x) satisfies: A(x,y) = xy + x*exp(A(x,y)) - x. Divides by n and shifts up under exponential transform.
T(n,k) = (n!/k!)*Stirling2(n-1, n-k). - Vladeta Jovovic, Jan 28 2004
T(n,k) = A055314(n,k)*(n-k) + A055314(n,k+1)*(k+1). The first term is the number of such trees with root degree > 1 while the second term is the number of such trees with root degree = 1. This simplifies to the above formula by Vladeta Jovovic. - Geoffrey Critzer, Dec 01 2012
E.g.f.: G(x,t) = log[1 + t * N(x*t,1/t)], where N(x,t) is the e.g.f. of A141618. Also, G(x*t,1/t)= log[1 + N(x,t)/t] is the comp. inverse in x of x / [1 + t * (e^x - 1)]. - Tom Copeland, Oct 26 2014
Previous Showing 21-30 of 102 results. Next