OFFSET
1,3
COMMENTS
p^2 divides a(2p-2) for prime p>3. a(2p-2)/p^2 = A061002(n) = A001008(p-1)/p^2 for prime p>2. - Alexander Adamchuk, Jul 07 2006
LINKS
Robert Israel, Table of n, a(n) for n = 1..2000
N. J. A. Sloane, Transforms
FORMULA
a(n) = Numerator[Sum[(-1)^(k+1)*Sum[(-1)^(i+1)*1/i,{i,1,k}],{k,1,n}]]. - Alexander Adamchuk, Jul 07 2006
a(n) = numerator((-1)^(n+1)*1/2*(log(2)+(-1)^(n+1)*(gamma+1/2*(psi(1+n/2)-psi(3/2+n/2))+psi(2+n)))), with gamma the Euler-Mascheroni constant. - - Gerry Martens, Apr 28 2011
MAPLE
S:= series(log(1-x)/(x^2-1), x, 101):
seq(numer(coeff(S, x, j)), j=1..100); # Robert Israel, Jun 02 2015
MATHEMATICA
Numerator[Table[Sum[(-1)^(k+1)*Sum[(-1)^(i+1)*1/i, {i, 1, k}], {k, 1, n}], {n, 1, 50}]] (* Alexander Adamchuk, Jul 07 2006 *)
PROG
(PARI) a(n)=numerator(polcoeff(log(1-x)/(x^2-1)+O(x^(n+1)), n))
CROSSREFS
KEYWORD
nonn,easy,frac
AUTHOR
STATUS
approved