Identifier
Identifier
Values
[2] generating graphics... => 0
[1,1] generating graphics... => 1
[3] generating graphics... => 0
[2,1] generating graphics... => 3
[1,1,1] generating graphics... => 4
[4] generating graphics... => 0
[3,1] generating graphics... => 13
[2,2] generating graphics... => 19
[2,1,1] generating graphics... => 22
[1,1,1,1] generating graphics... => 26
[5] generating graphics... => 0
[4,1] generating graphics... => 75
[3,2] generating graphics... => 141
[3,1,1] generating graphics... => 154
[2,2,1] generating graphics... => 188
[2,1,1,1] generating graphics... => 210
[1,1,1,1,1] generating graphics... => 236
[6] generating graphics... => 0
[5,1] generating graphics... => 541
[4,2] generating graphics... => 1231
[4,1,1] generating graphics... => 1306
[3,3] generating graphics... => 1543
[3,2,1] generating graphics... => 1864
[3,1,1,1] generating graphics... => 2018
[2,2,2] generating graphics... => 2118
[2,2,1,1] generating graphics... => 2306
[2,1,1,1,1] generating graphics... => 2516
[1,1,1,1,1,1] generating graphics... => 2752
[7] generating graphics... => 0
[6,1] generating graphics... => 4683
[5,2] generating graphics... => 12453
[5,1,1] generating graphics... => 12994
[4,3] generating graphics... => 18441
[4,2,1] generating graphics... => 21128
[4,1,1,1] generating graphics... => 22434
[3,3,1] generating graphics... => 24508
[3,2,2] generating graphics... => 26834
[3,2,1,1] generating graphics... => 28698
[3,1,1,1,1] generating graphics... => 30716
[2,2,2,1] generating graphics... => 31634
[2,2,1,1,1] generating graphics... => 33940
[2,1,1,1,1,1] generating graphics... => 36456
[1,1,1,1,1,1,1] generating graphics... => 39208
[8] generating graphics... => 0
[7,1] generating graphics... => 47293
[6,2] generating graphics... => 143599
[6,1,1] generating graphics... => 148282
[5,3] generating graphics... => 243343
[5,2,1] generating graphics... => 269872
[5,1,1,1] generating graphics... => 282866
[4,4] generating graphics... => 285811
[4,3,1] generating graphics... => 353188
[4,2,2] generating graphics... => 378234
[4,2,1,1] generating graphics... => 399362
[4,1,1,1,1] generating graphics... => 421796
[3,3,2] generating graphics... => 420982
[3,3,1,1] generating graphics... => 445490
[3,2,2,1] generating graphics... => 480242
[3,2,1,1,1] generating graphics... => 508940
[3,1,1,1,1,1] generating graphics... => 539656
[2,2,2,2] generating graphics... => 518794
[2,2,2,1,1] generating graphics... => 550428
[2,2,1,1,1,1] generating graphics... => 584368
[2,1,1,1,1,1,1] generating graphics... => 620824
[1,1,1,1,1,1,1,1] generating graphics... => 660032
[9] generating graphics... => 0
[8,1] generating graphics... => 545835
[7,2] generating graphics... => 1860981
[7,1,1] generating graphics... => 1908274
[6,3] generating graphics... => 3536121
[6,2,1] generating graphics... => 3837368
[6,1,1,1] generating graphics... => 3985650
[5,4] generating graphics... => 4717425
[5,3,1] generating graphics... => 5566564
[5,2,2] generating graphics... => 5875418
[5,2,1,1] generating graphics... => 6145290
[5,1,1,1,1] generating graphics... => 6428156
[4,4,1] generating graphics... => 6262244
[4,3,2] generating graphics... => 7175750
[4,3,1,1] generating graphics... => 7528938
[4,2,2,1] generating graphics... => 7995602
[4,2,1,1,1] generating graphics... => 8394964
[4,1,1,1,1,1] generating graphics... => 8816760
[3,3,3] generating graphics... => 7788946
[3,3,2,1] generating graphics... => 8704434
[3,3,1,1,1] generating graphics... => 9149924
[3,2,2,2] generating graphics... => 9270770
[3,2,2,1,1] generating graphics... => 9751012
[3,2,1,1,1,1] generating graphics... => 10259952
[3,1,1,1,1,1,1] generating graphics... => 10799608
[2,2,2,2,1] generating graphics... => 10403260
[2,2,2,1,1,1] generating graphics... => 10953688
[2,2,1,1,1,1,1] generating graphics... => 11538056
[2,1,1,1,1,1,1,1] generating graphics... => 12158880
[1,1,1,1,1,1,1,1,1] generating graphics... => 12818912
[10] generating graphics... => 0
[9,1] generating graphics... => 7087261
[8,2] generating graphics... => 26787871
[8,1,1] generating graphics... => 27333706
[7,3] generating graphics... => 56261743
[7,2,1] generating graphics... => 60125584
[7,1,1,1] generating graphics... => 62033858
[6,4] generating graphics... => 83516371
[6,3,1] generating graphics... => 95478004
[6,2,2] generating graphics... => 99760218
[6,2,1,1] generating graphics... => 103597586
[6,1,1,1,1] generating graphics... => 107583236
[5,5] generating graphics... => 94639351
[5,4,1] generating graphics... => 118577356
[5,3,2] generating graphics... => 132308254
[5,3,1,1] generating graphics... => 137874818
[5,2,2,1] generating graphics... => 144868706
[5,2,1,1,1] generating graphics... => 151013996
[5,1,1,1,1,1] generating graphics... => 157442152
[4,4,2] generating graphics... => 144918078
[4,4,1,1] generating graphics... => 151180322
[4,3,3] generating graphics... => 154713034
[4,3,2,1] generating graphics... => 170124098
[4,3,1,1,1] generating graphics... => 177653036
[4,2,2,2] generating graphics... => 179317786
[4,2,2,1,1] generating graphics... => 187313388
[4,2,1,1,1,1] generating graphics... => 195708352
[4,1,1,1,1,1,1] generating graphics... => 204525112
[3,3,3,1] generating graphics... => 182187314
[3,3,2,2] generating graphics... => 192228218
[3,3,2,1,1] generating graphics... => 200932652
[3,3,1,1,1,1] generating graphics... => 210082576
[3,2,2,2,1] generating graphics... => 212210484
[3,2,2,1,1,1] generating graphics... => 221961496
[3,2,1,1,1,1,1] generating graphics... => 232221448
[3,1,1,1,1,1,1,1] generating graphics... => 243021056
[2,2,2,2,2] generating graphics... => 224265028
[2,2,2,2,1,1] generating graphics... => 234668288
[2,2,2,1,1,1,1] generating graphics... => 245621976
[2,2,1,1,1,1,1,1] generating graphics... => 257160032
[2,1,1,1,1,1,1,1,1] generating graphics... => 269318912
[1,1,1,1,1,1,1,1,1,1] generating graphics... => 282137824
[11] generating graphics... => 0
[10,1] generating graphics... => 102247563
[9,2] generating graphics... => 424132773
[9,1,1] generating graphics... => 431220034
[8,3] generating graphics... => 974072121
[8,2,1] generating graphics... => 1029285368
[8,1,1,1] generating graphics... => 1056619074
[7,4] generating graphics... => 1586426721
[7,3,1] generating graphics... => 1772575588
[7,2,2] generating graphics... => 1838425994
[7,2,1,1] generating graphics... => 1898551578
[7,1,1,1,1] generating graphics... => 1960585436
[6,5] generating graphics... => 1993679073
[12] generating graphics... => 0
[11,1] generating graphics... => 1622632573
click to show generating function       
Description
The coefficient times the product of the factorials of the parts of the monomial symmetric function indexed by the partition in the formal group law for leaf labelled trees.
For a generating function $f$ the associated formal group law is the symmetric function $f(f^{(-1)}(x_1) + f^{(-1)}(x_2), \dots)$, see [1].
This statistic records the coefficient of the monomial symmetric function $m_\lambda$ times the product of the factorials of the parts of $\lambda$ in the formal group law for leaf labelled binary trees, whose generating function is the reversal of $f^{(-1)}(x) = 1+2x-\exp(x)$, see [1, sec. 3.2]
Fix a set of distinguishable vertices and a coloring of the vertices so that $\lambda_i$ are colored $i$. This statistic gives the number of rooted trees with leaves labeled with this set of vertices and internal vertices unlabeled so that no pair of 'twin' leaves have the same color.
References
[1] Taylor, J. Formal group laws and hypergraph colorings MathSciNet:3542357
Code
@cached_function
def data(n):
    R. = PowerSeriesRing(SR, default_prec=n+1)
    f_rev = 1+2*x-exp(x) # leaf labelled trees
    f = f_rev.reverse()
    f_coefficients = f.list()
    t = var('t')
    polynomials = (t*f_rev).exp().list()
    polynomials = [p.expand() for p in polynomials]
    return (f_coefficients, polynomials)

def statistic(P):
    f_coefficients, polynomials = data(P.size())
    p = SR(1)
    for i in P:
        p *= polynomials[i]
    p = p.expand()
    return (prod(factorial(e) for e in P)
	    *sum(p.coefficient(t,n) * f_coefficients[n] * factorial(n)
                 for n in range(p.degree(t)+1)).expand())

Created
Feb 02, 2018 at 20:16 by Martin Rubey
Updated
Feb 06, 2018 at 07:23 by Jair Taylor