A109534 a(0)=1, a(n)=n+a(n-1) if n mod 2=0, a(n)=3n-a(n-1) if n mod 2 = 1.

1, 2, 4, 5, 9, 6, 12, 9, 17, 10, 20, 13, 25, 14, 28, 17, 33, 18, 36, 21, 41, 22, 44, 25, 49, 26, 52, 29, 57, 30, 60, 33, 65, 34, 68, 37, 73, 38, 76, 41, 81, 42, 84, 45, 89, 46, 92, 49, 97, 50, 100, 53, 105, 54, 108, 57, 113, 58, 116, 61, 121, 62, 124, 65, 129, 66, 132, 69, 137

Offset

0,2

Links

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

Index entries for linear recurrences with constant coefficients, signature (0, 1, 0, 1, 0, -1).

Formula

G.f.: (1+2x+3x^2+3x^3+4x^4-x^5)/((1+x^2)(1-x)^2(1+x)^2). a(n)=a(n-2)+a(n-4)-a(n-6). a(2n)=A047461(n+1). a(2n+1)=A042963(n+1). [From R. J. Mathar, Oct 29 2008]

Maple

a:=proc(n) if n=0 then 1 elif n mod 2=0 then n+a(n-1) else 3*n-a(n-1) fi end: seq(a(n), n=0..68);

Mathematica

a[0] = 1; a[n_] := a[n] = If[Mod[n, 2] == 0, a[n - 1] + n, -a[n - 1] + 3*n] aa = Table[a[n], {n, 0, 100}]

Crossrefs

Sequence in context: A169624 A368797 A011175 * A276163 A011341 A116477

Adjacent sequences: A109531 A109532 A109533 * A109535 A109536 A109537

Keyword

nonn

Author

Roger L. Bagula, Jun 18 2005

Status

approved