[go: up one dir, main page]

login
Search: a115067 -id:a115067
     Sort: relevance | references | number | modified | created      Format: long | short | data
a(n) = 2^(n*(3n+5)/2) = 2^A115067(n+1).
+20
0
1, 16, 2048, 2097152, 17179869184, 1125899906842624, 590295810358705651712, 2475880078570760549798248448, 83076749736557242056487941267521536, 22300745198530623141535718272648361505980416
OFFSET
0,2
COMMENTS
Apparently this sequence is the Hankel transform of A098329 = [1, 1, 17, 49, 481, 2081, ...]; example : Det([1]) = 1, Det([1,1 ; 1,17]) = 16, Det([1,1,17 ; 1,17,49 ; 17,49,481]) = 2048, ...
MATHEMATICA
Table[2^(n(3n+5)/2), {n, 0, 9}] (* James C. McMahon, Oct 17 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Philippe Deléham, Sep 06 2006
STATUS
approved
Complement of A115067.
+20
0
1, 2, 3, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, 18, 19, 20, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86
OFFSET
1,2
FORMULA
(See the Mathematica code.)
MATHEMATICA
a=3/2; b=5/2;
F[n_]:=a*n^2+b*n;
R[n_]:=(n/a+((b-1)/(2a))^2)^(1/2);
G[n_]:=n-1+Ceiling[R[n]-(b-1)/(2a)];
Table[F[n], {n, 60}]
Table[G[n], {n, 100}]
CROSSREFS
Cf. A115067.
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jan 03 2011
STATUS
approved
The nonnegative integers.
+10
824
0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77
OFFSET
0,3
COMMENTS
Although this is a list, and lists normally have offset 1, it seems better to make an exception in this case. - N. J. A. Sloane, Mar 13 2010
The subsequence 0,1,2,3,4 gives the known values of n such that 2^(2^n)+1 is a prime (see A019434, the Fermat primes). - N. J. A. Sloane, Jun 16 2010
Also: The identity map, defined on the set of nonnegative integers. The restriction to the positive integers yields the sequence A000027. - M. F. Hasler, Nov 20 2013
The number of partitions of 2n into exactly 2 parts. - Colin Barker, Mar 22 2015
The number of orbits of Aut(Z^7) as function of the infinity norm n of the representative lattice point of the orbit, when the cardinality of the orbit is equal to 8960 or 168.- Philippe A.J.G. Chevalier, Dec 29 2015
Partial sums give A000217. - Omar E. Pol, Jul 26 2018
First differences are A000012 (the "all 1's" sequence). - M. F. Hasler, May 30 2020
See A061579 for the transposed infinite square matrix, or triangle with rows reversed. - M. F. Hasler, Nov 09 2021
This is the unique sequence (a(n)) that satisfies the inequality a(n+1) > a(a(n)) for all n in N. This simple and surprising result comes from the 6th problem proposed by Bulgaria during the second day of the 19th IMO (1977) in Belgrade (see link and reference). - Bernard Schott, Jan 25 2023
REFERENCES
Maurice Protat, Des Olympiades à l'Agrégation, suite vérifiant f(n+1) > f(f(n)), Problème 7, pp. 31-32, Ellipses, Paris 1997.
LINKS
Paul Barry, A Catalan Transform and Related Transformations on Integer Sequences, Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.5.
The IMO Compendium, Problem 6, 19th IMO 1977.
Tanya Khovanova, Recursive Sequences
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014-2015.
László Németh, The trinomial transform triangle, J. Int. Seqs., Vol. 21 (2018), Article 18.7.3. Also arXiv:1807.07109 [math.NT], 2018.
N. J. A. Sloane, "A Handbook of Integer Sequences" Fifty Years Later, arXiv:2301.03149 [math.NT], 2023, p. 12.
Eric Weisstein's World of Mathematics, Natural Number
Eric Weisstein's World of Mathematics, Nonnegative Integer
FORMULA
a(n) = n.
a(0) = 0, a(n) = a(n-1) + 1.
G.f.: x/(1-x)^2.
Multiplicative with a(p^e) = p^e. - David W. Wilson, Aug 01 2001
When seen as array: T(k, n) = n + (k+n)*(k+n+1)/2. Main diagonal is 2*n*(n+1) (A046092), antidiagonal sums are n*(n+1)*(n+2)/2 (A027480). - Ralf Stephan, Oct 17 2004
Dirichlet generating function: zeta(s-1). - Franklin T. Adams-Watters, Sep 11 2005
E.g.f.: x*e^x. - Franklin T. Adams-Watters, Sep 11 2005
a(0)=0, a(1)=1, a(n) = 2*a(n-1) - a(n-2). - Jaume Oliver Lafont, May 07 2008
Alternating partial sums give A001057 = A000217 - 2*(A008794). - Eric Desbiaux, Oct 28 2008
a(n) = 2*A080425(n) + 3*A008611(n-3), n>1. - Eric Desbiaux, Nov 15 2009
a(n) = A007966(n)*A007967(n). - Reinhard Zumkeller, Jun 18 2011
a(n) = Sum_{k>=0} A030308(n,k)*2^k. - Philippe Deléham, Oct 20 2011
a(n) = 2*A028242(n-1) + (-1)^n*A000034(n-1). - R. J. Mathar, Jul 20 2012
a(n+1) = det(C(i+1,j), 1 <= i, j <= n), where C(n,k) are binomial coefficients. - Mircea Merca, Apr 06 2013
a(n-1) = floor(n/e^(1/n)) for n > 0. - Richard R. Forberg, Jun 22 2013
a(n) = A000027(n) for all n>0.
a(n) = floor(cot(1/(n+1))). - Clark Kimberling, Oct 08 2014
a(0)=0, a(n>0) = 2*z(-1)^[( |z|/z + 3 )/2] + ( |z|/z - 1 )/2 for z = A130472(n>0); a 1 to 1 correspondence between integers and naturals. - Adriano Caroli, Mar 29 2015
EXAMPLE
Triangular view:
0
1 2
3 4 5
6 7 8 9
10 11 12 13 14
15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30 31 32 33 34 35
36 37 38 39 40 41 42 43 44
45 46 47 48 49 50 51 52 53 54
MAPLE
[ seq(n, n=0..100) ];
MATHEMATICA
Table[n, {n, 0, 100}] (* Stefan Steinerberger, Apr 08 2006 *)
LinearRecurrence[{2, -1}, {0, 1}, 77] (* Robert G. Wilson v, May 23 2013 *)
CoefficientList[ Series[x/(x - 1)^2, {x, 0, 76}], x] (* Robert G. Wilson v, May 23 2013 *)
PROG
(Magma) [ n : n in [0..100]];
(PARI) A001477(n)=n /* first term is a(0) */
(Haskell)
a001477 = id
a001477_list = [0..] -- Reinhard Zumkeller, May 07 2012
(Python)
def a(n): return n
print([a(n) for n in range(78)]) # Michael S. Branicky, Nov 13 2022
(Julia) print([n for n in 0:280]) # Paul Muljadi, Apr 15 2024
CROSSREFS
Cf. A000027 (n>=1).
Cf. A000012 (first differences).
Partial sums of A057427. - Jeremy Gardiner, Sep 08 2002
Cf. A038608 (alternating signs), A001787 (binomial transform).
Cf. A055112.
Cf. Boustrophedon transforms: A231179, A000737.
Cf. A245422.
Number of orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008586, A008585, A005843, A000217.
When written as an array, the rows/columns are A000217, A000124, A152948, A152950, A145018, A167499, A166136, A167487... and A000096, A034856, A055998, A046691, A052905, A055999... (with appropriate offsets); cf. analogous lists for A000027 in A185787.
Cf. A000290.
Cf. A061579 (transposed matrix / reversed triangle).
KEYWORD
core,nonn,easy,mult,tabl
STATUS
approved
Pentagonal numbers: a(n) = n*(3*n-1)/2.
(Formerly M3818 N1562)
+10
541
0, 1, 5, 12, 22, 35, 51, 70, 92, 117, 145, 176, 210, 247, 287, 330, 376, 425, 477, 532, 590, 651, 715, 782, 852, 925, 1001, 1080, 1162, 1247, 1335, 1426, 1520, 1617, 1717, 1820, 1926, 2035, 2147, 2262, 2380, 2501, 2625, 2752, 2882, 3015, 3151
OFFSET
0,3
COMMENTS
The average of the first n (n > 0) pentagonal numbers is the n-th triangular number. - Mario Catalani (mario.catalani(AT)unito.it), Apr 10 2003
a(n) is the sum of n integers starting from n, i.e., 1, 2 + 3, 3 + 4 + 5, 4 + 5 + 6 + 7, etc. - Jon Perry, Jan 15 2004
Partial sums of 1, 4, 7, 10, 13, 16, ... (1 mod 3), a(2k) = k(6k-1), a(2k-1) = (2k-1)(3k-2). - Jon Perry, Sep 10 2004
Starting with offset 1 = binomial transform of [1, 4, 3, 0, 0, 0, ...]. Also, A004736 * [1, 3, 3, 3, ...]. - Gary W. Adamson, Oct 25 2007
If Y is a 3-subset of an n-set X then, for n >= 4, a(n-3) is the number of 4-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
Solutions to the duplication formula 2*a(n) = a(k) are given by the index pairs (n, k) = (5,7), (5577, 7887), (6435661, 9101399), etc. The indices are integer solutions to the pair of equations 2(6n-1)^2 = 1 + y^2, k = (1+y)/6, so these n can be generated from the subset of numbers [1+A001653(i)]/6, any i, where these are integers, confined to the cases where the associated k=[1+A002315(i)]/6 are also integers. - R. J. Mathar, Feb 01 2008
a(n) is a binomial coefficient C(n,4) (A000332) if and only if n is a generalized pentagonal number (A001318). Also see A145920. - Matthew Vandermast, Oct 28 2008
Even octagonal numbers divided by 8. - Omar E. Pol, Aug 18 2011
Sequence found by reading the line from 0, in the direction 0, 5, ... and the line from 1, in the direction 1, 12, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Sep 08 2011
The hyper-Wiener index of the star-tree with n edges (see A196060, example). - Emeric Deutsch, Sep 30 2011
More generally the n-th k-gonal number is equal to n + (k-2)*A000217(n-1), n >= 1, k >= 3. In this case k = 5. - Omar E. Pol, Apr 06 2013
Note that both Euler's pentagonal theorem for the partition numbers and Euler's pentagonal theorem for the sum of divisors refer more exactly to the generalized pentagonal numbers, not this sequence. For more information see A001318, A175003, A238442. - Omar E. Pol, Mar 01 2014
The Fuss-Catalan numbers are Cat(d,k)= [1/(k*(d-1)+1)]*binomial(k*d,k) and enumerate the number of (d+1)-gon partitions of a (k*(d-1)+2)-gon (cf. Schuetz and Whieldon link). a(n)= Cat(n,3), so enumerates the number of (n+1)-gon partitions of a (3*(n-1)+2)-gon. Analogous sequences are A100157 (k=4) and A234043 (k=5). - Tom Copeland, Oct 05 2014
Binomial transform of (0, 1, 3, 0, 0, 0, ...) (A169585 with offset 1) and second partial sum of (0, 1, 3, 3, 3, ...). - Gary W. Adamson, Oct 05 2015
For n > 0, a(n) is the number of compositions of n+8 into n parts avoiding parts 2 and 3. - Milan Janjic, Jan 07 2016
a(n) is also the number of edges in the Mycielskian of the complete graph K[n]. Indeed, K[n] has n vertices and n(n-1)/2 edges. Then its Mycielskian has n + 3n(n-1)/2 = n(3n-1)/2. See p. 205 of the West reference. - Emeric Deutsch, Nov 04 2016
Sum of the numbers from n to 2n-1. - Wesley Ivan Hurt, Dec 03 2016
Also the number of maximal cliques in the n-Andrásfai graph. - Eric W. Weisstein, Dec 01 2017
Coefficients in the hypergeometric series identity 1 - 5*(x - 1)/(2*x + 1) + 12*(x - 1)*(x - 2)/((2*x + 1)*(2*x + 2)) - 22*(x - 1)*(x - 2)*(x - 3)/((2*x + 1)*(2*x + 2)*(2*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A002412 and A002418. Column 2 of A103450. - Peter Bala, Mar 14 2019
A generalization of the Comment dated Apr 10 2003 follows. (k-3)*A000292(n-2) plus the average of the first n (2k-1)-gonal numbers is the n-th k-gonal number. - Charlie Marion, Nov 01 2020
a(n+1) is the number of Dyck paths of size (3,3n+1); i.e., the number of NE lattice paths from (0,0) to (3,3n+1) which stay above the line connecting these points. - Harry Richman, Jul 13 2021
REFERENCES
Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, pages 2 and 311.
Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; p. 129.
Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
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. 1.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 284.
Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 64.
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).
André Weil, Number theory: an approach through history; from Hammurapi to Legendre, Birkhäuser, Boston, 1984; see p. 186.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, 1987, pp. 98-100.
Douglas B. West, Introduction to Graph Theory, 2nd ed., Prentice-Hall, NJ, 2001.
LINKS
Daniel Mondot, Table of n, a(n) for n = 0..10000 (first 1000 terms by T. D. Noe)
George E. Andrews, Euler's "De Partitio Numerorum", Bull. Amer. Math. Soc., 44 (No. 4, 2007), 561-573.
S. Barbero, U. Cerruti and N. Murru, Transforming Recurrent Sequences by Using the Binomial and Invert Operators, J. Int. Seq. 13 (2010) # 10.7.7, section 4.4.
Jordan Bell, Euler and the pentagonal number theorem, arXiv:math/0510054 [math.HO], 2005-2006.
Anicius Manlius Severinus Boethius, De institutione arithmetica libri duo, Book 2, sections 13-14.
C. K. Cook and M. R. Bacon, Some polygonal number summation formulas, Fib. Q., 52 (2014), 336-343.
Stephan Eberhart, Letter to N. J. A. Sloane, Jan 06 1978, also scanned copy of Mathematical-Physical Correspondence, No. 22, Christmas 1977.
Leonhard Euler, Observatio de summis divisorum p. 8.
Leonhard Euler, An observation on the sums of divisors, arXiv:math/0411587 [math.HO], 2004, 2009. See p. 8.
Leonhard Euler, On the remarkable properties of the pentagonal numbers, arXiv:math/0505373 [math.HO], 2005.
Rodney T. Hansen, Arithmetic of pentagonal numbers, Fib. Quart., 8 (1970), 83-87.
Milan Janjic and Boris Petkovic, A Counting Function, arXiv preprint arXiv:1301.4550 [math.CO], 2013. - From N. J. A. Sloane, Feb 13 2013
Jangwon Ju and Daejun Kim, The pentagonal theorem of sixty-three and generalizations of Cauchy's lemma, arXiv:2010.16123 [math.NT], 2020.
Bir Kafle, Florian Luca and Alain Togbé, Pentagonal and heptagonal repdigits, Annales Mathematicae et Informaticae. pp. 137-145.
Sameen Ahmed Khan, Sums of the powers of reciprocals of polygonal numbers, Int'l J. of Appl. Math. (2020) Vol. 33, No. 2, 265-282.
Hyun Kwang Kim, On Regular Polytope Numbers, Proc. Amer. Math. Soc., 131 (2002), 65-75.
Clark Kimberling, Complementary Equations, Journal of Integer Sequences, Vol. 10 (2007), Article 07.1.4.
Clark Kimberling and John E. Brown, Partial Complements and Transposable Dispersions, J. Integer Seqs., Vol. 7, 2004.
R. P. Loh, A. G. Shannon and A. F. Horadam, Divisibility Criteria and Sequence Generators Associated with Fermat Coefficients, Preprint, 1980.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Cliff Reiter, Polygonal Numbers and Fifty One Stars, Lafayette College, Easton, PA (2019).
Alison Schuetz and Gwyneth Whieldon, Polygonal Dissections and Reversions of Series, arXiv:1401.7194 [math.CO], 2014.
W. Sierpiński, Sur les nombres pentagonaux, Bull. Soc. Royale Sciences Liège, 33 (No. 9-10, 1964), 513-517. [Annotated scanned copy]
Michel Waldschmidt, Continued fractions, Ecole de recherche CIMPA-Oujda, Théorie des Nombres et ses Applications, 18 - 29 mai 2015: Oujda (Maroc).
Eric Weisstein's World of Mathematics, Andrásfai Graph.
Eric Weisstein's World of Mathematics, Maximal Clique.
Eric Weisstein's World of Mathematics, Pentagonal Number.
Wikipedia, Mycielskian.
FORMULA
Product_{m > 0} (1 - q^m) = Sum_{k} (-1)^k*x^a(k). - Paul Barry, Jul 20 2003
G.f.: x*(1+2*x)/(1-x)^3.
E.g.f.: exp(x)*(x+3*x^2/2).
a(n) = n*(3*n-1)/2.
a(-n) = A005449(n).
a(n) = binomial(3*n, 2)/3. - Paul Barry, Jul 20 2003
a(n) = A000290(n) + A000217(n-1). - Lekraj Beedassy, Jun 07 2004
a(0) = 0, a(1) = 1; for n >= 2, a(n) = 2*a(n-1) - a(n-2) + 3. - Miklos Kristof, Mar 09 2005
a(n) = Sum_{k=1..n} (2*n - k). - Paul Barry, Aug 19 2005
a(n) = 3*A000217(n) - 2*n. - Lekraj Beedassy, Sep 26 2006
a(n) = A126890(n, n-1) for n > 0. - Reinhard Zumkeller, Dec 30 2006
a(n) = A049452(n) - A022266(n) = A033991(n) - A005476(n). - Zerinvary Lajos, Jun 12 2007
Equals A034856(n) + (n - 1)^2. Also equals A051340 * [1,2,3,...]. - Gary W. Adamson, Jul 27 2007
a(n) = binomial(n+1, 2) + 2*binomial(n, 2).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0) = 0, a(1) = 1, a(2) = 5. - Jaume Oliver Lafont, Dec 02 2008
a(n) = a(n-1) + 3*n-2 with n > 0, a(0)=0. - Vincenzo Librandi, Nov 20 2010
a(n) = A000217(n) + 2*A000217(n-1). - Vincenzo Librandi, Nov 20 2010
a(n) = A014642(n)/8. - Omar E. Pol, Aug 18 2011
a(n) = A142150(n) + A191967(n). - Reinhard Zumkeller, Jul 07 2012
a(n) = (A000290(n) + A000384(n))/2 = (A000217(n) + A000566(n))/2 = A049450(n)/2. - Omar E. Pol, Jan 11 2013
a(n) = n*A000217(n) - (n-1)*A000217(n-1). - Bruno Berselli, Jan 18 2013
a(n) = A005449(n) - n. - Philippe Deléham, Mar 30 2013
From Oskar Wieland, Apr 10 2013: (Start)
a(n) = a(n+1) - A016777(n),
a(n) = a(n+2) - A016969(n),
a(n) = a(n+3) - A016777(n)*3 = a(n+3) - A017197(n),
a(n) = a(n+4) - A016969(n)*2 = a(n+4) - A017641(n),
a(n) = a(n+5) - A016777(n)*5,
a(n) = a(n+6) - A016969(n)*3,
a(n) = a(n+7) - A016777(n)*7,
a(n) = a(n+8) - A016969(n)*4,
a(n) = a(n+9) - A016777(n)*9. (End)
a(n) = A000217(2n-1) - A000217(n-1), for n > 0. - Ivan N. Ianakiev, Apr 17 2013
a(n) = A002411(n) - A002411(n-1). - J. M. Bergot, Jun 12 2013
Sum_{n>=1} a(n)/n! = 2.5*exp(1). - Richard R. Forberg, Jul 15 2013
a(n) = floor(n/(exp(2/(3*n)) - 1)), for n > 0. - Richard R. Forberg, Jul 27 2013
From Vladimir Shevelev, Jan 24 2014: (Start)
a(3*a(n) + 4*n + 1) = a(3*a(n) + 4*n) + a(3*n+1).
A generalization. Let {G_k(n)}_(n >= 0) be sequence of k-gonal numbers (k >= 3). Then the following identity holds: G_k((k-2)*G_k(n) + c(k-3)*n + 1) = G_k((k-2)*G_k(n) + c(k-3)*n) + G_k((k-2)*n + 1), where c = A000124. (End)
A242357(a(n)) = 1 for n > 0. - Reinhard Zumkeller, May 11 2014
Sum_{n>=1} 1/a(n)= (1/3)*(9*log(3) - sqrt(3)*Pi). - Enrique Pérez Herrero, Dec 02 2014. See the decimal expansion A244641.
a(n) = (A000292(6*n+k-1)-A000292(k))/(6*n-1)-A000217(3*n+k), for any k >= 0. - Manfred Arens, Apr 26 2015 [minor edits from Wolfdieter Lang, May 10 2015]
a(n) = A258708(3*n-1,1) for n > 0. - Reinhard Zumkeller, Jun 23 2015
a(n) = A007584(n) - A245301(n-1), for n > 0. - Manfred Arens, Jan 31 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(sqrt(3)*Pi - 6*log(2))/3 = 0.85501000622865446... - Ilya Gutkovskiy, Jul 28 2016
a(m+n) = a(m) + a(n) + 3*m*n. - Etienne Dupuis, Feb 16 2017
In general, let P(k,n) be the n-th k-gonal number. Then P(k,m+n) = P(k,m) + (k-2)mn + P(k,n). - Charlie Marion, Apr 16 2017
a(n) = A023855(2*n-1) - A023855(2*n-2). - Luc Rousseau, Feb 24 2018
a(n) = binomial(n,2) + n^2. - Pedro Caceres, Jul 28 2019
Product_{n>=2} (1 - 1/a(n)) = 3/5. - Amiram Eldar, Jan 21 2021
(n+1)*(a(n^2) + a(n^2+1) + ... + a(n^2+n)) = n*(a(n^2+n+1) + ... + a(n^2+2n)). - Charlie Marion, Apr 28 2024
a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * binomial(k, 2)*binomial(3*n+k-1, 2*k). - Peter Bala, Nov 04 2024
EXAMPLE
Illustration of initial terms:
.
. o
. o o
. o o o o
. o o o o o o
. o o o o o o o o o
. o o o o o o o o o o o
. o o o o o o o o o o o o o
. o o o o o o o o o o o o o o
. o o o o o o o o o o o o o o o
.
. 1 5 12 22 35
- Philippe Deléham, Mar 30 2013
MAPLE
A000326 := n->n*(3*n-1)/2: seq(A000326(n), n=0..100);
A000326:=-(1+2*z)/(z-1)**3; # Simon Plouffe in his 1992 dissertation
a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=2*a[n-1]-a[n-2]+3 od: seq(a[n], n=0..50); # Miklos Kristof, Zerinvary Lajos, Feb 18 2008
MATHEMATICA
Table[n (3 n - 1)/2, {n, 0, 60}] (* Stefan Steinerberger, Apr 01 2006 *)
Array[# (3 # - 1)/2 &, 47, 0] (* Zerinvary Lajos, Jul 10 2009 *)
LinearRecurrence[{3, -3, 1}, {0, 1, 5}, 61] (* Harvey P. Dale, Dec 27 2011 *)
pentQ[n_] := IntegerQ[(1 + Sqrt[24 n + 1])/6]; pentQ[0] = True; Select[Range[0, 3200], pentQ@# &] (* Robert G. Wilson v, Mar 31 2014 *)
Join[{0}, Accumulate[Range[1, 312, 3]]] (* Harvey P. Dale, Mar 26 2016 *)
(* For Mathematica 10.4+ *) Table[PolygonalNumber[RegularPolygon[5], n], {n, 0, 46}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
CoefficientList[Series[x (-1 - 2 x)/(-1 + x)^3, {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)
PolygonalNumber[5, Range[0, 20]] (* Eric W. Weisstein, Dec 01 2017 *)
PROG
(PARI) a(n)=n*(3*n-1)/2
(PARI) vector(100, n, n--; binomial(3*n, 2)/3) \\ Altug Alkan, Oct 06 2015
(PARI) is_a000326 = my(s); n==0 || (issquare (24*n+1, &s) && s%6==5) \\ Hugo Pfoertner, Aug 03 2023
(Haskell)
a000326 n = n * (3 * n - 1) `div` 2 -- Reinhard Zumkeller, Jul 07 2012
(Magma) [n*(3*n-1)/2 : n in [0..100]]; // Wesley Ivan Hurt, Oct 15 2015
(GAP) List([0..50], n->n*(3*n-1)/2); # Muniru A Asiru, Mar 18 2019
(Python) # Intended to compute the initial segment of the sequence, not isolated terms.
def aList():
x, y = 1, 1
yield 0
while True:
yield x
x, y = x + y + 3, y + 3
A000326 = aList()
print([next(A000326) for i in range(47)]) # Peter Luschny, Aug 04 2019
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.
Cf. A001318 (generalized pentagonal numbers), A049452, A033570, A010815, A034856, A051340, A004736, A033568, A049453, A002411 (partial sums), A033579.
See A220083 for a list of numbers of the form n*P(s,n)-(n-1)*P(s,n-1), where P(s,n) is the n-th polygonal number with s sides.
Cf. A240137: sum of n consecutive cubes starting from n^3.
Cf. similar sequences listed in A022288.
Partial sums of A016777.
KEYWORD
core,nonn,easy,nice,changed
EXTENSIONS
Incorrect example removed by Joerg Arndt, Mar 11 2010
STATUS
approved
Multiples of 4.
+10
193
0, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 60, 64, 68, 72, 76, 80, 84, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 128, 132, 136, 140, 144, 148, 152, 156, 160, 164, 168, 172, 176, 180, 184, 188, 192, 196, 200, 204, 208, 212, 216, 220, 224, 228
OFFSET
0,2
COMMENTS
Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 14 ).
A000466(n), a(n) and A053755(n) are Pythagorean triples. - Zak Seidov, Jan 16 2007
If X is an n-set and Y and Z disjoint 2-subsets of X then a(n-3) is equal to the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Aug 26 2007
Number of n-permutations (n>=1) of 5 objects u, v, z, x, y with repetition allowed, containing n-1 u's. Example: if n=1 then n-1 = zero (0) u, a(1)=4 because we have v, z, x, y. If n=2 then n-1 = one (1) u, a(2)=8 because we have vu, zu, xu, yu, uv, uz, ux, uy. A038231 formatted as a triangular array: diagonal: 4, 8, 12, 16, 20, 24, 28, 32, ... - Zerinvary Lajos, Aug 06 2008
For n > 0: numbers having more even than odd divisors: A048272(a(n)) < 0. - Reinhard Zumkeller, Jan 21 2012
A214546(a(n)) < 0 for n > 0. - Reinhard Zumkeller, Jul 20 2012
A090418(a(n)) = 0 for n > 0. - Reinhard Zumkeller, Aug 06 2012
Terms are the differences of consecutive centered square numbers (A001844). - Mihir Mathur, Apr 02 2013
a(n)*Pi = nonnegative zeros of the cycloid generated by a circle of radius 2 rolling along the positive x-axis from zero. - Wesley Ivan Hurt, Jul 01 2013
Apart from the initial term, number of vertices of minimal path on an n-dimensional cubic lattice (n>1) of side length 2, until a self-avoiding walk gets stuck. A004767 + 1. - Matthew Lehman, Dec 23 2013
The number of orbits of Aut(Z^7) as function of the infinity norm n of the representative lattice point of the orbit, when the cardinality of the orbit is equal to 2688. - Philippe A.J.G. Chevalier, Dec 29 2015
First differences of A001844. - Robert Price, May 13 2016
Numbers k such that Fibonacci(k) is a multiple of 3 (A033888). - Bruno Berselli, Oct 17 2017
LINKS
Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 3.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 316 [Broken link]
Tanya Khovanova, Recursive Sequences
Franck Ramaharo, Statistics on some classes of knot shadows, arXiv:1802.07701 [math.CO], 2018.
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014, 2015.
William A. Stein, The modular forms database
Eric Weisstein's World of Mathematics, Doubly Even Number
FORMULA
a(n) = A008574(n), n>0. - R. J. Mathar, Oct 28 2008
a(n) = Sum_{k>=0} A030308(n,k)*2^(k+2). - Philippe Deléham, Oct 17 2011
a(n+1) = A000290(n+2) - A000290(n). - Philippe Deléham, Mar 31 2013
G.f.: 4*x/(1-x)^2. - David Wilding, Jun 21 2014
E.g.f.: 4*x*exp(x). - Stefano Spezia, May 18 2021
MAPLE
A008586:=n->4*n; seq(A008586(n), n=0..100); # Wesley Ivan Hurt, Feb 24 2014
MATHEMATICA
Range[0, 500, 4] (* Vladimir Joseph Stephan Orlovsky, May 26 2011 *)
PROG
(PARI) a(n)=n<<2 \\ Charles R Greathouse IV, Oct 17 2011
(Haskell)
a008586 = (* 4)
a008586_list = [0, 4 ..] -- Reinhard Zumkeller, May 13 2014
CROSSREFS
Number of orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008585, A005843, A001477, A000217.
KEYWORD
nonn,easy
STATUS
approved
Second pentagonal numbers: a(n) = n*(3*n + 1)/2.
+10
142
0, 2, 7, 15, 26, 40, 57, 77, 100, 126, 155, 187, 222, 260, 301, 345, 392, 442, 495, 551, 610, 672, 737, 805, 876, 950, 1027, 1107, 1190, 1276, 1365, 1457, 1552, 1650, 1751, 1855, 1962, 2072, 2185, 2301, 2420, 2542, 2667, 2795, 2926, 3060, 3197, 3337, 3480
OFFSET
0,2
COMMENTS
Number of edges in the join of the complete graph and the cycle graph, both of order n, K_n * C_n. - Roberto E. Martinez II, Jan 07 2002
Also number of cards to build an n-tier house of cards. - Martin Wohlgemuth, Aug 11 2002
The modular form Delta(q) = q*Product_{n>=1} (1-q^n)^24 = q*(1 + Sum_{n>=1} (-1)^n*(q^(n*(3*n-1)/2)+q^(n*(3*n+1)/2))^24 = q*(1 + Sum_{n>=1} A033999(n)*(q^A000326(n)+q^a(n))^24. - Jonathan Vos Post, Mar 15 2006
Row sums of triangle A134403.
Bisection of A001318. - Omar E. Pol, Aug 22 2011
Sequence found by reading the line from 0 in the direction 0, 7, ... and the line from 2 in the direction 2, 15, ... in the square spiral whose vertices are the generalized pentagonal numbers, A001318. - Omar E. Pol, Sep 08 2011
A general formula for the n-th second k-gonal number is given by T(n, k) = n*((k-2)*n+k-4)/2, n>=0, k>=5. - Omar E. Pol, Aug 04 2012
Partial sums give A006002. - Denis Borris, Jan 07 2013
A002260 is the following array A read by antidiagonals:
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
0, 1, 2, 3, 4, 5, ...
and a(n) is the hook sum: Sum_{k=0..n} A(n,k) + Sum_{r=0..n-1} A(r,n). - R. J. Mathar, Jun 30 2013
From Klaus Purath, May 13 2021: (Start)
This sequence and A000326 provide all integers m such that 24*m + 1 is a square. The union of the two sequences is A001318.
If A is a sequence satisfying the recurrence t(n) = 3*t(n-1) - 2*t(n-2) with the initial values either A(0) = 1, A(1) = n + 2 or A(0) = -1, A(1) = n - 1, then a(n) = (A(i)^2 - A(i-1)*A(i+1))/2^i + n^2 for i>0. (End)
a(n+1) is the number of Dyck paths of size (3,3n+2), i.e., the number of NE lattice paths from (0,0) to (3,3n+2) which stay above the line connecting these points. - Harry Richman, Jul 13 2021
Binomial transform of [0, 2, 3, 0, 0, 0, ...], being a(n) = 2*binomial(n,1) + 3*binomial(n,2). a(3) = 15 = [0, 2, 3, 0] dot [1, 3, 3, 1] = [0 + 6 + 9 + 0]. - Gary W. Adamson, Dec 17 2022
a(n) is the sum of longest side length of all nondegenerate integer-sided triangles with shortest side length n and middle side length (n + 1), n > 0. - Torlach Rush, Feb 04 2024
REFERENCES
Henri Cohen, A Course in Computational Algebraic Number Theory, vol. 138 of Graduate Texts in Mathematics, Springer-Verlag, 2000.
LINKS
A. O. L. Atkin and F. Morain, Elliptic Curves and Primality Proving, Math. Comp., Vol. 61, No. 203 (1993), pp. 29-68.
Charles H. Conley and Valentin Ovsienko, Quiddities of polygon dissections and the Conway-Coxeter frieze equation, arXiv:2107.01234 [math.CO], 2021.
Leonhard Euler, De mirabilibus proprietatibus numerorum pentagonalium, Acta Academiae Scientiarum Imperialis Petropolitanae, Vol. 1780: I, pp. 56-75.
Leonhard Euler, Observatio de summis divisorum, Novi Commentarii academiae scientiarum Petropolitanae, Vol. 5, pp. 59-74.
Leonhard Euler, An observation on the sums of divisors, arXiv:math/0411587 [math.HO], 2004-2009, p. 8.
Leonhard Euler, On the remarkable properties of the pentagonal numbers, arXiv:math/0505373 [math.HO], 2005.
D. Suprijanto and I. W. Suwarno, Observation on Sums of Powers of Integers Divisible by 3k-1, Applied Mathematical Sciences, Vol. 8, No. 45 (2014), pp. 2211-2217.
FORMULA
a(n) = A110449(n, 1) for n>0.
G.f.: x*(2+x)/(1-x)^3. E.g.f.: exp(x)*(2*x + 3*x^2/2). a(n) = n*(3*n + 1)/2. a(-n) = A000326(n). - Michael Somos, Jul 18 2003
a(n) = A001844(n) - A000217(n+1) = A101164(n+2,2) for n>0. - Reinhard Zumkeller, Dec 03 2004
a(n) = Sum_{j=1..n} n+j. - Zerinvary Lajos, Sep 12 2006
a(n) = A126890(n,n). - Reinhard Zumkeller, Dec 30 2006
a(n) = 2*C(3*n,4)/C(3*n,2), n>=1. - Zerinvary Lajos, Jan 02 2007
a(n) = A000217(n) + A000290(n). - Zak Seidov, Apr 06 2008
a(n) = a(n-1) + 3*n - 1 for n>0, a(0)=0. - Vincenzo Librandi, Nov 18 2010
a(n) = A129267(n+5,n). - Philippe Deléham, Dec 21 2011
a(n) = 2*A000217(n) + A000217(n-1). - Philippe Deléham, Mar 25 2013
a(n) = A130518(3*n+1). - Philippe Deléham, Mar 26 2013
a(n) = (12/(n+2)!)*Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*j^(n+2). - Vladimir Kruchinin, Jun 04 2013
a(n) = floor(n/(1-exp(-2/(3*n)))) for n>0. - Richard R. Forberg, Jun 22 2013
a(n) = Sum_{i=1..n} (3*i - 1) for n >= 1. - Wesley Ivan Hurt, Oct 11 2013 [Corrected by Rémi Guillaume, Oct 24 2024]
a(n) = (A000292(6*n+k+1)-A000292(k))/(6*n+1) - A000217(3*n+k+1), for any k >= 0. - Manfred Arens, Apr 26 2015
Sum_{n>=1} 1/a(n) = 6 - Pi/sqrt(3) - 3*log(3) = 0.89036376976145307522... . - Vaclav Kotesovec, Apr 27 2016
a(n) = A000217(2*n) - A000217(n). - Bruno Berselli, Sep 21 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*Pi/sqrt(3) + 4*log(2) - 6. - Amiram Eldar, Jan 18 2021
From Klaus Purath, May 13 2021: (Start)
Partial sums of A016789 for n > 0.
a(n) = 3*n^2 - A000326(n).
a(n) = A000326(n) + n.
a(n) = A002378(n) + A000217(n-1) for n >= 1. [Corrected by Rémi Guillaume, Aug 14 2024] (End)
From Klaus Purath, Jul 14 2021: (Start)
b^2 = 24*a(n) + 1 we get by b^2 = (a(n+1) - a(n-1))^2 = (a(2*n)/n)^2.
a(2*n) = n*(a(n+1) - a(n-1)), n > 0.
a(2*n+1) = n*(a(n+1) - a(n)). (End)
A generalization of Lajos' formula, dated Sep 12 2006, follows. Let SP(k,n) = the n-th second k-gonal number. Then SP(2k+1,n) = Sum_{j=1..n} (k-1)*n+j+k-2. - Charlie Marion, Jul 13 2024
a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * binomial(k, 2) * binomial(3*n+k, 2*k). - Peter Bala, Nov 03 2024
EXAMPLE
From Omar E. Pol, Aug 22 2011: (Start)
Illustration of initial terms:
. O
. O O
. O O O O
. O O O O O O
. O O O O O O O O O
. O O O O O O O O O O O
. O O O O O O O O O O O O O
. O O O O O O O O O O O O O O
. O O O O O O O O O O O O O O O
. O O O O O O O O O O O O O O O
.
. 2 7 15 26 40
.
(End)
MAPLE
A005449:=n->n*(3*n + 1)/2; seq(A005449(k), k=0..100); # Wesley Ivan Hurt, Oct 11 2013
MATHEMATICA
Table[n (3 n + 1)/2, {n, 0, 100}] (* Zak Seidov, Jan 31 2012 *)
PROG
(PARI) {a(n) = n * (3*n + 1) / 2} /* Michael Somos, Jul 18 2003 */
(Magma) [n*(3*n + 1) / 2: n in [0..40]]; // Vincenzo Librandi, May 02 2011
CROSSREFS
Cf. A016789 (first differences), A006002 (partial sums).
The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, this sequence, A045943, A115067, A140090, A140091, A059845, A140672-A140675, A151542.
Cf. numbers of the form n*(n*k-k+4)/2 listed in A226488 (this sequence is the case k=3).
Cf. numbers of the form n*((2*k+1)*n+1)/2 listed in A022289 (this sequence is the case k=1).
KEYWORD
nonn,easy,changed
STATUS
approved
Triangular matchstick numbers: a(n) = 3*n*(n+1)/2.
+10
128
0, 3, 9, 18, 30, 45, 63, 84, 108, 135, 165, 198, 234, 273, 315, 360, 408, 459, 513, 570, 630, 693, 759, 828, 900, 975, 1053, 1134, 1218, 1305, 1395, 1488, 1584, 1683, 1785, 1890, 1998, 2109, 2223, 2340, 2460, 2583, 2709, 2838, 2970, 3105, 3243, 3384, 3528
OFFSET
0,2
COMMENTS
Also, 3 times triangular numbers, a(n) = 3*A000217(n).
In the 24-bit RGB color cube, the number of color-lattice-points in r+g+b = n planes at n < 256 equals the triangular numbers. For n = 256, ..., 765 the number of legitimate color partitions is less than A000217(n) because {r,g,b} components cannot exceed 255. For n = 256, ..., 511, the number of non-color partitions are computable with A045943(n-255), while for n = 512, ..., 765, the number of color points in r+g+b planes equals A000217(765-n). - Labos Elemer, Jun 20 2005
If a 3-set Y and an (n-3)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
a(n) is also the smallest number that may be written both as the sum of n-1 consecutive positive integers and n consecutive positive integers. - Claudio Meller, Oct 08 2010
For n >= 3, a(n) equals 4^(2+n)*Pi^(1 - n) times the coefficient of zeta(3) in the following integral with upper bound Pi/4 and lower bound 0: int x^(n+1) tan x dx. - John M. Campbell, Jul 17 2011
The difference a(n)-a(n-1) = 3*n, for n >= 1. - Stephen Balaban, Jul 25 2011 [Comment clarified by N. J. A. Sloane, Aug 01 2024]]
Sequence found by reading the line from 0, in the direction 0, 3, ..., and the same line from 0, in the direction 0, 9, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. This is one of the orthogonal axes of the spiral; the other is A032528. - Omar E. Pol, Sep 08 2011
A005449(a(n)) = A000332(3n + 3) = C(3n + 3, 4), a second pentagonal number of triangular matchstick number index number. Additionally, a(n) - 2n is a pentagonal number (A000326). - Raphie Frank, Dec 31 2012
Sum of the numbers from n to 2n. - Wesley Ivan Hurt, Nov 24 2015
Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 5376 or 17920 or 20160. - Philippe A.J.G. Chevalier, Dec 28 2015
Also the number of 4-cycles in the (n+4)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017
Number of terms less than 10^k, k=0,1,2,3,...: 1, 3, 8, 26, 82, 258, 816, 2582, 8165, 25820, 81650, 258199, 816497, 2581989, 8164966, ... - Muniru A Asiru, Jan 24 2018
Numbers of the form 3*m*(2*m + 1) for m = 0, -1, 1, -2, 2, -3, 3, ... - Bruno Berselli, Feb 26 2018
Partial sums of A008585. - Omar E. Pol, Jun 20 2018
Column 1 of A273464. (Number of ways to select a unit lozenge inside an isosceles triangle of side length n; all vertices on a hexagonal lattice.) - R. J. Mathar, Jul 10 2019
Total number of pips in the n-th suit of a double-n domino set. - Ivan N. Ianakiev, Aug 23 2020
REFERENCES
Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 543.
LINKS
Abderrahim Arabi, Hacène Belbachir, and Jean-Philippe Dubernard, Enumeration of parallelogram polycubes, arXiv:2105.00971 [cs.DM], 2021.
Francesco Brenti and Paolo Sentinelli, Wachs permutations, Bruhat order and weak order, arXiv:2212.04932 [math.CO], 2022.
Sela Fried, Counting r X s rectangles in nondecreasing and Smirnov words, arXiv:2406.18923 [math.CO], 2024. See p. 5.
Jose Manuel Garcia Calcines, Luis Javier Hernandez Paricio, and Maria Teresa Rivas Rodriguez, Semi-simplicial combinatorics of cyclinders and subdivisions, arXiv:2307.13749 [math.CO], 2023. See p. 29.
T. Aaron Gulliver, Sums of Powers of Integers Divisible by Three, Int. J. Contemp. Math. Sciences, Vol. 7, 2012, no. 38, 1895 - 1901.
Milan Janjic and Boris Petkovic, A Counting Function, arXiv preprint arXiv:1301.4550 [math.CO], 2013.
Milan Janjic and Boris Petkovic, A Counting Function Generalizing Binomial Coefficients and Some Other Classes of Integers, J. Int. Seq. 17 (2014) # 14.3.5.
E. Lábos, On the number of RGB-colors we can distinguish. Partition Spectra. Lecture at 7th Hungarian Conference on Biometry and Biomathematics. Budapest. Jul 06, 2005. Applied Ecology and Environmental Research 4(2): 159-169, 2006.
R. J. Mathar, Lozenge tilings of the equilateral triangle, arXiv:1909.06336 [math.CO], 2019.
Eric Weisstein's World of Mathematics, Graph Cycle.
FORMULA
a(n) is the sum of n+1 integers starting from n, i.e., 1+2, 2+3+4, 3+4+5+6, 4+5+6+7+8, etc. - Jon Perry, Jan 15 2004
a(n) = A126890(n+1,n-1) for n>1. - Reinhard Zumkeller, Dec 30 2006
a(n) + A145919(3*n+3) = 0. - Matthew Vandermast, Oct 28 2008
a(n) = A000217(2*n) - A000217(n-1); A179213(n) <= a(n). - Reinhard Zumkeller, Jul 05 2010
a(n) = a(n-1)+3*n, n>0. - Vincenzo Librandi, Nov 18 2010
G.f.: 3*x/(1-x)^3. - Bruno Berselli, Jan 21 2011
a(n) = A005448(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(n) = A001477(n)+A000290(n)+A000217(n). - J. M. Bergot, Dec 08 2012
a(n) = 3*a(n-1)-3*a(n-2)+a(n-3) for n>2. - Wesley Ivan Hurt, Nov 24 2015
a(n) = A027480(n)-A027480(n-1). - Peter M. Chema, Jan 18 2017.
2*a(n)+1 = A003215(n). - Miquel Cerda, Jan 22 2018
a(n) = T(2*n) - T(n-1), where T(n) = A000217(n). In general, T(k)*T(n) = Sum_{i=0..k-1} (-1)^i*T((k-i)*(n-i)). - Charlie Marion, Dec 06 2020
E.g.f.: 3*exp(x)*x*(2 + x)/2. - Stefano Spezia, May 19 2021
From Amiram Eldar, Jan 10 2022: (Start)
Sum_{n>=1} 1/a(n) = 2/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(2*log(2)-1)/3. (End)
Product_{n>=1} (1 - 1/a(n)) = -(3/(2*Pi))*cos(sqrt(11/3)*Pi/2). - Amiram Eldar, Feb 21 2023
EXAMPLE
From Stephen Balaban, Jul 25 2011: (Start)
T(n), the triangular numbers = number of nodes,
a(n-1) = number of edges in the T(n) graph:
o (T(1) = 1, a(0) = 0)
o
/ \ (T(2) = 3, a(1) = 3)
o - o
o
/ \
o - o (T(3) = 6, a(2) = 9)
/ \ / \
o - o - o
... [Corrected by N. J. A. Sloane, Aug 01 2024] (End)
MAPLE
seq(3*binomial(n+1, 2), n=0..49); # Zerinvary Lajos, Nov 24 2006
MATHEMATICA
Table[3 n (n + 1)/2, {n, 0, 50}] (* Vladimir Joseph Stephan Orlovsky, Oct 31 2008 *)
3 Accumulate@Range[0, 48] (* Arkadiusz Wesolowski, Oct 29 2012 *)
CoefficientList[Series[-3 x/(x - 1)^3, {x, 0, 47}], x] (* Robert G. Wilson v, Jan 29 2015 *)
LinearRecurrence[{3, -3, 1}, {0, 3, 9}, 50] (* Jean-François Alcover, Dec 12 2016 *)
PROG
(Common Lisp) (defun tri (i) (if (eq i 0) 0 (+ (* 3 (- i 1)) (tri (- i 1))))) // Stephen Balaban, Jul 25 2011
(Magma) [3*n*(n+1)/2: n in [0..50]]; // Vincenzo Librandi, May 02 2011
(PARI) a(n)=3*binomial(n+1, 2) \\ Charles R Greathouse IV, Jun 16 2011
(Haskell) a n = sum [x | x <- [n..2*n]] -- Peter Kagey, Jul 27 2015
(GAP) List([0..10^4], n -> 3*n*(n+1)/2); # Muniru A Asiru, Jan 24 2018
(Scala) (3 to 150 by 3).scanLeft(0)(_ + _) // Alonso del Arte, Sep 12 2019
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.
A diagonal of A010027.
Orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A115067, A008585, A005843, A001477, A000217.
Cf. A027480 (partial sums).
Cf. A002378 (3-cycles in triangular honeycomb acute knight graph), A028896 (5-cycles), A152773 (6-cycles).
This sequence: Sum_{k = n..2*n} k.
Cf. A304993: Sum_{k = n..2*n} k*(k+1)/2.
Cf. A050409: Sum_{k = n..2*n} k^2.
Similar sequences are listed in A316466.
KEYWORD
nonn,easy
AUTHOR
STATUS
approved
Hexagonal pyramidal numbers, or greengrocer's numbers.
(Formerly M4374 N1839)
+10
94
0, 1, 7, 22, 50, 95, 161, 252, 372, 525, 715, 946, 1222, 1547, 1925, 2360, 2856, 3417, 4047, 4750, 5530, 6391, 7337, 8372, 9500, 10725, 12051, 13482, 15022, 16675, 18445, 20336, 22352, 24497, 26775, 29190, 31746, 34447, 37297, 40300
OFFSET
0,3
COMMENTS
Binomial transform of (1, 6, 9, 4, 0, 0, 0, ...). - Gary W. Adamson, Oct 16 2007
a(n) is the sum of the maximum(m,n) over {(m,n): m,n in positive integers, m<=n}. - Geoffrey Critzer, Oct 11 2009
We obtain these numbers for d=2 in the identity n*(n*(d*n-d+2)/2)-sum(k*(d*k-d+2)/2, k=0..n-1) = n*(n+1)*(2*d*n-2*d+3)/6 (see Klaus Strassburger in Formula lines). - Bruno Berselli, Apr 21 2010, Nov 16 2010
q^a(n) is the Hankel transform of the q-Catalan numbers. - Paul Barry, Dec 15 2010
Row 1 of the convolution array A213835. - Clark Kimberling, Jul 04 2012
From Ant King, Oct 24 2012: (Start)
For n>0, the digital roots of this sequence A010888(A002412(n)) form the purely periodic 27-cycle {1,7,4,5,5,8,9,3,3,4,1,7,8,8,2,3,6,6,7,4,1,2,2,5,6,9,9}.
For n>0, the units' digits of this sequence A010879(A002412(n)) form the purely periodic 20-cycle {1,7,2,0,5,1,2,2,5,5,6,2,7,5,0,6,7,7,0,0}.
(End)
Partial sums of A000384. - Omar E. Pol, Jan 12 2013
Row sums of A094728. - J. M. Bergot, Jun 14 2013
Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 40320. - Philippe A.J.G. Chevalier, Dec 28 2015
Coefficients in the hypergeometric series identity 1 - 7*(x - 1)/(3*x + 1) + 22*(x - 1)*(x - 2)/((3*x + 1)*(3*x + 2)) - 50*(x - 1)*(x - 2)*(x - 3)/((3*x + 1)*(3*x + 2)*(3*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A000326 and A002418. Column 3 of A103450. - Peter Bala, Mar 14 2019
REFERENCES
A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
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.
T. A. Gulliver, Sequences from Arrays of Integers, Int. Math. Journal, Vol. 1, No. 4, pp. 323-332, 2002.
I. Siap, Linear codes over F_2 + u*F_2 and their complete weight enumerators, in Codes and Designs (Ohio State, May 18, 2000), pp. 259-271. De Gruyter, 2002.
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).
LINKS
T. D. Noe and William A. Tedeschi, Table of n, a(n) for n = 0..10000 (first 1000 terms computed by T. D. Noe)
Abdullah Atmaca and A. Yavuz Oruç, On the size of two families of unlabeled bipartite graphs, AKCE International Journal of Graphs and Combinatorics, Vo. 16, No. 2 (2019), pp. 222-229.
Bruno Berselli, A description of the transform in Comments lines: website Matem@ticamente (in Italian).
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
Luis Verde-Star, A Matrix Approach to Generalized Delannoy and Schröder Arrays, J. Int. Seq., Vol. 24 (2021), Article 21.4.1.
Eric Weisstein's World of Mathematics, Hexagonal Pyramidal Number.
FORMULA
a(n) = n(n + 1)(4n - 1)/6.
G.f.: x*(1+3*x)/(1-x)^4. - Simon Plouffe in his 1992 dissertation.
a(n) = n^3 - Sum_{i=1..n-1} i^2. - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de)
Partial sums of n odd-indexed triangular numbers, e.g., a(3) = t(1)+t(3)+t(5) = 1+6+15 = 22. - Jon Perry, Jul 23 2003
a(n) = Sum_{i=0..n-1} (n - i)*(n + i). - Jon Perry, Sep 26 2004
a(n) = n*A000292(n) - (n-1)*A000292(n-1) = n*binomial((n+2),3) - (n-1)*binomial((n+1),3); e.g., a(5) = 95 = 5*35 - 4*20. - Gary W. Adamson, Dec 28 2007
a(n) = Sum_{i=0..n} (2i^2 + 3i + 1), for n >= 0 (Omits the leading 0). - William A. Tedeschi, Aug 25 2010
a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4), with a(0)=0, a(1)=1, a(2)=7, a(3)=22. - Harvey P. Dale, Jul 16 2011
a(n) = sum a*b, where the summing is over all unordered partitions 2*n = a+b. - Vladimir Shevelev, May 11 2012
From Ant King, Oct 24 2012: (Start)
a(n) = a(n-1) + n*(2*n-1).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 4.
a(n) = (n+1)*(2*A000384(n) + n)/6 = (4*n-1)*A000217(n)/3.
a(n) = A000292(n) + 3*A000292(n-1) = A002411(n) + A000292(n-1).
a(n) = binomial(n+2,3) + 3*binomial(n+1,3) = (4*n-1)*binomial(n+1,2)/3.
Sum_{n>=1} 1/a(n) = 6*(12*log(2)-2*Pi-1)/5 = 1.2414...
(End)
a(n) = Sum_{i=1..n} Sum_{j=1..n} max(i,j) = Sum_{i=1..n} i*(2*n-i). - Enrique Pérez Herrero, Jan 15 2013
a(n) = A005900(n+1) - A000326(n+1) = Octahedral - Pentagonal Numbers. - Richard R. Forberg, Aug 07 2013
a(n) = n*A000217(n) + Sum_{i=0..n-1} A000217(i). - Bruno Berselli, Dec 18 2013
a(n) = 2n * A000217(n) - A000330(n). - J. M. Bergot, Apr 05 2014
a(n) = A080851(4,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: x*(6 + 15*x + 4*x^2)*exp(x)/6. - Ilya Gutkovskiy, May 12 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = 6*(1 + 2*sqrt(2)*Pi - 2*(3+sqrt(2))*log(2) + 4*sqrt(2)*log(2-sqrt(2)))/5. - Amiram Eldar, Jan 04 2022
EXAMPLE
Let n=5, 2*n=10. Since 10 = 1 + 9 = 2 + 8 = 3 + 7 = 4 + 6 = 5 + 5, a(5) = 1*9 + 2*8 + 3*7 + 4*6 + 5*5 = 95. - Vladimir Shevelev, May 11 2012
MAPLE
seq(sum(i*(2*k-i), i=1..k), k=0..100); # Wesley Ivan Hurt, Sep 25 2013
MATHEMATICA
Figurate[ ngon_, rank_, dim_] := Binomial[rank + dim - 2, dim - 1] ((rank - 1)*(ngon - 2) + dim)/dim; Table[ Figurate[6, r, 3], {r, 0, 40}] (* Robert G. Wilson v, Aug 22 2010 *)
Table[n(n+1)(4n-1)/6, {n, 0, 40}] (* or *) LinearRecurrence[{4, -6, 4, -1}, {0, 1, 7, 22}, 40] (* Harvey P. Dale, Jul 16 2011 *)
PROG
(PARI) v=vector(40, i, (i*(i+1))\2); s=0; print1(s", "); forstep(i=1, 40, 2, s+=v[i]; print1(s", "))
(Maxima) A002412(n):=n*(n+1)*(4*n-1)/6$ makelist(A002412(n), n, 0, 20); /* Martin Ettl, Dec 12 2012 */
(Magma) [n*(n+1)*(4*n-1)/6: n in [0..40]]; // Vincenzo Librandi, Nov 28 2015
(GAP) List([0..40], n->n*(n+1)*(4*n-1)/6); # Muniru A Asiru, Mar 18 2019
(Python) print([n*(n+1)*(4*n-1)//6 for n in range(40)]) # Michael S. Branicky, Mar 28 2022
CROSSREFS
Bisection of A002623. Equals A000578(n) - A000330(n-1).
a(n) = A093561(n+2, 3), (4, 1)-Pascal column.
Cf. A220084 for a list of numbers of the form n*P(k,n)-(n-1)*P(k,n-1), where P(k,n) is the n-th k-gonal pyramidal number (see Adamson's formula).
Cf. similar sequences listed in A237616.
Orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008585, A005843, A001477, A000217.
KEYWORD
nonn,easy,nice
STATUS
approved
a(n) = n*(3*n + 7)/2.
+10
35
0, 5, 13, 24, 38, 55, 75, 98, 124, 153, 185, 220, 258, 299, 343, 390, 440, 493, 549, 608, 670, 735, 803, 874, 948, 1025, 1105, 1188, 1274, 1363, 1455, 1550, 1648, 1749, 1853, 1960, 2070, 2183, 2299, 2418, 2540, 2665, 2793, 2924
OFFSET
0,2
COMMENTS
This sequence is mentioned in the Guo-Niu Han's paper, chapter 6: Dictionary of the standard puzzle sequences, p. 19 (see link). - Omar E. Pol, Oct 28 2011
Number of cards needed to build an n-tier house of cards with a flat, one-card-wide roof. - Tyler Busby, Dec 28 2022
LINKS
Guo-Niu Han, Enumeration of Standard Puzzles, 2011. [Cached copy]
Guo-Niu Han, Enumeration of Standard Puzzles, arXiv:2006.14070 [math.CO], 2020.
Sela Fried, Counting r X s rectangles in nondecreasing and Smirnov words, arXiv:2406.18923 [math.CO], 2024. See p. 5.
FORMULA
G.f.: x*(5 - 2*x)/(1 - x)^3. - Bruno Berselli, Feb 11 2011
a(n) = (3*n^2 + 7*n)/2.
a(n) = a(n-1) + 3*n + 2 (with a(0)=0). - Vincenzo Librandi, Nov 24 2010
E.g.f.: (1/2)*(3*x^2 + 10*x)*exp(x). - G. C. Greubel, Jul 17 2017
From Amiram Eldar, Feb 22 2022: (Start)
Sum_{n>=1} 1/a(n) = 117/98 - Pi/(7*sqrt(3)) - 3*log(3)/7.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*Pi/(7*sqrt(3)) + 4*log(2)/7 - 75/98. (End)
MATHEMATICA
LinearRecurrence[{3, -3, 1}, {0, 5, 13}, 50] (* Harvey P. Dale, Jan 17 2022 *)
PROG
(PARI) a(n)=n*(3*n+7)/2 \\ Charles R Greathouse IV, Sep 24 2015
CROSSREFS
The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, this sequence, A140091, A059845, A140672, A140673, A140674, A140675, A151542.
Cf. numbers of the form n*(d*n + 10 - d)/2: A008587, A056000, A028347, A014106, A028895, A045944, A186029, A007742, A022267, A033429, A022268, A049452, A186030, A135703, A152734, A139273.
KEYWORD
easy,nonn
AUTHOR
Omar E. Pol, May 22 2008
STATUS
approved
a(n) = 3*n*(n + 3)/2.
+10
27
0, 6, 15, 27, 42, 60, 81, 105, 132, 162, 195, 231, 270, 312, 357, 405, 456, 510, 567, 627, 690, 756, 825, 897, 972, 1050, 1131, 1215, 1302, 1392, 1485, 1581, 1680, 1782, 1887, 1995, 2106, 2220, 2337, 2457, 2580, 2706, 2835, 2967
OFFSET
0,2
COMMENTS
a(n) is also the dimension of the irreducible representation of the Lie algebra sl(3) with the highest weight 2*L_1+n*(L_1+L_2). - Leonid Bedratyuk, Jan 04 2010
Number of edges in the hexagonal triangle, T(n) (see the He et al. reference). - Emeric Deutsch, Nov 14 2014
a(n) = twice the area of a triangle having vertices at binomials (C(n,3),C(n+3,3)), (C(n+1,3),C(n+4,3)), and (C(n+2,3),C(n+5,3)) with n>=2. - J. M. Bergot, Mar 01 2018
REFERENCES
W. Fulton, J. Harris, Representation theory: a first course. (1991). page 224, Exercise 15.19. - Leonid Bedratyuk, Jan 04 2010
LINKS
Sela Fried, Counting r X s rectangles in nondecreasing and Smirnov words, arXiv:2406.18923 [math.CO], 2024. See p. 5.
Q. H. He, J. Z. Gu, S. J. Xu, and W. H. Chan, Hosoya polynomials of hexagonal triangles and trapeziums, MATCH, Commun. Math. Comput. Chem., Vol. 72, No. 3 (2014), pp. 835-843. - Emeric Deutsch, Nov 14 2014
FORMULA
a(n) = A000096(n)*3 = (3*n^2 + 9*n)/2 = n*(3*n+9)/2.
a(n) = a(n-1) + 3*n + 3 with n>0, a(0)=0. - Vincenzo Librandi, Nov 24 2010
G.f.: 3*x*(2 - x)/(1 - x)^3. - Arkadiusz Wesolowski, Dec 24 2011
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>2. - Harvey P. Dale, Aug 15 2015
E.g.f.: (1/2)*(3*x^2 + 12*x)*exp(x). - G. C. Greubel, Jul 17 2017
From Amiram Eldar, Feb 25 2022: (Start)
Sum_{n>=1} 1/a(n) = 11/27.
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(2)/9 - 5/27. (End)
MAPLE
A140091:=n->3*n*(n+3)/2: seq(A140091(n), n=0..50); # Wesley Ivan Hurt, Nov 14 2014
MATHEMATICA
Table[3 n (n + 3)/2, {n, 0, 50}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 6, 15}, 50] (* Harvey P. Dale, Aug 15 2015 *)
PROG
(Magma) [3*n*(n+3)/2 : n in [0..50]]; // Wesley Ivan Hurt, Nov 14 2014
(PARI) a(n)=3*n*(n+3)/2 \\ Charles R Greathouse IV, Sep 24 2015
CROSSREFS
The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, A140090, this sequence, A059845, A140672, A140673, A140674, A140675, A151542.
KEYWORD
easy,nonn
AUTHOR
Omar E. Pol, May 22 2008
STATUS
approved

Search completed in 0.052 seconds