cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A082687 Numerator of Sum_{k=1..n} 1/(n+k).

Original entry on oeis.org

1, 7, 37, 533, 1627, 18107, 237371, 95549, 1632341, 155685007, 156188887, 3602044091, 18051406831, 7751493599, 225175759291, 13981692518567, 14000078506967, 98115155543129, 3634060848592973, 3637485804655193
Offset: 1

Views

Author

Benoit Cloitre, Apr 12 2003

Keywords

Comments

Numerator of Sum_{k=0..n-1} 1/((k+1)(2k+1)) (denominator is A111876). - Paul Barry, Aug 19 2005
Numerator of the sum of all matrix elements of n X n Hilbert matrix M(i,j) = 1/(i+j-1) (i,j = 1..n). - Alexander Adamchuk, Apr 11 2006
Numerator of the 2n-th alternating harmonic number H'(2n) = Sum ((-1)^(k+1)/k, k=1..2n). H'(2n) = H(2n) - H(n), where H(n) = Sum_{k=1..n} 1/k is the n-th Harmonic Number. - Alexander Adamchuk, Apr 11 2006
a(n) almost always equals A117731(n) = numerator(n*Sum_{k=1..n} 1/(n+k)) = numerator(Sum_{j=1..n} Sum_{i=1..n} 1/(i+j-1)) but differs for n = 14, 53, 98, 105, 111, 114, 119, 164. - Alexander Adamchuk, Jul 16 2006
Sum_{k=1..n} 1/(n+k) = n!^2 *Sum_{j=1..n} (-1)^(j+1) /((n+j)!(n-j)!j). - Leroy Quet, May 20 2007
Seems to be the denominator of the harmonic mean of the first n hexagonal numbers. - Colin Barker, Nov 19 2014
Numerator of 2*n*binomial(2*n,n)*Sum_{k = 0..n-1} (-1)^k* binomial(n-1,k)/(n+k+1)^2. Cf. A049281. - Peter Bala, Feb 21 2017
From Peter Bala, Feb 16 2022: (Start)
2*Sum_{k = 1..n} 1/(n+k) = 1 + 1/(1*2)*(n-1)/(n+1) - 1/(2*3)*(n-1)*(n-2)/((n+1)*(n+2)) + 1/(3*4)*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) - 1/(4*5)*(n-1)*(n-2)*(n-3)*(n-4)/((n+1)*(n+2)*(n+3)*(n+4)) + - .... Cf. A101028.
2*Sum_{k = 1..n} 1/(n+k) = n - (1 + 1/2^2)*n*(n-1)/(n+1) + (1/2^2 + 1/3^2)*n*(n-1)*(n-2)/((n+1)*(n+2)) - (1/3^2 + 1/4^2)*n*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + (1/4^2 + 1/5^2)*n*(n-1)*(n-2)*(n-3)*(n-4)/((n+1)*(n+2)*(n+3)*(n+4)) - + .... Cf. A007406 and A120778.
These identities allow us to extend the definition of Sum_{k = 1..n} 1/(n+k) to non-integral values of n. (End)

Examples

			H'(2n) = H(2n) - H(n) = {1/2, 7/12, 37/60, 533/840, 1627/2520, 18107/27720, 237371/360360, 95549/144144, 1632341/2450448, 155685007/232792560, ...}, where H(n) = A001008/A002805.
n=2: HilbertMatrix(n,n)
   1  1/2
  1/2 1/3
so a(2) = Numerator(1 + 1/2 + 1/2 + 1/3) = Numerator(7/3) = 7.
The n X n Hilbert matrix begins:
   1   1/2  1/3  1/4  1/5  1/6  1/7  1/8  ...
  1/2  1/3  1/4  1/5  1/6  1/7  1/8  1/9  ...
  1/3  1/4  1/5  1/6  1/7  1/8  1/9  1/10 ...
  1/4  1/5  1/6  1/7  1/8  1/9  1/10 1/11 ...
  1/5  1/6  1/7  1/8  1/9  1/10 1/11 1/12 ...
  1/6  1/7  1/8  1/9  1/10 1/11 1/12 1/13 ...
		

Crossrefs

Bisection of A058313, A082688 (denominators).

Programs

  • Magma
    [Numerator((HarmonicNumber(2*n) -HarmonicNumber(n))): n in [1..40]]; // G. C. Greubel, Jul 24 2023
    
  • Maple
    a := n -> numer(harmonic(2*n) - harmonic(n)):
    seq(a(n), n=1..20); # Peter Luschny, Nov 02 2017
  • Mathematica
    Numerator[Sum[1/k,{k,1,2*n}] - Sum[1/k,{k,1,n}]] (* Alexander Adamchuk, Apr 11 2006 *)
    Table[Numerator[Sum[1/(i + j - 1), {i, n}, {j, n}]], {n, 20}] (* Alexander Adamchuk, Apr 11 2006 *)
    Table[HarmonicNumber[2 n] - HarmonicNumber[n], {n, 20}] // Numerator (* Eric W. Weisstein, Dec 14 2017 *)
  • PARI
    a(n) = numerator(sum(k=1, n, 1/(n+k))); \\ Michel Marcus, Dec 14 2017
    
  • SageMath
    [numerator(harmonic_number(2*n,1) - harmonic_number(n,1)) for n in range(1,41)] # G. C. Greubel, Jul 24 2023

