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

A000290 The squares: a(n) = n^2.

Original entry on oeis.org

0, 1, 4, 9, 16, 25, 36, 49, 64, 81, 100, 121, 144, 169, 196, 225, 256, 289, 324, 361, 400, 441, 484, 529, 576, 625, 676, 729, 784, 841, 900, 961, 1024, 1089, 1156, 1225, 1296, 1369, 1444, 1521, 1600, 1681, 1764, 1849, 1936, 2025, 2116, 2209, 2304, 2401, 2500
Offset: 0

Views

Author

Keywords

Comments

To test if a number is a square, see Cohen, p. 40. - N. J. A. Sloane, Jun 19 2011
Zero followed by partial sums of A005408 (odd numbers). - Jeremy Gardiner, Aug 13 2002
Begin with n, add the next number, subtract the previous number and so on ending with subtracting a 1: a(n) = n + (n+1) - (n-1) + (n+2) - (n-2) + (n+3) - (n-3) + ... + (2n-1) - 1 = n^2. - Amarnath Murthy, Mar 24 2004
Sum of two consecutive triangular numbers A000217. - Lekraj Beedassy, May 14 2004
Numbers with an odd number of divisors: {d(n^2) = A048691(n); for the first occurrence of 2n + 1 divisors, see A071571(n)}. - Lekraj Beedassy, Jun 30 2004
See also A000037.
First sequence ever computed by electronic computer, on EDSAC, May 06 1949 (see Renwick link). - Russ Cox, Apr 20 2006
Numbers k such that the imaginary quadratic field Q(sqrt(-k)) has four units. - Marc LeBrun, Apr 12 2006
For n > 0: number of divisors of (n-1)th power of any squarefree semiprime: a(n) = A000005(A006881(k)^(n-1)); a(n) = A000005(A000400(n-1)) = A000005(A011557(n-1)) = A000005(A001023(n-1)) = A000005(A001024(n-1)). - Reinhard Zumkeller, Mar 04 2007
If a 2-set Y and an (n-2)-set Z are disjoint subsets of an n-set X then a(n-2) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
Numbers a such that a^1/2 + b^1/2 = c^1/2 and a^2 + b = c. - Cino Hilliard, Feb 07 2008 (this comment needs clarification, Joerg Arndt, Sep 12 2013)
Numbers k such that the geometric mean of the divisors of k is an integer. - Ctibor O. Zizka, Jun 26 2008
Equals row sums of triangle A143470. Example: 36 = sum of row 6 terms: (23 + 7 + 3 + 1 + 1 + 1). - Gary W. Adamson, Aug 17 2008
Equals row sums of triangles A143595 and A056944. - Gary W. Adamson, Aug 26 2008
Number of divisors of 6^(n-1) for n > 0. - J. Lowell, Aug 30 2008
Denominators of Lyman spectrum of hydrogen atom. Numerators are A005563. A000290-A005563 = A000012. - Paul Curtz, Nov 06 2008
a(n) is the number of all partitions of the sum 2^2 + 2^2 + ... + 2^2, (n-1) times, into powers of 2. - Valentin Bakoev, Mar 03 2009
a(n) is the maximal number of squares that can be 'on' in an n X n board so that all the squares turn 'off' after applying the operation: in any 2 X 2 sub-board, a square turns from 'on' to 'off' if the other three are off. - Srikanth K S, Jun 25 2009
Zero together with the numbers k such that 2 is the number of perfect partitions of k. - Juri-Stepan Gerasimov, Sep 26 2009
Totally multiplicative sequence with a(p) = p^2 for prime p. - Jaroslav Krizek, Nov 01 2009
Satisfies A(x)/A(x^2), A(x) = A173277: (1, 4, 13, 32, 74, ...). - Gary W. Adamson, Feb 14 2010
Positive members are the integers with an odd number of odd divisors and an even number of even divisors. See also A120349, A120359, A181792, A181793, A181795. - Matthew Vandermast, Nov 14 2010
Besides the first term, this sequence is the denominator of Pi^2/6 = 1 + 1/4 + 1/9 + 1/16 + 1/25 + 1/36 + ... . - Mohammad K. Azarian, Nov 01 2011
Partial sums give A000330. - Omar E. Pol, Jan 12 2013
Drmota, Mauduit, and Rivat proved that the Thue-Morse sequence along the squares is normal; see A228039. - Jonathan Sondow, Sep 03 2013
a(n) can be decomposed into the sum of the four numbers [binomial(n, 1) + binomial(n, 2) + binomial(n-1, 1) + binomial(n-1, 2)] which form a "square" in Pascal's Triangle A007318, or the sum of the two numbers [binomial(n, 2) + binomial(n+1, 2)], or the difference of the two numbers [binomial(n+2, 3) - binomial(n, 3)]. - John Molokach, Sep 26 2013
In terms of triangular tiling, the number of equilateral triangles with side length 1 inside an equilateral triangle with side length n. - K. G. Stier, Oct 30 2013
Number of positive roots in the root systems of type B_n and C_n (when n > 1). - Tom Edgar, Nov 05 2013
Squares of squares (fourth powers) are also called biquadratic numbers: A000583. - M. F. Hasler, Dec 29 2013
For n > 0, a(n) is the largest integer k such that k^2 + n is a multiple of k + n. More generally, for m > 0 and n > 0, the largest integer k such that k^(2*m) + n is a multiple of k + n is given by k = n^(2*m). - Derek Orr, Sep 03 2014
For n > 0, a(n) is the number of compositions of n + 5 into n parts avoiding the part 2. - Milan Janjic, Jan 07 2016
a(n), for n >= 3, is also the number of all connected subtrees of a cycle graph, having n vertices. - Viktar Karatchenia, Mar 02 2016
On every sequence of natural continuous numbers with an even number of elements, the summatory of the second half of the sequence minus the summatory of the first half of the sequence is always a square. Example: Sequence from 61 to 70 has an even number of elements (10). Then 61 + 62 + 63 + 64 + 65 = 315; 66 + 67 + 68 + 69 + 70 = 340; 340 - 315 = 25. (n/2)^2 for n = number of elements. - César Aguilera, Jun 20 2016
On every sequence of natural continuous numbers from n^2 to (n+1)^2, the sum of the differences of pairs of elements of the two halves in every combination possible is always (n+1)^2. - César Aguilera, Jun 24 2016
Suppose two circles with radius 1 are tangent to each other as well as to a line not passing through the point of tangency. Create a third circle tangent to both circles as well as the line. If this process is continued, a(n) for n > 0 is the reciprocals of the radii of the circles, beginning with the largest circle. - Melvin Peralta, Aug 18 2016
Does not satisfy Benford's law [Ross, 2012]. - N. J. A. Sloane, Feb 08 2017
Numerators of the solution to the generalization of the Feynman triangle problem, with an offset of 2. If each vertex of a triangle is joined to the point (1/p) along the opposite side (measured say clockwise), then the area of the inner triangle formed by these lines is equal to (p - 2)^2/(p^2 - p + 1) times the area of the original triangle, p > 2. For example, when p = 3, the ratio of the areas is 1/7. The denominators of the ratio of the areas is given by A002061. [Cook & Wood, 2004] - Joe Marasco, Feb 20 2017
Equals row sums of triangle A004737, n >= 1. - Martin Michael Musatov, Nov 07 2017
Right-hand side of the binomial coefficient identity Sum_{k = 0..n} (-1)^(n+k+1)*binomial(n,k)*binomial(n + k,k)*(n - k) = n^2. - Peter Bala, Jan 12 2022
Conjecture: For n>0, min{k such that there exist subsets A,B of {0,1,2,...,a(n)-1} such that |A|=|B|=k and A+B contains {0,1,2,...,a(n)-1}} = n. - Michael Chu, Mar 09 2022
Number of 3-permutations of n elements avoiding the patterns 132, 213, 321. See Bonichon and Sun. - Michel Marcus, Aug 20 2022
Number of intercalates in cyclic Latin squares of order 2n (cyclic Latin squares of odd order do not have intercalates). - Eduard I. Vatutin, Feb 15 2024
a(n) is the number of ternary strings of length n with at most one 0, exactly one 1, and no restriction on the number of 2's. For example, a(3)=9, consisting of the 6 permutations of the string 102 and the 3 permutations of the string 122. - Enrique Navarrete, Mar 12 2025

