OFFSET
1,4
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..50
Eric Weisstein's World of Mathematics, Edge Cover
Eric Weisstein's World of Mathematics, Path Complement Graph
FORMULA
a(n) = Sum_{i=0..n} (Sum_{k=0, n-i} binomial(n-i, k)*A006129(k)*(2^i)^k*(2^i-1)^(n-i-k)) * (Sum_{j=floor(i/2}, i} (-1)^j*binomial(n-j, i-j)*binomial(j-1, 2*j-i)*2^binomial(i, 2)/2^j). - Andrew Howroyd, Apr 23 2018
MATHEMATICA
Table[Sum[Sum[Binomial[n - i, k] Sum[(-1)^(k - j) Binomial[k, j] 2^Binomial[j, 2], {j, 0, k}] (2^i)^k If[i == 0 && k == n, 1, (2^i - 1)^(n - i - k)], {k, 0, n - i}] Sum[(-1)^j Binomial[n - j, i - j] Binomial[j - 1, 2 j - i] 2^(Binomial[i, 2] - j), {j, Ceiling[i/2], i}], {i, 0, n}], {n, 10}] (* Eric W. Weisstein, Apr 24 2018 *)
PROG
(PARI)
a(n)={ my(p=serlaplace(sum(k=0, n, 2^binomial(k, 2)*x^k/k!)/exp(x+O(x*x^n))));
sum(i=0, n, sum(k=0, n-i, binomial(n-i, k)*polcoeff(p, k)*(2^i)^k*(2^i-1)^(n-i-k)) * sum(j=i\2, i, (-1)^j * binomial(n-j, i-j) * binomial(j-1, 2*j-i) * 2^binomial(i, 2)/2^j))} \\ Andrew Howroyd, Apr 23 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Eric W. Weisstein, Apr 12 2018
EXTENSIONS
Terms a(10) and beyond from Andrew Howroyd, Apr 23 2018
STATUS
approved