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

A259318 a(n) = A259109(n)*A006331(n) - A259108(n)^2.

Original entry on oeis.org

0, 0, 144, 6048, 85536, 679536, 3747744, 16039296, 56930688, 174978144, 479700144, 1198248480, 2770653600, 6002352720, 12298837824, 24014605824, 44957265408, 81097765056, 141549364944, 239891292576, 395928108576, 637992775728, 1005920381664, 1554840524160
Offset: 0

Views

Author

N. J. A. Sloane, Jun 24 2015

Keywords

Examples

			n=3: 6048 = 1588*28 - 196^2.
		

Crossrefs

Programs

  • PARI
    concat([0,0], Vec(-144*x^2*(x^6+31*x^5+187*x^4+330*x^3+187*x^2 +31*x+1)/(x-1)^11 + O(x^100))) \\ Colin Barker, Jun 29 2015

Formula

a(n) = (16*n^10+80*n^9+120*n^8-147*n^6-105*n^5+5*n^4+25*n^3+6*n^2)/525. - Colin Barker, Jun 29 2015
G.f.: -144*x^2*(x^6+31*x^5+187*x^4+330*x^3+187*x^2+31*x+1) / (x-1)^11. - Colin Barker, Jun 29 2015

A165453 Linear interpolation of the sequence that maps an entry of A002378 to the corresponding entry of A006331.

Original entry on oeis.org

0, 1, 2, 4, 6, 8, 10, 13, 16, 19, 22, 25, 28, 32, 36, 40, 44, 48, 52, 56, 60, 65, 70, 75, 80, 85, 90, 95, 100, 105, 110, 116, 122, 128, 134, 140, 146, 152, 158, 164, 170, 176, 182, 189, 196, 203, 210, 217, 224, 231, 238, 245, 252, 259, 266, 273, 280, 288, 296, 304
Offset: 0

Views

Author

Friedrich Regen (friedrich.regen(AT)tu-ilmenau.de), Sep 20 2009

Keywords

Crossrefs

Partial sums of A000194.

Programs

  • Mathematica
    f[n_]:=Round[Sqrt[n]]; a=0;lst={};Do[AppendTo[lst,a+=f[n]],{n,5!}];lst (* Vladimir Joseph Stephan Orlovsky, Oct 13 2009 *)
    Table[Sum[Floor[Sqrt[n + 1 - k] + 1/2], {k, n + 1}], {n, 0, 100}] (* Wesley Ivan Hurt, Dec 01 2020 *)
  • Python
    from math import isqrt
    def A165453(n): return (k:=(m:=isqrt(n))+(n-m*(m+1)>=1))*(3*n+1-k**2)//3 # Chai Wah Wu, Jun 19 2024

Formula

a(0)=0, a(n) = max(min(a(n-1)+x, n+a(n-x))), where the maximum is taken over all values for x from 1 to n.
a(n) = Sum_{k=1..n} floor(sqrt(k)+1/2). - Wesley Ivan Hurt, Dec 01 2020
a(n) = (1/3)*t*(3*n + 1 - t^2), where t = floor(sqrt(n)+1/2). - Ridouane Oudra, Feb 22 2021

A000330 Square pyramidal numbers: a(n) = 0^2 + 1^2 + 2^2 + ... + n^2 = n*(n+1)*(2*n+1)/6.

Original entry on oeis.org

0, 1, 5, 14, 30, 55, 91, 140, 204, 285, 385, 506, 650, 819, 1015, 1240, 1496, 1785, 2109, 2470, 2870, 3311, 3795, 4324, 4900, 5525, 6201, 6930, 7714, 8555, 9455, 10416, 11440, 12529, 13685, 14910, 16206, 17575, 19019, 20540, 22140, 23821, 25585, 27434, 29370
Offset: 0

Views

Author

Keywords

Comments