Examples

			For n = 8, a(8) = 8 * 15 - (1 + 3 + 5 + 7 + 9 + 11 + 13) - 7 = 8 * 15 - 49 - 7 = 64. - _Bruno Berselli_, May 04 2010
G.f. = x + 4*x^2 + 9*x^3 + 16*x^4 + 25*x^5 + 36*x^6 + 49*x^7 + 64*x^8 + 81*x^9 + ...
a(4) = 16. For n = 4 vertices, the cycle graph C4 is A-B-C-D-A. The subtrees are: 4 singles: A, B, C, D; 4 pairs: A-B, BC, C-D, A-D; 4 triples: A-B-C, B-C-D, C-D-A, D-A-B; 4 quads: A-B-C-D, B-C-D-A, C-D-A-B, D-A-B-C; 4 + 4 + 4 + 4 = 16. - _Viktar Karatchenia_, Mar 02 2016
		

References

  • G. L. Alexanderson et al., The William Lowell Putnam Mathematical Competition, Problems and Solutions: 1965-1984, "December 1967 Problem B4(a)", pp. 8(157) MAA Washington DC 1985.
  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 2.
  • Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See Chapter XV, pp. 135-167.
  • R. P. Burn & A. Chetwynd, A Cascade Of Numbers, "The prison door problem" Problem 4 pp. 5-7; 79-80 Arnold London 1996.
  • H. Cohen, A Course in Computational Algebraic Number Theory, Springer, 1996, p. 40.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 31, 36, 38, 63.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), p. 6.
  • M. Gardner, Time Travel and Other Mathematical Bewilderments, Chapter 6 pp. 71-2, W. H. Freeman NY 1988.
  • Granino A. Korn and Theresa M. Korn, Mathematical Handbook for Scientists and Engineers, McGraw-Hill Book Company, New York (1968), p. 982.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.1 Terminology and §8.6 Figurate Numbers, pp. 264, 290-291.
  • Alfred S. Posamentier, The Art of Problem Solving, Section 2.4 "The Long Cell Block" pp. 10-1; 12; 156-7 Corwin Press Thousand Oaks CA 1996.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 35, 52-53, 129-132, 244.
  • Michel Rigo, Formal Languages, Automata and Numeration Systems, 2 vols., Wiley, 2014. Mentions this sequence - see "List of Sequences" in Vol. 2.
  • 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).
  • J. K. Strayer, Elementary Number Theory, Exercise Set 3.3 Problems 32, 33, p. 88, PWS Publishing Co. Boston MA 1996.
  • C. W. Trigg, Mathematical Quickies, "The Lucky Prisoners" Problem 141 pp. 40, 141, Dover NY 1985.
  • R. Vakil, A Mathematical Mosaic, "The Painted Lockers" pp. 127;134 Brendan Kelly Burlington Ontario 1996.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 123.

Crossrefs

Cf. A092205, A128200, A005408, A128201, A002522, A005563, A008865, A059100, A143051, A143470, A143595, A056944, A001157 (inverse Möbius transform), A001788 (binomial transform), A228039, A001105, A004159, A159918, A173277, A095794, A162395, A186646 (Pisano periods), A028338 (2nd diagonal).
A row or column of A132191.
This sequence is related to partitions of 2^n into powers of 2, as it is shown in A002577. So A002577 connects the squares and A000447. - Valentin Bakoev, Mar 03 2009
Boustrophedon transforms: A000697, A000745.
Cf. A342819.
Cf. A013661.

