cp's OEIS Frontend

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

Showing 1-6 of 6 results.

A036039 Irregular triangle of multinomial coefficients of integer partitions read by rows (in Abramowitz and Stegun ordering) giving the coefficients of the cycle index polynomials for the symmetric groups S_n.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 6, 8, 3, 6, 1, 24, 30, 20, 20, 15, 10, 1, 120, 144, 90, 40, 90, 120, 15, 40, 45, 15, 1, 720, 840, 504, 420, 504, 630, 280, 210, 210, 420, 105, 70, 105, 21, 1, 5040, 5760, 3360, 2688, 1260, 3360, 4032, 3360, 1260, 1120, 1344, 2520, 1120, 1680, 105, 420, 1120, 420, 112, 210, 28, 1
Offset: 1

Views

Author

Keywords

Comments

The sequence of row lengths is A000041(n), n >= 1 (partition numbers).
Number of permutations whose cycle structure is the given partition. Row sums are factorials (A000142). - Franklin T. Adams-Watters, Jan 12 2006
A relation between partition polynomials formed from these "refined" Stirling numbers of the first kind and umbral operator trees and Lagrange inversion is presented in the link "Lagrange a la Lah".
These cycle index polynomials for the symmetric group S_n are also related to a raising operator / infinitesimal generator for fractional integro-derivatives, involving the digamma function and the Riemann zeta function values at positive integers, and to the characteristic polynomial for the adjacency matrix of complete n-graphs A055137 (cf. MathOverflow link). - Tom Copeland, Nov 03 2012
In the Lang link, replace all x(n) by t to obtain A132393. Furthermore replace x(1) by t and all other x(n) by 1 to obtain A008290. See A274760. - Tom Copeland, Nov 06 2012, Oct 29 2015 - corrected by Johannes W. Meijer, Jul 28 2016
The umbral compositional inverses of these polynomials are formed by negating the indeterminates x(n) for n>1, i.e., P(n,P(.,x(1),-x(2),-x(3),...),x(2),x(3),...) = x(1)^n (cf. A130561 for an example of umbral compositional inversion). The polynomials are an Appell sequence in x(1), i.e., dP(n,x(1))/dx(1) = n P(n-1, x(1)) and (P(.,x)+y)^n=P(n,x+y) umbrally, with P(0,x(1))=1. - Tom Copeland, Nov 14 2014
Regarded as the coefficients of the partition polynomials listed by Lang, a signed version of these polynomials IF(n,b1,b2,...,bn) (n! times polynomial on page 184 of Airault and Bouali) provides an inversion of the Faber polynomials F(n,b1,b2,...,bn) (page 52 of Bouali, A263916, and A115131). For example, F(3, IF(1,b1), IF(2,b1,b2)/2!, IF(3,b1,b2,b3)/3!) = b3 and IF(3, F(1,b1), F(2,b1,b2), F(3,b1,b2,b3))/3! = b3 with F(1,b1) = -b1. (Compare with A263634.) - Tom Copeland, Oct 28 2015; Sep 09 2016
The e.g.f. for the row partition polynomials is Sum_{n>=0} P_n(b_1,...,b_n) x^n/n! = exp[Sum_{n>=1} b_n x^n/n], or, exp[P.(b_1,...,b_n)x] = exp[-], expressed umbrally with <"power series"> denoting umbral evaluation (b.)^n = b_n within the power series. This e.g.f. is central to the paper by Maxim and Schuermannn on characteristic classes (cf. Friedrich and McKay also). - Tom Copeland, Nov 11 2015
The elementary Schur polynomials are given by S(n,x(1),x(2),...,x(n)) = P(n,x(1), 2*x(2),...,n*x(n)) / n!. See p. 12 of Carrell. - Tom Copeland, Feb 06 2016
These partition polynomials are also related to the Casimir invariants associated to quantum density states on p. 3 of Boya and Dixit and pp. 5 and 6 of Byrd and Khaneja. - Tom Copeland, Jul 24 2017
With the indeterminates (x_1,x_2,x_3,...) = (t,-c_2*t,-c_3*t,...) with c_n >0, umbrally P(n,a.) = P(n,t)|{t^n = a_n} = 0 and P(j,a.)P(k,a.) = P(j,t)P(k,t)|{t^n =a_n} = d_{j,k} >= 0 is the coefficient of x^j/j!*y^k/k! in the Taylor series expansion of the formal group law FGL(x,y) = f[f^{-1}(x)+f^{-1}(y)], where a_n are the inversion partition polynomials for calculating f(x) from the coefficients of the series expansion of f^{-1}(x) given in A133932. - Tom Copeland, Feb 09 2018
For relation to the Witt symmetric functions, as well as the basic power, elementary, and complete symmetric functions, see the Borger link p. 295. For relations to diverse zeta functions, determinants, and paths on graphs, see the MathOverflow question Cycling Through the Zeta Garden. - Tom Copeland, Mar 25 2018
Chmutov et al. identify the partition polynomials of this entry with the one-part Schur polynomials and assert that any linear combination with constant coefficients of these polynomials is a tau function for the KP hierarchy. - Tom Copeland, Apr 05 2018
With the indeterminates in the partition polynomials assigned as generalized harmonic numbers, i.e., as partial sums of the Dirichlet series for the Riemann zeta function, zeta(n), for integer n > 1, sums of simple normalizations of these polynomials give either unity or simple sums of consecutive zeta(n) (cf. Hoffman). Other identities involving these polynomials can be found in the Choi reference in Hoffman's paper. - Tom Copeland, Oct 05 2019
On p. 39 of Ma Luo's thesis is the e.g.f. of rational functions r_n obtained through the (umbral) formula 1/(1-r.T) = exp[log(1+P.T)], a differently signed e.g.f. of this entry, where (P.)^n = P_n are Eisenstein elliptic functions. P. 38 gives the example of 4! * r_4 as the signed 4th row partition polynomial of this entry. This series is equated through a simple proportionality factor to the Zagier Jacobi form on p. 25. Recurrence relations for the P_n are given on p. 24 involving the normalized k-weight Eisenstein series G_k introduced on p. 23 and related to the Bernoulli numbers. - Tom Copeland, Oct 16 2019
The Chern characteristic classes or forms of complex vector bundles and the characteristic polynomials of curvature forms for a smooth manifold can be expressed in terms of this entry's partition polynomials with the associated traces, or power sum polynomials, as the indeterminates. The Chern character is the e.g.f. of these traces and so its coefficients are given by the Faber polynomials with this entry's partition polynomials as the indeterminates. See the Mathoverflow question "A canonical reference for Chern characteristic classes". - Tom Copeland, Nov 04 2019
For an application to the physics of charged fermions in an external field, see Figueroa et al. - Tom Copeland, Dec 05 2019
Konopelchenko, in Proposition 5.2, p. 19, defines an operator P_k that is a differently signed operator version of the partition polynomials of this entry divided by a factorial. These operators give rise to bilinear Hirota equations for the KP hierarchy. These partition polynomials are also presented in Hopf algebras of symmetric functions by Cartier. - Tom Copeland, Dec 18 2019
For relationship of these partition polynomials to calculations of Pontryagin classes and the Riemann xi function, see A231846. - Tom Copeland, May 27 2020
Luest and Skliros summarize on p. 298 many of the properties of the cycle index polynomials given here; and Bianchi and Firrotta, a few on p. 6. - Tom Copeland, Oct 15 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

			The partition array T(n, k) begins (see the W. Lang link for rows 1..10):
  n\k   1    2    3    4    5    6    7    8    9   10   11  12   13  14 15 ...
  1:    1
  2:    1    1
  3:    2    3    1
  4:    6    8    3    6    1
  5:   24   30   20   20   15   10    1
  6:  120  144   90   40   90  120   15   40   45   15    1
  7:  720  840  504  420  504  630  280  210  210  420  105  70  105  21  1
