login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A033452 "STIRLING" transform of squares A000290. 5
0, 1, 5, 22, 99, 471, 2386, 12867, 73681, 446620, 2856457, 19217243, 135610448, 1001159901, 7714225057, 61904585510, 516347066551, 4468588592739, 40058673825258, 371421499686007, 3556976106133821 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

If an integer N is squarefree and has n+2 distinct prime factors then a(n) is the number of product signs needed to write the factorizations of N, so a(n)=A076277(N). - Floor van Lamoen, Oct 17 2002

Convolved with powers of 2 = A058681: (1, 7, 36, 171, 813,...). Cf. triangle A180338. - Gary W. Adamson, Aug 28 2010

LINKS

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

FORMULA

Representation as an infinite series, in Maple notation : a(n)=sum(k^n*k*(k-2)/k!, k=1..infinity)/exp(1), n=1, 2... . This is a Dobinski-type summation formula. - Karol A. Penson, Mar 21 2002

a(n) = A005493(n) - A000110(n). - Floor van Lamoen and _Christian Bower_, Oct 16 2002. (n^2 has e.g.f.: e^x * (x^2+x), a(n) thus has e.g.f: e^(e^x-1) * ( (e^x-1)^2 + (e^x-1) ) which simplifies to e^(e^x-1) * (e^2x - e^x). A005493 has e.g.f.: e^(e^x+2x-1), A000110 has e.g.f.: e^(e^x-1), A000110(n+1) has as e.g.f.: derivative of A000110 which is e^(e^x+x-1).)

a(n) = Bell(n+2) - 2*Bell(n+1). - Vladeta Jovovic, Jul 28 2003

G.f.: sum{k>=0, k^2*x^k/prod[l=1..k, 1-lx]}. - Ralf Stephan, Apr 18 2004

E.g.f.: exp( exp(x) - 1 + x) * (exp(x) - 1). - Michael Somos, Mar 28 2012

a(n) = A123158(n,3). - Philippe Deléham, Oct 06 2006

G.f.: G(0)/x -1/x, where G(k) = 1 - x^2*(k+1)/( x^2*(k+1) - (2*x+x*k-1)*(3*x+x*k-1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Feb 25 2014

EXAMPLE

G.f. = x + 5*x^2 + 22*x^3 + 99*x^4 + 471*x^5 + 2386*x^6 + 12867*x^7 + 73681*x^8 + ...

MAPLE

a:=n->(sum((j-1)*stirling2(n, j), j=2..n)): seq(a(n), n=2..21); # Zerinvary Lajos, Apr 18 2007

MATHEMATICA

max = 20; Clear[g]; g[max + 2] = 1; g[k_] := g[k] = 2 - 1/(1 - k*x)/(1 - x/(x - 1/g[k + 1])); gf = 1/x + 1/x^2 - g[0]/x^2; CoefficientList[ Series[gf, {x, 0, max}], x] (* Jean-François Alcover, Jan 24 2013, after Sergei N. Gladkovskii *)

PROG

(PARI) {a(n) = if( n<0, 0, n! * polcoeff( (exp(x + x * O(x^n)) - 1) * exp( exp(x + x * O(x^n)) - 1 + x), n))}; /* Michael Somos, Mar 28 2012 */

CROSSREFS

Partial sums of A005494.

Cf. A180338.

Sequence in context: A129164 A123347 A087439 * A295519 A179602 A262440

Adjacent sequences:  A033449 A033450 A033451 * A033453 A033454 A033455

KEYWORD

nonn

AUTHOR

N. J. A. Sloane

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 January 22 10:25 EST 2020. Contains 331144 sequences. (Running on oeis4.)