Programs

Formula

G.f.: x*(1 + x) / (1 - x)^3.
E.g.f.: exp(x)*(x + x^2).
Dirichlet g.f.: zeta(s-2).
a(n) = a(-n).
Multiplicative with a(p^e) = p^(2*e). - David W. Wilson, Aug 01 2001
Sum of all matrix elements M(i, j) = 2*i/(i+j) (i, j = 1..n). a(n) = Sum_{i = 1..n} Sum_{j = 1..n} 2*i/(i + j). - Alexander Adamchuk, Oct 24 2004
a(0) = 0, a(1) = 1, a(n) = 2*a(n-1) - a(n-2) + 2. - Miklos Kristof, Mar 09 2005
From Pierre CAMI, Oct 22 2006: (Start)
a(n) is the sum of the odd numbers from 1 to 2*n - 1.
a(0) = 0, a(1) = 1, then a(n) = a(n-1) + 2*n - 1. (End)
For n > 0: a(n) = A130064(n)*A130065(n). - Reinhard Zumkeller, May 05 2007
a(n) = Sum_{k = 1..n} A002024(n, k). - Reinhard Zumkeller, Jun 24 2007
Left edge of the triangle in A132111: a(n) = A132111(n, 0). - Reinhard Zumkeller, Aug 10 2007
Binomial transform of [1, 3, 2, 0, 0, 0, ...]. - Gary W. Adamson, Nov 21 2007
a(n) = binomial(n+1, 2) + binomial(n, 2).
This sequence could be derived from the following general formula (cf. A001286, A000330): n*(n+1)*...*(n+k)*(n + (n+1) + ... + (n+k))/((k+2)!*(k+1)/2) at k = 0. Indeed, using the formula for the sum of the arithmetic progression (n + (n+1) + ... + (n+k)) = (2*n + k)*(k + 1)/2 the general formula could be rewritten as: n*(n+1)*...*(n+k)*(2*n+k)/(k+2)! so for k = 0 above general formula degenerates to n*(2*n + 0)/(0 + 2) = n^2. - Alexander R. Povolotsky, May 18 2008
From a(4) recurrence formula a(n+3) = 3*a(n+2) - 3*a(n+1) + a(n) and a(1) = 1, a(2) = 4, a(3) = 9. - Artur Jasinski, Oct 21 2008
The recurrence a(n+3) = 3*a(n+2) - 3*a(n+1) + a(n) is satisfied by all k-gonal sequences from a(3), with a(0) = 0, a(1) = 1, a(2) = k. - Jaume Oliver Lafont, Nov 18 2008
a(n) = floor(n*(n+1)*(Sum_{i = 1..n} 1/(n*(n+1)))). - Ctibor O. Zizka, Mar 07 2009
Product_{i >= 2} 1 - 2/a(i) = -sin(A063448)/A063448. - R. J. Mathar, Mar 12 2009
a(n) = A002378(n-1) + n. - Jaroslav Krizek, Jun 14 2009
a(n) = n*A005408(n-1) - (Sum_{i = 1..n-2} A005408(i)) - (n-1) = n*A005408(n-1) - a(n-1) - (n-1). - Bruno Berselli, May 04 2010
a(n) == 1 (mod n+1). - Bruno Berselli, Jun 03 2010
a(n) = a(n-1) + a(n-2) - a(n-3) + 4, n > 2. - Gary Detlefs, Sep 07 2010
a(n+1) = Integral_{x >= 0} exp(-x)/( (Pn(x)*exp(-x)*Ei(x) - Qn(x))^2 +(Pi*exp(-x)*Pn(x))^2 ), with Pn the Laguerre polynomial of order n and Qn the secondary Laguerre polynomial defined by Qn(x) = Integral_{t >= 0} (Pn(x) - Pn(t))*exp(-t)/(x-t). - Groux Roland, Dec 08 2010
Euler transform of length-2 sequence [4, -1]. - Michael Somos, Feb 12 2011
A162395(n) = -(-1)^n * a(n). - Michael Somos, Mar 19 2011
a(n) = A004201(A000217(n)); A007606(a(n)) = A000384(n); A007607(a(n)) = A001105(n). - Reinhard Zumkeller, Feb 12 2011
Sum_{n >= 1} 1/a(n)^k = (2*Pi)^k*B_k/(2*k!) = zeta(2*k) with Bernoulli numbers B_k = -1, 1/6, 1/30, 1/42, ... for k >= 0. See A019673, A195055/10 etc. [Jolley eq 319].
Sum_{n>=1} (-1)^(n+1)/a(n)^k = 2^(k-1)*Pi^k*(1-1/2^(k-1))*B_k/k! [Jolley eq 320] with B_k as above.
A007968(a(n)) = 0. - Reinhard Zumkeller, Jun 18 2011
A071974(a(n)) = n; A071975(a(n)) = 1. - Reinhard Zumkeller, Jul 10 2011
a(n) = A199332(2*n - 1, n). - Reinhard Zumkeller, Nov 23 2011
For n >= 1, a(n) = Sum_{d|n} phi(d)*psi(d), where phi is A000010 and psi is A001615. - Enrique Pérez Herrero, Feb 29 2012
a(n) = A000217(n^2) - A000217(n^2 - 1), for n > 0. - Ivan N. Ianakiev, May 30 2012
a(n) = (A000217(n) + A000326(n))/2. - Omar E. Pol, Jan 11 2013
a(n) = A162610(n, n) = A209297(n, n) for n > 0. - Reinhard Zumkeller, Jan 19 2013
a(A000217(n)) = Sum_{i = 1..n} Sum_{j = 1..n} i*j, for n > 0. - Ivan N. Ianakiev, Apr 20 2013
a(n) = A133280(A000217(n)). - Ivan N. Ianakiev, Aug 13 2013
a(2*a(n)+2*n+1) = a(2*a(n)+2*n) + a(2*n+1). - Vladimir Shevelev, Jan 24 2014
a(n+1) = Sum_{t1+2*t2+...+n*tn = n} (-1)^(n+t1+t2+...+tn)*multinomial(t1+t2 +...+tn,t1,t2,...,tn)*4^(t1)*7^(t2)*8^(t3+...+tn). - Mircea Merca, Feb 27 2014
a(n) = floor(1/(1-cos(1/n)))/2 = floor(1/(1-n*sin(1/n)))/6, n > 0. - Clark Kimberling, Oct 08 2014
a(n) = ceiling(Sum_{k >= 1} log(k)/k^(1+1/n)) = -Zeta'[1+1/n]. Thus any exponent greater than 1 applied to k yields convergence. The fractional portion declines from A073002 = 0.93754... at n = 1 and converges slowly to 0.9271841545163232... for large n. - Richard R. Forberg, Dec 24 2014
a(n) = Sum_{j = 1..n} Sum_{i = 1..n} ceiling((i + j - n + 1)/3). - Wesley Ivan Hurt, Mar 12 2015
a(n) = Product_{j = 1..n-1} 2 - 2*cos(2*j*Pi/n). - Michel Marcus, Jul 24 2015
From Ilya Gutkovskiy, Jun 21 2016: (Start)
Product_{n >= 1} (1 + 1/a(n)) = sinh(Pi)/Pi = A156648.
Sum_{n >= 0} 1/a(n!) = BesselI(0, 2) = A070910. (End)
a(n) = A028338(n, n-1), n >= 1 (second diagonal). - Wolfdieter Lang, Jul 21 2017
For n >= 1, a(n) = Sum_{d|n} sigma_2(d)*mu(n/d) = Sum_{d|n} A001157(d)*A008683(n/d). - Ridouane Oudra, Apr 15 2021
a(n) = Sum_{i = 1..2*n-1} ceiling(n - i/2). - Stefano Spezia, Apr 16 2021
From Richard L. Ollerton, May 09 2021: (Start) For n >= 1,
a(n) = Sum_{k=1..n} psi(n/gcd(n,k)).
a(n) = Sum_{k=1..n} psi(gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)).
a(n) = Sum_{k=1..n} sigma_2(n/gcd(n,k))*mu(gcd(n,k))/phi(n/gcd(n,k)).
a(n) = Sum_{k=1..n} sigma_2(gcd(n,k))*mu(n/gcd(n,k))/phi(n/gcd(n,k)). (End)
a(n) = (A005449(n) + A000326(n))/3. - Klaus Purath, May 13 2021
Let T(n) = A000217(n), then a(T(n)) + a(T(n+1)) = T(a(n+1)). - Charlie Marion, Jun 27 2022
a(n) = Sum_{k=1..n} sigma_1(k) + Sum_{i=1..n} (n mod i). - Vadim Kataev, Dec 07 2022
a(n^2) + a(n^2+1) + ... + a(n^2+n) + 4*A000537(n) = a(n^2+n+1) + ... + a(n^2+2n). In general, if P(k,n) = the n-th k-gonal number, then P(2k,n^2) + P(2k,n^2+1) + ... + P(2k,n^2+n) + 4*(k-1)*A000537(n) = P(2k,n^2+n+1) + ... + P(2k,n^2+2n). - Charlie Marion, Apr 26 2024
Sum_{n>=1} 1/a(n) = A013661. - Alois P. Heinz, Oct 19 2024
a(n) = 1 + 3^3*((n-1)/(n+1))^2 + 5^3*((n-1)*(n-2)/((n+1)*(n+2)))^2 + 7^3*((n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)))^2 + ... for n >= 1. - Peter Bala, Dec 09 2024

