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

A120077 Denominators of row sums of rational triangle A120072/A120073.

Original entry on oeis.org

4, 36, 144, 3600, 3600, 176400, 705600, 6350400, 1270080, 153679680, 153679680, 25971865920, 25971865920, 129859329600, 519437318400, 150117385017600, 150117385017600, 54192375991353600, 2167695039654144, 1548353599752960, 221193371393280, 117011293467045120
Offset: 2

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

The first 19 terms coincide with A007407(n), for n>=2. However a(20) = 2167695039654144 and A007407(20) = 10838475198270720 = 5*a(20). Also a(21) = 1548353599752960 and A007407(21) = 221193371393280 = a(21)/7. From n = 22 up to at least n = 100 (checked) both sequences coincide again.
See the W. Lang link under A120072 for more details.
The corresponding numerators are given by A120076.
The n for which a(n) differs from A007407(n) are given by A309829. - Jeppe Stig Nielsen, Aug 18 2019

Examples

			The rationals A120076(m)/a(m), m>=2, begin with (3/4, 37/36, 169/144, 4549/3600, 4769/3600, ... ).
		

Crossrefs

Programs

  • Magma
    A120077:= func< n | Denominator( (&+[1/k^2: k in [1..n]]) -1/n) >;
    [A120077(n): n in [2..30]]; // G. C. Greubel, Apr 25 2023
    
  • Mathematica
    Table[Denominator[HarmonicNumber[n,2] -1/n], {n,2,40}] (* G. C. Greubel, Apr 25 2023 *)
  • PARI
    a(n) = denominator(sum(j=1,n-1,1/j^2-1/n^2)) \\ Jeppe Stig Nielsen, Aug 18 2019
    
  • PARI
    a(n) = denominator(sum(j=1,n,1/j^2) - 1/n) \\ Jeppe Stig Nielsen, Aug 18 2019
    
  • SageMath
    def A120077(n): return denominator(harmonic_number(n,2) - 1/n)
    [A120077(n) for n in range(2,31)] # G. C. Greubel, Apr 25 2023

Formula

a(n) = denominator(r(m)), with the rationals r(m) = Sum_{n=1..m-1} A120072(m,n)/A120073(m,n), m >= 2.
The rationals are r(m) = Zeta(2; m-1) - (m-1)/m^2, m>=2, with the partial sums Zeta(2; n) = Sum_{k=1..n} 1/k^2. See the W. Lang link under A103345.
O.g.f. for the rationals r(m), m>=2: log(1-x) + polylog(2,x)/(1-x).

Extensions

a(21)-a(23) from Jeppe Stig Nielsen, Aug 18 2019

A172370 Mirrored triangle A120072 read by rows.

Original entry on oeis.org

3, 5, 8, 7, 3, 15, 9, 16, 21, 24, 11, 5, 1, 2, 35, 13, 24, 33, 40, 45, 48, 15, 7, 39, 3, 55, 15, 63, 17, 32, 5, 56, 65, 8, 77, 80, 19, 9, 51, 4, 3, 21, 91, 6, 99, 21, 40, 57, 72, 85, 96, 105, 112, 117, 120, 23, 11, 7, 5, 95, 1, 119, 1, 5, 35, 143, 25, 48, 69, 88, 105, 120, 133, 144
Offset: 2

Views

Author

Paul Curtz, Feb 01 2010

Keywords

Comments

A table of numerators of 1/n^2 - 1/m^2 extended to negative m looks as follows, stacked such that values of common m are aligned
and the central column of -1 is defined for m=0:
.............................0..-1...0...3...8..15..24..35..48..63..80..99. A005563
.........................0..-3..-1..-3...0...5...3..21...2..45..15..77...6. A061037
.....................0..-5..-8..-1..-8..-5...0...7..16...1..40..55...8..91. A061039
.................0..-7..-3.-15..-1.-15..-3..-7...0...9...5..33...3..65..21. A061041
.............0..-9.-16.-21.-24..-1.-24.-21.-16..-9...0..11..24..39..56...3. A061043
.........0.-11..-5..-1..-2.-35..-1.-35..-2..-1..-5.-11...0..13...7...5...4. A061045
.....0.-13.-24.-33.-40.-45.-48..-1.-48.-45.-40.-33.-24.-13...0..15..32..51. A061047
.0.-15..-7.-39..-3.-55.-15.-63..-1.-63.-15.-55..-3.-39..-7.-15...0..17...9. A061049
The row-reversed variant of A120072 appears (negated) after the leftmost 0.
Equals A061035 with the first column removed. - Georg Fischer, Jul 26 2023

