A225825 a(2n)=A001896(n). a(2n+1)=(-1)^n*A110501(n+1).
1, 1, -1, -1, 7, 3, -31, -17, 127, 155, -2555, -2073, 1414477, 38227, -57337, -929569, 118518239, 28820619, -5749691557, -1109652905, 91546277357, 51943281731, -1792042792463, -2905151042481, 1982765468311237, 191329672483963, -286994504449393, -14655626154768697, 3187598676787461083, 1291885088448017715, -4625594554880206790555
Offset: 0
Keywords
Crossrefs
Cf. A083420.
Programs
-
Maple
A225825 := proc(n) local nhalf ; nhalf := floor(n/2) ; if type(n,'even') then A001896(nhalf) ; else (-1)^nhalf*A110501(nhalf+1) ; end if; end proc; # R. J. Mathar, Oct 28 2013
-
Mathematica
a[0] = 1; a[n_] := Numerator[BernoulliB[n, 1/2] - (n+1)*EulerE[n, 0]]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Aug 01 2013 *)
Extensions
More terms from Jean-François Alcover, Aug 01 2013
Definition corrected by R. J. Mathar, Oct 28 2013
Comments