... reformatted by _Wolfdieter Lang_, May 25 2019
		

References

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

Crossrefs

Cf. other versions based on different partition orderings: A102189 (rows reversed), A181897, A319192.
Cf. A133932.
Cf. A231846.
Cf. A127671.

Programs

  • Maple
    nmax:=7: 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 sA036039(n, m) := n!/ (mul((t)^q(t)*q(t)!, t=1..n)); od: od: seq(seq(A036039(n, m), m=1..numbpart(n)), n=1..nmax); # Johannes W. Meijer, Jul 14 2016
    # 2nd program:
    A036039 := proc(n,k)
        local a,prts,e,ai ;
        a := n! ;
        # ASPrts is implemented in A119441
        prts := ASPrts(n)[k] ;
        ai := 1;
        for e from 1 to nops(prts) do
            if e>1 then
                if op(e,prts) = op(e-1,prts) then
                    ai := ai+1 ;
                else
                    ai := 1;
                end if;
            end if;
            a := a/(op(e,prts)*ai) ;
        end do:
        a ;
    end proc:
    seq(seq(A036039(n,k),k=1..combinat[numbpart](n)),n=1..15) ; # R. J. Mathar, Dec 18 2016
  • Mathematica
    aspartitions[n_]:=Reverse/@Sort[Sort/@IntegerPartitions[n]];(* Abramowitz & Stegun ordering *);
    ascycleclasses[n_Integer]:=n!/(Times@@ #)&/@((#!
    Range[n]^#)&/@Function[par,Count[par,# ]&/@Range[n]]/@aspartitions[n])
    (* The function "ascycleclasses" is then identical with A&S multinomial M2. *)
    Table[ascycleclasses[n], {n, 1, 8}] // Flatten
    (* Wouter Meeussen, Jun 26 2009, Jun 27 2009 *)
  • Sage
    def PartAS(n):
        P = []
        for k in (1..n):
            Q = [p.to_list() for p in Partitions(n, length=k)]
            for q in Q: q.reverse()
            P = P + sorted(Q)
        return P
    def A036039_row(n):
        fn, C = factorial(n), []
        for q in PartAS(n):
            q.reverse()
            p = Partition(q)
            fp = 1; pf = 1
            for a, c in p.to_exp_dict().items():
                fp *= factorial(c)
                pf *= factorial(a)**c
            co = fn//(fp*pf)
            C.append(co*prod([factorial(i-1) for i in p]))
        return C
    for n in (1..10):
        print(A036039_row(n)) # Peter Luschny, Dec 18 2016

Formula

T(n,k) = n!/Product_{j=1..n} j^a(n,k,j)*a(n,k,j)!, with the k-th partition of n >= 1 in Abromowitz-Stegun order written as Product_{j=1..n} j^a(n,k,j) with nonnegative integers a(n,k,j) satisfying Sum_{j=1..n} j*a(n,k,j) = n, and the number of parts is Sum_{j=1..n} a(n,k,j) =: m(n,k). - Wolfdieter Lang, May 25 2019
Raising and lowering operators are given for the partition polynomials formed from this sequence in the link in "Lagrange a la Lah Part I" on p. 23. - Tom Copeland, Sep 18 2011
From Szabo p. 34, with b_n = q^n / (1-q^n)^2, the partition polynomials give an expansion of the MacMahon function M(q) = Product_{n>=1} 1/(1-q^n)^n = Sum_{n>=0} PL(n) q^n, the generating function for PL(n) = n! P_n(b_1,...,b_n), the number of plane partitions with sum n. - Tom Copeland, Nov 11 2015
From Tom Copeland, Nov 18 2015: (Start)
The partition polynomials of A036040 are obtained by substituting x[n]/(n-1)! for x[n] in the partition polynomials of this entry.
CIP_n(t-F(1,b1),-F(2,b1,b2),...,-F(n,b1,...,bn)) = P_n(b1,...,bn;t), where CIP_n are the partition polynomials of this entry; F(n,...), those of A263916; and P_n, those defined in my formula in A094587, e.g., P_2(b1,b2;t) = 2 b2 + 2 b1 t + t^2.
CIP_n(-F(1,b1),-F(2,b1,b2),...,-F(n,b1,...,bn)) = n! bn. (End)
From the relation to the elementary Schur polynomials given in A130561 and above, the partition polynomials of this array satisfy (d/d(x_m)) P(n,x_1,...,x_n) = (1/m) * (n!/(n-m)!) * P(n-m,x_1,...,x_(n-m)) with P(k,...) = 0 for k<0. - Tom Copeland, Sep 07 2016
Regarded as Appell polynomials in the indeterminate x(1)=u, the partition polynomials of this entry P_n(u) obey d/du P_n(u) = n * P_{n-1}(u), so the abscissas for the zeros of P_n(u) are the same as those of the extrema of P{n+1}(u). In addition, the coefficient of u^{n-1} in P_{n}(u) is zero since these polynomials are related to the characteristic polynomials of matrices with null main diagonals, and, therefore, the trace is zero, further implying the abscissa for any zero is the negative of the sum of the abscissas of the remaining zeros. This assumes all zeros are distinct and real. - Tom Copeland, Nov 10 2019

Extensions

More terms from David W. Wilson
Title expanded by Tom Copeland, Oct 15 2020

A059171 Size of largest conjugacy class in S_n, the symmetric group on n symbols.

Original entry on oeis.org

1, 1, 3, 8, 30, 144, 840, 5760, 45360, 403200, 3991680, 43545600, 518918400, 6706022400, 93405312000, 1394852659200, 22230464256000, 376610217984000, 6758061133824000, 128047474114560000, 2554547108585472000, 53523844179886080000, 1175091669949317120000
Offset: 1

Views

Author

Des MacHale, Feb 14 2001

Keywords

Comments

Apart from initial terms, same as A001048. The number a(n) is the maximum of row n in the triangle of refined rencontres numbers A181897. - Tilman Piesk, Apr 02 2012

Examples

			a(3) = 3 because the largest conjugacy class in S_3 consists of the three 2-cycles {(12),(13),(23)}.
G.f. = x + x^2 + 3*x^3 + 8*x^4 + 30*x^5 + 144*x^6 + 840*x^7 + 5760*x^8 + ...
		

Crossrefs

Programs

  • Magma
    [1, 1] cat [n*Factorial(n-2): n in [3..25]]; // Vincenzo Librandi, Oct 26 2011
    
  • Maple
    a := proc(n) if n<=2 then RETURN(1) else RETURN(n*(n-2)!) fi: end:for n from 1 to 40 do printf(`%d,`,a(n)) od:
  • Mathematica
    Join[{1,1},Table[n (n-2)!,{n,3,30}]] (* Harvey P. Dale, Oct 25 2011 *)
    a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ x - x^2/2 - x Log[1 - x], {x, 0, n}]]; (* Michael Somos, Aug 26 2015 *)
    a[ n_] := With[ {m = n - 1}, If[ m < 0, 0, m! SeriesCoefficient[ -Log[1 - x] - x + 1/(1 - x), {x, 0, m}]]]; (* Michael Somos, Aug 26 2015 *)
  • PARI
    Vec(1+x*serlaplace((1+x-x^2)/(1-x)^2+O(x^66))) \\ Joerg Arndt, Mar 28 2013
    
  • PARI
    a(n)=if(n<=1,1,n!/(n-1)); \\ Joerg Arndt, Mar 28 2013

