login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A267219 Expansion of exp( Sum_{n >= 1} A002895(n)*x^n/n ). 3
1, 4, 22, 152, 1241, 11444, 115390, 1243672, 14104480, 166460800, 2028202288, 25363355200, 324098616925, 4217387014948, 55737166570870, 746544123583928, 10116388473816503 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

LINKS

Table of n, a(n) for n=0..16.

FORMULA

n*a(n) = Sum_{k = 0..n-1} A002895(n-k)*a(k).

O.g.f. A(x) = exp( Sum_{n >= 1} A002895(n)*x^n/n ) = 1 + 4*x + 22*x^2 + 152*x^3 + 1241*x^4 + ....

The o.g.f. A(x) satisfies 1 + x* d/dx(log(A(x)) = Sum_{n >= 0} A002895(n)*x^n.

A(x)^(1/4) = 1 + x + 4*x*2 + 25*x^3 + 199*x^4 + 1837*x^5 + ... appears to have integer coefficients.

MAPLE

#A267219

#define the Domb numbers

A002895 := n -> add(binomial(n, k)^2*binomial(2*n-2*k, n-k)*binomial(2*k, k), k = 0..n):

A267219 := proc (n) option remember; if n = 0 then 1 else 1/n*add( A002895(n-k)*A267219(k), k = 0..n-1) end if; end proc:

seq(A267219(n), n = 0..20);

CROSSREFS

Cf. A002895, A267220.

Sequence in context: A307439 A189845 A039304 * A152404 A062817 A196275

Adjacent sequences:  A267216 A267217 A267218 * A267220 A267221 A267222

KEYWORD

nonn,easy

AUTHOR

Peter Bala, Jan 12 2016

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 17 18:51 EDT 2019. Contains 328127 sequences. (Running on oeis4.)