A099827 -id:A099827 - cp's OEIS frontend

A099828 Numerator of the generalized harmonic number H(n,5) = Sum_{k=1..n} 1/k^5.

1, 33, 8051, 257875, 806108207, 268736069, 4516906311683, 144545256245731, 105375212839937899, 105376229094957931, 16971048697474072945481, 16971114472329088045481, 6301272372663207205033976933

Offset: 1

Alexander Adamchuk, Oct 27 2004

From Alexander Adamchuk, Nov 07 2006: (Start)
a(n) is prime for n = {23, 25, 85, 147, 167, ...}.
There is a Wolstenholme-like theorem: p divides a(p-1) for prime p and p^2 divides a(p-1) for prime p > 7.
Also, p^3 divides a(p-1) for prime p = 5; p divides a((p-1)/2) for prime p = 37; p divides a((p-1)/3) for prime p = 37; p divides a((p-1)/4) for prime p = 37; p divides a((p-1)/5) for prime p = 11; p^2 divides a((p-1)/6) for prime p = 37; p divides a((p+1)/4) for prime p = 83; p divides a((p+1)/5) for prime p = 29; and p divides a((p+1)/6) for prime p = 11. (End)
See the Wolfdieter Lang link for information about Zeta(k, n) = H(n, k) with the rationals for k = 1..10, g.f.s, and polygamma formulas. - Wolfdieter Lang, Dec 03 2013

			H(n,5) = {1, 33/32, 8051/7776, 257875/248832, ... } = A099828/A069052.
For example, a(2) = numerator(1 + 1/2^5) = numerator(33/32) = 33 and a(3) = numerator(1 + 1/2^5 + 1/3^5) = numerator(8051/7776) = 8051. [Edited by _Petros Hadjicostas_, May 10 2020]

Alexander Adamchuk, Nov 07 2006, Table of n, a(n) for n = 1..100
Wolfdieter Lang, Rational Zeta(k,n) and more.
Eric Weisstein's World of Mathematics, Wolstenholme's Theorem.

Denominators are A069052.
A099827 = H(n,5) multiplied by (n!)^5.
Cf. A001008, A007406, A007408, A007410.

Numerator[Table[Sum[1/k^5, {k, 1, n}], {n, 1, 20}]]
Numerator[Table[HarmonicNumber[n, 5], {n, 1, 20}]]
Table[Numerator[Sum[1/k^5,{k,1,n}]],{n,1,100}] (* Alexander Adamchuk, Nov 07 2006 *)

a(n) = numerator(sum(k=1, n, 1/k^5)); \\ Michel Marcus, May 10 2020

a(n) = numerator(Sum_{k=1..n} 1/k^5) = numerator(HarmonicNumber[n, 5]).

A291556 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) = (n!)^k * Sum_{i=1..n} 1/i^k.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 1, 3, 3, 0, 1, 5, 11, 4, 0, 1, 9, 49, 50, 5, 0, 1, 17, 251, 820, 274, 6, 0, 1, 33, 1393, 16280, 21076, 1764, 7, 0, 1, 65, 8051, 357904, 2048824, 773136, 13068, 8, 0, 1, 129, 47449, 8252000, 224021776, 444273984, 38402064, 109584, 9

Offset: 0

Seiichi Manyama, Aug 26 2017

			Square array begins:
   0,  0,   0,     0,      0, ...
   1,  1,   1,     1,      1, ...
   2,  3,   5,     9,     17, ...
   3, 11,  49,   251,   1393, ...
   4, 50, 820, 16280, 357904, ...

Seiichi Manyama, Antidiagonals n = 0..59, flattened

Columns k=0-10 give: A001477, A000254, A001819, A066989, A203229, A099827, A291456, A291505, A291506, A291507, A291508.
Rows n=0-3 give: A000004, A000012, A000051, A074528.
Main diagonal gives A060943.

A:= (n, k)-> n!^k * add(1/i^k, i=1..n):
seq(seq(A(n, d-n), n=0..d), d=0..10);  # Alois P. Heinz, Aug 26 2017

A[0, ] = 0; A[1, ] = 1; A[n_, k_] := A[n, k] = ((n-1)^k + n^k) A[n-1, k] - (n-1)^(2k) A[n-2, k];
Table[A[n-k, k], {n, 0, 10}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, May 11 2019 *)

A(0, k) = 0, A(1, k) = 1, A(n+1, k) = (n^k+(n+1)^k)*A(n, k) - n^(2*k)*A(n-1, k).

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

Benoit Cloitre, Apr 03 2002

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

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

Wolfdieter Lang, Rational Zeta(k,n) and more.