Formula

a(1) = a(2) = 1; a(n) = n*(n-2)! = (n!)/(n-1) for n > 2. This is the number of (n-1)-cycles in S_n.
E.g.f.: -log(1-x) - x + 1/(1-x). [for a(n+1) - Michael Somos, Aug 26 2015]
E.g.f.: x - x^2/2 - x*log(1-x). - Michael Somos, Aug 26 2015
The sequence 1, 3, 8, ... has e.g.f. (1+x-x^2)/(1-x)^2 and a(n) = n!(n+2-0^n) = n!*A065475(n). - Paul Barry, May 14 2004
E.g.f.: E(0) - x, where E(k) = 1 + x/(k+1)/(1 - 1/(1 + 1/(k+1)/E(k+1))); (continued fraction). - Sergei N. Gladkovskii, Mar 27 2013
G.f.: 1 + x/Q(0), where Q(k)= 1 - x/(1+x) - x/(1+x)*(k+2)/(1 - x/(1+x)*(k+1)/Q(k+1)); (continued fraction). - Sergei N. Gladkovskii, Apr 22 2013
From Amiram Eldar, Jan 22 2023: (Start)
Sum_{n>=1} 1/a(n) = 5/2.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2/e - 1/2. (End)

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Fabian Rothelius and James A. Sellers, Feb 15 2001

