login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A081406 a(n) = (n+1)*a(n-3), a(0)=a(1)=a(2)=1 for n>1. 1
1, 1, 1, 4, 5, 6, 28, 40, 54, 280, 440, 648, 3640, 6160, 9720, 58240, 104720, 174960, 1106560, 2094400, 3674160, 24344320, 48171200, 88179840, 608608000, 1252451200, 2380855680, 17041024000, 36321084800, 71425670400, 528271744000, 1162274713600 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,4
LINKS
EXAMPLE
a(3n+2)=A034001[n]; while other subsequences are near(but not equal) to A001669, A000359.
MATHEMATICA
f[n_]:= (n+1)*f[n-3]; f[0]=1; f[1]=1; f[2]=1; Table[f[n], {n, 30}]
RecurrenceTable[{a[0]==a[1]==a[2]==1, a[n]==(n+1)a[n-3]}, a, {n, 30}] (* Harvey P. Dale, Mar 06 2019 *)
PROG
(PARI) a(n) = if(n<3, 1, (n+1)*a(n-3) );
vector(35, n, a(n-1)) \\ G. C. Greubel, Aug 24 2019
(Magma) a:= func< n | n le 2 select 1 else n in [3..5] select n+1 else (n+1)*Self(n-2) >;
[a(n): n in [0..35]]; // G. C. Greubel, Aug 24 2019
(Sage)
def a(n):
if n<3: return 1
elif 3<= n <= 5: return n+1
else: return (n+1)*a(n-3)
[a(n) for n in (0..35)] # G. C. Greubel, Aug 24 2019
(GAP)
a:= function(k)
if k<3 then return 1;
elif k<6 then return k+1;
else return (k+1)*a(k-3);
fi;
end;
List([0..35], n-> a(n) ); # G. C. Greubel, Aug 24 2019
CROSSREFS
Sequence in context: A048075 A048016 A287648 * A019067 A352865 A095212
KEYWORD
nonn
AUTHOR
Labos Elemer, Apr 01 2003
EXTENSIONS
Corrected and extended by Harvey P. Dale, Mar 06 2019
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 2 01:36 EDT 2023. Contains 363078 sequences. (Running on oeis4.)