Numerators are A099828.

Cf. A002805, A007480, A099827.

s = 0; lst = {}; Do[s += 1/n^5; AppendTo[lst, Denominator[s]], {n, 3 * 4!}]; lst (* Vladimir Joseph Stephan Orlovsky, Jan 24 2009 *)
Denominator[Table[Sum[1/i^5, {i, n}], {n, 20}]] (* Alonso del Arte, May 10 2020 *)

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

Alexander Adamchuk, Jun 13 2006

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.

			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.

Eric Weisstein's World of Mathematics, Wolstenholme's Theorem
Eric Weisstein's World of Mathematics, Harmonic Number

Cf. A099828, A099827, A001008, A007406, A007408, A007410.

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

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

A291456 a(n) = (n!)^6 * Sum_{i=1..n} 1/i^6.

Original entry on oeis.org

0, 1, 65, 47449, 194397760, 3037656102976, 141727869124448256, 16674281388691716870144, 4371079210518164503303028736, 2322975003299339366419974718488576, 2322977286679362958150790503464960000000

Offset: 0

Seiichi Manyama, Aug 24 2017

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..104

Column k=6 of A291556.
Cf. A000254 (k=1), A001819 (k=2), A066989 (k=3), A203229 (k=4), A099827 (k=5).
Cf. A008516, A103345, A103346.

Table[(n!)^6 * Sum[1/i^6, {i, 1, n}], {n, 0, 15}] (* Vaclav Kotesovec, Aug 27 2017 *)

a(0) = 0, a(1) = 1, a(n+1) = (n^6 + (n+1)^6)*a(n) - n^12*a(n-1) for n > 0.
a(n) ~ 8 * Pi^9 * n^(6*n+3) / (945 * exp(6*n)). - Vaclav Kotesovec, Aug 27 2017
a(n) = (n!)^6 * A103345(n)/A103346(n). - Petros Hadjicostas, May 10 2020
Sum_{n>=0} a(n) * x^n / (n!)^6 = polylog(6,x) / (1 - x). - Ilya Gutkovskiy, Jul 15 2020

A291505 a(n) = (n!)^7 * Sum_{i=1..n} 1/i^7.

Original entry on oeis.org

0, 1, 129, 282251, 4624680320, 361307736471424, 101143400834944548864, 83296040059942781485105152, 174684539610200377980575079727104, 835510910973061065615656036610946891776, 8355109938323553617123838798161699143680000000

Offset: 0

Seiichi Manyama, Aug 25 2017

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..92

Cf. A000254 (k=1), A001819 (k=2), A066989 (k=3), A099827 (k=5), A291456 (k=6), this sequence (k=7), A291506 (k=8), A291507 (k=9), A291508 (k=10).

Column k=7 of A291556.

Table[(n!)^7 * Sum[1/i^7, {i, 1, n}], {n, 0, 15}] (* Vaclav Kotesovec, Aug 27 2017 *)

a(n) = n!^7*sum(i=1, n, 1/i^7); \\ Michel Marcus, Aug 26 2017

a(0) = 0, a(1) = 1, a(n+1) = (n^7+(n+1)^7)*a(n) - n^14*a(n-1) for n > 0.
a(n) ~ zeta(7) * (2*Pi)^(7/2) * n^(7*n+7/2) / exp(7*n). - Vaclav Kotesovec, Aug 27 2017
Sum_{n>=0} a(n) * x^n / (n!)^7 = polylog(7,x) / (1 - x). - Ilya Gutkovskiy, Jul 15 2020

A291506 a(n) = (n!)^8 * Sum_{i=1..n} 1/i^8.

Original entry on oeis.org

0, 1, 257, 1686433, 110523752704, 43173450975314176, 72514862031522895036416, 418033821374598847702425993216, 7013444132843374500928464765799366656, 301905779820559925981495987360836056017534976

Offset: 0

Seiichi Manyama, Aug 25 2017

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..83

Cf. A000254 (k=1), A001819 (k=2), A066989 (k=3), A099827 (k=5), A291456 (k=6), A291505 (k=7), this sequence (k=8), A291507 (k=9), A291508 (k=10).

Column k=8 of A291556.

Table[(n!)^8 * Sum[1/i^8, {i, 1, n}], {n, 0, 15}] (* Vaclav Kotesovec, Aug 27 2017 *)

a(n) = n!^8*sum(i=1, n, 1/i^8); \\ Michel Marcus, Aug 26 2017