A073906 Number of conjugacy classes in the symmetric group S_n with distinct cardinality.

Original entry on oeis.org

1, 1, 3, 4, 6, 7, 11, 16, 23, 30, 40, 58, 69, 95, 119, 151, 184, 240, 297, 361, 452, 554, 663, 817, 980, 1177, 1402, 1665, 1995, 2346, 2774, 3259, 3837, 4466, 5222, 6057, 7061, 8159, 9450, 10917, 12533, 14408, 16570, 18958, 21623, 24681, 28123, 32000, 36232
Offset: 1

Views

Author

Vladeta Jovovic, Sep 03 2002

Keywords

Comments

a(n) is the number of distinct entries of row n in the triangle of refined rencontres numbers A181897. - Tilman Piesk, Apr 01 2012

Crossrefs

Programs

  • Mathematica
    Map[Length,Map[Union,Table[f[list_]:=Total[list]!/Apply[Times,list]/Apply[Times,Table[Count[list,n]!,{n,1,m}]];Map[f,Partitions[m]],{m,1,30}]]]  (* Geoffrey Critzer, Mar 25 2013 *)

Extensions

a(41)-a(49) from Alois P. Heinz, Feb 04 2014

A110141 Triangle, read by rows, where row n lists the denominators of unit fraction coefficients of the products of {c_k}, in ascending order by indices of {c_k}, in the coefficient of x^n in exp(Sum_{k>=1} c_k/k*x^k).