The sequence contains exactly one square greater than 1, namely 4900 (according to Gardner). - Jud McCranie, Mar 19 2001, Mar 22 2007 [This is a result from Watson. - Charles R Greathouse IV, Jun 21 2013] [See A351830 for further related comments and references.]
Number of rhombi in an n X n rhombus. - Matti De Craene (Matti.DeCraene(AT)rug.ac.be), May 14 2000
Number of acute triangles made from the vertices of a regular n-polygon when n is odd (cf. A007290). - Sen-Peng Eu, Apr 05 2001
Gives number of squares with sides parallel to the axes formed from an n X n square. In a 1 X 1 square, one is formed. In a 2 X 2 square, five squares are formed. In a 3 X 3 square, 14 squares are formed and so on. - Kristie Smith (kristie10spud(AT)hotmail.com), Apr 16 2002; edited by Eric W. Weisstein, Mar 05 2025
a(n-1) = B_3(n)/3, where B_3(x) = x(x-1)(x-1/2) is the third Bernoulli polynomial. - Michael Somos, Mar 13 2004
Number of permutations avoiding 13-2 that contain the pattern 32-1 exactly once.
Since 3*r = (r+1) + r + (r-1) = T(r+1) - T(r-2), where T(r) = r-th triangular number r*(r+1)/2, we have 3*r^2 = r*(T(r+1) - T(r-2)) = f(r+1) - f(r-1) ... (i), where f(r) = (r-1)*T(r) = (r+1)*T(r-1). Summing over n, the right hand side of relation (i) telescopes to f(n+1) + f(n) = T(n)*((n+2) + (n-1)), whence the result Sum_{r=1..n} r^2 = n*(n+1)*(2*n+1)/6 immediately follows. - Lekraj Beedassy, Aug 06 2004
Also as a(n) = (1/6)*(2*n^3 + 3*n^2 + n), n > 0: structured trigonal diamond numbers (vertex structure 5) (cf. A006003 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Number of triples of integers from {1, 2, ..., n} whose last component is greater than or equal to the others.
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Jun 12 2005
Sum of the first n positive squares. - Cino Hilliard, Jun 18 2007
Maximal number of cubes of side 1 in a right pyramid with a square base of side n and height n. - Pasquale CUTOLO (p.cutolo(AT)inwind.it), Jul 09 2007
If a 2-set Y and an (n-2)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 4-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
We also have the identity 1 + (1+4) + (1+4+9) + ... + (1+4+9+16+ ... + n^2) = n(n+1)(n+2)(n+(n+1)+(n+2))/36; ... and in general the k-fold nested sum of squares can be expressed as n(n+1)...(n+k)(n+(n+1)+...+(n+k))/((k+2)!(k+1)/2). - Alexander R. Povolotsky, Nov 21 2007
The terms of this sequence are coefficients of the Engel expansion of the following converging sum: 1/(1^2) + (1/1^2)*(1/(1^2+2^2)) + (1/1^2)*(1/(1^2+2^2))*(1/(1^2+2^2+3^2)) + ... - Alexander R. Povolotsky, Dec 10 2007
Convolution of A000290 with A000012. - Sergio Falcon, Feb 05 2008
Hankel transform of binomial(2*n-3, n-1) is -a(n). - Paul Barry, Feb 12 2008
Starting (1, 5, 14, 30, ...) = binomial transform of [1, 4, 5, 2, 0, 0, 0, ...]. - Gary W. Adamson, Jun 13 2008
Starting (1,5,14,30,...) = second partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = Sum_{i=0..n} binomial(n+2,i+2)*b(i), where b(i)=1,2,0,0,0,... - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
Convolution of A001477 with A005408: a(n) = Sum_{k=0..n} (2*k+1)*(n-k). - Reinhard Zumkeller, Mar 07 2009
Sequence of the absolute values of the z^1 coefficients of the polynomials in the GF1 denominators of A156921. See A157702 for background information. - Johannes W. Meijer, Mar 07 2009
The sequence is related to A000217 by a(n) = n*A000217(n) - Sum_{i=0..n-1} A000217(i) and this is the case d = 1 in the identity n^2*(d*n-d+2)/2 - Sum_{i=0..n-1} i*(d*i-d+2)/2 = n*(n+1)(2*d*n-2*d+3)/6, or also the case d = 0 in n^2*(n+2*d+1)/2 - Sum_{i=0..n-1} i*(i+2*d+1)/2 = n*(n+1)*(2*n+3*d+1)/6. - Bruno Berselli, Apr 21 2010, Apr 03 2012
a(n)/n = k^2 (k = integer) for n = 337; a(337) = 12814425, a(n)/n = 38025, k = 195, i.e., the number k = 195 is the quadratic mean (root mean square) of the first 337 positive integers. There are other such numbers -- see A084231 and A084232. - Jaroslav Krizek, May 23 2010
Also the number of moves to solve the "alternate coins game": given 2n+1 coins (n+1 Black, n White) set alternately in a row (BWBW...BWB) translate (not rotate) a pair of adjacent coins at a time (1 B and 1 W) so that at the end the arrangement shall be BBBBB..BW...WWWWW (Blacks separated by Whites). Isolated coins cannot be moved. - Carmine Suriano, Sep 10 2010
From J. M. Bergot, Aug 23 2011: (Start)
Using four consecutive numbers n, n+1, n+2, and n+3 take all possible pairs (n, n+1), (n, n+2), (n, n+3), (n+1, n+2), (n+1, n+3), (n+2, n+3) to create unreduced Pythagorean triangles. The sum of all six areas is 60*a(n+1).
Using three consecutive odd numbers j, k, m, (j+k+m)^3 - (j^3 + k^3 + m^3) equals 576*a(n) = 24^2*a(n) where n = (j+1)/2. (End)
From Ant King, Oct 17 2012: (Start)
For n > 0, the digital roots of this sequence A010888(a(n)) form the purely periodic 27-cycle {1, 5, 5, 3, 1, 1, 5, 6, 6, 7, 2, 2, 9, 7, 7, 2, 3, 3, 4, 8, 8, 6, 4, 4, 8, 9, 9}.
For n > 0, the units' digits of this sequence A010879(a(n)) form the purely periodic 20-cycle {1, 5, 4, 0, 5, 1, 0, 4, 5, 5, 6, 0, 9, 5, 0, 6, 5, 9, 0, 0}. (End)
Length of the Pisano period of this sequence mod n, n>=1: 1, 4, 9, 8, 5, 36, 7, 16, 27, 20, 11, 72, 13, 28, 45, 32, 17, 108, 19, 40, ... . - R. J. Mathar, Oct 17 2012
Sum of entries of n X n square matrix with elements min(i,j). - Enrique Pérez Herrero, Jan 16 2013
The number of intersections of diagonals in the interior of regular n-gon for odd n > 1 divided by n is a square pyramidal number; that is, A006561(2*n+1)/(2*n+1) = A000330(n-1) = (1/6)*n*(n-1)*(2*n-1). - Martin Renner, Mar 06 2013
For n > 1, a(n)/(2n+1) = A024702(m), for n such that 2n+1 = prime, which results in 2n+1 = A000040(m). For example, for n = 8, 2n+1 = 17 = A000040(7), a(8) = 204, 204/17 = 12 = A024702(7). - Richard R. Forberg, Aug 20 2013
A formula for the r-th successive summation of k^2, for k = 1 to n, is (2*n+r)*(n+r)!/((r+2)!*(n-1)!) (H. W. Gould). - Gary Detlefs, Jan 02 2014
The n-th square pyramidal number = the n-th triangular dipyramidal number (Johnson 12), which is the sum of the n-th + (n-1)-st tetrahedral numbers. E.g., the 3rd tetrahedral number is 10 = 1+3+6, the 2nd is 4 = 1+3. In triangular "dipyramidal form" these numbers can be written as 1+3+6+3+1 = 14. For "square pyramidal form", rebracket as 1+(1+3)+(3+6) = 14. - John F. Richardson, Mar 27 2014
Beukers and Top prove that no square pyramidal number > 1 equals a tetrahedral number A000292. - Jonathan Sondow, Jun 21 2014
Odd numbered entries are related to dissections of polygons through A100157. - Tom Copeland, Oct 05 2014
From Bui Quang Tuan, Apr 03 2015: (Start)
We construct a number triangle from the integers 1, 2, 3, ..., n as follows. The first column contains 2*n-1 integers 1. The second column contains 2*n-3 integers 2, ... The last column contains only one integer n. The sum of all the numbers in the triangle is a(n).
Here is an example with n = 5:
1
1 2
1 2 3
1 2 3 4
1 2 3 4 5
1 2 3 4
1 2 3
1 2
1
(End)
The Catalan number series A000108(n+3), offset 0, gives Hankel transform revealing the square pyramidal numbers starting at 5, A000330(n+2), offset 0 (empirical observation). - Tony Foster III, Sep 05 2016; see Dougherty et al. link p. 2. - Andrey Zabolotskiy, Oct 13 2016
Number of floating point additions in the factorization of an (n+1) X (n+1) real matrix by Gaussian elimination as e.g. implemented in LINPACK subroutines sgefa.f or dgefa.f. The number of multiplications is given by A007290. - Hugo Pfoertner, Mar 28 2018
The Jacobi polynomial P(n-1,-n+2,2,3) or equivalently the sum of dot products of vectors from the first n rows of Pascal's triangle (A007318) with the up-diagonal Chebyshev T coefficient vector (1,3,2,0,...) (A053120) or down-diagonal vector (1,-7,32,-120,400,...) (A001794). a(5) = 1 + (1,1).(1,3) + (1,2,1).(1,3,2) + (1,3,3,1).(1,3,2,0) + (1,4,6,4,1).(1,3,2,0,0) = (1 + (1,1).(1,-7) + (1,2,1).(1,-7,32) + (1,3,3,1).(1,-7,32,-120) + (1,4,6,4,1).(1,-7,32,-120,400))*(-1)^(n-1) = 55. - Richard Turk, Jul 03 2018
Coefficients in the terminating series identity 1 - 5*n/(n + 4) + 14*n*(n - 1)/((n + 4)*(n + 5)) - 30*n*(n - 1)*(n - 2)/((n + 4)*(n + 5)*(n + 6)) + ... = 0 for n = 1,2,3,.... Cf. A002415 and A108674. - Peter Bala, Feb 12 2019
n divides a(n) iff n == +- 1 (mod 6) (see A007310). (See De Koninck reference.) Examples: a(11) = 506 = 11 * 46, and a(13) = 819 = 13 * 63. - Bernard Schott, Jan 10 2020
For n > 0, a(n) is the number of ternary words of length n+2 having 3 letters equal to 2 and 0 only occurring as the last letter. For example, for n=2, the length 4 words are 2221,2212,2122,1222,2220. - Milan Janjic, Jan 28 2020
Conjecture: Every integer can be represented as a sum of three generalized square pyramidal numbers. A related conjecture is given in A336205 corresponding to pentagonal case. A stronger version of these conjectures is that every integer can be expressed as a sum of three generalized r-gonal pyramidal numbers for all r >= 3. In here "generalized" means negative indices are included. - Altug Alkan, Jul 30 2020
The natural number y is a term if and only if y = a(floor((3 * y)^(1/3))). - Robert Israel, Dec 04 2024
Also the number of directed bishop moves on an n X n chessboard, where two moves are considered the same if one can be obtained from the other by a rotation of the board. Reflections are ignored. Equivalently, number of directed bishop moves on an n X n chessboard, where two moves are considered the same if one can be obtained from the other by an axial reflection of the board (horizontal or vertical). Rotations and diagonal reflections are ignored. - Hilko Koning, Aug 22 2025

Examples

			G.f. = x + 5*x^2 + 14*x^3 + 30*x^4 + 55*x^5 + 91*x^6 + 140*x^7 + 204*x^8 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 813.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover Publications, NY, 1964, p. 194.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 215,223.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 122, see #19 (3(1)), I(n); p. 155.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 47-49.
  • H. S. M. Coxeter, Polyhedral numbers, pp. 25-35 of R. S. Cohen, J. J. Stachel and M. W. Wartofsky, eds., For Dirk Struik: Scientific, historical and political essays in honor of Dirk J. Struik, Reidel, Dordrecht, 1974.
  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p.165).
  • J. M. De Koninck and A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 310, pp. 46-196, Ellipses, Paris, 2004.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
  • 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. 2.
  • M. Gardner, Fractal Music, Hypercards and More, Freeman, NY, 1991, p. 293.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 293.
  • M. Holt, Math puzzles and games, Walker Publishing Company, 1977, p. 2 and p. 89.
  • Simon Singh, The Simpsons and Their Mathematical Secrets. London: Bloomsbury Publishing PLC (2013): 188.
  • 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).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 126.