a(0) = 0, a(1) = 1, a(n+1) = (n^8+(n+1)^8)*a(n) - n^16*a(n-1) for n > 0.
a(n) ~ 8 * Pi^12 * n^(8*n+4) / (4725 * exp(8*n)). - Vaclav Kotesovec, Aug 27 2017
Sum_{n>=0} a(n) * x^n / (n!)^8 = polylog(8,x) / (1 - x). - Ilya Gutkovskiy, Jul 15 2020

A291507 a(n) = (n!)^9 * Sum_{i=1..n} 1/i^9.

Original entry on oeis.org

0, 1, 513, 10097891, 2647111616000, 5170142516807540224, 52103129720841632885243904, 2102549272223560776918400601161728, 282199388424234851655058321255905292713984, 109329825340451764123791003609208862665771818418176

Offset: 0

Seiichi Manyama, Aug 25 2017

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..75

Cf. A000254 (k=1), A001819 (k=2), A066989 (k=3), A099827 (k=5), A291456 (k=6), A291505 (k=7), A291506 (k=8), this sequence (k=9), A291508 (k=10).

Column k=9 of A291556.

Table[(n!)^9 * Sum[1/i^9, {i, 1, n}], {n, 0, 12}] (* Vaclav Kotesovec, Aug 27 2017 *)

a(n) = n!^9*sum(i=1, n, 1/i^9); \\ Michel Marcus, Aug 26 2017

a(0) = 0, a(1) = 1, a(n+1) = (n^9+(n+1)^9)*a(n) - n^18*a(n-1) for n > 0.
a(n) ~ zeta(9) * (2*Pi)^(9/2) * n^(9*n+9/2) / exp(9*n). - Vaclav Kotesovec, Aug 27 2017
Sum_{n>=0} a(n) * x^n / (n!)^9 = polylog(9,x) / (1 - x). - Ilya Gutkovskiy, Jul 15 2020

A291508 a(n) = (n!)^10 * Sum_{i=1..n} 1/i^10.

Original entry on oeis.org

0, 1, 1025, 60526249, 63466432537600, 619789443653380965376, 37476298202061058687475122176, 10586126703664512292193022557971021824, 11366767006463449393869821987386636472445566976, 39633465899293694663690352980684333029782095493517541376

Offset: 0

Seiichi Manyama, Aug 25 2017

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..69

Cf. A000254 (k=1), A001819 (k=2), A066989 (k=3), A099827 (k=5), A291456 (k=6), A291505 (k=7), A291506 (k=8), A291507 (k=9), this sequence (k=10).

Column k=10 of A291556.

Table[(n!)^10 * Sum[1/i^10, {i, 1, n}], {n, 0, 12}] (* Vaclav Kotesovec, Aug 27 2017 *)

a(n) = n!^10*sum(i=1, n, 1/i^10); \\ Michel Marcus, Aug 26 2017

a(0) = 0, a(1) = 1, a(n+1) = (n^10+(n+1)^10)*a(n) - n^20*a(n-1) for n > 0.
a(n) ~ 32 * Pi^15 * n^(10*n+5) / (93555 * exp(10*n)). - Vaclav Kotesovec, Aug 27 2017
Sum_{n>=0} a(n) * x^n / (n!)^10 = polylog(10,x) / (1 - x). - Ilya Gutkovskiy, Jul 15 2020

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

Original entry on oeis.org

2479157521, 159936660724017234488561, 1119583852472161859174156302552583713828739479026834819554843860744244189

Offset: 3

Alexander Adamchuk, Jul 08 2006

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

			With prime(3) = 5, a(3) = numerator[ 1 + 1/2^10 + 1/3^10 + 1/4^10 ] / 5^2 = 61978938025 / 25 = 2479157521.

Alexander Adamchuk, First 5 terms.
Eric Weisstein's World of Mathematics, Harmonic Number.
Eric Weisstein's World of Mathematics, Wolstenholme's Theorem.

Cf. A119722, A099828, A099827, A001008, A007406, A007408, A007410.

Table[Numerator[Sum[1/k^(2*Prime[n]),{k,1,Prime[n]-1}]],{n,3,7}]/Table[Prime[n]^2,{n,3,7}]

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

cp's OEIS Frontend

A099828 Numerator of the generalized harmonic number H(n,5) = Sum_{k=1..n} 1/k^5.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A291556 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) = (n!)^k * Sum_{i=1..n} 1/i^k.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A291456 a(n) = (n!)^6 * Sum_{i=1..n} 1/i^6.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A291505 a(n) = (n!)^7 * Sum_{i=1..n} 1/i^7.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A291506 a(n) = (n!)^8 * Sum_{i=1..n} 1/i^8.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A291507 a(n) = (n!)^9 * Sum_{i=1..n} 1/i^9.

Views

Author

Keywords

Links