OFFSET
0,4
COMMENTS
There is more complicated way of defining the sequence: consider the sequence of modified Bernoulli numbers EVB(n) = A176327(n)/A176289(n) and its inverse binomial transform IEVB(n) = A176328(n)/A176591(n). Then a(n) is the numerator of the difference EVB(n)-IEVB(n). The denominator of the difference is 1 if n=0, else A040001(n-1).
A particularity of EVB(n) is: its (forward) binomial transform is 1, 1, 7/6, 3/2, 59/30,.. = (-1)^n*IEVB(n).
Note that A026741 is related to the Rydberg-Ritz spectrum of the hydrogen atom.
LINKS
R. J. Mathar, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,2,0,-1).
FORMULA
From R. J. Mathar, Dec 01 2010: (Start)
a(n) = (-1)^n*A026741(n) = n*(1-3*(-1)^n)/4.
G.f.: x*(1-x+x^2) / ( (x-1)^2*(1+x)^2 ).
a(n) = +2*a(n-2) -a(n-4). (End)
a(n) = -a(-n) for all n in Z. - Michael Somos, Jun 11 2013
From Michael Somos, Aug 30 2014: (Start)
Euler transform of length 6 sequence [ -1, 3, 1, 0, 0, -1].
0 = - 1 - a(n) - a(n+1) + a(n+2) + a(n+3) for all n in Z.
0 = 1 + a(n)*(-2 -a(n) + a(n+2)) - 2*a(n+1) - a(n+2) for all n in Z. (End)
From Michael Somos, May 04 2015: (Start)
a(n) is multiplicative with a(2^e) = -(2^(e-1)) if e>0, a(p^e) = p^e otherwise.
G.f.: (f(x) - 3 * f(-x)) / 4 where f(x) := x / (1 - x)^2.
G.f.: x * (1 - x) * (1 - x^6) / ((1 - x^2)^3 * (1 - x^3)). (End)
From Amiram Eldar, Sep 21 2023: (Start)
Dirichlet g.f.: zeta(s-1) * (1 - 3/2^s).
Sum_{k=0..n} a(k) = A008795(n-1), for n > 0.
Sum_{k=0..n} a(k) ~ n^2/8. (End)
EXAMPLE
G.f. = x - x^2 + 3*x^3 - 2*x^4 + 5*x^5 - 3*x^6 + 7*x^7 - 4*x^8 + 9*x^9 - 5*x^10 + ...
MATHEMATICA
a[n_?EvenQ]:=-(n/2); a[n_?OddQ]:=n; Table[a[n], {n, 100}] (* Alonso del Arte, Dec 01 2010 *)
a[ n_] := n / If[ Mod[ n, 2] == 1, 1, -2]; (* Michael Somos, Jun 11 2013 *)
CoefficientList[Series[x (1 - x + x^2)/((x - 1)^2*(1 + x)^2), {x, 0, 70}], x] (* Michael De Vlieger, Dec 10 2016 *)
LinearRecurrence[{0, 2, 0, -1}, {0, 1, -1, 3}, 80] (* Harvey P. Dale, Nov 01 2017 *)
PROG
(Magma) [n*(1-3*(-1)^n)/4: n in [0..60]]; // Vincenzo Librandi, Aug 04 2011
(PARI) {a(n) = n / if( n%2, 1, -2)}; /* Michael Somos, Jun 11 2013 */
CROSSREFS
KEYWORD
sign,easy,mult
AUTHOR
Paul Curtz, Apr 18 2010
STATUS
approved