A123751 -id:A123751 - cp's OEIS frontend

A007406 Wolstenholme numbers: numerator of Sum_{k=1..n} 1/k^2.

1, 5, 49, 205, 5269, 5369, 266681, 1077749, 9778141, 1968329, 239437889, 240505109, 40799043101, 40931552621, 205234915681, 822968714749, 238357395880861, 238820721143261, 86364397717734821, 17299975731542641, 353562301485889, 354019312583809, 187497409728228241

Offset: 1

N. J. A. Sloane, Mira Bernstein

By Wolstenholme's theorem, p divides a(p-1) for prime p > 3. - T. D. Noe, Sep 05 2002
Also p divides a( (p-1)/2 ) for prime p > 3. - Alexander Adamchuk, Jun 07 2006
The rationals a(n)/A007407(n) converge to Zeta(2) = (Pi^2)/6 = 1.6449340668... (see the decimal expansion A013661).
For the rationals a(n)/A007407(n), n >= 1, see the W. Lang link under A103345 (case k=2).
See the Wolfdieter Lang link under A103345 on Zeta(k, n) with the rationals for k=1..10, g.f.s and polygamma formulas. - Wolfdieter Lang, Dec 03 2013
Denominator of the harmonic mean of the first n squares. - Colin Barker, Nov 13 2014
Conjecture: for n > 3, gcd(n, a(n-1)) = A089026(n). Checked up to n = 10^5. - Amiram Eldar and Thomas Ordowski, Jul 28 2019
True if n is prime, by Wolstenholme's theorem. It remains to show that gcd(n, a(n-1)) = 1 if n > 3 is composite. - Jonathan Sondow, Jul 29 2019
From Peter Bala, Feb 16 2022: (Start)
Sum_{k = 1..n} 1/k^2 = 1 + (1 - 1/2^2)*(n-1)/(n+1) - (1/2^2 - 1/3^2)*(n-1)*(n-2)/((n+1)*(n+2)) + (1/3^2 - 1/4^2)*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) - (1/4^2 - 1/5^2)*(n-1)*(n-2)*(n-3)*(n-4)/((n+1)*(n+2)*(n+3)*(n+4)) + .... Cf. A082687 and A120778.
This identity allows us to extend the definition of Sum_{k = 1..n} 1/k^2 to non-integral values of n. (End)
Numerators of the Eulerian numbers T(-2,k) for k = 0,1..., if T(n,k) is extended to negative n by the recurrence T(n,k) = (k+1)*T(n-1,k) + (n-k)*T(n-1,k-1) (indexed as in A173018). - Michael J. Collins, Oct 10 2024

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Seiichi Manyama, Table of n, a(n) for n = 1..1152 (terms 1..200 from T. D. Noe)
Stephen Crowley, Two New Zeta Constants: Fractal String, Continued Fraction, and Hypergeometric Aspects of the Riemann Zeta Function, arXiv:1207.1126 [math.NT], 2012.
Romeo Mestrovic, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.
Hisanori Mishima, Factorizations of many number sequences
Hisanori Mishima, Factorizations of many number sequences
Hisanori Mishima, Factorizations of many number sequences
D. Y. Savio, E. A. Lamagna and S.-M. Liu, Summation of harmonic numbers, pp. 12-20 of E. Kaltofen and S. M. Watt, editors, Computers and Mathematics, Springer-Verlag, NY, 1989.
M. D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications , J. Int. Seq. 13 (2010), 10.6.7, Section 4.3.2.
Eric Weisstein's World of Mathematics, Wolstenholme's Theorem
Eric Weisstein's World of Mathematics, Wolstenholme Number

Cf. A001008, A007407 (denominators), A000290, A082687, A120778.

Numbers n such that a(n) is prime are listed in A111354. Primes in {a(n)} are listed in A123751. - Alexander Adamchuk, Oct 11 2006

import Data.Ratio ((%), numerator)
a007406 n = a007406_list !! (n-1)
a007406_list = map numerator $ scanl1 (+) $ map (1 %) $ tail a000290_list
-- Reinhard Zumkeller, Jul 06 2012