Examples

			The table starts
   3
   5   8
   7   3  15
   9  16  21  24
  11   5   1   2  35
  13  24  33  40  45  48
  15   7  39   3  55  15  63
  17  32   5  56  65   8  77  80
  19   9  51   4   3  21  91   6  99
		

Crossrefs

Lower diagonal gives: A070262, A061037(n+2).

Programs

  • Magma
    [[Numerator(1/(n-k)^2 -1/n^2): k in [1..n-1]]: n in [2..20]]; // G. C. Greubel, Sep 20 2018
  • Mathematica
    Table[Numerator[1/(n-k)^2 -1/n^2], {n, 2, 20}, {k, 1, n-1}]//Flatten (* G. C. Greubel, Sep 20 2018 *)
  • PARI
    for(n=2,20, for(k=1,n-1, print1(numerator(1/(n-k)^2 -1/n^2), ", "))) \\ G. C. Greubel, Sep 20 2018
    

Formula

T(n,m) = numerator of 1/(n-m)^2 - 1/n^2, n >= 2, 1 <= m < n. - R. J. Mathar, Nov 23 2010

Extensions

Comment rewritten and offset set to 2 by R. J. Mathar, Nov 23 2010

A120074 Row sums of triangle A120072 (numerator triangle for H atom spectrum).

Original entry on oeis.org

3, 13, 25, 70, 54, 203, 197, 340, 303, 825, 445, 1378, 892, 1221, 1565, 3128, 1545, 4389, 2427, 3592, 3688, 7843, 3589, 8420, 6191, 9097, 7135, 15834, 5774, 19375, 12493, 14814, 14147, 19647, 12264, 33078
Offset: 2

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Crossrefs

