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-5 of 5 results.

A058313 Numerator of the n-th alternating harmonic number, Sum_{k=1..n} (-1)^(k+1)/k.

Original entry on oeis.org

1, 1, 5, 7, 47, 37, 319, 533, 1879, 1627, 20417, 18107, 263111, 237371, 52279, 95549, 1768477, 1632341, 33464927, 155685007, 166770367, 156188887, 3825136961, 3602044091, 19081066231, 18051406831, 57128792093, 7751493599, 236266661971
Offset: 1

Views

Author

N. J. A. Sloane, Dec 09 2000

Keywords

Comments

A Wolstenholme-like theorem: for prime p > 3, if p = 6k-1, then p divides a(4k-1), otherwise if p = 6k+1, then p divides a(4k). - T. D. Noe, Apr 01 2004
For the limit n -> infinity of the partial sums of the alternating harmonic series see A002162. - Wolfdieter Lang, Sep 08 2015
a(n)/A058312(n) appears in the locker puzzle (see the links in A364317) as the probability of failures with the strategy used for n lockers and opening of up to floor(n/2) lockers. Note the alternative formula given below for a(n)/A058312(n) using only positive fractions. - Wolfdieter Lang, Aug 12 2023

Examples

			1, 1/2, 5/6, 7/12, 47/60, 37/60, 319/420, 533/840, 1879/2520, ...
For n=4: a(n)/A058312(n) = 7/12 because 1/1 - 1/2 + 1/3 - 1/4 = 7/12 = 1/4 + 1/3. - _Wolfdieter Lang_, Aug 12 2023
		

References

  • L. B. W. Jolley, Summation of Series, Dover Publications, 1961, page 14, #71.

Crossrefs

Denominators are A058312. Cf. A025530.
Apart from leading term, same as A075830.
Cf. A001008 (numerator of n-th harmonic number).
Bisections are A049281 and A082687.
Cf. A181983.

Programs

  • Haskell
    import Data.Ratio((%), numerator)
    a058313 n = a058313_list !! (n-1)
    a058313_list = map numerator $ scanl1 (+) $ map (1 %) $ tail a181983_list
    -- Reinhard Zumkeller, Mar 20 2013
  • Maple
    A058313 := n->numer(add((-1)^(k+1)/k,k=1..n));
    # Alternatively:
    a := n -> numer(harmonic(n) - harmonic((n-modp(n,2))/2)):
    seq(a(n), n=1..29); # Peter Luschny, May 03 2016
  • Mathematica
    Numerator[Table[Sum[(-1)^(k + 1)/k, {k, n}], {n, 30}]] (* Harvey P. Dale, Jul 18 2012 *)
    a[n_]:= (-1)^n (HarmonicNumber[n/2 - 1/2] - HarmonicNumber[n/2] + (-1)^n Log[4])/2; Table[a[n] // FullSimplify, {n, 29}] // Numerator (* Gerry Martens, Jul 05 2015 *)
    Rest[Numerator[CoefficientList[Series[Log[1 + x]/(1 - x), {x, 0, 33}], x]]] (* Vincenzo Librandi, Jul 06 2015 *)
    Table[Log[2] - (-1)^n LerchPhi[-1, 1, n + 1], {n, 20}] // Numerator (* Eric W. Weisstein, Aug 25 2023 *)
  • PARI
    a(n)=(-1)^n*numerator(polcoeff(log(1-x)/(x+1)+O(x^(n+1)), n))
    

Formula

G.f. for a(n)/A058312(n): log(1+x)/(1-x). - Benoit Cloitre, Jun 15 2003
a(n) = (n*a(n-1) + (-1)^(n+1)*A058312(n-1))/gcd(n*a(n-1) + (-1)^(n+1)*A058312(n-1), n*A058312(n-1)). - Robert Israel, Jul 06 2015
From Peter Luschny, May 03 2016: (Start)
Let H(n) denote the harmonic numbers, AH(n) denote the alternating harmonic numbers, Psi the polygamma function and euler(n,x) the Euler polynomials. Then:
AH(n) = H(n) - H((n - n mod 2)/2).
AH(z) = log(2)+(1/2)*cos(Pi*z)*(Psi(z/2+1/2)-Psi(z/2+1)).
AH(z) ~ log(2)+(1/2)*cos(Pi*z)*(-1/z+1/(2*z^2)-1/(4*z^4)+1/(2*z^6)-...).
AH(z) ~ log(2)-(1/2)*cos(Pi*z)*Sum_{n>=0} Euler(n,0)/z^(n+1). (End)
Sum_{k>=1} (-1)^(k+1)*AH(k)/k = Pi^2/12 + log(2)^2/2 (Boyadzhiev, 2013). - Amiram Eldar, Oct 04 2021
a(n)/A058312(n) = Sum_{j=0..ceiling(n/2) - 1} 1/(n-j), for n >= 1. (Proof by comparing the recurrences for even and odd n.) - Wolfdieter Lang, Aug 12 2023
For n >= 1, log(2) = a(n)/A058312(n) + (-1)^n*n!*Sum_{k >= 1} 1/(k*(k + 1)* ...*(k + n)*2^k). - Peter Bala, Dec 07 2023
a(n) = the (reduced) numerator of the continued fraction 1/(1 + 1^2/(1 + 2^2/(1 + 3^2/(1 + ... + (n-1)^2/(1))))). - Peter Bala, Feb 18 2024

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)).