Extensions

Incorrect comment and example removed by Joerg Arndt, Mar 11 2010

A115067 a(n) = (3*n^2 - n - 2)/2.

Original entry on oeis.org

0, 4, 11, 21, 34, 50, 69, 91, 116, 144, 175, 209, 246, 286, 329, 375, 424, 476, 531, 589, 650, 714, 781, 851, 924, 1000, 1079, 1161, 1246, 1334, 1425, 1519, 1616, 1716, 1819, 1925, 2034, 2146, 2261, 2379, 2500, 2624, 2751, 2881, 3014, 3150, 3289, 3431, 3576
Offset: 1

Views

Author

Roger L. Bagula, Mar 01 2006

Keywords

Comments

Number of orbits of Aut(Z^7) as function of the infinity norm n of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 6720. - Philippe A.J.G. Chevalier, Dec 28 2015
a(n) is the sum of the numerator and denominator of the reduced fraction resulting from the sum A000217(n-2)/A000217(n-1) + A000217(n-1)/A000217(n), n>1. - J. M. Bergot, Jun 10 2017
For n > 1, a(n) is also the number of (not necessarily maximal) cliques in the (n-1)-Andrásfai graph. - Eric W. Weisstein, Nov 29 2017
a(n+1) is the sum of the lengths of all the segments used to draw a square of side n representing the most classic pattern for walls made of 2 X 1 bricks, known as a 1-over-2 pattern, where each joint between neighboring bricks falls over the center of the brick below. - Stefano Spezia, Jun 05 2021

Examples

			Illustrations for n = 2..7 from _Stefano Spezia_, Jun 05 2021:
       _           _ _          _ _ _
      |_|         |_|_|        |_|_ _|
                  |_ _|        |_ _|_|
                               |_|_ _|
   a(2) = 4     a(3) = 11     a(4) = 21
    _ _ _ _     _ _ _ _ _    _ _ _ _ _ _
   |_ _|_ _|   |_ _|_ _|_|  |_ _|_ _|_ _|
   |_|_ _|_|   |_|_ _|_ _|  |_|_ _|_ _|_|
   |_ _|_ _|   |_ _|_ _|_|  |_ _|_ _|_ _|
   |_|_ _|_|   |_|_ _|_ _|  |_|_ _|_ _|_|
               |_ _|_ _|_|  |_ _|_ _|_ _|
                            |_|_ _|_ _|_|
   a(5) = 34    a(6) = 50     a(7) = 69
		

