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-10 of 23 results. Next

A103345 Numerator of Sum_{k=1..n} 1/k^6 = Zeta(6,n).

Original entry on oeis.org

1, 65, 47449, 3037465, 47463376609, 47464376609, 5584183099672241, 357389058474664049, 260537105518334091721, 52107472322919827957, 92311616995117182948130877, 92311647383100199924330877, 445570781131605573859221176881493, 445570839299219762020391212081493
Offset: 1

Views

Author

Wolfdieter Lang, Feb 15 2005

Keywords

Comments

For the rationals Zeta(k,n) for k = 1..10 and n = 1..20, see the W. Lang link.
a(n) gives the partial sum, Zeta(6,n), of Euler's (later Riemann's) Zeta(6). Zeta(k,n), k >= 2, is sometimes also called H(k,n) because for k = 1 these would be the harmonic numbers A001008/A002805. However, H(1,n) does not give partial sums of a convergent series.

Examples

			The first few fractions are 1, 65/64, 47449/46656, 3037465/2985984, 47463376609/46656000000, ... = A103345/A103346. - _Petros Hadjicostas_, May 10 2020
		

Crossrefs

Cf. A013664, A291456. For the denominators, see A103346.

Programs

Formula

a(n) = numerator(Sum_{k=1..n} 1/k^6) = numerator(A291456(n)/(n!)^6).
G.f. for rationals Zeta(6, n): polylogarithm(6, x)/(1-x).
Zeta(6, n) = (1/945)*Pi^6 - psi(5, n+1)/5!, see eq. 6.4.3 with m = 5, p. 260, of the Abramowitz-Stegun reference. - Wolfdieter Lang, Dec 03 2013

A099827 Generalized harmonic number H(n,5) = Sum_{k=1..n} 1/k^5 multiplied by (n!)^5.

Original entry on oeis.org

0, 1, 33, 8051, 8252000, 25795462624, 200610400564224, 3371852494046112768, 110492114540967125581824, 6524555433591956305924325376, 652461835742417609568446054400000, 105080260346474296336209157187174400000
Offset: 0

Views

Author

Alexander Adamchuk, Oct 27 2004

Keywords

Comments

Note that a(n) is divisible by n, except when n is prime. Also, a(n+1) is divisible by n, except when n is prime or n = 0.

Examples

			a(2) = (2!)^5 * (1 + 1/2^5) = 2^5 + 1 = 33,
a(3) = (3!)^5 * (1 + 1/2^5 + 1/3^5) = 6^5 + 3^5 + 1 = 8051.
		

Crossrefs

Column k = 5 of A291556.

Programs

  • Mathematica
    Table[(n!)^5*Sum[1/k^5, {k, 1, n}], {n, 0, 13}] or Table[(n!)^5*HarmonicNumber[n, 5], {n, 0, 13}]

Formula

a(n) = (n!)^5 * Sum_{k=1..n} 1/k^5 = (n!)^5 * HarmonicNumber[n, 5] = (n!)^5 * A099828(n)/A069052(n).
a(0) = 0, a(1) = 1, a(n+1) = (n^5 + (n+1)^5)*a(n) - n^10*a(n-1) for n > 0. - Seiichi Manyama, Aug 24 2017
a(n) ~ Zeta(5) * (2*Pi)^(5/2) * n^(5*n+5/2) / exp(5*n). - Vaclav Kotesovec, Aug 27 2017
Sum_{n>=0} a(n) * x^n / (n!)^5 = polylog(5,x) / (1 - x). - Ilya Gutkovskiy, Jul 14 2020

Extensions

a(0) = 0 prepended by Seiichi Manyama, Aug 23 2017
Name edited by Petros Hadjicostas, May 10 2020

A136675 Numerator of Sum_{k=1..n} (-1)^(k+1)/k^3.

Original entry on oeis.org

1, 7, 197, 1549, 195353, 194353, 66879079, 533875007, 14436577189, 14420574181, 19209787242911, 19197460851911, 42198121495296467, 6025866788581781, 6027847576222613, 48209723660000029, 236907853607882606477
Offset: 1

Views

Author

Alexander Adamchuk, Jan 16 2008

Keywords

Comments

a(n) is prime for n in A136683.
Lim_{n -> infinity} a(n)/A334582(n) = A197070. - Petros Hadjicostas, May 07 2020

Examples

			The first few fractions are 1, 7/8, 197/216, 1549/1728, 195353/216000, 194353/216000, 66879079/74088000, 533875007/592704000, ... = a(n)/A334582(n). - _Petros Hadjicostas_, May 06 2020
		

Crossrefs