Crossrefs

Sums of 2 consecutive terms give A005900.
Column 0 of triangle A094414.
Column 1 of triangle A008955.
Right side of triangle A082652.
Row 2 of array A103438.
Partial sums of A000290.
Cf. similar sequences listed in A237616 and A254142.
Cf. |A084930(n, 1)|.
Cf. A253903 (characteristic function).
Cf. A034705 (differences of any two terms).

Programs

  • GAP
    List([0..30], n-> n*(n+1)*(2*n+1)/6); # G. C. Greubel, Dec 31 2019
  • Haskell
    a000330 n = n * (n + 1) * (2 * n + 1) `div` 6
    a000330_list = scanl1 (+) a000290_list
    -- Reinhard Zumkeller, Nov 11 2012, Feb 03 2012
    
  • Magma
    [n*(n+1)*(2*n+1)/6: n in [0..50]]; // Wesley Ivan Hurt, Jun 28 2014
    
  • Magma
    [0] cat [((2*n+3)*Binomial(n+2,2))/3: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014
    
  • Maple
    A000330 := n -> n*(n+1)*(2*n+1)/6;
    a := n->(1/6)*n*(n+1)*(2*n+1): seq(a(n),n=0..53); # Emeric Deutsch
    with(combstruct): ZL:=[st, {st=Prod(left, right), left=Set(U, card=r), right=Set(U, card=r), U=Sequence(Z, card>=1)}, unlabeled]: subs(r=1, stack): seq(count(subs(r=2, ZL), size=m*2), m=1..45) ; # Zerinvary Lajos, Jan 02 2008
    nmax := 44; for n from 0 to nmax do fz(n) := product( (1-(2*m-1)*z)^(n+1-m) , m=1..n); c(n) := abs(coeff(fz(n),z,1)); end do: a := n-> c(n): seq(a(n), n=0..nmax); # Johannes W. Meijer, Mar 07 2009
  • Mathematica
    Table[Binomial[w+2, 3] + Binomial[w+1, 3], {w, 0, 30}]
    CoefficientList[Series[x(1+x)/(1-x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 30 2014 *)
    Accumulate[Range[0,50]^2] (* Harvey P. Dale, Sep 25 2014 *)
  • Maxima
    A000330(n):=binomial(n+2,3)+binomial(n+1,3)$
    makelist(A000330(n),n,0,20); /* Martin Ettl, Nov 12 2012 */
    
  • PARI
    {a(n) = n * (n+1) * (2*n+1) / 6};
    
  • PARI
    upto(n) = [x*(x+1)*(2*x+1)/6 | x<-[0..n]] \\ Cino Hilliard, Jun 18 2007, edited by M. F. Hasler, Jan 02 2024
    
  • Python
    a=lambda n: (n*(n+1)*(2*n+1))//6 # Indranil Ghosh, Jan 04 2017
    
  • Sage
    [n*(n+1)*(2*n+1)/6 for n in (0..30)] # G. C. Greubel, Dec 31 2019
    

Formula

G.f.: x*(1+x)/(1-x)^4. - Simon Plouffe (in his 1992 dissertation: generating function for sequence starting at a(1))
E.g.f.: (x + 3*x^2/2 + x^3/3)*exp(x).
a(n) = n*(n+1)*(2*n+1)/6 = binomial(n+2, 3) + binomial(n+1, 3).
2*a(n) = A006331(n). - N. J. A. Sloane, Dec 11 1999
Can be extended to Z with a(n) = -a(-1-n) for all n in Z.
a(n) = A002492(n)/4. - Paul Barry, Jul 19 2003
a(n) = (((n+1)^4 - n^4) - ((n+1)^2 - n^2))/12. - Xavier Acloque, Oct 16 2003
From Alexander Adamchuk, Oct 26 2004: (Start)
a(n) = sqrt(A271535(n)).
a(n) = (Sum_{k=1..n} Sum_{j=1..n} Sum_{i=1..n} (i*j*k)^2)^(1/3). (End)
a(n) = Sum_{i=1..n} i*(2*n-2*i+1); sum of squares gives 1 + (1+3) + (1+3+5) + ... - Jon Perry, Dec 08 2004
a(n+1) = A000217(n+1) + 2*A000292(n). - Creighton Dement, Mar 10 2005
Sum_{n>=1} 1/a(n) = 6*(3-4*log(2)); Sum_{n>=1} (-1)^(n+1)*1/a(n) = 6*(Pi-3). - Philippe Deléham, May 31 2005
Sum of two consecutive tetrahedral (or pyramidal) numbers a(n) = A000292(n-1) + A000292(n). - Alexander Adamchuk, May 17 2006
Euler transform of length-2 sequence [ 5, -1 ]. - Michael Somos, Sep 04 2006
a(n) = a(n-1) + n^2. - Rolf Pleisch, Jul 22 2007
a(n) = A132121(n,0). - Reinhard Zumkeller, Aug 12 2007
a(n) = binomial(n, 2) + 2*binomial(n, 3). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009, corrected by M. F. Hasler, Jan 02 2024
a(n) = A168559(n) + 1 for n > 0. - Reinhard Zumkeller, Feb 03 2012
a(n) = Sum_{i=1..n} J_2(i)*floor(n/i), where J_2 is A007434. - Enrique Pérez Herrero, Feb 26 2012
a(n) = s(n+1, n)^2 - 2*s(n+1, n-1), where s(n, k) are Stirling numbers of the first kind, A048994. - Mircea Merca, Apr 03 2012
a(n) = A001477(n) + A000217(n) + A007290(n+2) + 1. - J. M. Bergot, May 31 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 2. - Ant King, Oct 17 2012
a(n) = Sum_{i = 1..n} Sum_{j = 1..n} min(i,j). - Enrique Pérez Herrero, Jan 15 2013
a(n) = A000217(n) + A007290(n+1). - Ivan N. Ianakiev, May 10 2013
a(n) = (A047486(n+2)^3 - A047486(n+2))/24. - Richard R. Forberg, Dec 25 2013
a(n) = Sum_{i=0..n-1} (n-i)*(2*i+1), with a(0) = 0. After 0, row sums of the triangle in A101447. - Bruno Berselli, Feb 10 2014
a(n) = n + 1 + Sum_{i=1..n+1} (i^2 - 2i). - Wesley Ivan Hurt, Feb 25 2014
a(n) = A000578(n+1) - A002412(n+1). - Wesley Ivan Hurt, Jun 28 2014
a(n) = Sum_{i = 1..n} Sum_{j = i..n} max(i,j). - Enrique Pérez Herrero, Dec 03 2014
a(n) = A055112(n)/6, see Singh (2013). - Alonso del Arte, Feb 20 2015
For n >= 2, a(n) = A028347(n+1) + A101986(n-2). - Bui Quang Tuan, Apr 03 2015
For n > 0: a(n) = A258708(n+3,n-1). - Reinhard Zumkeller, Jun 23 2015
a(n) = A175254(n) + A072481(n), n >= 1. - Omar E. Pol, Aug 12 2015
a(n) = A000332(n+3) - A000332(n+1). - Antal Pinter, Dec 27 2015
Dirichlet g.f.: zeta(s-3)/3 + zeta(s-2)/2 + zeta(s-1)/6. - Ilya Gutkovskiy, Jun 26 2016
a(n) = A080851(2,n-1). - R. J. Mathar, Jul 28 2016
a(n) = (A005408(n) * A046092(n))/12 = (2*n+1)*(2*n*(n+1))/12. - Bruce J. Nicholson, May 18 2017
12*a(n) = (n+1)*A001105(n) + n*A001105(n+1). - Bruno Berselli, Jul 03 2017
a(n) = binomial(n-1, 1) + binomial(n-1, 2) + binomial(n, 3) + binomial(n+1, 2) + binomial(n+1, 3). - Tony Foster III, Aug 24 2018
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Nathan Fox, Dec 04 2019
Let T(n) = A000217(n), the n-th triangular number. Then a(n) = (T(n)+1)^2 + (T(n)+2)^2 + ... + (T(n)+n)^2 - (n+2)*T(n)^2. - Charlie Marion, Dec 31 2019
a(n) = 2*n - 1 - a(n-2) + 2*a(n-1). - Boštjan Gec, Nov 09 2023
a(n) = 2/(2*n)! * Sum_{j = 1..n} (-1)^(n+j) * j^(2*n+2) * binomial(2*n, n-j). Cf. A060493. - Peter Bala, Mar 31 2025

Extensions

Partially edited by Joerg Arndt, Mar 11 2010

A001105 a(n) = 2*n^2.

Original entry on oeis.org

0, 2, 8, 18, 32, 50, 72, 98, 128, 162, 200, 242, 288, 338, 392, 450, 512, 578, 648, 722, 800, 882, 968, 1058, 1152, 1250, 1352, 1458, 1568, 1682, 1800, 1922, 2048, 2178, 2312, 2450, 2592, 2738, 2888, 3042, 3200, 3362, 3528, 3698, 3872, 4050, 4232, 4418
Offset: 0

Views

Author

Bernd.Walter(AT)frankfurt.netsurf.de

Keywords

Comments

Number of edges of the complete bipartite graph of order 3n, K_{n,2n}. - Roberto E. Martinez II, Jan 07 2002
"If each period in the periodic system ends in a rare gas ..., the number of elements in a period can be found from the ordinal number n of the period by the formula: L = ((2n+3+(-1)^n)^2)/8..." - Nature, Jun 09 1951; Nature 411 (Jun 07 2001), p. 648. This produces the present sequence doubled up.
Let z(1) = i = sqrt(-1), z(k+1) = 1/(z(k)+2i); then a(n) = (-1)*Imag(z(n+1))/Real(z(n+1)). - Benoit Cloitre, Aug 06 2002
Maximum number of electrons in an atomic shell with total quantum number n. Partial sums of A016825. - Jeremy Gardiner, Dec 19 2004
Arithmetic mean of triangular numbers in pairs: (1+3)/2, (6+10)/2, (15+21)/2, ... . - Amarnath Murthy, Aug 05 2005
These numbers form a pattern on the Ulam spiral similar to that of the triangular numbers. - G. Roda, Oct 20 2010
Integral areas of isosceles right triangles with rational legs (legs are 2n and triangles are nondegenerate for n > 0). - Rick L. Shepherd, Sep 29 2009
Even squares divided by 2. - Omar E. Pol, Aug 18 2011
Number of stars when distributed as in the U.S.A. flag: n rows with n+1 stars and, between each pair of these, one row with n stars (i.e., n-1 of these), i.e., n*(n+1)+(n-1)*n = 2*n^2 = A001105(n). - César Eliud Lozada, Sep 17 2012
Apparently the number of Dyck paths with semilength n+3 and an odd number of peaks and the central peak having height n-3. - David Scambler, Apr 29 2013
Sum of the partition parts of 2n into exactly two parts. - Wesley Ivan Hurt, Jun 01 2013
Consider primitive Pythagorean triangles (a^2 + b^2 = c^2, gcd(a, b) = 1) with hypotenuse c (A020882) and respective odd leg a (A180620); sequence gives values c-a, sorted with duplicates removed. - K. G. Stier, Nov 04 2013
Number of roots in the root systems of type B_n and C_n (for n > 1). - Tom Edgar, Nov 05 2013
Area of a square with diagonal 2n. - Wesley Ivan Hurt, Jun 18 2014
This sequence appears also as the first and second member of the quartet [a(n), a(n), p(n), p(n)] of the square of [n, n, n+1, n+1] in the Clifford algebra Cl_2 for n >= 0. p(n) = A046092(n). See an Oct 15 2014 comment on A147973 where also a reference is given. - Wolfdieter Lang, Oct 16 2014
a(n) are the only integers m where (A000005(m) + A000203(m)) = (number of divisors of m + sum of divisors of m) is an odd number. - Richard R. Forberg, Jan 09 2015
a(n) represents the first term in a sum of consecutive integers running to a(n+1)-1 that equals (2n+1)^3. - Patrick J. McNab, Dec 24 2016
Also the number of 3-cycles in the (n+4)-triangular honeycomb obtuse knight graph. - Eric W. Weisstein, Jul 29 2017
Also the Wiener index of the n-cocktail party graph for n > 1. - Eric W. Weisstein, Sep 07 2017
Numbers represented as the palindrome 242 in number base B including B=2 (binary), 3 (ternary) and 4: 242(2)=18, 242(3)=32, 242(4)=50, ... 242(9)=200, 242(10)=242, ... - Ron Knott, Nov 14 2017
a(n) is the square of the hypotenuse of an isosceles right triangle whose sides are equal to n. - Thomas M. Green, Aug 20 2019
The sequence contains all odd powers of 2 (A004171) but no even power of 2 (A000302). - Torlach Rush, Oct 10 2019
From Bernard Schott, Aug 31 2021 and Sep 16 2021: (Start)
Apart from 0, integers such that the number of even divisors (A183063) is odd.
Proof: every n = 2^q * (2k+1), q, k >= 0, then 2*n^2 = 2^(2q+1) * (2k+1)^2; now, gcd(2, 2k+1) = 1, tau(2^(2q+1)) = 2q+2 and tau((2k+1)^2) = 2u+1 because (2k+1)^2 is square, so, tau(2*n^2) = (2q+2) * (2u+1).
The 2q+2 divisors of 2^(2q+1) are {1, 2, 2^2, 2^3, ..., 2^(2q+1)}, so 2^(2q+1) has 2q+1 even divisors {2^1, 2^2, 2^3, ..., 2^(2q+1)}.
Conclusion: these 2q+1 even divisors create with the 2u+1 odd divisors of (2k+1)^2 exactly (2q+1)*(2u+1) even divisors of 2*n^2, and (2q+1)*(2u+1) is odd. (End)
a(n) with n>0 are the numbers with period length 2 for Bulgarian and Mancala solitaire. - Paul Weisenhorn, Jan 29 2022
Number of points at L1 distance = 2 from any given point in Z^n. - Shel Kaphan, Feb 25 2023
Integer that multiplies (h^2)/(m*L^2) to give the energy of a 1-D quantum mechanical particle in a box whenever it is an integer multiple of (h^2)/(m*L^2), where h = Planck's constant, m = mass of particle, and L = length of box. - A. Timothy Royappa, Mar 14 2025

Examples

			a(3) = 18; since 2(3) = 6 has 3 partitions with exactly two parts: (5,1), (4,2), (3,3).  Adding all the parts, we get: 1 + 2 + 3 + 3 + 4 + 5 = 18. - _Wesley Ivan Hurt_, Jun 01 2013
		

References

  • Peter Atkins, Julio De Paula, and James Keeler, "Atkins' Physical Chemistry," Oxford University Press, 2023, p. 31.
  • Arthur Beiser, Concepts of Modern Physics, 2nd Ed., McGraw-Hill, 1973.
  • Martin Gardner, The Colossal Book of Mathematics, Classic Puzzles, Paradoxes and Problems, Chapter 2 entitled "The Calculus of Finite Differences," W. W. Norton and Company, New York, 2001, pages 12-13.
  • L. B. W. Jolley, "Summation of Series", Dover Publications, 1961, p. 44.
  • Alain M. Robert, A Course in p-adic Analysis, Springer-Verlag, 2000, p. 213.

Crossrefs

Cf. numbers of the form n*(n*k-k+4)/2 listed in A226488.
Cf. A058331 and A247375. - Bruno Berselli, Sep 16 2014
Cf. A194715 (4-cycles in the triangular honeycomb obtuse knight graph), A290391 (5-cycles), A290392 (6-cycles). - Eric W. Weisstein, Jul 29 2017
Integers such that: this sequence (the number of even divisors is odd), A028982 (the number of odd divisors is odd), A028983 (the number of odd divisors is even), A183300 (the number of even divisors is even).

Programs

Formula

a(n) = (-1)^(n+1) * A053120(2*n, 2).
G.f.: 2*x*(1+x)/(1-x)^3.
a(n) = A100345(n, n).
Sum_{n>=1} 1/a(n) = Pi^2/12 =A072691. [Jolley eq. 319]. - Gary W. Adamson, Dec 21 2006
a(n) = A049452(n) - A033991(n). - Zerinvary Lajos, Jun 12 2007
a(n) = A016742(n)/2. - Zerinvary Lajos, Jun 20 2008
a(n) = 2 * A000290(n). - Omar E. Pol, May 14 2008
a(n) = 4*n + a(n-1) - 2, n > 0. - Vincenzo Librandi
a(n) = A002378(n-1) + A002378(n). - Joerg M. Schuetze (joerg(AT)cyberheim.de), Mar 08 2010 [Corrected by Klaus Purath, Jun 18 2020]
a(n) = A176271(n,k) + A176271(n,n-k+1), 1 <= k <= n. - Reinhard Zumkeller, Apr 13 2010
a(n) = A007607(A000290(n)). - Reinhard Zumkeller, Feb 12 2011
For n > 0, a(n) = 1/coefficient of x^2 in the Maclaurin expansion of 1/(cos(x)+n-1). - Francesco Daddi, Aug 04 2011
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Artur Jasinski, Nov 24 2011
a(n) = A070216(n,n) for n > 0. - Reinhard Zumkeller, Nov 11 2012
a(n) = A014132(2*n-1,n) for n > 0. - Reinhard Zumkeller, Dec 12 2012
a(n) = A000217(n) + A000326(n). - Omar E. Pol, Jan 11 2013
(a(n) - A000217(k))^2 = A000217(2*n-1-k)*A000217(2*n+k) + n^2, for all k. - Charlie Marion, May 04 2013
a(n) = floor(1/(1-cos(1/n))), n > 0. - Clark Kimberling, Oct 08 2014
a(n) = A251599(3*n-1) for n > 0. - Reinhard Zumkeller, Dec 13 2014
a(n) = Sum_{j=1..n} Sum_{i=1..n} ceiling((i+j-n+4)/3). - Wesley Ivan Hurt, Mar 12 2015
a(n) = A002061(n+1) + A165900(n). - Torlach Rush, Feb 21 2019
E.g.f.: 2*exp(x)*x*(1 + x). - Stefano Spezia, Oct 12 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/24 (A222171). - Amiram Eldar, Jul 03 2020
From Amiram Eldar, Feb 03 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = sqrt(2)*sinh(Pi/sqrt(2))/Pi.
Product_{n>=1} (1 - 1/a(n)) = sqrt(2)*sin(Pi/sqrt(2))/Pi. (End)

A005902 Centered icosahedral (or cuboctahedral) numbers, also crystal ball sequence for f.c.c. lattice.

Original entry on oeis.org

1, 13, 55, 147, 309, 561, 923, 1415, 2057, 2869, 3871, 5083, 6525, 8217, 10179, 12431, 14993, 17885, 21127, 24739, 28741, 33153, 37995, 43287, 49049, 55301, 62063, 69355, 77197, 85609, 94611, 104223, 114465, 125357, 136919, 149171, 162133, 175825, 190267, 205479
Offset: 0

Views

Author

Keywords

Comments

Called "magic numbers" in some chemical contexts.
Partial sums of A005901(n). - Lekraj Beedassy, Oct 30 2003
Equals binomial transform of [1, 12, 30, 20, 0, 0, 0, ...]. - Gary W. Adamson, Aug 01 2008
Crystal ball sequence for A_3 lattice. - Michael Somos, Jun 03 2012

Examples

			a(4) = 147 = (1, 3, 3, 1) dot (1, 12, 30, 20) = (1 + 36 + 90 + 20). - _Gary W. Adamson_, Aug 01 2008
G.f. = 1 + 13*x + 55*x^2 + 147*x^3 + 309*x^4 + 561*x^5 + 923*x^6 + 1415*x^7 + ...
		

References

  • H. S. M. Coxeter, Polyhedral numbers, pp. 25-35 of R. S. Cohen, J. J. Stachel and M. W. Wartofsky, eds., For Dirk Struik: Scientific, historical and political essays in honor of Dirk J. Struik, Reidel, Dordrecht, 1974.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

(1/12)*t*(2*n^3-3*n^2+n)+2*n-1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.
The 28 uniform 3D tilings: cab: A299266, A299267; crs: A299268, A299269; fcu: A005901, A005902; fee: A299259, A299265; flu-e: A299272, A299273; fst: A299258, A299264; hal: A299274, A299275; hcp: A007899, A007202; hex: A005897, A005898; kag: A299256, A299262; lta: A008137, A299276; pcu: A005899, A001845; pcu-i: A299277, A299278; reo: A299279, A299280; reo-e: A299281, A299282; rho: A008137, A299276; sod: A005893, A005894; sve: A299255, A299261; svh: A299283, A299284; svj: A299254, A299260; svk: A010001, A063489; tca: A299285, A299286; tcd: A299287, A299288; tfs: A005899, A001845; tsi: A299289, A299290; ttw: A299257, A299263; ubt: A299291, A299292; bnn: A007899, A007202. See the Proserpio link in A299266 for overview.

Programs

  • Magma
    [(2*n+1)*(5*n^2+5*n+3)/3: n in [0..30]]; // G. C. Greubel, Dec 01 2017
    
  • Maple
    A005902 := n -> (2*n+1)*(5*n^2+5*n+3)/3;
    A005902:=(z+1)*(z**2+8*z+1)/(z-1)**4; # Simon Plouffe in his 1992 dissertation
  • Mathematica
    f[n_] := (2n + 1)(5n^2 + 5n + 3)/3; Array[f, 36, 0] (* Robert G. Wilson v, Feb 02 2011 *)
    LinearRecurrence[{4,-6,4,-1},{1,13,55,147},50] (* Harvey P. Dale, Oct 08 2015 *)
    CoefficientList[Series[(x^3 + 9*x^2 + 9*x + 1)/(x - 1)^4, {x, 0, 50}], x] (* Indranil Ghosh, Apr 08 2017 *)
  • PARI
    {a(n) = (2*n + 1) * (5*n^2 + 5*n + 3) / 3}; /* Michael Somos, Jun 03 2012 */
    
  • PARI
    x='x+O('x^50); Vec((x^3 + 9*x^2 + 9*x + 1)/(x - 1)^4) \\ Indranil Ghosh, Apr 08 2017
    
  • Python
    def a(n): return (2*n+1)*(5*n**2+5*n+3)//3
    print([a(n) for n in range(40)]) # Michael S. Branicky, Jan 13 2021

Formula

a(n) = (2*n+1)*(5*n^2+5*n+3)/3.
For n > 0, n*a(n) = (Sum_{i=0..n-1} a(i)) + 2*A005891(n)*A000217(n). - Bruno Berselli, Feb 02 2011
a(-1 - n) = -a(n). - Michael Somos, Jun 03 2012
From Indranil Ghosh, Apr 08 2017: (Start)
G.f.: (x^3 + 9x^2 + 9x + 1)/(x - 1)^4.
E.g.f.: (1/3)*exp(x)*(10x^3 + 45x^2 + 36x + 3).
(End)
a(n) = A100171(n+1) - A008778(n-1) = A100174(n+1) - A000290(n) = A005917(n+1) - A006331(n) = A051673(n+1) + A000578(n). - Bruce J. Nicholson, Jul 05 2018

A212959 Number of (w,x,y) such that w,x,y are all in {0,...,n} and |w-x| = |x-y|.

Original entry on oeis.org

1, 4, 11, 20, 33, 48, 67, 88, 113, 140, 171, 204, 241, 280, 323, 368, 417, 468, 523, 580, 641, 704, 771, 840, 913, 988, 1067, 1148, 1233, 1320, 1411, 1504, 1601, 1700, 1803, 1908, 2017, 2128, 2243, 2360, 2481, 2604, 2731, 2860, 2993, 3128, 3267
Offset: 0

Views

Author

Clark Kimberling, Jun 01 2012

Keywords

Comments

In the following guide to related sequences: M=max(x,y,z), m=min(x,y,z), and R=range=M-m. In some cases, it is an offset of the listed sequence which fits the conditions shown for w,x,y. Each sequence satisfies a linear recurrence relation, some of which are identified in the list by the following code (signature):
A: 2, 0, -2, 1, i.e., a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4);
B: 3, -2, -2, 3, -1;
C: 4, -6, 4, -1;
D: 1, 2, -2, -1, 1;
E: 2, 1, -4, 1, 2, -1;
F: 2, -1, 1, -2, 1;
G: 2, -1, 0, 1, -2, 1;
H: 2, -1, 2, -4, 2, -1, 2, -1;
I: 3, -3, 2, -3, 3, -1;
J: 4, -7, 8, -7, 4, -1.
...
A212959 ... |w-x|=|x-y| ...... recurrence type A
A212960 ... |w-x| != |x-y| ................... B
A212683 ... |w-x| < |x-y| .................... B
A212684 ... |w-x| >= |x-y| ................... B
A212963 ... see entry for definition ......... B
A212964 ... |w-x| < |x-y| < |y-w| ............ B
A006331 ... |w-x| < y ........................ C
A005900 ... |w-x| <= y ....................... C
A212965 ... w = R ............................ D
A212966 ... 2*w = R
A212967 ... w < R ............................ E
A212968 ... w >= R ........................... E
A077043 ... w = x > R ........................ A
A212969 ... w != x and x > R ................. E
A212970 ... w != x and x < R ................. E
A055998 ... w = x + y - 1
A011934 ... w < floor((x+y)/2) ............... B
A182260 ... w > floor((x+y)/2) ............... B
A055232 ... w <= floor((x+y)/2) .............. B
A011934 ... w >= floor((x+y)/2) .............. B
A212971 ... w < floor((x+y)/3) ............... B
A212972 ... w >= floor((x+y)/3) .............. B
A212973 ... w <= floor((x+y)/3) .............. B
A212974 ... w > floor((x+y)/3) ............... B
A212975 ... R is even ........................ E
A212976 ... R is odd ......................... E
A212978 ... R = 2*n - w - x
A212979 ... R = average{w,x,y}
A212980 ... w < x + y and x < y .............. B
A212981 ... w <= x+y and x < y ............... B
A212982 ... w < x + y and x <= y ............. B
A212983 ... w <= x + y and x <= y ............ B
A002623 ... w >= x + y and x <= y ............ B
A087811 ... w = 2*x + y ...................... A
A008805 ... w = 2*x + 2*y .................... D
A000982 ... 2*w = x + y ...................... F
A001318 ... 2*w = 2*x + y .................... F
A001840 ... w = 3*x + y
A212984 ... 3*w = x + y
A212985 ... 3*w = 3*x + y
A001399 ... w = 2*x + 3*y
A212986 ... 2*w = 3*x + y
A008810 ... 3*x = 2*x + y .................... F
A212987 ... 3*w = 2*x + 2*y
A001972 ... w = 4*x + y ...................... G
A212988 ... 4*w = x + y ...................... G
A212989 ... 4*w = 4*x + y
A008812 ... 5*w = 2*x + 3*y
A016061 ... n < w + x + y <= 2*n ............. C
A000292 ... w + x + y <=n .................... C
A000292 ... 2*n < w + x + y <= 3*n ........... C
A212977 ... n/2 < w + x + y <= n
A143785 ... w < R < x ........................ E
A005996 ... w < R <= x ....................... E
A128624 ... w <= R <= x ...................... E
A213041 ... R = 2*|w - x| .................... A
A213045 ... R < 2*|w - x| .................... B
A087035 ... R >= 2*|w - x| ................... B
A213388 ... R <= 2*|w - x| ................... B
A171218 ... M < 2*m .......................... B
A213389 ... R < 2|w - x| ..................... E
A213390 ... M >= 2*m ......................... E
A213391 ... 2*M < 3*m ........................ H
A213392 ... 2*M >= 3*m ....................... H
A213393 ... 2*M > 3*m ........................ H
A213391 ... 2*M <= 3*m ....................... H
A047838 ... w = |x + y - w| .................. A
A213396 ... 2*w < |x + y - w| ................ I
A213397 ... 2*w >= |x + y - w| ............... I
A213400 ... w < R < 2*w
A069894 ... min(|w-x|,|x-y|) = 1
A000384 ... max(|w-x|,|x-y|) = |w-y|
A213395 ... max(|w-x|,|x-y|) = w
A213398 ... min(|w-x|,|x-y|) = x ............. A
A213399 ... max(|w-x|,|x-y|) = x ............. D
A213479 ... max(|w-x|,|x-y|) = w+x+y ......... D
A213480 ... max(|w-x|,|x-y|) != w+x+y ........ E
A006918 ... |w-x| + |x-y| > w+x+y ............ E
A213481 ... |w-x| + |x-y| <= w+x+y ........... E
A213482 ... |w-x| + |x-y| < w+x+y ............ E
A213483 ... |w-x| + |x-y| >= w+x+y ........... E
A213484 ... |w-x|+|x-y|+|y-w| = w+x+y
A213485 ... |w-x|+|x-y|+|y-w| != w+x+y ....... J
A213486 ... |w-x|+|x-y|+|y-w| > w+x+y ........ J
A213487 ... |w-x|+|x-y|+|y-w| >= w+x+y ....... J
A213488 ... |w-x|+|x-y|+|y-w| < w+x+y ........ J
A213489 ... |w-x|+|x-y|+|y-w| <= w+x+y ....... J
A213490 ... w,x,y,|w-x|,|x-y| distinct
A213491 ... w,x,y,|w-x|,|x-y| not distinct
A213493 ... w,x,y,|w-x|,|x-y|,|w-y| distinct
A213495 ... w = min(|w-x|,|x-y|,|w-y|)
A213492 ... w != min(|w-x|,|x-y|,|w-y|)
A213496 ... x != max(|w-x|,|x-y|)
A213498 ... w != max(|w-x|,|x-y|,|w-y|)
A213497 ... w = min(|w-x|,|x-y|)
A213499 ... w != min(|w-x|,|x-y|)
A213501 ... w != max(|w-x|,|x-y|)
A213502 ... x != min(|w-x|,|x-y|)
...
A211795 includes a guide for sequences that count 4-tuples (w,x,y,z) having all terms in {0,...,n} and satisfying selected properties. Some of the sequences indexed at A211795 satisfy recurrences that are represented in the above list.
Partial sums of the numbers congruent to {1,3} mod 6 (see A047241). - Philippe Deléham, Mar 16 2014

Examples

			a(1)=4 counts these (x,y,z): (0,0,0), (1,1,1), (0,1,0), (1,0,1).
Numbers congruent to {1, 3} mod 6: 1, 3, 7, 9, 13, 15, 19, ...
a(0) = 1;
a(1) = 1 + 3 = 4;
a(2) = 1 + 3 + 7 = 11;
a(3) = 1 + 3 + 7 + 9 = 20;
a(4) = 1 + 3 + 7 + 9 + 13 = 33;
a(5) = 1 + 3 + 7 + 9 + 13 + 15 = 48; etc. - _Philippe Deléham_, Mar 16 2014
		

References

  • A. Barvinok, Lattice Points and Lattice Polytopes, Chapter 7 in Handbook of Discrete and Computational Geometry, CRC Press, 1997, 133-152.
  • P. Gritzmann and J. M. Wills, Lattice Points, Chapter 3.2 in Handbook of Convex Geometry, vol. B, North-Holland, 1993, 765-797.

Crossrefs

Programs

  • Mathematica
    t = Compile[{{n, _Integer}}, Module[{s = 0},
    (Do[If[Abs[w - x] == Abs[x - y], s = s + 1],
    {w, 0, n}, {x, 0, n}, {y, 0, n}]; s)]];
    m = Map[t[#] &, Range[0, 50]]   (* A212959 *)
  • PARI
    a(n)=(6*n^2+8*n+3)\/4 \\ Charles R Greathouse IV, Jul 28 2015

Formula

a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4).
G.f.: (1+2*x+3*x^2)/((1+x)*(1-x)^3).
a(n) + A212960(n) = (n+1)^3.
a(n) = (6*n^2 + 8*n + 3 + (-1)^n)/4. - Luce ETIENNE, Apr 05 2014
a(n) = 2*A069905(3*(n+1)+2) - 3*(n+1). - Ayoub Saber Rguez, Aug 31 2021

A002492 Sum of the first n even squares: 2*n*(n+1)*(2*n+1)/3.

Original entry on oeis.org

0, 4, 20, 56, 120, 220, 364, 560, 816, 1140, 1540, 2024, 2600, 3276, 4060, 4960, 5984, 7140, 8436, 9880, 11480, 13244, 15180, 17296, 19600, 22100, 24804, 27720, 30856, 34220, 37820, 41664, 45760, 50116, 54740, 59640, 64824, 70300, 76076, 82160
Offset: 0

Views

Author

Keywords

Comments

Total number of possible bishop moves on an n+1 X n+1 chessboard, if the bishop is placed anywhere. E.g., on a 3 X 3-Board: bishop has 8 X 2 moves and 1 X 4 moves, so a(2)=20. - Ulrich Schimke (ulrschimke(AT)aol.com)
Let M_n denote the n X n matrix M_n(i,j)=(i+j)^2; then the characteristic polynomial of M_n is x^n - a(n)x^(n-1) - .... - Michael Somos, Nov 14 2002
Partial sums of A016742. - Lekraj Beedassy, Jun 19 2004
0,4,20,56,120 gives the number of electrons in closed shells in the double shell periodic system of elements. This is a new interpretation of the periodic system of the elements. The factor 4 in the formula 4*n(n+1)(2n+1)/6 plays a significant role, since it designates the degeneracy of electronic states in this system. Closed shells with more than 120 electrons are not expected to exist. - Karl-Dietrich Neubert (kdn(AT)neubert.net)
Inverse binomial transform of A240434. - Wesley Ivan Hurt, Apr 13 2014
Atomic number of alkaline-earth metals of period 2n. - Natan Arie Consigli, Jul 03 2016
a(n) are the negative cubic coefficients in the expansion of sin(kx) into powers of sin(x) for the odd k: sin(kx) = k sin(x) - c(k) sin^3(x) + O(sin^5(x)); a(n) = c(2n+1) = A000292(2n). - Mathias Zechmeister, Jul 24 2022
Also the number of distinct series-parallel networks under series-parallel reduction on three unlabeled edges of n element kinds. - Michael R. Hayashi, Aug 02 2023

References

  • A. O. Barut, Group Structure of the Periodic System, in Wybourne, Ed., The Structure of Matter, University of Canterbury Press, Christchurch, 1972, p. 126.
  • Edward G. Mazur, Graphic Representation of the Periodic System during One Hundred Years, University of Alabama Press, Alabama, 1974.
  • W. Permans and J. Kemperman, "Nummeringspribleem van S. Dockx, Mathematisch Centrum. Amsterdam," Rapport ZW; 1949-005, 4 leaves, 19.8 X 34 cm.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A033586 (King), A035005 (Queen), A035006 (Rook), A035008 (Knight) and A049450 (Pawn).

Programs

  • Magma
    [2*n*(n+1)*(2*n+1)/3: n in [0..40]]; // Vincenzo Librandi, Jun 16 2011
  • Maple
    A002492:=n->2*n*(n+1)*(2*n+1)/3; seq(A002492(n), n=0..50); # Wesley Ivan Hurt, Apr 04 2014
  • Mathematica
    Table[2n(n+1)(2n+1)/3, {n,0,40}] (* or *) Binomial[2*Range[0,40]+2,3] (* or *) LinearRecurrence[{4,-6,4,-1}, {0,4,20,56},40] (* Harvey P. Dale, Aug 15 2012 *)
    Accumulate[(2*Range[0,40])^2] (* Harvey P. Dale, Jun 04 2019 *)
  • PARI
    a(n)=2*n*(n+1)*(2*n+1)/3
    

Formula

G.f.: 4*x*(1+x)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
a(-1-n) = -a(n).
a(n) = 4*A000330(n) = 2*A006331(n) = A000292(2*n).
a(n) = (-1)^(n+1)*A053120(2*n+1,3) (fourth unsigned column of Chebyshev T-triangle, zeros omitted).
a(n) = binomial(2*n+2, 3). - Lekraj Beedassy, Jun 19 2004
A035005(n+1) = a(n) + A035006(n+1) since Queen = Bishop + Rook. - Johannes W. Meijer, Feb 04 2010
a(n) - a(n-1) = 4*n^2. - Joerg Arndt, Jun 16 2011
a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4) for n>3. - Harvey P. Dale, Aug 15 2012
a(n) = Sum_{k=0..3} C(n-2+k,n-2)*C(n+3-k,n), for n>2. - J. M. Bergot, Jun 14 2014
a(n) = 2*A006331(n). - R. J. Mathar, May 28 2016
From Natan Arie Consigli Jul 03 2016: (Start)
a(n) = A166464(n) - 1.
a(n) = A168380(2*n). (End)
a(n) = Sum_{i=0..n} A005408(i)*A005408(i-1)+1 with A005408(-1):=-1. - Bruno Berselli, Jan 09 2017
a(n) = A002412(n) + A016061(n). - Bruce J. Nicholson, Nov 12 2017
From Amiram Eldar, Jan 04 2022: (Start)
Sum_{n>=1} 1/a(n) = 9/2 - 6*log(2).
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*Pi/2 - 9/2. (End)
a(n) = A081277(3, n-1) = (1+2*n)*binomial(n+1, n-2)*2^2/(n-1) for n > 0. - Mathias Zechmeister, Jul 26 2022
E.g.f.: 2*exp(x)*x*(6 + 9*x + 2*x^2)/3. - Stefano Spezia, Jul 31 2022

Extensions

Minor errors corrected and edited by Johannes W. Meijer, Feb 04 2010
Title modified by Charles R Greathouse IV at the suggestion of J. M. Bergot, Apr 05 2014

A059270 a(n) is both the sum of n+1 consecutive integers and the sum of the n immediately higher consecutive integers.

Original entry on oeis.org

0, 3, 15, 42, 90, 165, 273, 420, 612, 855, 1155, 1518, 1950, 2457, 3045, 3720, 4488, 5355, 6327, 7410, 8610, 9933, 11385, 12972, 14700, 16575, 18603, 20790, 23142, 25665, 28365, 31248, 34320, 37587, 41055, 44730, 48618, 52725, 57057, 61620
Offset: 0

Views

Author

Henry Bottomley, Jan 24 2001

Keywords

Comments

Group the non-multiples of n as follows, e.g., for n = 4: (1,2,3), (5,6,7), (9,10,11), (13,14,15), ... Then a(n) is the sum of the members of the n-th group. Or, the sum of (n-1)successive numbers preceding n^2. - Amarnath Murthy, Jan 19 2004
Convolution of odds (A005408) and multiples of three (A008585). G.f. is the product of the g.f. of A005408 by the g.f. of A008585. - Graeme McRae, Jun 06 2006
Sums of rows of the triangle in A126890. - Reinhard Zumkeller, Dec 30 2006
Corresponds to the Wiener indices of C_{2n+1} i.e., the cycle on 2n+1 vertices (n > 0). - K.V.Iyer, Mar 16 2009
Also the product of the three numbers from A005843(n) up to A163300(n), divided by 8. - Juri-Stepan Gerasimov, Jul 26 2009
Partial sums of A033428. - Charlie Marion, Dec 08 2013
For n > 0, sum of multiples of n and (n+1) from 1 to n*(n+1). - Zak Seidov, Aug 07 2016
A generalization of Ianakiev's formula, a(n) = A005408(n)*A000217(n), follows. A005408(n+k)*A000217(n) is the sum of n+1 consecutive integers and, after skipping k integers, the sum of the n immediately higher consecutive integers. For example, for n = 3 and k = 2, 9*6 = 54 = 12+13+14+15 = 17+18+19. - Charlie Marion, Jan 25 2022

Examples

			a(5) = 25 + 26 + 27 + 28 + 29 + 30 = 31 + 32 + 33 + 34 + 35 = 165.
		

Crossrefs

Cf. A059255 for analog for sum of squares.
Cf. A222716 for the analogous sum of triangular numbers.
Cf. A234319 for nonexistence of analogs for sums of n-th powers, n > 2. - Jonathan Sondow, Apr 23 2014
Cf. A098737 (first subdiagonal).
Bisection of A109900.

Programs

  • Magma
    I:=[0, 3, 15, 42]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // Vincenzo Librandi, Jun 23 2012
    
  • Maple
    A059270 := proc(n) n*(n+1)*(2*n+1)/2 ; end proc: # R. J. Mathar, Jul 10 2011
  • Mathematica
    # (#+1)(2#+1)/2 &/@ Range[0,39] (* Ant King, Jan 03 2011 *)
    CoefficientList[Series[3 x (1 + x)/(x - 1)^4, {x, 0, 39}], x]
    LinearRecurrence[{4,-6,4,-1},{0,3,15,42},50] (* Vincenzo Librandi, Jun 23 2012 *)
  • PARI
    a(n) = n*(n+1)*(2*n+1)/2 \\ Charles R Greathouse IV, Mar 08 2013
  • Sage
    [bernoulli_polynomial(n+1,3) for n in range(0, 41)] # Zerinvary Lajos, May 17 2009
    

Formula

a(n) = n*(n+1)*(2*n+1)/2.
a(n) = A000330(n)*3 = A006331(n)*3/2 = A055112(n)/2 = A000217(A002378(n)) - A000217(A005563(n-1)) = A000217(A005563(n)) - A000217(A002378(n)).
a(n) = A110449(n+1, n-1) for n > 1.
a(n) = Sum_{k=A000290(n) .. A002378(n)} k = Sum_{k=n^2..n^2+n} k.
a(n) = Sum_{k=n^2+n+1 .. n^2+2*n} k = Sum_{k=A002061(n+1) .. A005563(n)} k.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 6 = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Ant King, Jan 03 2011
G.f.: 3*x*(1+x)/(1-x)^4. - Ant King, Jan 03 2011
a(n) = A000578(n+1) - A000326(n+1). - Ivan N. Ianakiev, Nov 29 2012
a(n) = A005408(n)*A000217(n) = a(n-1) + 3*A000290(n). -Ivan N. Ianakiev, Mar 08 2013
a(n) = n^3 + n^2 + A000217(n). - Charlie Marion, Dec 04 2013
From Ilya Gutkovskiy, Aug 08 2016: (Start)
E.g.f.: x*(6 + 9*x + 2*x^2)*exp(x)/2.
Sum_{n>=1} 1/a(n) = 2*(3 - 4*log(2)) = 0.4548225555204375246621... (End)
a(n) = Sum_{k=0..2*n} A001318(k). - Jacob Szlachetka, Dec 20 2021
a(n) = Sum_{k=0..n} A000326(k) + A005449(k). - Jacob Szlachetka, Dec 21 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(Pi-3). - Amiram Eldar, Sep 17 2022

A036464 Number of ways to place two nonattacking queens on an n X n board.

Original entry on oeis.org

0, 0, 8, 44, 140, 340, 700, 1288, 2184, 3480, 5280, 7700, 10868, 14924, 20020, 26320, 34000, 43248, 54264, 67260, 82460, 100100, 120428, 143704, 170200, 200200, 234000, 271908, 314244, 361340, 413540, 471200, 534688, 604384
Offset: 1

Views

Author

Robert G. Wilson v, Raymond Bush (c17h21no4(AT)hotmail.com), Kirk Conely, N. J. A. Sloane

Keywords

Crossrefs

Column k=2 of A348129.

Programs

  • Maple
    f:=n->n^4/2 - 5*n^3/3 + 3*n^2/2 - n/3; [seq(f(n),n=1..200)]; # N. J. A. Sloane, Feb 16 2013
  • Mathematica
    f[k_] := 2 k; t[n_] := Table[f[k], {k, 1, n}]
    a[n_] := SymmetricPolynomial[2, t[n]]
    Table[a[n], {n, 2, 50}]   (* A036464 *)
    Table[a[n]/4, {n, 2, 50}] (* A000914 *)
    (* Clark Kimberling, Dec 31 2011 *)
    CoefficientList[Series[4 x^2 (2 + x) / (1-x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, May 02 2013 *)
    LinearRecurrence[{5,-10,10,-5,1},{0,0,8,44,140},50] (* Harvey P. Dale, Mar 26 2015 *)

Formula

a(n) = C(n, 3)*(3*n-1).
G.f.: 4*x^3*(2+x)/(1-x)^5. - Colin Barker, May 02 2012
a(n) = 2*sum_{i=1..n-2} i(i + 1)^2. - Wesley Ivan Hurt, Mar 18 2014
E.g.f.: (exp(x) * x^3 * (8 + 3*x))/6. - Vaclav Kotesovec, Feb 15 2015
For n>0, a(n) = A163102(n-1) - A006331(n-1). - Antal Pinter, Sep 20 2015

A323540 a(n) = Product_{k=0..n} (k^2 + (n-k)^2).

Original entry on oeis.org

0, 1, 32, 2025, 204800, 30525625, 6307891200, 1727713080625, 606076928000000, 265058191985900625, 141409376995328000000, 90403125002859606705625, 68229510086445571768320000, 60026603304487418050791015625, 60893916244529680380723200000000
Offset: 0

Views

Author

Vaclav Kotesovec, Jan 17 2019

Keywords

Crossrefs

Cf. 2*A000330 and A006331 (with sum instead of product).

Programs

  • Magma
    m:=2; [(&*[k^m + (n-k)^m: k in [0..n]]): n in [0..20]]; // G. C. Greubel, Jan 18 2019
    
  • Mathematica
    Table[Product[k^2+(n-k)^2, {k, 0, n}], {n, 0, 20}]
  • PARI
    m=2; vector(20, n, n--; prod(k=0,n, k^m + (n-k)^m)) \\ G. C. Greubel, Jan 18 2019
    
  • Sage
    m=2; [product(k^m +(n-k)^m for k in (0..n)) for n in (0..20)] # G. C. Greubel, Jan 18 2019

Formula

a(n) ~ exp((Pi-4)*n/2) * n^(2*n+2).
Showing 1-10 of 47 results. Next