Crossrefs

The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542.
Orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A008585, A005843, A001477, A000217.

Programs

Formula

a(n) = (3*n+2)*(n-1)/2.
a(n+1) = n*(3*n + 5)/2. - Omar E. Pol, May 21 2008
a(n) = 3*n + a(n-1) - 2 for n>1, a(1)=0. - Vincenzo Librandi, Nov 13 2010
a(n) = A095794(-n). - Bruno Berselli, Sep 02 2011
G.f.: x^2*(4-x) / (1-x)^3. - R. J. Mathar, Sep 02 2011
a(n) = A055998(2*n-2) - A055998(n-1). - Bruno Berselli, Sep 23 2016
E.g.f.: exp(x)*x*(8 + 3*x)/2. - Stefano Spezia, May 19 2021
From Amiram Eldar, Feb 22 2022: (Start)
Sum_{n>=2} 1/a(n) = Pi/(5*sqrt(3)) - 3*log(3)/5 + 21/25.
Sum_{n>=2} (-1)^n/a(n) = 4*log(2)/5 - 2*Pi/(5*sqrt(3)) + 9/25. (End)
a(n) = Sum_{j=0..n-2} (2*n-j) = Sum_{j=0..n-2} (n+2+j), for n>=1. See the trapezoid link. - Leo Tavares, May 20 2022

Extensions

Edited by N. J. A. Sloane, Mar 05 2006

A126890 Triangle read by rows: T(n,k) = n*(n+2*k+1)/2, 0 <= k <= n.

Original entry on oeis.org

0, 1, 2, 3, 5, 7, 6, 9, 12, 15, 10, 14, 18, 22, 26, 15, 20, 25, 30, 35, 40, 21, 27, 33, 39, 45, 51, 57, 28, 35, 42, 49, 56, 63, 70, 77, 36, 44, 52, 60, 68, 76, 84, 92, 100, 45, 54, 63, 72, 81, 90, 99, 108, 117, 126, 55, 65, 75, 85, 95, 105, 115, 125, 135, 145, 155, 66, 77, 88
Offset: 0

Views

Author

Reinhard Zumkeller, Dec 30 2006

Keywords

Comments

T(n,k) + T(n,n-k) = A014105(n);
row sums give A059270; Sum_{k=0..n-1} T(n,k) = A000578(n);
central terms give A007742; T(2*n+1,n) = A016754(n);
T(n,0) = A000217(n);
T(n,1) = A000096(n) for n > 0;
T(n,2) = A055998(n) for n > 1;
T(n,3) = A055999(n) for n > 2;
T(n,4) = A056000(n) for n > 3;
T(n,5) = A056115(n) for n > 4;
T(n,6) = A056119(n) for n > 5;
T(n,7) = A056121(n) for n > 6;
T(n,8) = A056126(n) for n > 7;
T(n,10) = A101859(n-1) for n > 9;
T(n,n-3) = A095794(n-1) for n > 2;
T(n,n-2) = A045943(n-1) for n > 1;
T(n,n-1) = A000326(n) for n > 0;
T(n,n) = A005449(n).

Examples

			From _Philippe Deléham_, Oct 03 2011: (Start)
Triangle begins:
   0;
   1,  2;
   3,  5,  7;
   6,  9, 12, 15;
  10, 14, 18, 22, 26;
  15, 20, 25, 30, 35, 40;
  21, 27, 33, 39, 45, 51, 57;
  28, 35, 42, 49, 56, 63, 70, 77; (End)
		

References

  • Léonard Euler, Introduction à l'analyse infinitésimale, tome premier, ACL-Editions, Paris, 1987, p. 353-354.

Crossrefs

Cf. A110449.

Programs

  • Haskell
    a126890 n k = a126890_tabl !! n !! k
    a126890_row n = a126890_tabl !! n
    a126890_tabl = map fst $ iterate
       (\(xs@(x:_), i) -> (zipWith (+) ((x-i):xs) [2*i+1 ..], i+1)) ([0], 0)
    -- Reinhard Zumkeller, Nov 10 2013
  • Mathematica
    Flatten[Table[(n(n+2k+1))/2,{n,0,20},{k,0,n}]] (* Harvey P. Dale, Jun 21 2013 *)

Formula

T(n,k) = T(n,k-1) + n, for k <= n. - Philippe Deléham, Oct 03 2011

A049777 Triangular array read by rows: T(m,n) = n + n+1 + ... + m = (m+n)(m-n+1)/2.

Original entry on oeis.org

1, 3, 2, 6, 5, 3, 10, 9, 7, 4, 15, 14, 12, 9, 5, 21, 20, 18, 15, 11, 6, 28, 27, 25, 22, 18, 13, 7, 36, 35, 33, 30, 26, 21, 15, 8, 45, 44, 42, 39, 35, 30, 24, 17, 9, 55, 54, 52, 49, 45, 40, 34, 27, 19, 10, 66, 65, 63, 60, 56, 51, 45, 38, 30, 21, 11, 78, 77, 75, 72, 68, 63, 57, 50
Offset: 1

Views

Author

Keywords

Comments

