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

A001108 a(n)-th triangular number is a square: a(n+1) = 6*a(n) - a(n-1) + 2, with a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 8, 49, 288, 1681, 9800, 57121, 332928, 1940449, 11309768, 65918161, 384199200, 2239277041, 13051463048, 76069501249, 443365544448, 2584123765441, 15061377048200, 87784138523761, 511643454094368, 2982076586042449, 17380816062160328, 101302819786919521
Offset: 0

Views

Author

Keywords

Comments

b(0)=0, c(0)=1, b(i+1)=b(i)+c(i), c(i+1)=b(i+1)+b(i); then a(i) (the number in the sequence) is 2b(i)^2 if i is even, c(i)^2 if i is odd and b(n)=A000129(n) and c(n)=A001333(n). - Darin Stephenson (stephenson(AT)cs.hope.edu) and Alan Koch
For n > 1 gives solutions to A007913(2x) = A007913(x+1). - Benoit Cloitre, Apr 07 2002
If (X,X+1,Z) is a Pythagorean triple, then Z-X-1 and Z+X are in the sequence.
For n >= 2, a(n) gives exactly the positive integers m such that 1,2,...,m has a perfect median. The sequence of associated perfect medians is A001109. Let a_1,...,a_m be an (ordered) sequence of real numbers, then a term a_k is a perfect median if Sum_{j=1..k-1} a_j = Sum_{j=k+1..m} a_j. See Puzzle 1 in MSRI Emissary, Fall 2005. - Asher Auel, Jan 12 2006
This is the r=8 member of the r-family of sequences S_r(n) defined in A092184 where more information can be found.
Also, 1^3 + 2^3 + 3^3 + ... + a(n)^3 = k(n)^4 where k(n) is A001109. - Anton Vrba (antonvrba(AT)yahoo.com), Nov 18 2006
If T_x = y^2 is a triangular number which is also a square, the least number which is both triangular and square and greater than T_x is T_(3*x + 4*y + 1) = (2*x + 3*y + 1)^2 (W. Sierpiński 1961). - Richard Choulet, Apr 28 2009
If (a,b) is a solution of the Diophantine equation 0 + 1 + 2 + ... + x = y^2, then a or (a+1) is a perfect square. If (a,b) is a solution of the Diophantine equation 0 + 1 + 2 + ... + x = y^2, then a or a/8 is a perfect square. If (a,b) and (c,d) are two consecutive solutions of the Diophantine equation 0 + 1 + 2 + ... + x = y^2 with a < c, then a+b = c-d and ((d+b)^2, d^2-b^2) is a solution, too. If (a,b), (c,d) and (e,f) are three consecutive solutions of the Diophantine equation 0 + 1 + 2 + ... + x = y^2 with a < c < e, then (8*d^2, d*(f-b)) is a solution, too. - Mohamed Bouhamida, Aug 29 2009
If (p,q) and (r,s) are two consecutive solutions of the Diophantine equation 0 + 1 + 2 + ... + x = y^2 with p < r, then r = 3p + 4q + 1 and s = 2p + 3q + 1. - Mohamed Bouhamida, Sep 02 2009
Also numbers k such that (ceiling(sqrt(k*(k+1)/2)))^2 - k*(k+1)/2 = 0. - Ctibor O. Zizka, Nov 10 2009
From Lekraj Beedassy, Mar 04 2011: (Start)
Let x=a(n) be the index of the associated triangular number T_x=1+2+3+...+x and y=A001109(n) be the base of the associated perfect square S_y=y^2. Now using the identity S_y = T_y + T_{y-1}, the defining T_x = S_y may be rewritten as T_y = T_x - T_{y-1}, or 1+2+3+...+y = y+(y+1)+...+x. This solves the Strand Magazine House Number problem mentioned in A001109 in references from Poo-Sung Park and John C. Butcher. In a variant of the problem, solving the equation 1+3+5+...+(2*x+1) = (2*x+1)+(2*x+3)+...+(2*y-1) implies S_(x+1) = S_y - S_x, i.e., with (x,x+1,y) forming a Pythagorean triple, the solutions are given by pairs of x=A001652(n), y=A001653(n). (End)
If P = 8*n +- 1 is a prime, then P divides a((P-1)/2); e.g., 7 divides a(3) and 41 divides a(20). Also, if P = 8*n +- 3 is prime, then 4*P divides (a((P-1)/2) + a((P+1)/2) + 3). - Kenneth J Ramsey, Mar 05 2012
Starting at a(2), a(n) gives all the dimensions of Euclidean k-space in which the ratio of outer to inner Soddy hyperspheres' radii for k+1 identical kissing hyperspheres is rational. The formula for this ratio is (1+3k+2*sqrt(2k*(k+1)))/(k-1) where k is the dimension. So for a(3) = 49, the ratio is 6 in the 49th dimension. See comment for A010502. - Frank M Jackson, Feb 09 2013
Conjecture: For n>1 a(n) is the index of the first occurrence of -n in sequence A123737. - Vaclav Kotesovec, Jun 02 2015
For n=2*k, k>0, a(n) is divisible by 8 (deficient), so since all proper divisors of deficient numbers are deficient, then a(n) is deficient. For n=2*k+1, k>0, a(n) is odd. If a(n) is a prime number, it is deficient; otherwise a(n) has one or two distinct prime factors and is therefore deficient again. sigma(a(5)) = 1723 < 3362 = 2*a(5). In either case, a(n) is deficient. - Muniru A Asiru, Apr 14 2016
The squares of NSW numbers (A008843) interleaved with twice squares from A084703, where A008843(n) = A002315(n)^2 and A084703(n) = A001542(n)^2. Conjecture: Also numbers n such that sigma(n) = A000203(n) and sigma(n-th triangular number) = A074285(n) are both odd numbers. - Jaroslav Krizek, Aug 05 2016
For n > 0, numbers for which the number of odd divisors of both n and of n + 1 is odd. - Gionata Neri, Apr 30 2018
a(n) will be solutions to some (A000217(k) + A000217(k+1))/2. - Art Baker, Jul 16 2019
For n >= 2, a(n) is the base for which A058331(A001109(n)) is a length-3 repunit. Example: for n=2, A001109(2)=6 and A058331(6)=73 and 73 in base a(2)=8 is 111. See Grantham and Graves. - Michel Marcus, Sep 11 2020