A074029 Number of binary Lyndon words of length n with trace 1 and subtrace 0 over Z_2.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 4, 8, 15, 27, 48, 85, 155, 288, 541, 1024, 1935, 3654, 6912, 13107, 24940, 47616, 91136, 174760, 335626, 645435, 1242904, 2396745, 4627915, 8947294, 17317888, 33554432, 65076240, 126324495, 245428574, 477218560, 928638035, 1808400384, 3524068955
Offset: 1

Views

Author

Frank Ruskey and Nate Kube, Aug 21 2002

Keywords

Comments

Same as the number of binary Lyndon words of length n with trace 1 and subtrace 0 over GF(2).

Examples

			a(3;1,0)=1 since the one binary Lyndon word of trace 1, subtrace 0 and length 3 is { 001 }.
		

Crossrefs

Formula

a(2n) = A042982(2n), a(2n+1) = A049281(2n+1). This follows from Cattell et al. (see A042979), Main Theorem on p. 33 and Theorem 4 on p. 44.

Extensions

Terms a(33) onward from Max Alekseyev, Apr 09 2013

A354894 a(n) is the numerator of the n-th hyperharmonic number of order n.

Original entry on oeis.org

1, 5, 47, 319, 1879, 20417, 263111, 261395, 8842385, 33464927, 166770367, 3825136961, 19081066231, 57128792093, 236266661971, 7313175618421, 14606816124167, 102126365345729, 3774664307989373, 3771059091081773, 154479849447926113, 6637417807457499259, 6632660439700528339
Offset: 1

Views

Author

Ilya Gutkovskiy, Jun 10 2022

Keywords

Examples

			1, 5/2, 47/6, 319/12, 1879/20, 20417/60, 263111/210, 261395/56, 8842385/504, ...
		

References

  • J. H. Conway and R. K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996, p. 258.

Crossrefs

Differs from A049281.

Programs

  • Mathematica
    Table[SeriesCoefficient[-Log[1 - x]/(1 - x)^n, {x, 0, n}], {n, 1, 23}] // Numerator
    Table[Binomial[2 n - 1, n - 1] (HarmonicNumber[2 n - 1] - HarmonicNumber[n - 1]), {n, 1, 23}] // Numerator
  • PARI
    H(n) = sum(i=1, n, 1/i);
    a(n) = numerator(binomial(2*n-1,n-1) * (H(2*n-1) - H(n-1))); \\ Michel Marcus, Jun 10 2022
    
  • Python
    from math import comb
    from sympy import harmonic
    def A354894(n): return (comb(2*n-1,n-1)*(harmonic(2*n-1)-harmonic(n-1))).p # Chai Wah Wu, Jun 18 2022

Formula

a(n) is the numerator of the coefficient of x^n in the expansion of -log(1 - x) / (1 - x)^n.
a(n) is the numerator of binomial(2*n-1,n-1) * (H(2*n-1) - H(n-1)), where H(n) is the n-th harmonic number.
a(n) / A354895(n) ~ log(2) * 2^(2*n-1) / sqrt(Pi * n).

A069685 Denominators of coefficients in -log(1+x)log(1-x) power series.

Original entry on oeis.org

1, 12, 180, 1680, 12600, 166320, 2522520, 576576, 22054032, 465585120, 2560718160, 64250746560, 348024877200, 1124388064800, 4990906206000, 165032631878400, 350694342741600, 2599263952084800, 101515697684200800
Offset: 1

Views

Author

Benoit Cloitre, May 03 2002

Keywords

Crossrefs

Programs

  • PARI
    my(x='x+O('x^50), v=Vec(-log(1+x)*log(1-x))); apply(denominator, vector(#v\2, k, v[2*k-1])) \\ Michel Marcus, Feb 16 2021

Formula

a(n) = n*A076100(n) (conjectured). - F. Chapoton, Nov 09 2009
Showing 1-5 of 5 results.