OFFSET
1,5
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..1275
Andrew Howroyd, Additional PARI Programs
FORMULA
T(n,k) = T(k,k) for n > k.
EXAMPLE
Array begins:
===========================================================
n\k| 1 2 3 4 5 6 7 8 9 10 11 12
---+-------------------------------------------------------
1 | 1 1 1 1 1 1 1 1 1 1 1 1 ...
2 | 1 4 5 11 14 24 30 45 55 76 91 119 ...
3 | 1 4 10 24 51 114 219 424 768 1352 2278 3759 ...
4 | 1 4 10 33 78 224 549 1403 3292 7677 16934 36581 ...
5 | 1 4 10 33 91 277 792 2341 6654 18802 51508 138147 ...
6 | 1 4 10 33 91 298 881 2825 8791 27947 87410 272991 ...
7 | 1 4 10 33 91 298 910 2974 9655 32287 108274 367489 ...
8 | 1 4 10 33 91 298 910 3017 9886 33767 116325 410298 ...
9 | 1 4 10 33 91 298 910 3017 9945 34124 118729 424498 ...
...
MATHEMATICA
permcount[v_List] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
c[p_List, q_List, k_] := SeriesCoefficient[1/Product[(1 - x^LCM[p[[i]], q[[j]]])^GCD[p[[i]], q[[j]]], {j, 1, Length[q]}, {i, 1, Length[p]}], {x, 0, k}];
M[m_, n_, k_] := Module[{s=0}, Do[Do[s += permcount[p]*permcount[q]*c[p, q, k], {q, IntegerPartitions[n]}], {p, IntegerPartitions[m]}]; s/(m!*n!)];
Table[M[n-k+1, n-k+1, k], {n, 1, 12}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Sep 12 2018, after Andrew Howroyd *)
PROG
(PARI) \\ see also link.
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}
K(q, t, k)={1/prod(j=1, #q, (1-y^lcm(t, q[j]) + O(y*y^k))^gcd(t, q[j]))}
M(m, n, k)={my(s=0); forpart(q=m, s+=permcount(q)*polcoef(polcoef(exp(sum(t=1, n, K(q, t, k)/t*x^t) + O(x*x^n)), n), k)); s/m!}
for(n=1, 10, for(k=1, 12, print1(M(n, n, k), ", ")); print); \\ updated Andrew Howroyd, Mar 29 2020
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Andrew Howroyd, Sep 03 2018
STATUS
approved