OFFSET
1,3
COMMENTS
Up to n = 18, a(n) is the same as A058313(n) = numerator of the n-th alternating harmonic number, Sum_{k=1..n} (-1)^(k+1)/k. a(n) differs from A058313(n) only for n = 18, 28, 87, 99.
Up to n = 100 the ratio a(n)/A058313(n) = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 11, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1}.
A Wolstenholme-like theorem: for prime p > 3, if p = 6*k - 1, then p divides a(4*k-1), otherwise if p = 6*k + 1, then p divides a(4*k).
Sum_{j=1..n} Sum_{i=1..n} (-1)^(i+j)*i/j = -1/4 * (2*(-1)^n*n + (-1)^n - 1) * Sum_{k=1..n} (-1)^(k+1)/k.
FORMULA
a(n) = abs(numerator(Sum_{j=1..n} Sum_{i=1..n} (-1)^(i+j)*i/j)).
EXAMPLE
The absolute values of the first few fractions are 1, 1/2, 5/3, 7/6, 47/20, 37/20, 319/105, 533/210, 1879/504, ... = A120301/A334724. - Petros Hadjicostas, May 09 2020
MATHEMATICA
Abs[Numerator[Table[Sum[Sum[(-1)^(i+j)*i/j, {i, 1, n}], {j, 1, n}], {n, 1, 50}]]]
PROG
(PARI) a(n) = abs(numerator(sum(j=1, n, sum(i=1, n, (-1)^(i+j)*i/j)))); \\ Michel Marcus, May 09 2020
CROSSREFS
KEYWORD
frac,nonn
AUTHOR
Alexander Adamchuk, Jul 12 2006
STATUS
approved