Triangle read by rows, T(n,k) = A000217(n) - A000217(k), 0 <= k < n. - Philippe Deléham, Mar 07 2013
Subtriangle of triangle in A049780. - Philippe Deléham, Mar 07 2013
No primes and all composite numbers (except 2^x) are generated after the first two columns of the square array for this sequence. In other words, no primes and all composites except 2^x are generated when m-n >= 2. - Bob Selcoe, Jun 18 2013
Diagonal sums in the square array equal partial sums of squares (A000330). - Bob Selcoe, Feb 14 2014
From Bob Selcoe, Oct 27 2014: (Start)
The following apply to the triangle as a square array read by rows unless otherwise specified (see Table link);
Conjecture: There is at least one prime in interval [T(n,k), T(n,k+1)]. Since T(n,k+1)/T(n,k) decreases to (k+1)/k as n increases, this is true for k=1 ("Bertrand's Postulate", first proved by P. Chebyshev), k=2 (proved by El Bachraoui) and k=3 (proved by Loo).
Starting with T(1,1), The falling diagonal of the first 2 numbers in each column (read by column) are the generalized pentagonal numbers (A001318). That is, the coefficients of T(1,1), T(2,1), T(2,2), T(3,2), T(3,3), T(4,3), T(4,4) etc. are the generalized pentagonal numbers. These are A000326 and A005449 (Pentagonal and Second pentagonal numbers: n*(3*n+1)/2, respectively), interweaved.
Let D(n,k) denote falling diagonals starting with T(n,k):
Treating n as constant: pentagonal numbers of the form n*k + 3*k*(k-1)/2 are D(n,1); sequences A000326, 005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542 are formed by n = 1 through 12, respectively.
Treating k as constant: D(1,k) are (3*n^2 + (4k-5)*n + (k-1)*(k-2))/2. When k = 2(mod3), D(1,k), is same as D(k+1,1) omitting the first (k-2)/3 numbers in the sequences. So D(1,2) is same as D(3,1); D(1,5) is same as D(6,1) omitting the 6; D(1,8) is same as D(9,1) omitting the 9 and 21; etc.
D(1,3) and D(1,4) are sequences A095794 and A140229, respectively.
(End)

Examples

			Rows: {1}; {3,2}; {6,5,3}; ...
Triangle begins:
   1;
   3,  2;
   6,  5,  3;
  10,  9,  7,  4;
  15, 14, 12,  9,  5;
  21, 20, 18, 15, 11,  6;
  28, 27, 25, 22, 18, 13,  7;
  36, 35, 33, 30, 26, 21, 15,  8;
  45, 44, 42, 39, 35, 30, 24, 17,  9;
  55, 54, 52, 49, 45, 40, 34, 27, 19, 10; ...
		

Crossrefs

Row sums = A000330.
Cf. A001318 (generalized pentagonal numbers).
Cf. A000326, 005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542 (pentagonal numbers of form n*k + 3*k*(k-1)/2).

