reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
editing
approved
Main diagonal gives A188388.
approved
editing
Andrew Howroyd, <a href="/A188392/b188392_1.txt">Table of n, a(n) for n = 1..181</a> (terms 1..69 from R. H. Hardin)
proposed
approved
editing
proposed
R. H. Hardin, Andrew Howroyd, <a href="/A188392/b188392_1.txt">Table of n, a(n) for n = 1..181</a> (terms 1..69</a> from R. H. Hardin)
permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); WeighT(v)[n]^k/prod(i=1, #v, i^v[i]*v[i]!)}
T(n, k)={my(m=n*k, q=Vec(exp(O(x*x^m) + intformal((x^n-1)/(1-x)))/(1-x))); if(n==0, 1, sum(j=0, m, my(s=0); forpart(p=n*k, j, s+=permcount(p)*polcoef(prodD(i=1, #p, 1 + x^p[i] + O(x*x^, n), k), [1, n)^k]); s/(n*kq[#q-j]))!} \\ Andrew Howroyd, Dec 10 12 2018