Programs

  • Magma
    A120072:= func< n,k | Numerator(1/k^2 - 1/n^2) >;
    [(&+[A120072(n,k): k in [1..n-1]]): n in [2..50]]; // G. C. Greubel, Apr 24 2023
    
  • Mathematica
    Table[Sum[1/n^2 - 1/m^2 //Numerator, {n,m-1}], {m,2,40}]  (* Jean-François Alcover, Sep 16 2013 *)
  • SageMath
    def A120072(n,k): return numerator(1/k^2 - 1/n^2)
    [sum(A120072(n,k) for k in range(1,n)) for n in range(2,51)] # G. C. Greubel, Apr 24 2023

Formula

a(n) = Sum_{k=1..n-1} A120072(n,k) for n >= 2.

A120076 Numerators of row sums of rational triangle A120072/A120073.

Original entry on oeis.org

3, 37, 169, 4549, 4769, 241481, 989549, 9072541, 1841321, 225467009, 227698469, 38801207261, 39076419341, 196577627041, 790503882349, 229526961468061, 230480866420061, 83512167402400421, 3351610394325821
Offset: 2

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

The corresponding denominators are given by A120077.
See the W. Lang link under A120072 for more details.

Examples

			The rationals a(m)/A120077(m), m>=2, begin with (3/4, 37/36, 169/144, 4549/3600, 4769/3600, ...).
		

Crossrefs

Programs

  • Magma
    A120076:= func< n | Numerator( (&+[1/k^2: k in [1..n]]) -1/n) >;
    [A120076(n): n in [2..30]]; // G. C. Greubel, Apr 24 2023
    
  • Mathematica
    Table[Numerator[HarmonicNumber[n,2] -1/n], {n,2,40}] (* G. C. Greubel, Apr 24 2023 *)
  • SageMath
    def A120076(n): return numerator(harmonic_number(n,2) - 1/n)
    [A120076(n) for n in range(2,31)] # G. C. Greubel, Apr 24 2023

Formula

a(n) = numerator(r(m)), with the rationals r(m) = Sum_{n=1..m-1} A120072(m,n)/A120073(m,n), m >= 2.
The rationals are r(m) = Zeta(2; m-1) - (m-1)/m^2, m >= 2, with the partial sums Zeta(2; n) = Sum_{k=1..n} 1/k^2. See the W. Lang link in A103345.
O.g.f. for the rationals r(m), m>=2: log(1-x) + polylog(2,x)/(1-x).

A164561 Triangle with elements A120070(m,n)/A120072(m,n) read by rows, m>=2, 1<=n

Original entry on oeis.org

1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 16, 27, 4, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 16, 1, 4, 1, 1, 1, 9, 1, 1, 9, 1, 1, 1, 16, 1, 4, 25, 16, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 27, 128, 1, 108, 1, 16, 9, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 16, 1, 4, 1, 16, 49, 4, 1, 16, 1, 4, 1, 1, 1, 27, 1
Offset: 2

Views

Author

Paul Curtz, Aug 16 2009

Keywords

Examples

			The triangle starts in row m=2 as
1;
1, 1;
1, 4, 1;
1, 1, 1, 1;
1, 16, 27, 4, 1;
1, 1, 1, 1, 1, 1;
1, 4, 1, 16, 1, 4, 1;
		

Programs

Extensions

Edited and extended by R. J. Mathar, Aug 19 2009

A005563 a(n) = n*(n+2) = (n+1)^2 - 1.

Original entry on oeis.org

0, 3, 8, 15, 24, 35, 48, 63, 80, 99, 120, 143, 168, 195, 224, 255, 288, 323, 360, 399, 440, 483, 528, 575, 624, 675, 728, 783, 840, 899, 960, 1023, 1088, 1155, 1224, 1295, 1368, 1443, 1520, 1599, 1680, 1763, 1848, 1935, 2024, 2115, 2208, 2303, 2400, 2499, 2600
Offset: 0

Views

Author

Keywords

Comments

Erdős conjectured that n^2 - 1 = k! has a solution if and only if n is 5, 11 or 71 (when k is 4, 5 or 7).
Second-order linear recurrences y(m) = 2y(m-1) + a(n)*y(m-2), y(0) = y(1) = 1, have closed form solutions involving only powers of integers. - Len Smiley, Dec 08 2001
Number of edges in the join of two cycle graphs, both of order n, C_n * C_n. - Roberto E. Martinez II, Jan 07 2002
Let k be a positive integer, M_n be the n X n matrix m_(i,j) = k^abs(i-j) then det(M_n) = (-1)^(n-1)*a(k-1)^(n-1). - Benoit Cloitre, May 28 2002
Also numbers k such that 4*k + 4 is a square. - Cino Hilliard, Dec 18 2003
For each term k, the function sqrt(x^2 + 1), starting with 1, produces an integer after k iterations. - Gerald McGarvey, Aug 19 2004
a(n) mod 3 = 0 if and only if n mod 3 > 0: a(A008585(n)) = 2; a(A001651(n)) = 0; a(n) mod 3 = 2*(1-A079978(n)). - Reinhard Zumkeller, Oct 16 2006
a(n) is the number of divisors of a(n+1) that are not greater than n. - Reinhard Zumkeller, Apr 09 2007
Nonnegative X values of solutions to the equation X^3 + X^2 = Y^2. To find Y values: b(n) = n(n+1)(n+2). - Mohamed Bouhamida, Nov 06 2007
Sequence allows us to find X values of the equation: X + (X + 1)^2 + (X + 2)^3 = Y^2. To prove that X = n^2 + 2n: Y^2 = X + (X + 1)^2 + (X + 2)^3 = X^3 + 7*X^2 + 15X + 9 = (X + 1)(X^2 + 6X + 9) = (X + 1)*(X + 3)^2 it means: (X + 1) must be a perfect square, so X = k^2 - 1 with k>=1. we can put: k = n + 1, which gives: X = n^2 + 2n and Y = (n + 1)(n^2 + 2n + 3). - Mohamed Bouhamida, Nov 12 2007
From R. K. Guy, Feb 01 2008: (Start)
Toads and Frogs puzzle:
This is also the number of moves that it takes n frogs to swap places with n toads on a strip of 2n + 1 squares (or positions, or lily pads) where a move is a single slide or jump, illustrated for n = 2, a(n) = 8 by
T T - F F
T - T F F
T F T - F
T F T F -
T F - F T
- F T F T
F - T F T
F F T - T
F F - T T
I was alerted to this by the Holton article, but on consulting Singmaster's sources, I find that the puzzle goes back at least to 1867.
Probably the first to publish the number of moves for n of each animal was Edouard Lucas in 1883. (End)
a(n+1) = terms of rank 0, 1, 3, 6, 10 = A000217 of A120072 (3, 8, 5, 15). - Paul Curtz, Oct 28 2008
Row 3 of array A163280, n >= 1. - Omar E. Pol, Aug 08 2009
Final digit belongs to a periodic sequence: 0, 3, 8, 5, 4, 5, 8, 3, 0, 9. - Mohamed Bouhamida, Sep 04 2009 [Comment edited by N. J. A. Sloane, Sep 24 2009]
Let f(x) be a polynomial in x. Then f(x + n*f(x)) is congruent to 0 (mod f(x)); here n belongs to N. There is nothing interesting in the quotients f(x + n*f(x))/f(x) when x belongs to Z. However, when x is irrational these quotients consist of two parts, a) rational integers and b) integer multiples of x. The present sequence represents the non-integer part when the polynomial is x^2 + x + 1 and x = sqrt(2), f(x+n*f(x))/f(x) = A056108(n) + a(n)*sqrt(2). - A.K. Devaraj, Sep 18 2009
For n >= 1, a(n) is the number for which 1/a(n) = 0.0101... (A000035) in base (n+1). - Rick L. Shepherd, Sep 27 2009
For n > 0, continued fraction [n, 1, n] = (n+1)/a(n); e.g., [6, 1, 6] = 7/48. - Gary W. Adamson, Jul 15 2010
Starting (3, 8, 15, ...) = binomial transform of [3, 5, 2, 0, 0, 0, ...]; e.g., a(3) = 15 = (1*3 + 2*5 +1*2) = (3 + 10 + 2). - Gary W. Adamson, Jul 30 2010
a(n) is essentially the case 0 of the polygonal numbers. The polygonal numbers are defined as P_k(n) = Sum_{i=1..n} ((k-2)*i-(k-3)). Thus P_0(n) = 2*n-n^2 and a(n) = -P_0(n+2). See also A067998 and for the case k=1 A080956. - Peter Luschny, Jul 08 2011
a(n) is the maximal determinant of a 2 X 2 matrix with integer elements from {1, ..., n+1}, so the maximum determinant of a 2x2 matrix with integer elements from {1, ..., 5} = 5^2 - 1 = a(4) = 24. - Aldo González Lorenzo, Oct 12 2011
Using four consecutive triangular numbers t1, t2, t3 and t4, plot the points (0, 0), (t1, t2), and (t3, t4) to create a triangle. Twice the area of this triangle are the numbers in this sequence beginning with n = 1 to give 8. - J. M. Bergot, May 03 2012
Given a particle with spin S = n/2 (always a half-integer value), the quantum-mechanical expectation value of the square of the magnitude of its spin vector evaluates to = S(S+1) = n(n+2)/4, i.e., one quarter of a(n) with n = 2S. This plays an important role in the theory of magnetism and magnetic resonance. - Stanislav Sykora, May 26 2012
Twice the harmonic mean [H(x, y) = (2*x*y)/(x + y)] of consecutive triangular numbers A000217(n) and A000217(n+1). - Raphie Frank, Sep 28 2012
Number m such that floor(sqrt(m)) = floor(m/floor(sqrt(m))) - 2 for m > 0. - Takumi Sato, Oct 10 2012
The solutions of equation 1/(i - sqrt(j)) = i + sqrt(j), when i = (n+1), j = a(n). For n = 1, 2 + sqrt(3) = 3.732050.. = A019973. For n = 2, 3 + sqrt(8) = 5.828427... = A156035. - Kival Ngaokrajang, Sep 07 2013
The integers in the closed form solution of a(n) = 2*a(n-1) + a(m-2)*a(n-2), n >= 2, a(0) = 0, a(1) = 1 mentioned by Len Smiley, Dec 08 2001, are m and -m + 2 where m >= 3 is a positive integer. - Felix P. Muga II, Mar 18 2014
Let m >= 3 be a positive integer. If a(n) = 2*a(n-1) + a(m-2) * a(n-2), n >= 2, a(0) = 0, a(1) = 1, then lim_{n->oo} a(n+1)/a(n) = m. - Felix P. Muga II, Mar 18 2014
For n >= 4 the Szeged index of the wheel graph W_n (with n + 1 vertices). In the Sarma et al. reference, Theorem 2.7 is incorrect. - Emeric Deutsch, Aug 07 2014
If P_{k}(n) is the n-th k-gonal number, then a(n) = t*P_{s}(n+2) - s*P_{t}(n+2) for s=t+1. - Bruno Berselli, Sep 04 2014
For n >= 1, a(n) is the dimension of the simple Lie algebra A_n. - Wolfdieter Lang, Oct 21 2015
Finding all positive integers (n, k) such that n^2 - 1 = k! is known as Brocard's problem, (see A085692). - David Covert, Jan 15 2016
For n > 0, a(n) mod (n+1) = a(n) / (n+1) = n. - Torlach Rush, Apr 04 2016
Conjecture: When using the Sieve of Eratosthenes and sieving (n+1..a(n)), with divisors (1..n) and n>0, there will be no more than a(n-1) composite numbers. - Fred Daniel Kline, Apr 08 2016
a(n) mod 8 is periodic with period 4 repeating (0,3,0,7), that is a(n) mod 8 = 5/2 - (5/2) cos(n*Pi) - sin(n*Pi/2) + sin(3*n*Pi/2). - Andres Cicuttin, Jun 02 2016
Also for n > 0, a(n) is the number of times that n-1 occurs among the first (n+1)! terms of A055881. - R. J. Cano, Dec 21 2016
The second diagonal of composites (the only prime is number 3) from the right on the Klauber triangle (see Kival Ngaokrajang link), which is formed by taking the positive integers and taking the first 1, the next 3, the following 5, and so on, each centered below the last. - Charles Kusniec, Jul 03 2017
Also the number of independent vertex sets in the n-barbell graph. - Eric W. Weisstein, Aug 16 2017
Interleaving of A000466 and A033996. - Bruce J. Nicholson, Nov 08 2019
a(n) is the number of degrees of freedom in a triangular cell for a Raviart-Thomas or Nédélec first kind finite element space of order n. - Matthew Scroggs, Apr 22 2020
From Muge Olucoglu, Jan 19 2021: (Start)
For n > 1, a(n-2) is the maximum number of elements in the second stage of the Quine-McCluskey algorithm whose minterms are not covered by the functions of n bits. At n=3, we have a(3-2) = a(1) = 1*(1+2) = 3 and f(A,B,C) = sigma(0,1,2,5,6,7).
.
0 1 2 5 6 7
+---------------
*(0,1)| X X
(0,2)| X X
(1,5)| X X
*(2,6)| X X
*(5,7)| X X
(6,7)| X X
.
*: represents the elements that are covered. (End)
1/a(n) is the ratio of the sum of the first k odd numbers and the sum of the next n*k odd numbers. - Melvin Peralta, Jul 15 2021
For n >= 1, the continued fraction expansion of sqrt(a(n)) is [n; {1, 2n}]. - Magus K. Chu, Sep 09 2022
Number of diagonals parallel to an edge in a regular (2*n+4)-gon (cf. A367204). - Paolo Xausa, Nov 21 2023
For n >= 1, also the number of minimum cyclic edge cuts in the (n+2)-trapezohedron graph. - Eric W. Weisstein, Nov 21 2024
For n >= 1, a(n) is the sum of the interior angles of a polygon with n+2 sides, in radians, multiplied by (n+2)/Pi. - Stuart E Anderson, Aug 06 2025