Programs

  • Magma
    /* As triangle */ [[(m+n)*(m-n+1) div 2: n in [1..m]]: m in [1.. 15]]; // Vincenzo Librandi, Oct 27 2014
  • Mathematica
    Flatten[Table[(n+k) (n-k+1)/2,{n,15},{k,n}]] (* Harvey P. Dale, Feb 27 2012 *)
  • PARI
    {T(n,k) = if( k<1 || nMichael Somos, Oct 06 2007 */
    

Formula

Partial sums of A002260 row terms, starting from the right; e.g., row 3 of A002260 = (1, 2, 3), giving (6, 5, 3). - Gary W. Adamson, Oct 23 2007
Sum_{k=0..n-1} (-1)^k*(2*k+1)*A000203(T(n,k)) = (-1)^(n-1)*A000330(n). - Philippe Deléham, Mar 07 2013
Read as a square array: T(n,k) = k*(k+2n-1)/2. - Bob Selcoe, Oct 27 2014

A254963 a(n) = n*(11*n + 3)/2.

Original entry on oeis.org

0, 7, 25, 54, 94, 145, 207, 280, 364, 459, 565, 682, 810, 949, 1099, 1260, 1432, 1615, 1809, 2014, 2230, 2457, 2695, 2944, 3204, 3475, 3757, 4050, 4354, 4669, 4995, 5332, 5680, 6039, 6409, 6790, 7182, 7585, 7999, 8424, 8860, 9307, 9765, 10234, 10714, 11205, 11707
Offset: 0

Views

Author

Bruno Berselli, Feb 11 2015

Keywords

Comments

This sequence provides the first differences of A254407 and the partial sums of A017473.
Also:
a(n) - n = A022269(n);
a(n) + n = n*(11*n+5)/2: 0, 8, 27, 57, 98, 150, 213, 287, ...;
a(n) - 2*n = A022268(n);
a(n) + 2*n = n*(11*n+7)/2: 0, 9, 29, 60, 102, 155, 219, 294, ...;
a(n) - 3*n = n*(11*n-3)/2: 0, 4, 19, 45, 82, 130, 189, 259, ...;
a(n) + 3*n = A211013(n);
a(n) - 4*n = A226492(n);
a(n) + 4*n = A152740(n);
a(n) - 5*n = A180223(n);
a(n) + 5*n = n*(11*n+13)/2: 0, 12, 35, 69, 114, 170, 237, 315, ...;
a(n) - 6*n = A051865(n);
a(n) + 6*n = n*(11*n+15)/2: 0, 13, 37, 72, 118, 175, 243, 322, ...;
a(n) - 7*n = A152740(n-1) with A152740(-1) = 0;
a(n) + 7*n = n*(11*n+17)/2: 0, 14, 39, 75, 122, 180, 249, 329, ...;
a(n) - n*(n-1)/2 = A168668(n);
a(n) + n*(n-1)/2 = A049453(n);
a(n) - n*(n+1)/2 = A202803(n);
a(n) + n*(n+1)/2 = A033580(n).

Crossrefs

Cf. A008729 and A218530 (seventh column); A017473, A254407.
Cf. similar sequences of the type 4*n^2 + k*n*(n+1)/2: A055999 (k=-7, n>6), A028552 (k=-6, n>2), A095794 (k=-5, n>1), A046092 (k=-4, n>0), A000566 (k=-3), A049450 (k=-2), A022264 (k=-1), A016742 (k=0), A022267 (k=1), A202803 (k=2), this sequence (k=3), A033580 (k=4).
Cf. A069125: (2*n+1)^2 + 3*n*(n+1)/2; A147875: n^2 + 3*n*(n+1)/2.

Programs

  • Magma
    [n*(11*n+3)/2: n in [0..50]];
    
  • Mathematica
    Table[n (11 n + 3)/2, {n, 0, 50}]
    LinearRecurrence[{3,-3,1},{0,7,25},50] (* Harvey P. Dale, Mar 25 2018 *)
  • Maxima
    makelist(n*(11*n+3)/2, n, 0, 50);
  • PARI
    vector(50, n, n--; n*(11*n+3)/2)
    
  • Sage
    [n*(11*n+3)/2 for n in (0..50)]
    

Formula

G.f.: x*(7 + 4*x)/(1 - x)^3.
From Elmo R. Oliveira, Dec 15 2024: (Start)
E.g.f.: exp(x)*x*(14 + 11*x)/2.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2. (End)

A327618 Number A(n,k) of parts in all k-times partitions of n; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 5, 6, 1, 0, 1, 7, 14, 12, 1, 0, 1, 9, 25, 44, 20, 1, 0, 1, 11, 39, 109, 100, 35, 1, 0, 1, 13, 56, 219, 315, 274, 54, 1, 0, 1, 15, 76, 386, 769, 1179, 581, 86, 1, 0, 1, 17, 99, 622, 1596, 3643, 3234, 1417, 128, 1, 0, 1, 19, 125, 939, 2960, 9135, 12336, 10789, 2978, 192, 1
Offset: 0

Views

Author

Alois P. Heinz, Sep 19 2019

Keywords

Comments

Row n is binomial transform of the n-th row of triangle A327631.

Examples

			A(2,2) = 5 = 1+2+2 counting the parts in 2, 11, 1|1.
Square array A(n,k) begins:
  0,  0,   0,    0,     0,     0,     0,      0, ...
  1,  1,   1,    1,     1,     1,     1,      1, ...
  1,  3,   5,    7,     9,    11,    13,     15, ...
  1,  6,  14,   25,    39,    56,    76,     99, ...
  1, 12,  44,  109,   219,   386,   622,    939, ...
  1, 20, 100,  315,   769,  1596,  2960,   5055, ...
  1, 35, 274, 1179,  3643,  9135, 19844,  38823, ...
  1, 54, 581, 3234, 12336, 36911, 93302, 208377, ...
		

Crossrefs

Columns k=0-3 give: A057427, A006128, A327594, A327627.
Rows n=0-3 give: A000004, A000012, A005408, A095794(k+1).
Main diagonal gives A327619.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, [1, 0],
         `if`(k=0, [1, 1], `if`(i<2, 0, b(n, i-1, k))+
             (h-> (f-> f +[0, f[1]*h[2]/h[1]])(h[1]*
            b(n-i, min(n-i, i), k)))(b(i$2, k-1))))
        end:
    A:= (n, k)-> b(n$2, k)[2]:
    seq(seq(A(n, d-n), n=0..d), d=0..14);
  • Mathematica
    b[n_, i_, k_] := b[n, i, k] = If[n == 0, {1, 0}, If[k == 0, {1, 1}, If[i < 2, 0, b[n, i - 1, k]] + Function[h, Function[f, f + {0, f[[1]] h[[2]]/ h[[1]]}][h[[1]] b[n - i, Min[n - i, i], k]]][b[i, i, k - 1]]]];
    A[n_, k_] := b[n, n, k][[2]];
    Table[A[n, d-n], {d, 0, 14}, {n, 0, d}] // Flatten (* Jean-François Alcover, Apr 30 2020, after Alois P. Heinz *)

Formula

A(n,k) = Sum_{i=0..k} binomial(k,i) * A327631(n,i).

A003777 a(n) = n^3 + n^2 - 1.

Original entry on oeis.org

1, 11, 35, 79, 149, 251, 391, 575, 809, 1099, 1451, 1871, 2365, 2939, 3599, 4351, 5201, 6155, 7219, 8399, 9701, 11131, 12695, 14399, 16249, 18251, 20411, 22735, 25229, 27899, 30751, 33791, 37025, 40459, 44099, 47951, 52021, 56315, 60839, 65599, 70601, 75851
Offset: 1

Views

Author

N. J. A. Sloane, Jun 14 1998

Keywords

Comments

This sequence in related to A095794 by a(n) = n*A095794(n) - Sum_{i=1..n-1} A095794(i), for n > 1. - Bruno Berselli, Dec 28 2010
a(n) is the area of a triangle with vertices at points (n-1,(n-1)^2), (n,n^2), and ((n+1)^2,n+1). - J. M. Bergot, Jun 03 2014
Old name was: "Number of stacks of n pikelets, distance 3 flips from a well-ordered stack".

Crossrefs

Programs

Formula

G.f.: x*(1+7*x-3*x^2+x^3)/(1-x)^4. Also, a(n) = 2*A002411(n) - 1 = A000578(n-1) + A001107(n). - Bruno Berselli, Dec 28 2010
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n > 4. - Wesley Ivan Hurt, Oct 08 2017
E.g.f.: 1 + (-1 + 2*x + 4*x^2 + x^3)*exp(x). - G. C. Greubel, Jan 03 2024

Extensions

More terms from Wesley Ivan Hurt, Jun 04 2014
Entry revised by N. J. A. Sloane, Jun 15 2014

A283394 a(n) = 3*n*(3*n + 7)/2 + 4.

Original entry on oeis.org

4, 19, 43, 76, 118, 169, 229, 298, 376, 463, 559, 664, 778, 901, 1033, 1174, 1324, 1483, 1651, 1828, 2014, 2209, 2413, 2626, 2848, 3079, 3319, 3568, 3826, 4093, 4369, 4654, 4948, 5251, 5563, 5884, 6214, 6553, 6901, 7258, 7624, 7999, 8383, 8776, 9178, 9589, 10009
Offset: 0

Views

Author

Bruno Berselli, Mar 23 2017

Keywords

Comments

Sum_{k = 0..n} (3*k + r)^3 is divisible by 3*n*(3*n + 2*r + 3)/2 + r^2: the sequence corresponds to the case r = 2 of this formula (other cases are listed in Crossrefs section).
Also, Sum_{k = 0..n} (3*k + 2)^3 / a(n) gives 2, 7, 15, 26, 40, 57, 77, 100, 126, 155, 187, 222, ... (A005449).
a(n) is even if n belongs to A014601. No term is divisible by 3, 5, 7 and 11.

Crossrefs

Sequences with formula 3*n*(3*n + 2*r + 3)/2 + r^2: A038764 (r=-1), A027468 (r=0), A081271 (r=1), this sequence (r=2), A027468 (r=3; offset: -1), A080855 (r=4; offset: -2).

Programs

  • Magma
    [3*n*(3*n+7)/2+4: n in [0..50]];
    
  • Mathematica
    Table[3 n (3 n + 7)/2 + 4, {n, 0, 50}]
    LinearRecurrence[{3,-3,1},{4,19,43},50] (* Harvey P. Dale, Mar 02 2019 *)
  • Maxima
    makelist(3*n*(3*n+7)/2+4, n, 0, 50);
    
  • PARI
    a(n) = 3*n*(3*n + 7)/2 + 4; \\ Indranil Ghosh, Mar 24 2017
  • Python
    [3*n*(3*n+7)/2+4 for n in range(50)]
    
  • Sage
    [3*n*(3*n+7)/2+4 for n in range(50)]
    

Formula

O.g.f.: (4 + 7*x - 2*x^2)/(1 - x)^3.
E.g.f.: (8 + 30*x + 9*x^2)*exp(x)/2.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = A081271(-n-2).
a(n) = 3*A095794(n+1) + 1.
a(n) = A034856(3*n+2) = A101881(6*n+2) = A165157(6*n+3) = A186349(6*n+3).
The inverse binomial transform yields 4, 15, 9, 0 (0 continued), therefore:
a(n) = 4*binomial(n,0) + 15*binomial(n,1) + 9*binomial(n,2).

A096038 Triangle T(n,m) = (3*n^2-3*m^2+5*m-4+n)/2 read by rows.

Original entry on oeis.org

1, 6, 4, 14, 12, 7, 25, 23, 18, 10, 39, 37, 32, 24, 13, 56, 54, 49, 41, 30, 16, 76, 74, 69, 61, 50, 36, 19, 99, 97, 92, 84, 73, 59, 42, 22, 125, 123, 118, 110, 99, 85, 68, 48, 25, 154, 152, 147, 139, 128, 114, 97, 77, 54, 28, 186, 184, 179, 171, 160, 146, 129, 109, 86, 60, 31
Offset: 1

Views

Author

Gary W. Adamson, Jun 17 2004

Keywords

Comments

The triangle is obtained by subtracting the triangle A094930 from
its square root (also described in A094930) and then dividing each element of column m through 3*m-1.
For the first three rows n=1 to 3 this yields for example:
4;.................2;............2......................1;
14,25;......minus..2,5;.......=..12,20;......->.divide..6,4;
30,65,64;..........2,5,8;........28,60,56;..............14;12,7;

Crossrefs

Programs

  • Python
    def A096038(n,m):
        return (3*n**2-3*m**2+5*m-4+n)//2
    print( [A096038(n,m) for n in range(20) for m in range(1,n+1)] )
    # R. J. Mathar, Oct 11 2009

Formula

T(n,1) = A095794(n).
T(n,n) = 3*n-2.
T(n,m) = A094930(n,m)/(3*m-1)-1.

Extensions

Edited, T(3,2) corrected, and extended by R. J. Mathar, Oct 11 2009

A116689 Partial sums of dodecahedral numbers (A006566).

Original entry on oeis.org

0, 1, 21, 105, 325, 780, 1596, 2926, 4950, 7875, 11935, 17391, 24531, 33670, 45150, 59340, 76636, 97461, 122265, 151525, 185745, 225456, 271216, 323610, 383250, 450775, 526851, 612171, 707455, 813450, 930930, 1060696, 1203576, 1360425
Offset: 0

Views

Author

Jonathan Vos Post, Mar 15 2006

Keywords

Comments

Geometrically, the partial sums of dodecahedral numbers may be interpreted as 4-dimensional dodecahedral hyperpyramidal numbers. It is somewhat surprising that this is (with proper offset) the triangular number of the "second pentagonal numbers, minus 1."
Also, the sequence is related to A004188 by a(n) = n*A004188(n)-sum(A004188(i), i=0..n-1). [Bruno Berselli, Apr 05 2012]

Crossrefs

Programs

Formula

a(n) = Sum_{i = 0..n} A006566(i).
a(n) = Sum_{i = 0..n} i*(3*i-1)*(3*i-2)/2.
a(n+1) = A000217(A095794(n)).
a(n+1) = A000217(A005449(n) - 1).
a(n+1) = A000217(n*(3n+1)/2-1).
a(n+1) = A000217(A001844(n) - A000217(n+1) - 1).
a(n) = n*(9*n^3+6*n^2-5*n-2)/8. G.f.: x*(1+16*x+10*x^2)/(1-x)^5. [Colin Barker, Apr 04 2012]
a(n) = binomial(A005449(n), 2). - Wesley Ivan Hurt, Oct 06 2013
From Peter Bala, Sep 03 2023: (Start)
a(n) = n*(n + 1)*(3*n + 1)*(3*n - 2)/8.
a(n) = Sum_{0 <= i <= j <= n-1} (3*i + 1)*(3*j + 1). Cf. A024212 (End)
Showing 1-10 of 25 results. Next