Formula

Limit_{n -> oo} Sum_{k=1..n} 1/(n+k) = log(2).
Numerator of Psi(2*n+1) - Psi(n+1). - Vladeta Jovovic, Aug 24 2003
a(n) = numerator((Sum_{k=1..2*n} 1/k) - Sum_{k=1..n} 1/k). - Alexander Adamchuk, Apr 11 2006
a(n) = numerator(Sum_{j=1..n} (Sum_{i=1..n} 1/(i+j-1))). - Alexander Adamchuk, Apr 11 2006
The o.g.f for Sum_{k=1..n} 1/(n+k) is f(x) = (sqrt(x)*log((1+sqrt(x))/(1-sqrt(x))) + log(1-x))/(2*x*(1-x)).

A117664 Denominator of the sum of all elements in the n X n Hilbert matrix M(i,j) = 1/(i+j-1), where i,j = 1..n.

Original entry on oeis.org

1, 3, 10, 105, 252, 2310, 25740, 9009, 136136, 11639628, 10581480, 223092870, 1029659400, 2868336900, 11090902680, 644658718275, 606737617200, 4011209802600, 140603459396400, 133573286426580, 5215718803323600
Offset: 1

Views

Author

Alexander Adamchuk, Apr 11 2006

Keywords

Comments

Sum_{j=1..n} Sum_{i=1..n} 1/(i+j-1) = A117731(n) / A117664(n) = 2n * H'(2n) = 2n * A058313(2n) / A058312(2n), where H'(2n) is 2n-th alternating sign Harmonic Number. H'(2n) = H(2n) - H(n), where H(n) is n-th Harmonic Number. - Alexander Adamchuk, Apr 23 2006

Examples

			For n=2, the 2 X 2 Hilbert matrix is [1, 1/2; 1/2, 1/3], so a(2) = denominator(1 + 1/2 + 1/2 + 1/3) = denominator(7/3) = 3.
The n X n Hilbert matrix begins:
    1 1/2 1/3 1/4  1/5  1/6  1/7  1/8 ...
  1/2 1/3 1/4 1/5  1/6  1/7  1/8  1/9 ...
  1/3 1/4 1/5 1/6  1/7  1/8  1/9 1/10 ...
  1/4 1/5 1/6 1/7  1/8  1/9 1/10 1/11 ...
  1/5 1/6 1/7 1/8  1/9 1/10 1/11 1/12 ...
  1/6 1/7 1/8 1/9 1/10 1/11 1/12 1/13 ...
  ...
		

Crossrefs

Programs

  • Mathematica
    Table[Denominator[Sum[1/(i + j - 1), {i, n}, {j, n}]], {n, 30}]

Formula

a(n) = A111876(n-1)/n.
a(n) = denominator( Sum_{j=1..n} Sum_{i=1..n} 1/(i+j-1) ). Numerator is A117731(n). - Alexander Adamchuk, Apr 23 2006
a(n) = denominator( Sum_{k=1..n} (2*k)/(n+k) ). - Peter Bala, Oct 10 2021

A296519 Denominator of n*Sum_{k=1..n} 1/(n+k).

Original entry on oeis.org

2, 6, 20, 210, 504, 4620, 51480, 18018, 272272, 23279256, 21162960, 446185740, 2059318800, 5736673800, 22181805360, 1289317436550, 1213475234400, 8022419605200, 281206918792800, 267146572853160, 10431437606647200, 428163098127382800, 409547311252279200
Offset: 1

Views

Author

Eric W. Weisstein, Dec 14 2017

Keywords

Comments

a(n) is divisible by all primes p such that the numerator of Sum_{n < k*p <= n} 1/k is not divisible by p, in particular by all primes from n+1 to 2*n-1. - Robert Israel, May 21 2020

Examples

			The first few fractions are 1/2, 7/6, 37/20, 533/210, 1627/504, 18107/4620, 237371/51480, ... = A117731/a(n).
		

Crossrefs

Cf. A111876, A117731 (numerators), A117664.

Programs

  • Magma
    [Denominator(n*(HarmonicNumber(2*n) -HarmonicNumber(n))): n in [1..40]]; // G. C. Greubel, Jul 24 2023
    
  • Maple
    N:= 30: # for a(1)..a(N)
    H:= ListTools:-PartialSums([seq(1/i,i=1..2*N)]):
    map(n -> denom(n*(H[2*n]-H[n])), [$1..N]); # Robert Israel, May 21 2020
  • Mathematica
    Table[n (HarmonicNumber[2 n] - HarmonicNumber[n]), {n, 30}] // Denominator
  • PARI
    a(n) = denominator(n*sum(k=1, n, 1/(n+k))); \\ Michel Marcus, Dec 14 2017
    
  • SageMath
    [denominator(n*(harmonic_number(2*n,1) - harmonic_number(n,1))) for n in range(1,41)] # G. C. Greubel, Jul 24 2023

Formula

From G. C. Greubel, Jul 24 2023: (Start)
a(n) = 2*A117664(n).
a(n) = 2*A111876(n-1)/n. (End)
Showing 1-3 of 3 results.