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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A156626 a(0)=1; a(1)=2; a(2)=6; a(n+1) = 2*(n+1)*a(n) - n^2*a(n-1), n > 1. 1
1, 2, 6, 28, 170, 1252, 10774, 105764, 1164298, 14188468, 189461222, 2749300084, 43058394154, 723619035908, 12984464393398, 247704600763972, 5005042735932554, 106759075226130004, 2396869357456172038, 56491095210068416148, 1394373970361058540202 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..442

MAPLE

A156626 := proc(n)

        if n<=1 then

                n+1 ;

        elif n = 2 then

                6 ;

        else

                2*n*procname(n-1)-(n-1)^2*procname(n-2) ;

        end if;

end proc:

seq(A156626(n), n=0..20) ; # R. J. Mathar, Sep 27 2011

MATHEMATICA

Join[{1}, RecurrenceTable[{a[n] == 2*n*a[n-1] - (n-1)^2*a[n-2], a[1] == 2, a[2] == 6}, a, {n, 1, 50}]] (* G. C. Greubel, Sep 01 2018 *)

PROG

(PARI) m=30; v=concat([2, 6], vector(m-2)); for(n=3, m, v[n] = 2*n*v[n-1]-(n-1)^2*v[n-2]); concat([1], v) \\ G. C. Greubel, Sep 01 2018

(MAGMA) I:=[2, 6]; [1] cat [n le 2 select I[n] else 2*n*Self(n-1) - (n-1)^2*Self(n-2): n in [1..30]]; // G. C. Greubel, Sep 01 2018

CROSSREFS

Cf. A001620, A002720.

Sequence in context: A109570 A262002 A245633 * A207386 A088501 A140092

Adjacent sequences:  A156623 A156624 A156625 * A156627 A156628 A156629

KEYWORD

easy,nonn

AUTHOR

Jonathan Vos Post, Feb 12 2009

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 July 24 05:05 EDT 2021. Contains 346273 sequences. (Running on oeis4.)