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 26 results. Next

A002805 Denominators of harmonic numbers H(n) = Sum_{i=1..n} 1/i.

Original entry on oeis.org

1, 2, 6, 12, 60, 20, 140, 280, 2520, 2520, 27720, 27720, 360360, 360360, 360360, 720720, 12252240, 4084080, 77597520, 15519504, 5173168, 5173168, 118982864, 356948592, 8923714800, 8923714800, 80313433200, 80313433200, 2329089562800, 2329089562800, 72201776446800
Offset: 1

Views

Author

Keywords

Comments

H(n)/2 is the maximal distance that a stack of n cards can project beyond the edge of a table without toppling.
If n is not in {1, 2, 6} then a(n) has at least one prime factor other than 2 or 5. E.g., a(5) = 60 has a prime factor 3 and a(7) = 140 has a prime factor 7. This implies that every H(n) = A001008(n)/A002805(n), n not from {1, 2, 6}, has an infinite decimal representation. For a proof see the J. Havil reference. - Wolfdieter Lang, Jun 29 2007
a(n) = A213999(n,n-1). - Reinhard Zumkeller, Jul 03 2012
From Wolfdieter Lang, Apr 16 2015: (Start)
a(n)/A001008(n) = 1/H(n) is the solution of the following version of the classical cistern and pipes problem. A cistern is connected to n different pipes of water. For the k-th pipe it takes k time units (say, days) to fill the empty cistern, for k = 1, 2, ..., n. How long does it take for the n pipes together to fill the empty cistern? 1/H(n) gives the answer as a fraction of the time unit.
In general, if the k-th pipe needs d(k) days to fill the empty cistern then all pipes together need 1/Sum_{k=1..n} 1/d(k) = HM(d(1), ..., d(n))/n days, where HM denotes the harmonic mean HM. For the described problem, HM(1, 2, ..., n)/n = A102928(n)/(n*A175441(n)) = 1/H(n).
For a classical cistern and pipes problem see, e.g., the Hunger-Vogel reference (in Greek and German) given in A256101, problem 27, p. 29, where n = 3, and d(1), d(2) and d(3) are 6, 4 and 3 days. On p. 97 of this reference one finds remarks on the history of such problems (called in German 'Brunnenaufgabe'). (End)
From Wolfdieter Lang, Apr 17 2015: (Start)
An example of the above mentioned cistern and pipes problems appears in Chiu Chang Suan Shu (nine books on arithmetic) in book VI, problem 26. The numbers are there 1/2, 1, 5/2, 3 and 5 (days) and the result is 15/75 (day). See the reference (in German) on p. 68.
A historical account on such cistern problems is found in the Johannes Tropfke reference, given in A256101, section 4.2.1.2 Zisternenprobleme (Leistungsprobleme), pp. 578-579.
In Fibonacci's Liber Abaci such problems appear on p. 281 and p. 284 of L. E. Sigler's translation. (End)
All terms > 1 are even while corresponding numerators (A001008) are all odd (proof in Pólya and Szegő). - Bernard Schott, Dec 24 2021

Examples

			H(n) = [ 1, 3/2, 11/6, 25/12, 137/60, 49/20, 363/140, 761/280, 7129/2520, ... ] = A001008/A002805.
		

References

  • Chiu Chang Suan Shu, Neun Bücher arithmetischer Technik, translated and commented by Kurt Vogel, Ostwalds Klassiker der exakten Wissenschaften, Band 4, Friedr. Vieweg & Sohn, Braunschweig, 1968, p. 68.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 258-261.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 259.
  • J. Havil, Gamma, (in German), Springer, 2007, p. 35-6; Gamma: Exploring Euler's Constant, Princeton Univ. Press, 2003.
  • D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 1, p. 615.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, volume II, Springer, reprint of the 1976 edition, 1998, problem 251, p. 154.
  • L. E. Sigler, Fibonacci's Liber Abaci, Springer, 2003, pp. 281, 284.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A001008 (numerators), A075135, A025529, A203810, A203811, A203812.
Partial sums: A027612/A027611 = 1, 5/2, 13/3, 77/12, 87/10, 223/20,...
The following fractions are all related to each other: Sum 1/n: A001008/A002805, Sum 1/prime(n): A024451/A002110 and A106830/A034386, Sum 1/nonprime(n): A282511/A282512, Sum 1/composite(n): A250133/A296358, Sum 1/n^2: A007406/A007407, Sum 1/n^3: A007408/A007409.