Original entry on oeis.org

1, 1, 2, 2, 6, 2, 3, 24, 4, 3, 8, 4, 120, 12, 6, 8, 4, 6, 5, 720, 48, 18, 16, 8, 6, 5, 48, 8, 18, 6, 5040, 240, 72, 48, 24, 12, 10, 48, 8, 18, 6, 24, 10, 12, 7, 40320, 1440, 360, 192, 96, 36, 30, 96, 16, 36, 12, 24, 10, 12, 7, 384, 32, 36, 12, 15, 32, 8, 362880, 10080, 2160, 960
Offset: 0

Views

Author

Paul D. Hanna, Jul 13 2005

Keywords

Comments

Row n starts with n!, after which the following pattern holds. When terms of row n are divided by a list of factorials, with (n-j-1)! repeated A002865(j+1) times in the list as j=1..n-1, the result is the initial terms of A110142. E.g., row 6 is: {720,48,18,16,8,6,5,48,8,18,6}; divide by respective factorials: {6!,4!,3!,2!,2!,1!,1!,0!,0!,0!,0!} with {4!,3!,2!,1!,0!} respectively occurring {1,1,2,2,4} times (A002865), yields the initial terms of A110142: {1,2,3,8,4,6,5,48,8,18,6}.
The term of the sequence corresponding to the product c_1^{n_1}c_2^{n_2}...c_k^{n_k} is equal to the number of elements in the centralizer of a permutation of n_1+2n_2+...+kn_k elements whose cycle type is 1^{n_1}2^{n_2}...k^{n^k}. (This fact is very standard, in particular, for the theory of symmetric functions.) - Vladimir Dotsenko, Apr 19 2009
Multiplying the values of row n by the corresponding values in row n of A102189, one obtains n!. - Jaimal Ichharam, Aug 06 2015
a(n,k) is the number of permutations in S_n that commute with a permutation having cycle type "k". Here, the cycle type of an n-permutation pi is the vector (i_1,...,i_n) where i_j is the number of cycles in pi of length j. These A000041(n) vectors can be ordered in reverse lexicographic order. The k-th cycle type is the k-th vector in this ordering. - Geoffrey Critzer, Jan 18 2019

