login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A120285 Numerator of harmonic number H(p-1) = Sum_{k=1..p-1} 1/k for prime p. 2
1, 3, 25, 49, 7381, 86021, 2436559, 14274301, 19093197, 315404588903, 9304682830147, 54801925434709, 2078178381193813, 12309312989335019, 5943339269060627227, 14063600165435720745359, 254381445831833111660789 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Prime(n)^2 divides a(n) for n>2.

LINKS

Robert Israel, Table of n, a(n) for n = 1..342

R. Mestrovic, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.

Eric Weisstein's World of Mathematics, Wolstenholme's Theorem.

FORMULA

a(n) = numerator(Sum_{k=1..prime(n)-1} 1/k).

a(n) = A001008(prime(n)-1).

a(n) = A061002(n)*prime(n)^2 for n > 2.

MAPLE

f3:=proc(n) local p;

p:=ithprime(n);

numer(add(1/i, i=1..p-1));

end proc;

[seq(f3(n), n=1..20)];

MATHEMATICA

Numerator[Table[Sum[1/k, {k, 1, Prime[n]-1}], {n, 1, 20}]]

PROG

(PARI) a(n) = my(p=prime(n)); numerator(sum(k=1, p-1, 1/k)); \\ Michel Marcus, Dec 25 2018

CROSSREFS

Cf. A001008, A061002, A185399.

Sequence in context: A051280 A145609 A259923 * A041897 A242974 A006222

Adjacent sequences:  A120282 A120283 A120284 * A120286 A120287 A120288

KEYWORD

frac,nonn

AUTHOR

Alexander Adamchuk, Jul 07 2006

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 24 19:14 EDT 2020. Contains 334580 sequences. (Running on oeis4.)