Programs

  • GAP
    List([1..30],n->DenominatorRat(Sum([1..n],i->1/i))); # Muniru A Asiru, Dec 20 2018
    
  • Haskell
    import Data.Ratio ((%), denominator)
    a002805 = denominator . sum . map (1 %) . enumFromTo 1
    a002805_list = map denominator $ scanl1 (+) $ map (1 %) [1..]
    -- Reinhard Zumkeller, Jul 03 2012
    
  • Magma
    [Denominator(HarmonicNumber(n)): n in [1..40]]; // Vincenzo Librandi, Apr 16 2015
    
  • Maple
    seq(denom(sum((2*k-1)/k, k=1..n), n=1..30); # Gary Detlefs, Jul 18 2011
    f:=n->denom(add(1/k, k=1..n)); # N. J. A. Sloane, Nov 15 2013
  • Mathematica
    Denominator[ Drop[ FoldList[ #1 + 1/#2 &, 0, Range[ 30 ] ], 1 ] ] (* Harvey P. Dale, Feb 09 2000 *)
    Table[Denominator[HarmonicNumber[n]], {n, 1, 40}] (* Stefan Steinerberger, Apr 20 2006 *)
    Denominator[Accumulate[1/Range[25]]] (* Alonso del Arte, Nov 21 2018 *)
  • PARI
    a(n)=denominator(sum(k=2,n,1/k)) \\ Charles R Greathouse IV, Feb 11 2011
    
  • Python
    from fractions import Fraction
    def a(n): return sum(Fraction(1, i) for i in range(1, n+1)).denominator
    print([a(n) for n in range(1, 30)]) # Michael S. Branicky, Dec 24 2021
  • Sage
    def harmonic(a, b): # See the F. Johansson link.
        if b - a == 1 : return 1, a
        m = (a+b)//2
        p, q = harmonic(a,m)
        r, s = harmonic(m,b)
        return p*s+q*r, q*s
    def A002805(n) : H = harmonic(1,n+1); return denominator(H[0]/H[1])
    [A002805(n) for n in (1..29)] # Peter Luschny, Sep 01 2012
    

Formula

a(n) = Denominator(Sum_{k=1..n} (2*k-1)/k). - Gary Detlefs, Jul 18 2011
a(n) = n! / gcd(Stirling1(n+1, 2), n!) = n! / gcd(A000254(n),n!). - Max Alekseyev, Mar 01 2018
a(n) = the (reduced) denominator of the continued fraction 1/(1 - 1^2/(3 - 2^2/(5 - 3^2/(7 - ... - (n-1)^2/(2*n-1))))). - Peter Bala, Feb 18 2024

Extensions

Definition edited by Daniel Forgues, May 19 2010

A098849 a(n) = n*(n + 16).

Original entry on oeis.org

0, 17, 36, 57, 80, 105, 132, 161, 192, 225, 260, 297, 336, 377, 420, 465, 512, 561, 612, 665, 720, 777, 836, 897, 960, 1025, 1092, 1161, 1232, 1305, 1380, 1457, 1536, 1617, 1700, 1785, 1872, 1961, 2052, 2145, 2240, 2337, 2436, 2537, 2640, 2745, 2852, 2961
Offset: 0

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Nov 04 2004

Keywords

Crossrefs

a(n-8), n>=9, eighth column (used for the n=8 series of the hydrogen atom) of triangle A120070.

Programs

Formula

a(n) = (n+8)^2 - 8^2 = n*(n + 16), n>=0.
G.f.: x*(17 - 15*x)/(1-x)^3.
a(n) = a(n-1) + 2*n + 15 (with a(0)=0). - Vincenzo Librandi, Nov 17 2010
From G. C. Greubel, Jul 29 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
E.g.f.: x*(17 + x)*exp(x). (End)
From Amiram Eldar, Jan 15 2021: (Start)
Sum_{n>=1} 1/a(n) = H(16)/16 = A001008(16)/A102928(16) = 2436559/11531520, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 95549/2306304. (End)

Extensions

More terms from Emeric Deutsch, Mar 26 2005

A132761 a(n) = n*(n+17).

Original entry on oeis.org

0, 18, 38, 60, 84, 110, 138, 168, 200, 234, 270, 308, 348, 390, 434, 480, 528, 578, 630, 684, 740, 798, 858, 920, 984, 1050, 1118, 1188, 1260, 1334, 1410, 1488, 1568, 1650, 1734, 1820, 1908, 1998, 2090, 2184, 2280, 2378, 2478, 2580, 2684, 2790, 2898, 3008, 3120
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Comments

a(n) is the first Zagreb index of the helm graph H[n] (n>=3). - Emeric Deutsch, Nov 05 2016
From Emeric Deutsch, Nov 07 2016: (Start)
a(n) is the first Zagreb index of the graph obtained by joining one vertex of the cycle graph C[n] with each vertex of a second cycle graph C[n].
The first Zagreb index of a simple connected graph is the sum of the squared degrees of its vertices. Alternately, it is the sum of the degree sums d(i) + d(j) over all edges ij of the graph. (End)
From Emeric Deutsch, May 11 2018: (Start)
The M-polynomial of the Helm graph H[n] is M(H[n];x,y) = n*x*y^4 + n*x^4*y^4 + n*x^4*y^n.
The helm graph H[n] is the graph obtained from an n-wheel graph by adjoining a pendant edge at each node of the cycle. (End)

Crossrefs

Programs

Formula

a(n) = n*(n + 17).
a(n) = A132760(n) + 2*n = A132765(n) - 6*n = A098849(n) + n = A120071(n) - 3*n. - Zerinvary Lajos, Feb 17 2008
a(n) = 2*n + a(n-1) + 16 for n > 0, a(0)=0. - Vincenzo Librandi, Aug 03 2010
G.f.: 2*x*(9 - 8*x)/(1 - x)^3. - Emeric Deutsch, Nov 07 2016
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(17)/17 = A001008(17)/A102928(17) = 42142223/208288080, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/17 - 1768477/41657616. (End)
From Elmo R. Oliveira, Dec 12 2024: (Start)
E.g.f.: exp(x)*x*(18 + x).
a(n) = 2*A056126(n).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2. (End)

A098850 a(n) = n*(n + 18).

Original entry on oeis.org

0, 19, 40, 63, 88, 115, 144, 175, 208, 243, 280, 319, 360, 403, 448, 495, 544, 595, 648, 703, 760, 819, 880, 943, 1008, 1075, 1144, 1215, 1288, 1363, 1440, 1519, 1600, 1683, 1768, 1855, 1944, 2035, 2128, 2223, 2320, 2419, 2520, 2623, 2728, 2835, 2944, 3055
Offset: 0

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Nov 04 2004

Keywords

Crossrefs

a(n-9), n>=10, ninth column (used for the n=9 series of the hydrogen atom) of triangle A120070.

Programs

Formula

a(n) = (n+9)^2 - 9^2 = n*(n + 18), n>=0.
G.f.: x*(19 - 17*x)/(1-x)^3.
a(n) = 2*n + a(n-1) + 17 (with a(0)=0). - Vincenzo Librandi, Nov 17 2010
From G. C. Greubel, Jul 29 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
E.g.f.: x*(19 + x)*exp(x). (End)
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(18)/18 = A001008(18)/A102928(18) = 14274301/73513440, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 1632341/44108064. (End)

Extensions

More terms from Emeric Deutsch, Mar 06 2005

A120071 a(n) = n*(n+20).

Original entry on oeis.org

0, 21, 44, 69, 96, 125, 156, 189, 224, 261, 300, 341, 384, 429, 476, 525, 576, 629, 684, 741, 800, 861, 924, 989, 1056, 1125, 1196, 1269, 1344, 1421, 1500, 1581, 1664, 1749, 1836, 1925, 2016, 2109, 2204, 2301, 2400, 2501, 2604, 2709, 2816, 2925, 3036, 3149, 3264
Offset: 0

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Crossrefs

a(n-10), n >= 11, tenth column (used for the n=10 series of the hydrogen atom) of triangle A120070.
For n*(n+18) see A098850.

Programs

Formula

a(n) = (n+10)^2 - 10^2 = n*(n+20), n >= 0.
G.f.: x*(21-19*x)/(1-x)^3.
a(n) = 2*n + a(n-1) + 19 (with a(0)=0). - Vincenzo Librandi, Nov 13 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(20)/20 = A001008(20)/A102928(20) = 11167027/62078016, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 155685007/4655851200. (End)
From Elmo R. Oliveira, Jan 12 2025: (Start)
E.g.f.: exp(x)*x*(21 + x).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n >= 3. (End)

A132765 a(n) = n*(n + 23).

Original entry on oeis.org

0, 24, 50, 78, 108, 140, 174, 210, 248, 288, 330, 374, 420, 468, 518, 570, 624, 680, 738, 798, 860, 924, 990, 1058, 1128, 1200, 1274, 1350, 1428, 1508, 1590, 1674, 1760, 1848, 1938, 2030, 2124, 2220, 2318, 2418, 2520, 2624, 2730, 2838, 2948, 3060, 3174, 3290, 3408
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

Formula

a(n) = n*(n + 23).
a(n) = 2*n + a(n-1) + 22 for n>0, a(0)=0. - Vincenzo Librandi, Aug 03 2010
From Chai Wah Wu, Dec 17 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
G.f.: 2*x*(12 - 11*x)/(1-x)^3. (End)
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(23)/23 = A001008(23)/A102928(23) = 444316699/2736605872, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/23 - 3825136961/123147264240. (End)
E.g.f.: x*(24 + x)*exp(x). - G. C. Greubel, Mar 14 2022

A132762 a(n) = n*(n + 19).

Original entry on oeis.org

0, 20, 42, 66, 92, 120, 150, 182, 216, 252, 290, 330, 372, 416, 462, 510, 560, 612, 666, 722, 780, 840, 902, 966, 1032, 1100, 1170, 1242, 1316, 1392, 1470, 1550, 1632, 1716, 1802, 1890, 1980, 2072, 2166, 2262, 2360, 2460, 2562, 2666, 2772, 2880, 2990, 3102, 3216
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

Formula

a(n) = 2*n + a(n-1) + 18 for n > 0, a(0) = 0. - Vincenzo Librandi, Aug 03 2010
From Chai Wah Wu, Dec 17 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
G.f.: 2*x*(10 - 9*x)/(1-x)^3. (End)
a(n) = 2*A051942(n+9). - R. J. Mathar, Sep 05 2018
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(19)/19 = A001008(19)/A102928(19) = 275295799/1474352880, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/19 - 33464927/884611728. (End)
E.g.f.: x*(20 + x)*exp(x). - G. C. Greubel, Mar 14 2022

A132764 a(n) = n*(n+22).

Original entry on oeis.org

0, 23, 48, 75, 104, 135, 168, 203, 240, 279, 320, 363, 408, 455, 504, 555, 608, 663, 720, 779, 840, 903, 968, 1035, 1104, 1175, 1248, 1323, 1400, 1479, 1560, 1643, 1728, 1815, 1904, 1995, 2088, 2183, 2280, 2379, 2480, 2583, 2688, 2795, 2904, 3015, 3128, 3243, 3360
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Examples

			a(1)=2*1+0+21=23; a(2)=2*2+23+21=48; a(3)=2*3+48+21=75. - _Vincenzo Librandi_, Aug 03 2010
		

Crossrefs

Programs

Formula

a(n) = n*(n + 22).
a(n) = 2*n + a(n-1) + 21 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=23, a(2)=48, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, May 02 2012
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(22)/22 = A001008(22)/A102928(22) = 19093197/113809696, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 156188887/5121436320. (End)
From G. C. Greubel, Mar 14 2022: (Start)
G.f.: x*(23 - 21*x)/(1-x)^3.
E.g.f.: x*(23 + x)*exp(x). (End)

A132763 a(n) = n*(n+21).

Original entry on oeis.org

0, 22, 46, 72, 100, 130, 162, 196, 232, 270, 310, 352, 396, 442, 490, 540, 592, 646, 702, 760, 820, 882, 946, 1012, 1080, 1150, 1222, 1296, 1372, 1450, 1530, 1612, 1696, 1782, 1870, 1960, 2052, 2146, 2242, 2340, 2440, 2542, 2646, 2752, 2860, 2970, 3082, 3196, 3312
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

Formula

a(n) = n*(n + 21).
a(n) = 2*n + a(n-1) + 20 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=22, a(2)=46, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, May 25 2014
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(21)/21 = A001008(21)/A102928(21) = 18858053/108636528, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/21 - 166770367/4888643760. (End)
From Stefano Spezia, Jan 30 2021: (Start)
O.g.f.: 2*x*(11 - 10*x)/(1 - x)^3.
E.g.f.: x*(22 + x)*exp(x). (End)

A132766 a(n) = n*(n+24).

Original entry on oeis.org

0, 25, 52, 81, 112, 145, 180, 217, 256, 297, 340, 385, 432, 481, 532, 585, 640, 697, 756, 817, 880, 945, 1012, 1081, 1152, 1225, 1300, 1377, 1456, 1537, 1620, 1705, 1792, 1881, 1972, 2065, 2160, 2257, 2356, 2457, 2560, 2665, 2772, 2881, 2992, 3105, 3220, 3337
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

  • Mathematica
    Table[n (n + 24), {n, 0, 50}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 25, 52}, 50] (* Harvey P. Dale, Feb 11 2016 *)
  • PARI
    a(n)=n*(n+24) \\ Charles R Greathouse IV, Jun 17 2017
    
  • Sage
    [n*(n+24) for n in (0..50)] # G. C. Greubel, Mar 14 2022

Formula

a(n) = n*(n + 24).
a(n) = 2*n + a(n-1) + 23 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=25, a(2)=52; for n>2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Feb 11 2016
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(24)/24 = A001008(24)/A102928(24) = 1347822955/8566766208, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 3602044091/128501493120. (End)
From G. C. Greubel, Mar 14 2022: (Start)
G.f.: 2*x*(13 - 12*x)/(1-x)^3.
E.g.f.: x*(26 + x)*exp(x). (End)
Showing 1-10 of 26 results. Next