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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A003048 a(n+1)=n*a(n)-(-1)^n. 1
1, 2, 3, 10, 39, 196, 1175, 8226, 65807, 592264, 5922639, 65149030, 781788359, 10163248668, 142285481351, 2134282220266, 34148515524255, 580524763912336, 10449445750422047, 198539469258018894 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

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

FORMULA

E.g.f.: (2-exp(-x))/(1-x) (if offset 0).

a(n+2)=n(a(n)+a(n+1)) for n>0. - Amarnath Murthy (amarnath_murthy(AT)yahoo.com), Oct 05 2002

a(n) = 2*(n-1)! - floor(((n-1)!+1)/e), n>1 [From Gary Detlefs, Apr 11 2010]

MAPLE

a:= proc(p) option remember; p*a(p-1)-(-1)^p end proc: a(0):= 1: seq(a(p), p=0..19); - Robert Israel, Jan 05 2008

MATHEMATICA

a[0] = 1; a[p_] := p*a[p - 1] -(-1)^p; a /@ Range[0, 19] - Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Mar 29 2007

FoldList[#1*#2 - (-1)^#2 &, 1, Range[19]] (* Robert G. Wilson v, Jul 07 2012 *)

PROG

(PARI) a(n)=if(n<2, n>0, (n-1)*a(n-1)+(-1)^n)

(PARI) a(n)=if(n<1, 0, (n-1)!*polcoeff((2-exp(-x+O(x^n)))/(1-x), n-1))

CROSSREFS

a(n)=(n-1)!+A002467(n), n>0.

Sequence in context: A141102 A144720 A164933 * A008980 A064183 A050381

Adjacent sequences:  A003045 A003046 A003047 * A003049 A003050 A003051

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 | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified May 22 09:31 EDT 2013. Contains 225519 sequences.