Examples

			G.f. = 3*x + 8*x^2 + 15*x^3 + 24*x^4 + 35*x^5 + 48*x^6 + 63*x^7 + 80*x^8 + ...
		

References

  • E. R. Berlekamp, J. H. Conway and R. K. Guy, Winning Ways, Academic Press, NY, 2 vols., 1982, see index under Toads and Frogs Puzzle.
  • Martin Gardner, Perplexing Puzzles and Tantalizing Teasers, p. 21 (for "The Dime and Penny Switcheroo").
  • R. K. Guy, Unsolved Problems in Theory of Numbers, Section D25.
  • Derek Holton, Math in School, 37 #1 (Jan 2008) 20-22.
  • Edouard Lucas, Récréations Mathématiques, Gauthier-Villars, Vol. 2 (1883) 141-143.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

Formula

G.f.: x*(3-x)/(1-x)^3. - Simon Plouffe in his 1992 dissertation
a(n) = A000290(n+1) - 1.
A002378(a(n)) = A002378(n)*A002378(n+1); e.g., A002378(15)=240=12*20. - Charlie Marion, Dec 29 2003
a(n) = A067725(n)/3. - Zerinvary Lajos, Mar 06 2007
a(n) = Sum_{k=1..n} A144396(k). - Zerinvary Lajos, May 11 2007
a(n) = A134582(n+1)/4. - Zerinvary Lajos, Feb 01 2008
A143053(a(n)) = A000290(n+1), for n > 0. - Reinhard Zumkeller, Jul 20 2008
a(n) = Real((n+1+i)^2). - Gerald Hillier, Oct 12 2008
A053186(a(n)) = 2*n. - Reinhard Zumkeller, May 20 2009
a(n) = (n! + (n+1)!)/(n-1)!, n > 0. - Gary Detlefs, Aug 10 2009
a(n) = floor(n^5/(n^3+1)) with offset 1 (a(1)=0). - Gary Detlefs, Feb 11 2010
a(n) = a(n-1) + 2*n + 1 (with a(0)=0). - Vincenzo Librandi, Nov 18 2010
Sum_{n>=1} 1/a(n) = 3/4. - Mohammad K. Azarian, Dec 29 2010
a(n) = 2/(Integral_{x=0..Pi/2} (sin(x))^(n-1)*(cos(x))^3), for n > 0. - Francesco Daddi, Aug 02 2011
a(n) = A002378(n) + floor(sqrt(A002378(n))); pronic number + its root. - Fred Daniel Kline, Sep 16 2011
a(n-1) = A008833(n) * A068310(n) for n > 1. - Reinhard Zumkeller, Nov 26 2011
G.f.: U(0) where U(k) = -1 + (k+1)^2/(1 - x/(x + (k+1)^2/U(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 19 2012
a(n) = 15*C(n+4,3)*C(n+4,5)/(C(n+4,2)*C(n+4,4)). - Gary Detlefs, Aug 05 2013
a(n) = (n+2)!/((n-1)! + n!), n > 0. - Ivan N. Ianakiev, Nov 11 2013
a(n) = 3*C(n+1,2) - C(n,2) for n >= 0. - Felix P. Muga II, Mar 11 2014
a(n) = (A016742(n+1) - 4)/4 for n >= 0. - Felix P. Muga II, Mar 11 2014
a(-2 - n) = a(n) for all n in Z. - Michael Somos, Aug 07 2014
A253607(a(n)) = 1. - Reinhard Zumkeller, Jan 05 2015
E.g.f.: x*(x + 3)*exp(x). - Ilya Gutkovskiy, Jun 03 2016
For n >= 1, a(n^2 + n - 2) = a(n-1) * a(n). - Miko Labalan, Oct 15 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = 1/4. - Amiram Eldar, Nov 04 2020
From Amiram Eldar, Feb 17 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = 2.
Product_{n>=1} (1 - 1/a(n)) = -sqrt(2)*sin(sqrt(2)*Pi)/Pi. (End)
a(n) = A000290(n+2) - n*2. See Bounded Squares illustration. - Leo Tavares, Oct 05 2021
From Leo Tavares, Oct 10 2021: (Start)
a(n) = A008585(n) + 2*A000217(n-1). See Trapezoids illustration.
2*A005563 = A054000(n+1). See Trapagons illustration.
a(n) = 2*A000217(n) + n. (End)
a(n) = (n+2)!!/(n-2)!! for n > 1. - Jacob Szlachetka, Jan 02 2022

Extensions

Partially edited by Joerg Arndt, Mar 11 2010
More terms from N. J. A. Sloane, Aug 01 2010

A120070 Triangle of numbers used to compute the frequencies of the spectral lines of the hydrogen atom.

Original entry on oeis.org

3, 8, 5, 15, 12, 7, 24, 21, 16, 9, 35, 32, 27, 20, 11, 48, 45, 40, 33, 24, 13, 63, 60, 55, 48, 39, 28, 15, 80, 77, 72, 65, 56, 45, 32, 17, 99, 96, 91, 84, 75, 64, 51, 36, 19, 120, 117, 112, 105, 96, 85, 72, 57, 40, 21
Offset: 2

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

The rationals r(m,n):=a(m,n)/(m^2*n^2), for m-1 >= n, else 0, are used to compute the frequencies of the spectral lines of the H-atom according to quantum theory: nu(m,n) = r(m,n)*c*R' with c*R'=3.287*10^15 s^(-1) an approximation for the Rydberg frequency. R' indicates, that the correction factor 1/(1+m_e/m_p), approximately 0.9995, with the masses for the electron and proton, has been used for the Rydberg constant R_infinity. c:=299792458 m/s is, per definition, the velocity of light in vacuo (see A003678).
In order to compute the wave length of the spectral lines approximately one uses the reciprocal rationals: lambda(m,n):= c/nu(m,n) = (1/r(m,n))*91.1961 nm. 1 nm = 10^{-9} m. For the corresponding energies one uses approximately E(m,n) = r(m,n)*13.599 eV (electron Volts).
The author was inspired by Dewdney's book to compile this table and related ones.
For the approximate frequencies, energies and wavelengths of the first members of the Lyman (n=1, m>=2), Balmer (n=2, m>=3), Paschen (n=3, m>=4), Brackett (n=4, m>=5) and Pfund (n=5, m>=6) series see the W. Lang link under A120072.
Frenicle wrote this as a(n+1) = A140978(n) - A133819(n-1). - Paul Curtz, Aug 19 2008
This triangle also has an interpretation related to particle spin. For proper offset such that T(0,0) = 3, then, where h-bar = h/(2*Pi) = A003676/A019692 (= The Dirac constant, also known as Planck's reduced constant) and Spin(n/2) = h-bar/2*sqrt(n(n+2)), it follows that: h-bar/2*sqrt(T(r,k)) = h-bar/2*sqrt(T(r,0) - T(k-1,0)) = sqrt((Spin((r+1)/2))^2 - (Spin(k/2))^2). For example, for r = k = 4, then h-bar/2*sqrt(11) = h-bar/2*sqrt(T(4,4)) = h-bar/2*sqrt(T(4,0) - T(3,0)) = sqrt(h-bar^2/4*T(4,0) - h-bar^2/4*T(3,0)) = sqrt(h-bar^2/4*35 - h-bar^2/4*24) = sqrt((Spin((4+1)/2))^2 - (Spin(4/2))^2); 35 = 5*(5+2) & 24 = 4*(4+2). - Raphie Frank, Dec 30 2012

Examples

			Triangle begins
  [ 3];
  [ 8, 5];
  [15, 12,  7];
  [24, 21, 16,  9];
  [35, 32, 27, 20, 11];
  ...
		

References

  • A. K. Dewdney, Reise in das Innere der Mathematik, Birkhäuser, Basel, 2000, pp. 148-154; engl.: A Mathematical Mystery Tour, John Wiley & Sons, N.Y., 1999.

Crossrefs

Row sums give A016061(n-1), n>=2.
Cf. A120072/A120073 numerator and denominator tables for rationals r(m, n).

Programs

  • Mathematica
    ColumnForm[Table[n^2 - k^2, {n, 2, 13}, {k, n - 1}], Center] (* Alonso del Arte, Oct 26 2011 *)
  • PARI
    nmax=400;a=vector(1+nmax*(nmax-1)\2);idx=1;for(n=2,nmax,for(k=1,n-1,a[idx]=n*n-k*k;idx++)) \\ Stanislav Sykora, Feb 17 2014
    
  • PARI
    T(n,k)=n^2-k^2;
    for (n=1,10,for(k=1,n-1, print1(T(n,k),", ")));
    \\ Joerg Arndt, Feb 23 2014

Formula

a(m,n) = m^2 - n^2 for m-1 >= n, otherwise 0.
G.f. for column n=1,2,...: x^(n+1)*((2*n+1)- (2*n-1)*x)/(1-x)^3.
G.f. for rationals r(m,n), n=1,2,...,10 see W. Lang link.
T(r,k) = T(r,0) - T(k-1,0), T(0,0) = 3. - Raphie Frank, Dec 27 2012

A120073 Denominator triangle for hydrogen spectrum rationals.

Original entry on oeis.org

4, 9, 36, 16, 16, 144, 25, 100, 225, 400, 36, 9, 12, 144, 900, 49, 196, 441, 784, 1225, 1764, 64, 64, 576, 64, 1600, 576, 3136, 81, 324, 81, 1296, 2025, 324, 3969, 5184, 100, 25, 900, 400, 100, 225, 4900, 1600, 8100, 121, 484, 1089, 1936, 3025, 4356, 5929, 7744, 9801, 12100
Offset: 2

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

The corresponding numerator triangle is A120072.
See A120072 and A120070 for more details.

Examples

			For the rational triangle see W. Lang link.
Denominator triangle begins as:
    4;
    9,  36;
   16,  16, 144;
   25, 100, 225,  400;
   36,   9,  12,  144,  900;
   49, 196, 441,  784, 1225, 1764;
   64,  64, 576,   64, 1600,  576, 3136;
   81, 324,  81, 1296, 2025,  324, 3969, 5184;
  100,  25, 900,  400,  100,  225, 4900, 1600, 8100;
		

Crossrefs

Programs

  • Magma
    [Denominator(1/k^2 - 1/n^2): k in [1..n-1], n in [2..18]]; // G. C. Greubel, Apr 24 2023
    
  • Mathematica
    Table[(1/n^2 - 1/m^2)//Denominator, {m,2,15}, {n,m-1}]//Flatten (* Jean-François Alcover, Sep 16 2013 *)
  • SageMath
    def A120073(n,k): return denominator(1/k^2 - 1/n^2)
    flatten([[A120073(n,k) for k in range(1,n)] for n in range(2,19)]) # G. C. Greubel, Apr 24 2023

Formula

a(m,n) = denominator(r(m,n)) with r(m,n) = 1/n^2 - 1/m^2, m>=2, n=1..m-1.

A119937 Triangle of numbers related to the spectrum of the hydrogen (H) atom.

Original entry on oeis.org

3, 32, 5, 135, 27, 7, 3456, 756, 256, 81, 3500, 800, 300, 125, 44, 172800, 40500, 16000, 7425, 3456, 1300, 694575, 165375, 67375, 33075, 17199, 8575, 3375, 6272000, 1509200, 627200, 318500, 175616, 98000, 51200
Offset: 2

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

The rational number triangle r(m,n):=A120072(m,n)/A120073(m,n), used to compute the spectral series of the hydrogen atom, is mapped to this nonnegative number triangle by multiplying the least common multiples (LCM) for each row m.

Examples

			[3]; [32,5]; [135,27,7]; [3456,756,256,81]; [3500,800,300,125,44]; ...
		

Crossrefs

The LCM sequence which has been used here is [4, 36, 144, 3600, 3600, 176400, 705600, 6350400, 6350400, 768398400, ...] = A051418(m) = (A003418(m))^2 = (2*A025555(m-1))^2, m >= 2.
The row sums give A119938.

Formula

a(m,n) = r(m,n)*lcm_{k=1..m-1} seq(r(m,k)) with r(m,n) = 1/n^2 - 1/m^2 = A120072(m,n)/A120073(m,n), m >= 2, n = 1..m-1.

A172157 Triangle T(n,m) = numerator of 1/n^2 - 1/m^2, read by rows, with T(n,0) = -1.

Original entry on oeis.org

-1, -1, -3, -1, -8, -5, -1, -15, -3, -7, -1, -24, -21, -16, -9, -1, -35, -2, -1, -5, -11, -1, -48, -45, -40, -33, -24, -13, -1, -63, -15, -55, -3, -39, -7, -15, -1, -80, -77, -8, -65, -56, -5, -32, -17, -1, -99, -6, -91, -21, -3, -4, -51, -9, -19, -1, -120, -117, -112
Offset: 1

Views

Author

Paul Curtz, Jan 27 2010

Keywords

Comments

The triangle obtained by negating the values of the triangle A120072 and adding a row T(n,0) = -1.

Examples

			The full array of numerators starts in row n=1 with columns m>=0 as:
-1...0...3...8..15..24..35..48..63..80..99. A005563
-1..-3...0...5...3..21...2..45..15..77...6. A061037, A070262
-1..-8..-5...0...7..16...1..40..55...8..91. A061039
-1.-15..-3..-7...0...9...5..33...3..65..21. A061041
-1.-24.-21.-16..-9...0..11..24..39..56...3. A061043
-1.-35..-2..-1..-5.-11...0..13...7...5...4. A061045
-1.-48.-45.-40.-33.-24.-13...0..15..32..51. A061047
-1.-63.-15.-55..-3.-39..-7.-15...0..17...9. A061049
The triangle is the portion below the main diagonal, left from the zeros, 0<=m<n.
		

Crossrefs

Programs

  • Mathematica
    T[n_, 0] := -1; T[n_, k_] := 1/n^2 - 1/k^2; Table[Numerator[T[n, k]], {n, 1, 100}, {k, 0, n - 1}] // Flatten (* G. C. Greubel, Sep 19 2018 *)
Showing 1-10 of 18 results. Next