OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..50
Eric Weisstein's World of Mathematics, Rook Graph
Eric Weisstein's World of Mathematics, Total Dominating Set
FORMULA
a(n) = Sum_{k=0..n} (-1)^k*binomial(n,k)^2*k!*A287065(n-k). - Andrew Howroyd, Apr 20 2018
a(n) ~ 2^(n^2). - Vaclav Kotesovec, Apr 20 2018
MATHEMATICA
b[0] = 1; b[n_] := (2^n - 1)^n + Sum[Binomial[n, i] Sum[(-1)^j (-1 + 2^(n - j))^i Binomial[n, j], {j, 0, n}], {i, n - 1}]; Table[Sum[(-1)^k Binomial[n, k]^2 k! b[n - k], {k, 0, n}], {n, 10}]
PROG
(PARI) \\ here c(n) is A287065.
b(m, n)=sum(j=0, m, (-1)^j*binomial(m, j)*(2^(m - j) - 1)^n);
c(n)=(2^n-1)^n + sum(i=1, n-1, b(n, i)*binomial(n, i));
a(n) = {sum(k=0, n, (-1)^k*binomial(n, k)^2*k!*c(n-k))} \\ Andrew Howroyd, Apr 20 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Eric W. Weisstein, Apr 19 2018
EXTENSIONS
Terms a(6) and beyond from Andrew Howroyd, Apr 20 2018
STATUS
approved