Identifier
Identifier
Values
[2] generating graphics... => 0
[1,1] generating graphics... => 2
[3] generating graphics... => 0
[2,1] generating graphics... => 5
[1,1,1] generating graphics... => 9
[4] generating graphics... => 0
[3,1] generating graphics... => 16
[2,2] generating graphics... => 36
[2,1,1] generating graphics... => 46
[1,1,1,1] generating graphics... => 64
[5] generating graphics... => 0
[4,1] generating graphics... => 65
[3,2] generating graphics... => 236
[3,1,1] generating graphics... => 268
[2,2,1] generating graphics... => 405
[2,1,1,1] generating graphics... => 497
[1,1,1,1,1] generating graphics... => 625
[6] generating graphics... => 0
[5,1] generating graphics... => 326
[4,2] generating graphics... => 1646
[4,1,1] generating graphics... => 1776
[3,3] generating graphics... => 2658
[3,2,1] generating graphics... => 3682
[3,1,1,1] generating graphics... => 4218
[2,2,2] generating graphics... => 4722
[2,2,1,1] generating graphics... => 5532
[2,1,1,1,1] generating graphics... => 6526
[1,1,1,1,1,1] generating graphics... => 7776
[7] generating graphics... => 0
[6,1] generating graphics... => 1957
[5,2] generating graphics... => 12652
[5,1,1] generating graphics... => 13304
[4,3] generating graphics... => 28620
[4,2,1] generating graphics... => 35529
[4,1,1,1] generating graphics... => 39081
[3,3,1] generating graphics... => 48364
[3,2,2] generating graphics... => 57068
[3,2,1,1] generating graphics... => 64432
[3,1,1,1,1] generating graphics... => 72868
[2,2,2,1] generating graphics... => 77981
[2,2,1,1,1] generating graphics... => 89045
[2,1,1,1,1,1] generating graphics... => 102097
[1,1,1,1,1,1,1] generating graphics... => 117649
[8] generating graphics... => 0
[7,1] generating graphics... => 13700
[6,2] generating graphics... => 107814
[6,1,1] generating graphics... => 111728
[5,3] generating graphics... => 315486
[5,2,1] generating graphics... => 367724
[5,1,1,1] generating graphics... => 394332
[4,4] generating graphics... => 442880
[4,3,1] generating graphics... => 640330
[4,2,2] generating graphics... => 720268
[4,2,1,1] generating graphics... => 791326
[4,1,1,1,1] generating graphics... => 869488
[3,3,2] generating graphics... => 893304
[3,3,1,1] generating graphics... => 990032
[3,2,2,1] generating graphics... => 1139986
[3,2,1,1,1] generating graphics... => 1268850
[3,1,1,1,1,1] generating graphics... => 1414586
[2,2,2,2] generating graphics... => 1323608
[2,2,2,1,1] generating graphics... => 1479570
[2,2,1,1,1,1] generating graphics... => 1657660
[2,1,1,1,1,1,1] generating graphics... => 1861854
[1,1,1,1,1,1,1,1] generating graphics... => 2097152
[9] generating graphics... => 0
[8,1] generating graphics... => 109601
[7,2] generating graphics... => 1015352
[7,1,1] generating graphics... => 1042752
[6,3] generating graphics... => 3654000
[6,2,1] generating graphics... => 4095041
[6,1,1,1] generating graphics... => 4318497
[5,4] generating graphics... => 6659144
[5,3,1] generating graphics... => 8747056
[5,2,2] generating graphics... => 9549024
[5,2,1,1] generating graphics... => 10284472
[5,1,1,1,1] generating graphics... => 11073136
[4,4,1] generating graphics... => 11170353
[4,3,2] generating graphics... => 14293024
[4,3,1,1] generating graphics... => 15573684
[4,2,2,1] generating graphics... => 17351741
[4,2,1,1,1] generating graphics... => 18934393
[4,1,1,1,1,1] generating graphics... => 20673369
[3,3,3] generating graphics... => 16752744
[3,3,2,1] generating graphics... => 20567780
[3,3,1,1,1] generating graphics... => 22547844
[3,2,2,2] generating graphics... => 23169912
[3,2,2,1,1] generating graphics... => 25449884
[3,2,1,1,1,1] generating graphics... => 27987584
[3,1,1,1,1,1,1] generating graphics... => 30816756
[2,2,2,2,1] generating graphics... => 28854249
[2,2,2,1,1,1] generating graphics... => 31813389
[2,2,1,1,1,1,1] generating graphics... => 35128709
[2,1,1,1,1,1,1,1] generating graphics... => 38852417
[1,1,1,1,1,1,1,1,1] generating graphics... => 43046721
[10] generating graphics... => 0
[9,1] generating graphics... => 986410
[8,2] generating graphics... => 10506174
[8,1,1] generating graphics... => 10725376
[7,3] generating graphics... => 44918754
[7,2,1] generating graphics... => 49048662
[7,1,1,1] generating graphics... => 51134166
[6,4] generating graphics... => 101098560
[6,3,1] generating graphics... => 124671082
[6,2,2] generating graphics... => 133419804
[6,2,1,1] generating graphics... => 141609886
[6,1,1,1,1] generating graphics... => 150246880
[5,5] generating graphics... => 131400690
[5,4,1] generating graphics... => 194969340
[5,3,2] generating graphics... => 235686288
[5,3,1,1] generating graphics... => 253180400
[5,2,2,1] generating graphics... => 275721004
[5,2,1,1,1] generating graphics... => 296289948
[5,1,1,1,1,1] generating graphics... => 318436220
[4,4,2] generating graphics... => 283590654
[4,4,1,1] generating graphics... => 305931360
[4,3,3] generating graphics... => 320347536
[4,3,2,1] generating graphics... => 380721282
[4,3,1,1,1] generating graphics... => 411868650
[4,2,2,2] generating graphics... => 419381394
[4,2,2,1,1] generating graphics... => 454084876
[4,2,1,1,1,1] generating graphics... => 491953662
[4,1,1,1,1,1,1] generating graphics... => 533300400
[3,3,3,1] generating graphics... => 435037384
[3,3,2,2] generating graphics... => 481123104
[3,3,2,1,1] generating graphics... => 522258664
[3,3,1,1,1,1] generating graphics... => 567354352
[3,2,2,2,1] generating graphics... => 579502682
[3,2,2,1,1,1] generating graphics... => 630402450
[3,2,1,1,1,1,1] generating graphics... => 686377618
[3,1,1,1,1,1,1,1] generating graphics... => 748011130
[2,2,2,2,2] generating graphics... => 644609030
[2,2,2,2,1,1] generating graphics... => 702317528
[2,2,2,1,1,1,1] generating graphics... => 765944306
[2,2,1,1,1,1,1,1] generating graphics... => 836201724
[2,1,1,1,1,1,1,1,1] generating graphics... => 913906558
[1,1,1,1,1,1,1,1,1,1] generating graphics... => 1000000000
[11] generating graphics... => 0
[10,1] generating graphics... => 9864101
[9,2] generating graphics... => 118687532
[9,1,1] generating graphics... => 120660352
[8,3] generating graphics... => 588005676
[8,2,1] generating graphics... => 630578377
[8,1,1,1] generating graphics... => 652029129
[7,4] generating graphics... => 1579007720
[7,3,1] generating graphics... => 1863969724
[7,2,2] generating graphics... => 1967280808
[7,2,1,1] generating graphics... => 2065378132
[12] generating graphics... => 0
[11,1] generating graphics... => 108505112
[10,2] generating graphics... => 1455009206
[10,1,1] generating graphics... => 1474737408
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 vertex 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 vertex labelled trees, whose reversal of the generating function $f^{(-1)}(x) = x\exp(-x)$, see [1, sec. 3.3]
Fix a set of distinguishable vertices and a coloring of the vertices so that $\lambda_i$ are colored $i$. Then this statistic gives the number of ways of putting a rooted tree on this set of colored vertices so that no leaf is the same color as its parent.
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 = x*exp(-x) # 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:09 by Martin Rubey
Updated
Feb 04, 2018 at 21:51 by Jair Taylor