Examples

			Coefficients [x^n] exp(c1*x + (c2/2)*x^2 + (c3/3)*x^3 + ...) begin:
[x^0]: 1;
[x^1]: 1*c1;
[x^2]: (1/2)*c1^2 + (1/2)*c2;
[x^3]: (1/6)*c1^3 + (1/2)*c1*c2 + (1/3)*c3;
[x^4]: (1/24)*c1^4 + (1/4)*c1^2*c2 + (1/3)*c1*c3 + (1/8)*c2^2 + (1/4)*c4;
[x^5]: (1/120)*c1^5 + (1/12)*c1^3*c2 + (1/6)*c1^2*c3 + (1/8)*c1*c2^2 + (1/4)*c1*c4 + (1/6)*c2*c3 + (1/5)*c5;
[x^6]: (1/720)*c1^6 + (1/48)*c1^4*c2 + (1/18)*c1^3*c3 + (1/16)*c1^2*c2^2 + (1/8)*c1^2*c4 + (1/6)*c1*c2*c3 + (1/5)*c1*c5 + (1/48)*c2^3 + (1/8)*c2*c4 + (1/18)*c3^2 + (1/6)*c6;
forming this triangle of unit fraction coefficients:
1;
1;
2,2;
6,2,3;
24,4,3,8,4;
120,12,6,8,4,6,5;
720,48,18,16,8,6,5,48,8,18,6;
5040,240,72,48,24,12,10,48,8,18,6,24,10,12,7;
40320,1440,360,192,96,36,30,96,16,36,12,24,10,12,7,384,32,36,12,15,32,8;
362880,10080,2160,960,480,144,120,288,48,108,36,48,20,24,14,384,32,36,12,15,32,8,144,40,24,14,162,18,20,9; ...
		

References

  • Macdonald, I. G. Symmetric functions and Hall polynomials. Oxford University Press, 1995. [From Vladimir Dotsenko, Apr 19 2009]

Crossrefs

Cf. A000041, A002865, A102189, A110142, A110143 (row sums).
First, second and third entries of each row are given (up to an offset) by A000142, A052849, and A052560 respectively. - Vladimir Dotsenko, Apr 19 2009

Programs

  • Mathematica
    Table[n!/CoefficientRules[n! CycleIndex[SymmetricGroup[n], s]][[All, 2]], {n, 1, 8}] // Grid (* Geoffrey Critzer, Jan 18 2019 *)

Formula

Number of terms in row n is A000041(n) (partition numbers). The unit fractions of each row sum to unity: Sum_{k=1..A000041(n)} 1/T(n, k) = 1.
a(n,k) = n!/A181897(n,k). - Geoffrey Critzer, Jan 18 2019

A198380 Cycle type of the n-th finite permutation represented by index number of A194602.

Original entry on oeis.org

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

Views

Author

Tilman Piesk, Oct 23 2011

Keywords

Comments

This sequence shows the cycle type of each finite permutation (A195663) as the index number of the corresponding partition. (When a permutation has a 3-cycle and a 2-cycle, this corresponds to the partition 3+2, etc.) Partitions can be ordered, so each partition can be denoted by its index in this order, e.g. 6 for the partition 3+2. Compare A194602.
From the properties of A194602 follows:
Entries 1,2,4,6,10,14,21... ( A000041(n)-1 from n=2 ) correspond to permutations with exactly one n-cycle (and no other cycles).
Entries 1,3,7,15,30,56,101... ( A000041(2n-1) from n=1 ) correspond to permutations with exactly n 2-cycles (and no other cycles), so these are the symmetric permutations.
Entries n = 1,3,4,7,9,10,12... ( A194602(n) has an even binary digit sum ) correspond to even permutations. This goes along with the fact, that a permutation is even when its partition contains an even number of even addends.
(Compare "Table for A194602" in section LINKS. Concerning the first two properties see especially the end of this file.)

