OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..181 (terms 1..69 from R. H. Hardin)
FORMULA
A(n,k) = 0 for n > 2^(k-1). - Andrew Howroyd, Jan 24 2020
EXAMPLE
Array begins:
============================================================================
n\k| 1 2 3 4 5 6 7 8 9
---+------------------------------------------------------------------------
1 | 1 2 5 15 52 203 877 4140 21147
2 | 0 1 8 80 1088 19232 424400 11361786 361058000
3 | 0 0 5 205 11301 904580 101173251 15207243828 2975725761202
4 | 0 0 1 278 67198 24537905 13744869502 11385203921707 ...
5 | 0 0 0 205 250735 425677958 1184910460297 ...
6 | 0 0 0 80 621348 5064948309 ...
7 | 0 0 0 15 1058139 ...
8 | 0 0 0 1 ...
...
Some solutions for 16 X 4:
1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1
1 0 1 1 1 1 0 1 1 1 0 0 1 0 1 1 1 1 0 0
1 0 1 0 1 0 1 1 1 0 1 1 1 0 0 1 1 0 1 1
1 0 0 1 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0
0 1 1 1 0 1 1 0 0 1 1 1 0 1 1 0 0 1 1 1
0 1 0 1 0 1 0 0 0 1 0 0 0 1 0 1 0 1 0 0
0 1 0 0 0 0 1 1 0 0 1 1 0 1 0 0 0 0 1 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
PROG
(PARI)
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+1, q=Vec(exp(intformal(O(x^m) - x^n/(1-x)))/(1+x))); if(n==0, 1, (-1)^m*sum(j=0, m, my(s=0); forpart(p=j, s+=(-1)^#p*D(p, n, k), [1, n]); s*q[#q-j])/2)} \\ Andrew Howroyd, Dec 16 2018
CROSSREFS
Column sums are A319190.
KEYWORD
nonn,tabl
AUTHOR
R. H. Hardin, Mar 31 2011
STATUS
approved