Examples

			a(1) = ((3 + 2*sqrt(2)) + (3 - 2*sqrt(2)) - 2) / 4 = (3 + 3 - 2) / 4 = 4 / 4 = 1;
a(2) = ((3 + 2*sqrt(2))^2 + (3 - 2*sqrt(2))^2 - 2) / 4 = (9 + 4*sqrt(2) + 8 + 9 - 4*sqrt(2) + 8 - 2) / 4 = (18 + 16 - 2) / 4 = (34 - 2) / 4 = 32 / 4 = 8, etc.
		

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 193.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 204.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 10.
  • M. S. Klamkin, "International Mathematical Olympiads 1978-1985," (Supplementary problem N.T.6)
  • W. Sierpiński, Pythagorean triangles, Dover Publications, Inc., Mineola, NY, 2003, pp. 21-22 MR2002669
  • 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).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 257-258.

Crossrefs

Partial sums of A002315. A000129, A005319.
a(n) = A115598(n), n > 0. - Hermann Stamm-Wilbrandt, Jul 27 2014

Programs

  • Haskell
    a001108 n = a001108_list !! n
    a001108_list = 0 : 1 : map (+ 2)
       (zipWith (-) (map (* 6) (tail a001108_list)) a001108_list)
    -- Reinhard Zumkeller, Jan 10 2012
    
  • Magma
    m:=30; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(x*(1+x)/((1-x)*(1-6*x+x^2)))); // G. C. Greubel, Jul 15 2018
  • Maple
    A001108:=-(1+z)/(z-1)/(z**2-6*z+1); # Simon Plouffe in his 1992 dissertation, without the leading 0
  • Mathematica
    Table[(1/2)(-1 + Sqrt[1 + Expand[8(((3 + 2Sqrt[2])^n - (3 - 2Sqrt[2])^n)/(4Sqrt[2]))^2]]), {n, 0, 100}] (* Artur Jasinski, Dec 10 2006 *)
    Transpose[NestList[{#[[2]],#[[3]],6#[[3]]-#[[2]]+2}&,{0,1,8},20]][[1]] (* Harvey P. Dale, Sep 04 2011 *)
    LinearRecurrence[{7, -7, 1}, {0, 1, 8}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 12 2012 *)
  • PARI
    a(n)=(real((3+quadgen(32))^n)-1)/2
    
  • PARI
    a(n)=(subst(poltchebi(abs(n)),x,3)-1)/2
    
  • PARI
    a(n)=if(n<0,a(-n),(polsym(1-6*x+x^2,n)[n+1]-2)/4)
    
  • PARI
    x='x+O('x^99); concat(0, Vec(x*(1+x)/((1-x)*(1-6*x+x^2)))) \\ Altug Alkan, May 01 2018
    

Formula

a(0) = 0, a(n+1) = 3*a(n) + 1 + 2*sqrt(2*a(n)*(a(n)+1)). - Jim Nastos, Jun 18 2002
a(n) = floor( (1/4) * (3+2*sqrt(2))^n ). - Benoit Cloitre, Sep 04 2002
a(n) = A001653(k)*A001653(k+n) - A001652(k)*A001652(k+n) - A046090(k)*A046090(k+n). - Charlie Marion, Jul 01 2003
a(n) = A001652(n-1) + A001653(n-1) = A001653(n) - A046090(n) = (A001541(n)-1)/2 = a(-n). - Michael Somos, Mar 03 2004
a(n) = 7*a(n-1) - 7*a(n-2) + a(n-3). - Antonio Alberto Olivares, Oct 23 2003
a(n) = Sum_{r=1..n} 2^(r-1)*binomial(2n, 2r). - Lekraj Beedassy, Aug 21 2004
If n > 1, then both A000203(n) and A000203(n+1) are odd numbers: n is either a square or twice a square. - Labos Elemer, Aug 23 2004
a(n) = (T(n, 3)-1)/2 with Chebyshev's polynomials of the first kind evaluated at x=3: T(n, 3) = A001541(n). - Wolfdieter Lang, Oct 18 2004
G.f.: x*(1+x)/((1-x)*(1-6*x+x^2)). Binet form: a(n) = ((3+2*sqrt(2))^n + (3-2*sqrt(2))^n - 2)/4. - Bruce Corrigan (scentman(AT)myfamily.com), Oct 26 2002
a(n) = floor(sqrt(2*A001110(n))) = floor(A001109(n)*sqrt(2)) = 2*(A000129(n)^2) - (n mod 2) = A001333(n)^2 - 1 + (n mod 2). - Henry Bottomley, Apr 19 2000, corrected by Eric Rowland, Jun 23 2017
A072221(n) = 3*a(n) + 1. - David Scheers, Dec 25 2006
A028982(a(n)) + 1 = A028982(a(n) + 1). - Juri-Stepan Gerasimov, Mar 28 2011
a(n+1)^2 + a(n)^2 + 1 = 6*a(n+1)*a(n) + 2*a(n+1) + 2*a(n). - Charlie Marion, Sep 28 2011
a(n) = 2*A001653(m)*A053141(n-m-1) + A002315(m)*A046090(n-m-1) + a(m) with m < n; otherwise, a(n) = 2*A001653(m)*A053141(m-n) - A002315(m)*A001652(m-n) + a(m). See Link to Generalized Proof re Square Triangular Numbers. - Kenneth J Ramsey, Oct 13 2011
a(n) = A048739(2n-2), n > 0. - Richard R. Forberg, Aug 31 2013
From Peter Bala, Jan 28 2014: (Start)
A divisibility sequence: that is, a(n) divides a(n*m) for all n and m. Case P1 = 8, P2 = 12, Q = 1 of the 3-parameter family of linear divisibility sequences found by Williams and Guy.
a(2*n+1) = A002315(n)^2 = Sum_{k = 0..4*n + 1} Pell(n), where Pell(n) = A000129(n).
a(2*n) = (1/2)*A005319(n)^2 = 8*A001109(n)^2.
(2,1) entry of the 2 X 2 matrix T(n,M), where M = [0, -3; 1, 4] and T(n,x) is the Chebyshev polynomial of the first kind. (End)
E.g.f.: exp(x)*(exp(2*x)*cosh(2*sqrt(2)*x) - 1)/2. - Stefano Spezia, Oct 25 2024

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Apr 19 2000
More terms from Lekraj Beedassy, Aug 21 2004

A346868 Sum of divisors of the numbers with no middle divisors.

Original entry on oeis.org

4, 6, 8, 18, 12, 14, 24, 18, 20, 32, 36, 24, 42, 40, 30, 32, 48, 54, 38, 60, 56, 42, 44, 84, 72, 48, 72, 98, 54, 72, 80, 90, 60, 62, 96, 84, 68, 126, 96, 72, 74, 114, 124, 140, 168, 80, 126, 84, 108, 132, 120, 90, 168, 128, 144, 120, 98, 102, 216, 104, 192, 162, 108, 110
Offset: 1

Views

Author

Omar E. Pol, Aug 18 2021

Keywords

Comments

The characteristic shape of the symmetric representation of a(n) consists in that in the main diagonal of the diagram the width is equal to zero.
So knowing this characteristic shape we can know if a number has middle divisors (or not) just by looking at the diagram, even ignoring the concept of middle divisors.
Therefore we can see a geometric pattern of the distribution of the numbers with no middle divisors in the stepped pyramid described in A245092.
For the definition of "width" see A249351.
All terms are even numbers.

Examples

			a(4) = 18 because the sum of divisors of the fourth number with no middle divisors (i.e., 10) is 1 + 2 + 5 + 10 = 18.
On the other hand we can see that in the main diagonal of every diagram the width is equal to zero as shown below.
Illustration of initial terms:
m(n) = A071561(n).
.
   n   m(n) a(n)   Diagram
.                      _   _   _     _ _   _ _     _   _   _ _ _     _
                      | | | | | |   | | | | | |   | | | | | | | |   | |
                   _ _|_| | | | |   | | | | | |   | | | | | | | |   | |
   1    3    4    |_ _|  _|_| | |   | | | | | |   | | | | | | | |   | |
                   _ _ _|    _|_|   | | | | | |   | | | | | | | |   | |
   2    5    6    |_ _ _|  _|    _ _| | | | | |   | | | | | | | |   | |
                   _ _ _ _|     |  _ _|_| | | |   | | | | | | | |   | |
   3    7    8    |_ _ _ _|  _ _|_|    _ _|_| |   | | | | | | | |   | |
                            |  _|     |  _ _ _|   | | | | | | | |   | |
                   _ _ _ _ _| |      _|_|    _ _ _|_| | | | | | |   | |
   4   10   18    |_ _ _ _ _ _|  _ _|       |    _ _ _|_| | | | |   | |
   5   11   12    |_ _ _ _ _ _| |  _|      _|   |  _ _ _ _|_| | |   | |
                   _ _ _ _ _ _ _| |      _|  _ _| | |  _ _ _ _|_|   | |
   6   13   14    |_ _ _ _ _ _ _| |  _ _|  _|    _| | |    _ _ _ _ _| |
   7   14   24    |_ _ _ _ _ _ _ _| |     |     |  _|_|   |  _ _ _ _ _|
                                    |  _ _|  _ _|_|       | |
                   _ _ _ _ _ _ _ _ _| |  _ _|  _|        _|_|
   8   17   18    |_ _ _ _ _ _ _ _ _| | |_ _ _|         |
                   _ _ _ _ _ _ _ _ _ _| |  _ _|        _|
   9   19   20    |_ _ _ _ _ _ _ _ _ _| | |        _ _|
                   _ _ _ _ _ _ _ _ _ _ _| |  _ _ _|
  10   21   32    |_ _ _ _ _ _ _ _ _ _ _| | |  _ _|
  11   22   36    |_ _ _ _ _ _ _ _ _ _ _ _| | |
  12   23   24    |_ _ _ _ _ _ _ _ _ _ _ _| | |
                                            | |
                   _ _ _ _ _ _ _ _ _ _ _ _ _| |
  13   26   42    |_ _ _ _ _ _ _ _ _ _ _ _ _ _|
.
		

Crossrefs

Some sequences that gives sum of divisors: A000225 (of powers of 2), A008864 (of prime numbers), A065764 (of squares), A073255 (of composites), A074285 (of triangular numbers, also of generalized hexagonal numbers), A139256 (of perfect numbers), A175926 (of cubes), A224613 (of multiples of 6), A346865 (of hexagonal numbers), A346866 (of second hexagonal numbers), A346867 (of numbers with middle divisors).

Programs

  • Mathematica
    s[n_] := Module[{d = Divisors[n]}, If[AnyTrue[d, Sqrt[n/2] <= # < Sqrt[n*2] &], 0, Plus @@ d]]; Select[Array[s, 110], # > 0 &] (* Amiram Eldar, Aug 19 2021 *)
  • PARI
    is(n) = fordiv(n, d, if(sqrt(n/2) <= d && d < sqrt(2*n), return(0))); 1; \\ A071561 apply(sigma, select(is, [1..150])) \\ Michel Marcus, Aug 19 2021

Formula

a(n) = A000203(A071561(n)).

A346873 Triangle read by rows in which row n lists the row A000217(n) of A237591, n >= 1.

Original entry on oeis.org

1, 2, 1, 4, 1, 1, 6, 2, 1, 1, 8, 3, 2, 1, 1, 11, 4, 3, 1, 1, 1, 15, 5, 3, 2, 1, 1, 1, 19, 6, 4, 2, 2, 1, 1, 1, 23, 8, 5, 2, 2, 2, 1, 1, 1, 28, 10, 5, 3, 3, 2, 1, 1, 1, 1, 34, 11, 6, 4, 3, 2, 2, 1, 1, 1, 1, 40, 13, 7, 5, 3, 2, 2, 2, 1, 1, 1, 1, 46, 16, 8, 5, 4, 2, 3
Offset: 1

Views

Author

Omar E. Pol, Aug 06 2021

Keywords

Comments

The characteristic shape of the symmetric representation of sigma(A000217(n)) consists in that in the main diagonal of the diagram the smallest Dyck path has a valley and the largest Dyck path has a peak, or vice versa, the smallest Dyck path has a peak and the largest Dyck path has valley.
So knowing this characteristic shape we can know if a number is a triangular number (or not) just by looking at the diagram, even ignoring the concept of triangular number.
Therefore we can see a geometric pattern of the distribution of the triangular numbers in the stepped pyramid described in A245092.
T(n,k) is also the length of the k-th line segment of the largest Dyck path of the symmetric representation of sigma(A000217(n)), from the border to the center, hence the sum of the n-th row of triangle is equal to A000217(n).
T(n,k) is also the difference between the total number of partitions of all positive integers <= n-th triangular number into exactly k consecutive parts, and the total number of partitions of all positive integers <= n-th triangular number into exactly k + 1 consecutive parts.

Examples

			Triangle begins:
   1;
   2,  1;
   4,  1, 1;
   6,  2, 1, 1;
   8,  3, 2, 1, 1;
  11,  4, 3, 1, 1, 1;
  15,  5, 3, 2, 1, 1, 1;
  19,  6, 4, 2, 2, 1, 1, 1;
  23,  8, 5, 2, 2, 2, 1, 1, 1;
  28, 10, 5, 3, 3, 2, 1, 1, 1, 1;
  34, 11, 6, 4, 3, 2, 2, 1, 1, 1, 1;
  40, 13, 7, 5, 3, 2, 2, 2, 1, 1, 1, 1;
  46, 16, 8, 5, 4, 2, 3, 1, 2, 1, 1, 1, 1;
...
Illustration of initial terms:
Column T gives the triangular numbers (A000217).
Column S gives A074285, the sum of the divisors of the triangular numbers which equals the area (and the number of cells) of the associated diagram.
-------------------------------------------------------------------------
  n    T    S   Diagram
-------------------------------------------------------------------------
                 _   _     _       _         _           _             _
  1    1    1   |_| | |   | |     | |       | |         | |           | |
               1 _ _|_|   | |     | |       | |         | |           | |
  2    3    4   |_ _|  _ _| |     | |       | |         | |           | |
                  2  1|    _|     | |       | |         | |           | |
                 _ _ _|  _|    _ _| |       | |         | |           | |
  3    6   12   |_ _ _ _| 1   |  _ _|       | |         | |           | |
                    4    1 _ _|_|           | |         | |           | |
                          |  _|1       _ _ _|_|         | |           | |
                 _ _ _ _ _| | 1    _ _| |               | |           | |
  4   10   18   |_ _ _ _ _ _|2    |    _|               | |           | |
                      6          _|  _|          _ _ _ _|_|           | |
                                |_ _|1 1        | |                   | |
                                | 2            _| |                   | |
                 _ _ _ _ _ _ _ _|4            |  _|          _ _ _ _ _| |
  3   15   24   |_ _ _ _ _ _ _ _|          _ _|_|           |  _ _ _ _ _|
                        8              _ _|  _|1            | |
                                      |_ _ _|1 1         _ _| |
                                      |  3           _ _|  _ _|
                                      |4            |    _|
                 _ _ _ _ _ _ _ _ _ _ _|            _|  _|
  4   21   32   |_ _ _ _ _ _ _ _ _ _ _|      _ _ _|  _|1 1
                          11                |  _ _ _|2
                                            | |  3
                                            | |
                                            | |5
                 _ _ _ _ _ _ _ _ _ _ _ _ _ _| |
  5   28   56   |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _|
                              15
.
		

Crossrefs

Row sums give A000217, n >= 1.
Column 1 gives A039823.
For the characteristic shape of sigma(A000040(n)) see A346871.
For the characteristic shape of sigma(A000079(n)) see A346872.
For the visualization of Mersenne numbers A000225 see A346874.
For the characteristic shape of sigma(A000384(n)) see A346875.
For the characteristic shape of sigma(A000396(n)) see A346876.
For the characteristic shape of sigma(A008588(n)) see A224613.

Formula

T(n,k) = A237591(A000217(n),k). - Omar E. Pol, Feb 06 2023

Extensions

Name corrected by Omar E. Pol, Feb 06 2023

A256151 Triangular numbers n such that sigma(n) is a square number.

Original entry on oeis.org

1, 3, 66, 210, 820, 2346, 4278, 22578, 27966, 32131, 35511, 51681, 53956, 102378, 169653, 173755, 177906, 223446, 241860, 256686, 306153, 310866, 349866, 431056, 434778, 470935, 491536, 512578, 567645, 579426, 688551, 799480, 845650, 893116, 963966, 1031766, 1110795, 1200475, 1613706, 1719585
Offset: 1

Views

Author

Antonio Roldán, Mar 16 2015

Keywords

Comments

This sequence is the intersection of A000217 and A006532.
The corresponding triangular indices are in A116990. - Michel Marcus, Mar 17 2015

Examples

			3 is in the sequence because 3=2*3/2 is triangular, and sigma(3)=1+3=4=2^2 is square.
		

Crossrefs

Programs

  • Magma
    [n*(n+1) div 2: n in [1..2000] | IsSquare(SumOfDivisors(n*(n+1) div 2))]; // Vincenzo Librandi, Mar 17 2015
  • Mathematica
    Select[Accumulate[Range[0, 2000]], IntegerQ@Sqrt@DivisorSigma[1, #] &] (* Michael De Vlieger, Mar 17 2015 *)
  • PARI
    {for(i=1,2*10^3,n=i*(i+1)/2;if(issquare(sigma(n)),print1(n,", ")))}
    

A346865 Sum of divisors of the n-th hexagonal number.

Original entry on oeis.org

1, 12, 24, 56, 78, 144, 112, 360, 234, 360, 384, 672, 434, 960, 720, 992, 864, 1872, 760, 2352, 1344, 1584, 1872, 2880, 1767, 3024, 2160, 4032, 2400, 4320, 1984, 6552, 4032, 3672, 4608, 6552, 2812, 7440, 5376, 7200, 5082, 8064, 4752, 10080, 7020, 8064, 6144
Offset: 1

Views

Author

Omar E. Pol, Aug 17 2021

Keywords

Comments

The characteristic shape of the symmetric representation of a(n) consists in that in the main diagonal of the diagram the smallest Dyck path has a valley and the largest Dyck path has a peak.
So knowing this characteristic shape we can know if a number is an hexagonal number (or not) just by looking at the diagram, even ignoring the concept of hexagonal number.
Therefore we can see a geometric pattern of the distribution of the hexagonal numbers in the stepped pyramid described in A245092.

Examples

			a(3) = 24 because the sum of divisors of the third hexagonal number (i.e., 15) is 1 + 3 + 5 + 15 = 24.
On the other hand we can see that in the main diagonal of every diagram the smallest Dyck path has a valley and the largest Dyck path has a peak as shown below.
Illustration of initial terms:
-------------------------------------------------------------------------
  n  H(n)  a(n)  Diagram
-------------------------------------------------------------------------
                 _         _                 _                         _
  1    1    1   |_|       | |               | |                       | |
                          | |               | |                       | |
                       _ _| |               | |                       | |
                      |    _|               | |                       | |
                 _ _ _|  _|                 | |                       | |
  2    6   12   |_ _ _ _|                   | |                       | |
                                            | |                       | |
                                       _ _ _|_|                       | |
                                   _ _| |                             | |
                                  |    _|                             | |
                                 _|  _|                               | |
                                |_ _|                                 | |
                                |                                     | |
                 _ _ _ _ _ _ _ _|                            _ _ _ _ _| |
  3   15   24   |_ _ _ _ _ _ _ _|                           |  _ _ _ _ _|
                                                            | |
                                                         _ _| |
                                                     _ _|  _ _|
                                                    |    _|
                                                   _|  _|
                                                  |  _|
                                             _ _ _| |
                                            |  _ _ _|
                                            | |
                                            | |
                                            | |
                 _ _ _ _ _ _ _ _ _ _ _ _ _ _| |
  4   28   56   |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _|
.
Column H gives the nonzero hexagonal numbers (A000384).
a(n) is also the area (and the number of cells) of the n-th diagram.
For n = 3 the sum of the regions (or parts) of the third diagram is 8 + 8 + 8 = 24, so a(3) = 24.
For more information see A237593.
		

Crossrefs

Bisection of A074285.
Some sequences that gives sum of divisors: A000225 (of powers of 2), A008864 (of prime numbers), A065764 (of squares), A073255 (of composites), A074285 (of triangular numbers, also of generalized hexagonal numbers), A139256 (of perfect numbers), A175926 (of cubes), A224613 (of multiples of 6), A346866 (of second hexagonal numbers), A346867 (of numbers with middle divisors), A346868 (of numbers with no middle divisors).

Programs

  • Mathematica
    a[n_] := DivisorSigma[1, n*(2*n - 1)]; Array[a, 50] (* Amiram Eldar, Aug 18 2021 *)
  • PARI
    a(n) = sigma(n*(2*n-1)); \\ Michel Marcus, Aug 18 2021
    
  • Python
    from sympy import divisors
    def a(n): return sum(divisors(n*(2*n - 1)))
    print([a(n) for n in range(1, 48)]) # Michael S. Branicky, Aug 20 2021

Formula

a(n) = A000203(A000384(n)).
Sum_{k=1..n} a(k) ~ 4*n^3/3. - Vaclav Kotesovec, Aug 18 2021

A346866 Sum of divisors of the n-th second hexagonal number.

Original entry on oeis.org

4, 18, 32, 91, 72, 168, 192, 270, 260, 576, 288, 868, 560, 720, 768, 1488, 864, 1482, 1120, 1764, 1408, 2808, 1152, 3420, 2232, 2268, 2880, 4480, 1800, 4464, 3328, 5292, 3264, 5184, 3456, 6734, 4712, 5760, 4480, 10890, 3528, 10368, 5280, 7560, 8736, 9216, 5760, 12152
Offset: 1

Views

Author

Omar E. Pol, Aug 17 2021

Keywords

Comments

The characteristic shape of the symmetric representation of a(n) consists in that in the main diagonal of the diagram the smallest Dyck path has a peak and the largest Dyck path has a valley.
So knowing this characteristic shape we can know if a number is a second hexagonal number (or not) just by looking at the diagram, even ignoring the concept of second hexagonal number.
Therefore we can see a geometric pattern of the distribution of the second hexagonal numbers in the stepped pyramid described in A245092.

Examples

			a(3) = 32 because the sum of divisors of the third second hexagonal number (i.e., 21) is 1 + 3 + 7 + 21 = 32.
On the other hand we can see that in the main diagonal of every diagram the smallest Dyck path has a peak and the largest Dyck path has a valley as shown below.
Illustration of initial terms:
---------------------------------------------------------------------------------------
  n  h(n)  a(n)  Diagram
---------------------------------------------------------------------------------------
                    _             _                     _                            _
                   | |           | |                   | |                          | |
                _ _|_|           | |                   | |                          | |
  1    3    4  |_ _|             | |                   | |                          | |
                                 | |                   | |                          | |
                              _ _| |                   | |                          | |
                             |  _ _|                   | |                          | |
                          _ _|_|                       | |                          | |
                         |  _|                         | |                          | |
                _ _ _ _ _| |                           | |                          | |
  2   10   18  |_ _ _ _ _ _|                           | |                          | |
                                                _ _ _ _|_|                          | |
                                               | |                                  | |
                                              _| |                                  | |
                                             |  _|                                  | |
                                          _ _|_|                                    | |
                                      _ _|  _|                                      | |
                                     |_ _ _|                                        | |
                                     |                                 _ _ _ _ _ _ _| |
                                     |                                |    _ _ _ _ _ _|
                _ _ _ _ _ _ _ _ _ _ _|                                |   |
  3   21   32  |_ _ _ _ _ _ _ _ _ _ _|                             _ _|   |
                                                                  |       |
                                                                 _|    _ _|
                                                                |     |
                                                             _ _|    _|
                                                         _ _|      _|
                                                        |        _|
                                                   _ _ _|    _ _|
                                                  |         |
                                                  |  _ _ _ _|
                                                  | |
                                                  | |
                                                  | |
                                                  | |
               _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _| |
  4   36   91 |_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _|
.
Column h gives the n-th second hexagonal number (A014105).
The widths of the main diagonal of the diagrams are [0, 0, 0, 1] respectively.
a(n) is the area (and the number of cells) of the n-th diagram.
For n = 3 the sum of the regions (or parts) of the third diagram is 11 + 5 + 5 + 11 = 32, so a(3) = 32.
For n = 4 there is only one region (or part) of size 91 in the fourth diagram so a(4) = 91.
		

Crossrefs

Bisection of A074285.
Some sequences that gives sum of divisors: A000225 (of powers of 2), A008864 (of prime numbers), A065764 (of squares), A073255 (of composites), A074285 (of triangular numbers, also of generalized hexagonal numbers), A139256 (of perfect numbers), A175926 (of cubes), A224613 (of multiples of 6), A346865 (of hexagonal numbers), A346867 (of numbers with middle divisors), A346868 (of numbers with no middle divisors), A347155 (of nontriangular numbers).

Programs

  • Mathematica
    a[n_] := DivisorSigma[1, n*(2*n + 1)]; Array[a, 50] (* Amiram Eldar, Aug 18 2021 *)
  • PARI
    a(n) = sigma(n*(2*n + 1)); \\ Michel Marcus, Aug 18 2021

Formula

a(n) = A000203(A014105(n)).
Sum_{k=1..n} a(k) ~ 4*n^3/3. - Amiram Eldar, Dec 31 2024

A116990 Indices of triangular numbers whose sum of divisors is square.

Original entry on oeis.org

1, 2, 11, 20, 40, 68, 92, 212, 236, 253, 266, 321, 328, 452, 582, 589, 596, 668, 695, 716, 782, 788, 836, 928, 932, 970, 991, 1012, 1065, 1076, 1173, 1264, 1300, 1336, 1388, 1436, 1490, 1549, 1796, 1854, 1927, 1995, 2159, 2228, 2252, 2468, 2545, 2588
Offset: 1

Views

Author

Jonathan Vos Post, Apr 04 2006

Keywords

Comments

From Zak Seidov, Oct 19 2010: (Start)
A074285(n) = A000203(A000217(n)) = s^2.
Corresponding values of s begin: 1,2,12,24,42,72,96,216,240,192,240,288,336,456,504, 480,600,672,840,720,720,792,960,930,936,756,992,936,1008,1080,... (are most values of s multiples of 3?).
(End)

Examples

			a(1) = 1 because sigma(1*2/2) = sigma(1) = 1 = 1^2,
a(2) = 2 because sigma(2*3/2) = sigma(3) = 2^2,
a(3) = 11 because sigma(11*12/2) = sigma(66) = 144 = 12^2.
		

Crossrefs

See also: A000217 Triangular numbers: a(n) = C(n+1,2) = n(n+1)/2 = 0+1+2+...+n. A074285 Sum of the divisors of n-th triangular number. A083675 Triangular number for which the sum of the proper divisors is also a triangular number. A000203 sigma(n) = sum of divisors of n. Also called sigma_1(n).

Programs

  • Maple
    with(numtheory): a:=proc(n) if type(sqrt(sigma(n*(n+1)/2)),integer)=true then n else fi end: seq(a(n),n=0..3100); # Emeric Deutsch, Apr 06 2006
  • Mathematica
    Flatten@ Position[Accumulate[Range@ 2600], n_ /; IntegerQ@ Sqrt@ DivisorSigma[1, n] == True] (* Michael De Vlieger, Mar 17 2015 *)
    Select[Range[2600],IntegerQ[Sqrt[DivisorSigma[1,(#(#+1))/2]]]&] (* Harvey P. Dale, Nov 19 2022 *)
  • PARI
    for(n=1,1000,if(issquare(sigma(n*(n+1)/2)),print1(n","))) \\ Zak Seidov, Mar 21 2015

Formula

n such that A074285(n) is in A000290.
n such that sum( d | A000217(n), d ) is in A000290.
n such that A000203(A000217(n)) is in A000290.
n such that sum( d | n*(n+1)/2, d ) = k^2 for integer k.

Extensions

More terms from Emeric Deutsch, Apr 06 2006
Incorrect term 0 removed by Michel Marcus, Mar 17 2015

A346867 Sum of divisors of the numbers that have middle divisors.

Original entry on oeis.org

1, 3, 7, 12, 15, 13, 28, 24, 31, 39, 42, 60, 31, 56, 72, 63, 48, 91, 90, 96, 78, 124, 57, 93, 120, 120, 168, 104, 127, 144, 144, 195, 96, 186, 121, 224, 180, 234, 112, 252, 171, 156, 217, 210, 280, 216, 248, 182, 360, 133, 312, 255, 252, 336, 240, 336, 168, 403, 372, 234
Offset: 1

Views

Author

Omar E. Pol, Aug 18 2021

Keywords

Comments

The characteristic shape of the symmetric representation of a(n) consists in that in the main diagonal of the diagram the width is >= 1.
Also the width on the main diagonal equals the number of middle divisors.
So knowing this characteristic shape we can know if a number has middle divisors (or not) and the number of them just by looking at the diagram, even ignoring the concept of middle divisors.
Therefore we can see a geometric pattern of the distribution of the numbers with middle divisors in the stepped pyramid described in A245092.
For the definition of "width" see A249351.

Examples

			a(4) = 12 because the sum of divisors of the fourth number that has middle divisors (i.e., 6) is 1 + 2 + 3 + 6 = 12.
On the other hand we can see that in the main diagonal of every diagram the width is >= 1 as shown below.
Illustration of initial terms:
m(n) = A071562(n).
.
   n   m(n) a(n)   Diagram
.                  _ _   _   _   _ _     _     _ _   _   _       _
   1    1    1    |_| | | | | | | | |   | |   | | | | | | |     | |
   2    2    3    |_ _|_| | | | | | |   | |   | | | | | | |     | |
                   _ _|  _|_| | | | |   | |   | | | | | | |     | |
   3    4    7    |_ _ _|    _|_| | |   | |   | | | | | | |     | |
                   _ _ _|  _|  _ _|_|   | |   | | | | | | |     | |
   4    6   12    |_ _ _ _|  _| |  _ _ _| |   | | | | | | |     | |
                   _ _ _ _| |_ _|_|    _ _|   | | | | | | |     | |
   5    8   15    |_ _ _ _ _|  _|     |  _ _ _|_| | | | | |     | |
   6    9   13    |_ _ _ _ _| |      _|_| |  _ _ _|_| | | |     | |
                              |  _ _|    _| |    _ _ _|_| |     | |
                   _ _ _ _ _ _| |  _|  _|  _|   |  _ _ _ _|     | |
   7   12   28    |_ _ _ _ _ _ _| |_ _|  _|  _ _| |    _ _ _ _ _| |
                                  |  _ _|  _|    _|   |    _ _ _ _|
                   _ _ _ _ _ _ _ _| |     |     |  _ _|   |
   8   15   24    |_ _ _ _ _ _ _ _| |  _ _|  _ _|_|       |
   9   16   31    |_ _ _ _ _ _ _ _ _| |  _ _|  _|      _ _|
                   _ _ _ _ _ _ _ _ _| | |     |      _|
  10   18   39    |_ _ _ _ _ _ _ _ _ _| |  _ _|    _|
                   _ _ _ _ _ _ _ _ _ _| | |       |
  11   20   42    |_ _ _ _ _ _ _ _ _ _ _| |  _ _ _|
                                          | |
                                          | |
                   _ _ _ _ _ _ _ _ _ _ _ _| |
  12   24   60    |_ _ _ _ _ _ _ _ _ _ _ _ _|
.
The n-th diagram has the property that at least it shares a vertex with the (n+1)-st diagram.
		

Crossrefs

Some sequences that gives sum of divisors: A000225 (of powers of 2), A008864 (of prime numbers), A065764 (of squares), A073255 (of composites), A074285 (of triangular numbers, also of generalized hexagonal numbers), A139256 (of perfect numbers), A175926 (of cubes), A224613 (of multiples of 6), A346865 (of hexagonal numbers), A346866 (of second hexagonal numbers), A346868 (of numbers with no middle divisors).

Programs

  • Mathematica
    s[n_] := Module[{d = Divisors[n]}, If[AnyTrue[d, Sqrt[n/2] <= # < Sqrt[n*2] &], Plus @@ d, 0]]; Select[Array[s, 150], # > 0 &] (* Amiram Eldar, Aug 19 2021 *)
  • PARI
    is(n) = fordiv(n, d, if(d^2>=n/2 && d^2<2*n, return(1))); 0 ; \\ A071562
    apply(sigma, select(is, [1..200])) \\ Michel Marcus, Aug 19 2021

Formula

a(n) = A000203(A071562(n)).

A117948 Sum of the divisors of pentagonal numbers.

Original entry on oeis.org

1, 6, 28, 36, 48, 72, 144, 168, 182, 180, 372, 576, 280, 336, 864, 720, 558, 702, 1120, 1080, 1024, 1008, 1296, 2016, 1178, 1344, 3600, 2016, 1320, 2160, 2304, 3720, 2736, 1836, 4704, 4212, 2736, 2280, 5040, 6048, 2604, 4992, 5588, 4752, 5304, 3312
Offset: 1

Views

Author

Jonathan Vos Post, Apr 04 2006

Keywords

Examples

			a(5) = 48 because the divisors of 5*(3*5-1)/2 = 35 are 1, 5, 7 and 35, having sum 48.
		

Crossrefs

Programs

  • Maple
    with(numtheory): seq(sigma(n*(3*n-1)/2),n=1..46); # Emeric Deutsch, Apr 15 2006
  • Mathematica
    Table[DivisorSigma[1,(n(3n-1))/2],{n,50}] (* Harvey P. Dale, Jul 05 2015 *)
  • PARI
    a(n) = sigma(n*(3*n-1)/2); \\ Amiram Eldar, Dec 14 2024

Formula

a(n) = sigma(n(3n-1)/2), where sigma is the sum of divisors function. - Emeric Deutsch, Apr 15 2006
Sum_{k=1..n} a(k) ~ (9/10) * n^3. - Amiram Eldar, Dec 14 2024

Extensions

a(0) removed by Amiram Eldar, Aug 17 2019

A347108 a(n) = Sum_{k=1..n} sigma(k)*sigma(2*k), where sigma(n) = A000203(n) is the sum of the divisors of n.

Original entry on oeis.org

3, 24, 72, 177, 285, 621, 813, 1278, 1785, 2541, 2973, 4653, 5241, 6585, 8313, 10266, 11238, 14787, 15987, 19767, 22839, 25863, 27591, 35031, 37914, 42030, 46830, 53550, 56250, 68346, 71418, 79419, 86331, 93135, 100047, 117792, 122124, 130524, 139932, 156672
Offset: 1

Views

Author

Vaclav Kotesovec, Aug 18 2021

Keywords

Crossrefs

Programs

  • Mathematica
    Accumulate[Table[DivisorSigma[1,k] * DivisorSigma[1,2*k], {k, 1, 100}]]
  • PARI
    a(n) = sum(k=1, n, sigma(k)*sigma(2*k)); \\ Michel Marcus, Aug 18 2021

Formula

a(n) ~ 2*zeta(3)*n^3.
Showing 1-10 of 17 results. Next