A280166 a(2*n) = 4*n if n>0, a(2*n + 1) = -(2*n + 1), a(0) = 1.

1, -1, 4, -3, 8, -5, 12, -7, 16, -9, 20, -11, 24, -13, 28, -15, 32, -17, 36, -19, 40, -21, 44, -23, 48, -25, 52, -27, 56, -29, 60, -31, 64, -33, 68, -35, 72, -37, 76, -39, 80, -41, 84, -43, 88, -45, 92, -47, 96, -49, 100, -51, 104, -53, 108, -55, 112, -57, 116

Offset

0,3

Links

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

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

Formula

Euler transform of length 6 sequence [-1, 4, 1, -1, 0, -1].

a(n) = (-1)^n * A257088(n), with A257088 multiplicative (see there).

a(n) = n * A168361(n+1) if n>0.

a(2*n) = A008574(n). a(2*n + 1) = - A005408(n).

G.f.: (1 - x + x^2) * (1 + x^2) / (1 - x^2)^2.

Example

G.f. = 1 - x + 4*x^2 - 3*x^3 + 8*x^4 - 5*x^5 + 12*x^6 - 7*x^7 + 16*x^8 + ...

Mathematica

a[ n_] := Which[ n < 1, Boole[n == 0], OddQ[n], -n, True, 2 n];
a[ n_] := SeriesCoefficient[ (1 - x + x^2) (1 + x^2) / (1 - x^2)^2, {x, 0, n}];

Prog

(PARI) {a(n) = if( n<1, n==0, n%2, -n, 2*n)};
(PARI) x='x+O('x^50); Vec((1-x+x^2)*(1+x^2)/(1-x^2)^2) \\ G. C. Greubel, Aug 04 2018
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!((1 - x+x^2)*(1+x^2)/(1-x^2)^2)); // G. C. Greubel, Aug 04 2018

Crossrefs

Cf. A005408, A008574, A168361, A257088.

Sequence in context: A156028 A021232 A263616 * A257088 A022998 A082895

Adjacent sequences: A280163 A280164 A280165 * A280167 A280168 A280169

Keyword

sign,easy

Author

Michael Somos, Dec 27 2016

Extensions

Edited by M. F. Hasler, May 08 2018

Status

approved