Crossrefs

Cf. A195663, A195664, A055089 (ordered finite permutations).
Cf. A194602 (ordered partitions interpreted as binary numbers).
Cf. A181897 (number of n-permutations with cycle type k).

Extensions

Changed offset to 0 by Tilman Piesk, Jan 25 2012

A385081 Irregular triangle T(n,k) of refined derangement counts in the symmetric group S_(n+1), refined per cycle type.

Original entry on oeis.org

1, 2, 3, 6, 20, 24, 15, 90, 40, 120, 210, 504, 420, 720, 105, 1260, 1120, 3360, 2688, 1260, 5040, 2520, 9072, 15120, 25920, 2240, 20160, 18144, 40320, 945, 18900, 25200, 75600, 120960, 56700, 226800, 50400, 172800, 151200, 72576, 362880
Offset: 1

Views

Author

Gregory Gerard Wojnar, Jun 16 2025

Keywords

Comments

The triangle consists of selected entries from A181897.
Each permutation in the symmetric group S_N has a cycle type specified by an integer partition of N. We encode a partition of N as an N-tuple of multiplicities of j=1..N; e.g., the partition 8 = 1+1+2+4 is encoded as (2,1,0,1,0,0,0,0), abbreviated (2,1,0,1). In general (m_j: j=1..N), and this partition corresponds to a cycle type (a)(b)(cd)(efgh) in S_8. Derangements have no fixed points, hence correspond to partitions with no addends of 1; i.e., m_1=0. A181897 counts all permutations via cycle type (not just derangements), with the partitions ordered reverse lexicographically in terms of their multiplicities N-tuples. E.g., partitions of 4 are ordered: (4), (2,1), (1,0,1), (0,0,0,1). In fact, the columns of A181897 (as a triangular array) are scaled columns of binomial coefficients binomial(p,q) with q fixed for each column, scaled by the entries in this listing (see new comment in A181897 for details).
The number of terms in the n-th row of this irregular triangle is given by p(n+1)-p(n) where p(n) = A000041(n) is the partitions counting function. The row sum of row n is A000166(n+1). The number of 'parts' of a partition P is K(P) := Sum_{j=1..N} m_j; if this sequence is signed by (-1)^(1+N+K), then the signed sum of row n is equal to n. The last entry in row n is n!. The first entry in odd rows n is equal to n!!. The first entry in even rows n is equal to n*(n+1)!!/3.

Examples

			The triangle begins:
    1
    2
    3,   6
   20,  24
   15,  90,  40, 120
  210, 504, 420, 720
		

Crossrefs

Programs

  • Mathematica
    partitionMultiplicities[aPartn_]:=Table[Count[aPartn,m],{m,Total[aPartn]}]
    partitionBase[aPartn_]:=Sum[m*aPartn[[m]],{m,Length[aPartn]}]
    partitionFactorial[aPartn_]:=Product[m^aPartn[[m]],{m,partitionBase[aPartn]}]
    partitionParts[aPartn_]:=Sum[aPartn[[m]],{m,Length[aPartn]}]
    A385081[aPartn_]:=Multinomial@@aPartn*partitionBase[aPartn]!/(partitionFactorial[aPartn]*partitionParts[aPartn]!)
    Grid[Table[Map[A385081,Select[ReverseSort[Map[partitionMultiplicities,Partitions[n]],LexicographicOrder],#[[1]]==0&]],{n,2,12}]]

Formula

Given a partition P of N encoded as its multiplicities N-tuple (m_j: j=1..N) with K(P) = Sum_{j=1..N} m_j 'parts'. Define P! := Product_{j=1..N} j^m_j. Then the number of permutations in the symmetric group S_N of cycles type labelled by P is #P = multinomial(K(P); (m_j){j=1..N}) * N! / (P! * K(P)!), as stated in A181897 by Carlos Mafra.
Showing 1-6 of 6 results.