[Numerator(&+[1/k^2:k in [1..n]]):n in [1..23]]; // Marius A. Burtea, Aug 02 2019

a:= n-> numer(add(1/i^2, i=1..n)): seq(a(n), n=1..24);  # Zerinvary Lajos, Mar 28 2007

a[n_] := If[ n<1, 0, Numerator[HarmonicNumber[n, 2]]]; Table[a[n], {n, 100}]
Numerator[HarmonicNumber[Range[20],2]] (* Harvey P. Dale, Jul 06 2014 *)

{a(n) = if( n<1, 0, numerator( sum( k=1, n, 1 / k^2 ) ) )} /* Michael Somos, Jan 16 2011 */

Sum_{k=1..n} 1/k^2 = sqrt(Sum_{j=1..n} Sum_{i=1..n} 1/(i*j)^2). - Alexander Adamchuk, Oct 26 2004
G.f. for rationals a(n)/A007407(n), n >= 1: polylog(2,x)/(1-x).
a(n) = Numerator of (Pi^2)/6 - Zeta(2,n). - Artur Jasinski, Mar 03 2010

A111354 Numbers n such that the numerator of Sum_{i=1..n} (1/i^2), in reduced form, is prime.

Original entry on oeis.org

2, 7, 13, 19, 121, 188, 252, 368, 605, 745, 1085, 1127, 1406, 1743, 1774, 2042, 2087, 2936, 3196, 3207, 3457, 4045, 7584, 10307, 12603, 12632, 14438, 14526, 14641, 15662, 15950, 16261, 18084, 18937, 19676, 40984, 45531, 46009, 48292, 48590

Offset: 1

Ryan Propper, Nov 05 2005

nonn

Numbers n such that A007406(n) is prime.
Some of the larger entries may only correspond to probable primes.
A007406(n) are the Wolstenholme numbers: numerator of Sum 1/k^2, k = 1..n. Primes in A007406(n) are listed in A123751(n) = A007406(a(n)) = {5,266681,40799043101,86364397717734821,...}.
For prime p>3, Wolstenholme's theorem says that p divides A007406(p-1). Hence n+1 cannot be prime for any n>2 in this sequence. - 12 more terms from T. D. Noe, Nov 11 2005
No other n<50000. All n<=1406 yield provable primes. - T. D. Noe, Mar 08 2006

			A007406(n) begins {1, 5, 49, 205, 5269, 5369, 266681, 1077749, 9778141,...}.
Thus a(1) = 2 because A007406(2) = 5 is prime but A007406(1) = 1 is not prime.
a(2) = 7 because A007406(7) = 266681 is prime but all A007406(k) are composite for 2 < k < 7.

Carlos M. da Fonseca, M. Lawrence Glasser, Victor Kowalenko, Generalized cosecant numbers and trigonometric inverse power sums, Applicable Analysis and Discrete Mathematics, Vol. 12, No. 1 (2018), 70-109.
Eric Weisstein's World of Mathematics, Wolstenholme's Theorem
Eric Weisstein's World of Mathematics, Harmonic Number.
Eric Weisstein's World of Mathematics, Wolstenholme Number

Cf. A007406 (numerator of Sum_{i=1..n} (1/i^2)).

Cf. A123751, A001008, A007407, A067567, A056903.

s = 0; Do[s += 1/n^2; If[PrimeQ[Numerator[s]], Print[n]], {n, 1, 10^4}]
Module[{nn=10400,t},t=Accumulate[1/Range[nn]^2];Select[Thread[{Range[nn],Numerator[t]}],PrimeQ[#[[2]]]&]][[;;,1]] (* The program generates the first 24 terms of the sequence. *) (* Harvey P. Dale, May 18 2025 *)

12 more terms from T. D. Noe, Nov 11 2005
More terms from T. D. Noe, Mar 08 2006
Additional comments from Alexander Adamchuk, Oct 11 2006
Edited by N. J. A. Sloane, Nov 11 2006

cp's OEIS Frontend

A007406 Wolstenholme numbers: numerator of Sum_{k=1..n} 1/k^2.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Formula