Programs

  • Maple
    map(numer,ListTools:-PartialSums([seq((-1)^(k+1)/k^3, k=1..100)])); # Robert Israel, Nov 09 2023
  • Mathematica
    (* Program #1 *) Table[Numerator[Sum[(-1)^(k+1)/k^3, {k,1,n}]], {n,1,50}]
    (* Program #2 *) Numerator[Accumulate[Table[(-1)^(k+1) 1/k^3, {k,50}]]] (* Harvey P. Dale, Feb 12 2013 *)
  • PARI
    a(n) = numerator(sum(k=1, n, (-1)^(k+1)/k^3)); \\ Michel Marcus, May 07 2020

A136677 Numerator of Sum_{k=1..n} (-1)^(k+1)/k^6.

Original entry on oeis.org

1, 63, 45991, 2942695, 45982595359, 5109066151, 601081707598999, 38469080386820311, 252396118308232060471, 252395862211967012407, 447134922152359540530757327, 447134770212444455649757327, 2158234586764514215343657417779543, 308319185132349039219686748825649
Offset: 1

Views

Author

Alexander Adamchuk, Jan 16 2008

Keywords

Comments

p divides a(p-1) for prime p > 2. a(n) is prime for n = {19, 47, 164, ...} = A136686.
Lim_{n -> infinity} a(n)/A334605(n) = A275703 = (31/32)*A013664. - Petros Hadjicostas, May 07 2020

Examples

			The first few fractions are 1, 63/64, 45991/46656, 2942695/2985984, 45982595359/46656000000, 5109066151/5184000000, ... = a(n)/A334605(n). - _Petros Hadjicostas_, May 07 2020
		

Crossrefs

Programs

  • Mathematica
    Table[ Numerator[ Sum[ (-1)^(k+1)/k^6, {k,1,n} ] ], {n,1,30} ]
    Accumulate[Table[(-1)^(k+1)/k^6,{k,20}]]//Numerator (* Harvey P. Dale, Aug 21 2023 *)

A136676 Numerator of Sum_{k=1..n} (-1)^(k+1)/k^5.

Original entry on oeis.org

1, 31, 7565, 241837, 755989457, 755889457, 12705011703799, 406547611705943, 98792790681344149, 98791774426324117, 15910615688635928566967, 15910549913780913466967, 5907492176026179821253778331
Offset: 1

Views

Author

Alexander Adamchuk, Jan 16 2008

Keywords

Comments

a(n) is prime for n in A136685.
Lim_{n -> infinity} a(n)/A334604(n) = A267316 = (15/16)*A013663. - Petros Hadjicostas, May 07 2020

Examples

			The first few fractions are 1, 31/32, 7565/7776, 241837/248832, 755989457/777600000, 755889457/777600000, ... = a(n)/A334604(n). - _Petros Hadjicostas_, May 07 2020
		

Crossrefs

Programs

  • Mathematica
    Table[ Numerator[ Sum[ (-1)^(k+1)/k^5, {k,1,n} ] ], {n,1,30} ]
  • PARI
    a(n) = numerator(sum(k=1, n, (-1)^(k+1)/k^5)); \\ Michel Marcus, May 07 2020

A069052 Denominator of Sum_{i = 1..n} 1/i^5.

Original entry on oeis.org

1, 32, 7776, 248832, 777600000, 259200000, 4356374400000, 139403980800000, 101625502003200000, 101625502003200000, 16366888723117363200000, 16366888723117363200000, 6076911214672415134617600000
Offset: 1

Views

Author

Benoit Cloitre, Apr 03 2002

Keywords

Comments

If 1 <= n <= 19, a(n) = A007480(n) * A002805(n) = denominator(Sum_{i = 1..n} 1/i^4) * denominator(Sum_{i = 1..n} 1/i).

Examples

			The first few fractions are 1, 33/32, 8051/7776, 257875/248832, ... = A099828/A069052. - _Petros Hadjicostas_, May 10 2020
		

Crossrefs

Numerators are A099828.

Programs

Formula

a(n) = denominator(Sum_{k=1..n} 1/k^5) = denominator(A099827(n)/(n!)^5). - Petros Hadjicostas, May 10 2020

A119722 Numerator of generalized harmonic number H(p-1,p)= Sum[ 1/k^p, {k,1,p-1}] divided by p^3 for prime p>3.

Original entry on oeis.org

2063, 2743174627, 19563315706517008974432827112201617, 2597378078067393746941400113704449589199274012223316613, 777478358612529699991463948563778410644748121498526065585976638854277886379480749840301120148933
Offset: 3

Views

Author

Alexander Adamchuk, Jun 13 2006

Keywords

Comments

Generalized harmonic number is H(n,m)= Sum[ 1/k^m, {k,1,n} ]. The numerator of generalized harmonic number H(p-1,p) is divisible by p^3 for prime p>3.

Examples

			Prime[3] = 5.
a(3) = numerator[ 1 + 1/2^5 + 1/3^5 + 1/4^5 ] / 5^3 = 257875/125 = 2063.
Prime[4] = 7
a(4) = numerator[ 1 + 1/2^7 + 1/3^7 + 1/4^7 + 1/5^7 + 1/6^7 ] / 7^3 = 2743174627.
		

Crossrefs

Programs

  • Mathematica
    Numerator[Table[Sum[1/k^Prime[n],{k,1,Prime[n]-1}],{n,3,9}]]/Table[Prime[n]^3,{n,3,9}]

Formula

a(n) = numerator[ Sum[ 1/k^Prime[n], {k,1,Prime[n]-1} ]] / Prime[n]^3 for n>2.

A136681 Numbers k such that A058313(k) is prime.

Original entry on oeis.org

3, 4, 5, 6, 9, 10, 13, 16, 17, 18, 37, 43, 58, 121, 124, 126, 137, 203, 247, 283, 285, 286, 289, 317, 424, 508, 751, 790, 937, 958, 1066, 1097, 1151, 1166, 1194, 1199, 1235, 1414, 1418, 1460, 1498, 1573, 2090, 2122, 2691, 2718, 3030, 3426, 3600, 3653, 3737
Offset: 1

Views

Author

Alexander Adamchuk, Jan 16 2008

Keywords

Comments

A058313(k) = Numerator of Sum_{j=1..k} (-1)^(j+1)/j.

Crossrefs

Programs

  • Mathematica
    Do[ f=Numerator[ Sum[ (-1)^(k+1)*1/k, {k,1,n} ] ]; If[ PrimeQ[f], Print[ {n,f} ] ], {n,1,317} ]
  • PARI
    isok(n) = isprime(numerator(sum(k=1, n, (-1)^(k+1)/k))); \\ Michel Marcus, Mar 14 2019

Extensions

a(25)-a(30) from James R. Buddenhagen, Sep 22 2015
a(31)-a(51) from Amiram Eldar, Mar 14 2019

A136682 Numbers k such that A119682(k) is prime.

Original entry on oeis.org

2, 3, 5, 8, 23, 41, 47, 48, 49, 95, 125, 203, 209, 284, 323, 395, 504, 553, 655, 781, 954, 1022, 1474, 1797, 1869, 1923, 1934, 1968, 2043, 2678, 3413, 3439, 4032, 4142, 4540, 4895, 5018, 5110, 5194, 5357, 6591, 11504, 11949, 14084, 20365
Offset: 1

Views

Author

Alexander Adamchuk, Jan 16 2008

Keywords

Comments

A119682(k) = Numerator of Sum_{j=1..k} (-1)^(j+1)/j^2.

Crossrefs

Programs

  • Mathematica
    Do[ f=Numerator[ Sum[ (-1)^(k+1)*1/k^2, {k,1,n} ] ]; If[ PrimeQ[f], Print[ {n,f} ] ], {n,1,125} ]

Extensions

a(12)-a(17) from Alexander Adamchuk, Apr 28 2008
a(18)-a(31) from Amiram Eldar, Mar 14 2019
a(32)-a(45) from Robert Price, Apr 14 2019

A136683 Numbers k such that A136675(k) is prime.

Original entry on oeis.org

2, 3, 4, 5, 6, 9, 20, 21, 29, 119, 132, 151, 351, 434, 457, 462, 572, 611, 930, 1107, 1157, 1452, 1515, 2838, 3997, 5346, 6463, 6725, 7664, 10234, 14168, 14299
Offset: 1

Views

Author

Alexander Adamchuk, Jan 16 2008

Keywords

Comments

A136675(k) = Numerator of Sum_{j=1..k} (-1)^(j+1)/j^3.

Crossrefs

Programs

  • Mathematica
    Do[ f=Numerator[ Sum[ (-1)^(k+1)*1/k^3, {k,1,n} ] ]; If[ PrimeQ[f], Print[ {n,f} ] ], {n,1,151} ]
    Flatten[Position[Numerator[Accumulate[Table[(-1)^(k+1) 1/k^3,{k,3000}]]],?PrimeQ] ] (* _Harvey P. Dale, Feb 12 2013 *)
  • PARI
    isok(n) = ispseudoprime(numerator(sum(k=1, n, (-1)^(k+1) / k^3))); \\ Daniel Suteu, Mar 15 2019

Extensions

More terms from Harvey P. Dale, Feb 12 2013
a(25)-a(28) from Amiram Eldar, Mar 15 2019
a(29)-a(32) from Robert Price, Apr 22 2019
Showing 1-10 of 23 results. Next