A120286 Numerator of 1/n^2 + 2/(n-1)^2 + 3/(n-2)^2 +...+ (n-1)/2^2 + n.

1, 9, 65, 725, 3899, 28763, 419017, 864669, 7981633, 3586319, 200763407, 2649665993, 34899471137, 176508049513, 356606957297, 12234391348253, 209672027529221, 4012917216669239, 15350275129353301, 15443118015171841

Offset

1,2

Comments

p^2 divides a(p-1) for prime p>2. p divides a(p-2) for prime p>3.

Links

Table of n, a(n) for n=1..20.

Formula

a(n) = numerator[Sum[Sum[1/i^2,{i,1,k}],{k,1,n}]].

Mathematica

Numerator[Table[Sum[Sum[1/i^2, {i, 1, k}], {k, 1, n}], {n, 1, 30}]]
Table[-EulerGamma + HarmonicNumber[1 + n, 2] + n*HarmonicNumber[1 + n, 2] - PolyGamma[0, 2 + n], {n, 1, 20}] // Numerator (* Vaclav Kotesovec, May 02 2024 *)

Prog

(Python)
from fractions import Fraction
def A120286(n): return sum(Fraction(n-i+1, i**2) for i in range(1, n+1)).numerator # Chai Wah Wu, May 01 2024

Crossrefs

Cf. A027612, A001008, A007406, A007407, A370774 (denom.).

Sequence in context: A339688 A100311 A259242 * A152581 A122733 A118465

Adjacent sequences: A120283 A120284 A120285 * A120287 A120288 A120289

Keyword

frac,nonn

Author

Alexander Adamchuk, Jul 07 2006

Status

approved