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

A001477 The nonnegative integers.

Original entry on oeis.org

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

Views

Author

Keywords

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

Examples

			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
		

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.

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).

Programs

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
G.f. as triangle: x*(1 + (x^2 - 5*x + 2)*y + x*(2*x - 1)*y^2)/((1 - x)^3*(1 - x*y)^3). - Stefano Spezia, Jul 22 2025

A008586 Multiples of 4.

Original entry on oeis.org

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

Views

Author

Keywords

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

Crossrefs

Number of orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008585, A005843, A001477, A000217.

Programs

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

A045943 Triangular matchstick numbers: a(n) = 3*n*(n+1)/2.

Original entry on oeis.org

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

Views

Author

Keywords

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

Examples

			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)
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 543.

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.

Programs

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

A002412 Hexagonal pyramidal numbers, or greengrocer's numbers.

Original entry on oeis.org

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

Views

Author

Keywords

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

Examples

			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
		

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).

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.

Programs

  • GAP
    List([0..40],n->n*(n+1)*(4*n-1)/6); # Muniru A Asiru, Mar 18 2019
    
  • Magma
    [n*(n+1)*(4*n-1)/6: n in [0..40]]; // Vincenzo Librandi, Nov 28 2015
    
  • 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 *)
  • Maxima
    A002412(n):=n*(n+1)*(4*n-1)/6$ makelist(A002412(n),n,0,20); /* Martin Ettl, Dec 12 2012 */
    
  • PARI
    v=vector(40,i,(i*(i+1))\2); s=0; print1(s","); forstep(i=1,40,2,s+=v[i]; print1(s","))
    
  • Python
    print([n*(n+1)*(4*n-1)//6 for n in range(40)]) # Michael S. Branicky, Mar 28 2022

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

A153641 Nonzero coefficients of the Swiss-Knife polynomials for the computation of Euler, tangent, and Bernoulli numbers (triangle read by rows).

Original entry on oeis.org

1, 1, 1, -1, 1, -3, 1, -6, 5, 1, -10, 25, 1, -15, 75, -61, 1, -21, 175, -427, 1, -28, 350, -1708, 1385, 1, -36, 630, -5124, 12465, 1, -45, 1050, -12810, 62325, -50521, 1, -55, 1650, -28182, 228525, -555731, 1, -66, 2475, -56364, 685575, -3334386, 2702765, 1
Offset: 0

Views

Author

Peter Luschny, Dec 29 2008

Keywords

Comments

In the following the expression [n odd] is 1 if n is odd, 0 otherwise.
(+) W_n(0) = E_n are the Euler (or secant) numbers A122045.
(+) W_n(1) = T_n are the signed tangent numbers, see A009006.
(+) W_{n-1}(1) n / (4^n - 2^n) = B_n gives for n > 1 the Bernoulli number A027641/A027642.
(+) W_n(-1) 2^{-n}(n+1) = G_n the Genocchi number A036968.
(+) W_n(1/2) 2^{n} are the signed generalized Euler (Springer) number, see A001586.
(+) | W_n([n odd]) | the number of alternating permutations A000111.
(+) | W_n([n odd]) / n! | for 0<=n the Euler zeta number A099612/A099617 (see Wikipedia on Bernoulli number). - Peter Luschny, Dec 29 2008
The diagonals in the full triangle (with zero coefficients) of the polynomials have the general form E(k)*binomial(n+k,k) (k>=0 fixed, n=0,1,...) where E(n) are the Euler numbers in the enumeration A122045. For k=2 we find the triangular numbers A000217 and for k=4 A154286. - Peter Luschny, Jan 06 2009
From Peter Bala, Jun 10 2009: (Start)
The Swiss-Knife polynomials W_n(x) may be expressed in terms of the Bernoulli polynomials B(n,x) as
... W_n(x) = 4^(n+1)/(2*n+2)*[B(n+1,(x+3)/4) - B(n+1,(x+1)/4)].
The Swiss-Knife polynomials are, apart from a multiplying factor, examples of generalized Bernoulli polynomials.
Let X be the Dirichlet character modulus 4 defined by X(4*n+1) = 1, X(4*n+3) = -1 and X(2*n) = 0. The generalized Bernoulli polynomials B(X;n,x), n = 1,2,..., associated with the character X are defined by means of the generating function
... t*exp(x*t)*(exp(t)-exp(3*t))/(exp(4*t)-1) = sum {n = 1..inf} B(X;n,x)*t^n/n!.
The first few values are B(X;1,x) = -1/2, B(X;2,x) = -x, B(X,3,x) = -3/2*(x^2-1) and B(X;4,x) = -2*(x^3-3*x).
In general, W_n(x) = -2/(n+1)*B(X;n+1,x).
For the theory of generalized Bernoulli polynomials associated to a periodic arithmetical function see [Cohen, Section 9.4].
The generalized Bernoulli polynomials may be used to evaluate twisted sums of k-th powers. For the present case the result is
sum{n = 0..4*N-1} X(n)*n^k = 1^k - 3^k + 5^k - 7^k + ... - (4*N-1)^k
= [B(X;k+1,4*N) - B(X;k+1,0)]/(k+1) = [W_k(0) - W_k(4*N)]/2.
For the proof apply [Cohen, Corollary 9.4.17 with m = 4 and x = 0].
The generalized Bernoulli polynomials and the Swiss-Knife polynomials are also related to infinite sums of powers through their Fourier series - see the formula section below. For a table of the coefficients of generalized Bernoulli polynomials attached to a Dirichlet character modulus 8 see A151751.
(End)
The Swiss-Knife polynomials provide a general formula for alternating sums of powers similar to the formula which are provided by the Bernoulli polynomials for non-alternating sums of powers (see the Luschny link). Sequences covered by this formula include A001057, A062393, A062392, A011934, A144129, A077221, A137501, A046092. - Peter Luschny, Jul 12 2009
The greatest common divisor of the nonzero coefficients of the decapitated Swiss-Knife polynomials is exp(Lambda(n)), where Lambda(n) is the von Mangoldt function for odd primes, symbolically:
gcd(coeffs(SKP_{n}(x) - x^n)) = A155457(n) (n>1). - Peter Luschny, Dec 16 2009
Another version is at A119879. - Philippe Deléham, Oct 26 2013

Examples

			1
x
x^2  -1
x^3  -3x
x^4  -6x^2   +5
x^5 -10x^3  +25x
x^6 -15x^4  +75x^2  -61
x^7 -21x^5 +175x^3 -427x
		

References

  • H. Cohen, Number Theory - Volume II: Analytic and Modern Tools, Graduate Texts in Mathematics. Springer-Verlag. [From Peter Bala, Jun 10 2009]

Crossrefs

W_n(k), k=0,1,...
W_0: 1, 1, 1, 1, 1, 1, ........ A000012
W_1: 0, 1, 2, 3, 4, 5, ........ A001477
W_2: -1, 0, 3, 8, 15, 24, ........ A067998
W_3: 0, -2, 2, 18, 52, 110, ........ A121670
W_4: 5, 0, -3, 32, 165, 480, ........
W_n(k), n=0,1,...
k=0: 1, 0, -1, 0, 5, 0, -61, ... A122045
k=1: 1, 1, 0, -2, 0, 16, 0, ... A155585
k=2: 1, 2, 3, 2, -3, 2, 63, ... A119880
k=3: 1, 3, 8, 18, 32, 48, 128, ... A119881
k=4: 1, 4, 15, 52, 165, 484, ........ [Peter Luschny, Jul 07 2009]

Programs

  • Maple
    w := proc(n,x) local v,k,pow,chen; pow := (a,b) -> if a = 0 and b = 0 then 1 else a^b fi; chen := proc(m) if irem(m+1,4) = 0 then RETURN(0) fi; 1/((-1)^iquo(m+1,4) *2^iquo(m,2)) end; add(add((-1)^v*binomial(k,v)*pow(v+x+1,n)*chen(k),v=0..k), k=0..n) end:
    # Coefficients with zeros:
    seq(print(seq(coeff(i!*coeff(series(exp(x*t)*sech(t),t,16),t,i),x,i-n),n=0..i)), i=0..8);
    # Recursion
    W := proc(n,z) option remember; local k,p;
    if n = 0 then 1 else p := irem(n+1,2);
    z^n - p + add(`if`(irem(k,2)=1,0,
    W(k,0)*binomial(n,k)*(power(z,n-k)-p)),k=2..n-1) fi end:
    # Peter Luschny, edited and additions Jul 07 2009, May 13 2010, Oct 24 2011
  • Mathematica
    max = 9; rows = (Reverse[ CoefficientList[ #, x]] & ) /@ CoefficientList[ Series[ Exp[x*t]*Sech[t], {t, 0, max}], t]*Range[0, max]!; par[coefs_] := (p = Partition[ coefs, 2][[All, 1]]; If[ EvenQ[ Length[ coefs]], p, Append[ p, Last[ coefs]]]); Flatten[ par /@ rows] (* Jean-François Alcover, Oct 03 2011, after g.f. *)
    sk[n_, x_] := Sum[Binomial[n, k]*EulerE[k]*x^(n-k), {k, 0, n}]; Table[CoefficientList[sk[n, x], x] // Reverse // Select[#, # =!= 0 &] &, {n, 0, 13}] // Flatten (* Jean-François Alcover, May 21 2013 *)
    Flatten@Table[Binomial[n, 2k] EulerE[2k], {n, 0, 12}, {k, 0, n/2}](* Oliver Seipel, Jan 14 2025 *)
  • Sage
    def A046978(k):
        if k % 4 == 0:
            return 0
        return (-1)**(k // 4)
    def A153641_poly(n, x):
        return expand(add(2**(-(k // 2))*A046978(k+1)*add((-1)**v*binomial(k,v)*(v+x+1)**n for v in (0..k)) for k in (0..n)))
    for n in (0..7): print(A153641_poly(n, x))  # Peter Luschny, Oct 24 2011

Formula

W_n(x) = Sum_{k=0..n}{v=0..k} (-1)^v binomial(k,v)*c_k*(x+v+1)^n where c_k = frac((-1)^(floor(k/4))/2^(floor(k/2))) [4 not div k] (Iverson notation).
From Peter Bala, Jun 10 2009: (Start)
E.g.f.: 2*exp(x*t)*(exp(t)-exp(3*t))/(1-exp(4*t))= 1 + x*t + (x^2-1)*t^2/2! + (x^3-3*x)*t^3/3! + ....
W_n(x) = 1/(2*n+2)*Sum_{k=0..n+1} 1/(k+1)*Sum_{i=0..k} (-1)^i*binomial(k,i)*((x+4*i+3)^(n+1) - (x+4*i+1)^(n+1)).
Fourier series expansion for the generalized Bernoulli polynomials:
B(X;2*n,x) = (-1)^n*(2/Pi)^(2*n)*(2*n)! * {sin(Pi*x/2)/1^(2*n) - sin(3*Pi*x/2)/3^(2*n) + sin(5*Pi*x/2)/5^(2*n) - ...}, valid for 0 <= x <= 1 when n >= 1.
B(X;2*n+1,x) = (-1)^(n+1)*(2/Pi)^(2*n+1)*(2*n+1)! * {cos(Pi*x/2)/1^(2*n+1) - cos(3*Pi*x/2)/3^(2*n+1) + cos(5*Pi*x/2)/5^(2*n+1) - ...}, valid for 0 <= x <= 1 when n >= 1 and for 0 <= x < 1 when n = 0.
(End)
E.g.f.: exp(x*t) * sech(t). - Peter Luschny, Jul 07 2009
O.g.f. as a J-fraction: z/(1-x*z+z^2/(1-x*z+4*z^2/(1-x*z+9*z^2/(1-x*z+...)))) = z + x*z^2 + (x^2-1)*z^3 + (x^3-3*x)*z^4 + .... - Peter Bala, Mar 11 2012
Conjectural o.g.f.: Sum_{n >= 0} (1/2^((n-1)/2))*cos((n+1)*Pi/4)*( Sum_{k = 0..n} (-1)^k*binomial(n,k)/(1 - (k + x)*t) ) = 1 + x*t + (x^2 - 1)*t^2 + (x^3 - 3*x)*t^3 + ... (checked up to O(t^13)), which leads to W_n(x) = Sum_{k = 0..n} 1/2^((k - 1)/2)*cos((k + 1)*Pi/4)*( Sum_{j = 0..k} (-1)^j*binomial(k, j)*(j + x)^n ). - Peter Bala, Oct 03 2016

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

Original entry on oeis.org

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

Views

Author

Roger L. Bagula, Mar 01 2006

Keywords

Comments

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

Examples

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

Crossrefs

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

Programs

Formula

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

Extensions

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

A277950 Triangle read by rows, in which row n gives coefficients in expansion of ((x^n - 1)/(x - 1))^5.

Original entry on oeis.org

1, 1, 5, 10, 10, 5, 1, 1, 5, 15, 30, 45, 51, 45, 30, 15, 5, 1, 1, 5, 15, 35, 65, 101, 135, 155, 155, 135, 101, 65, 35, 15, 5, 1, 1, 5, 15, 35, 70, 121, 185, 255, 320, 365, 381, 365, 320, 255, 185, 121, 70, 35, 15, 5, 1
Offset: 1

Views

Author

Juan Pablo Herrera P., Nov 05 2016

Keywords

Comments

Sum of n-th row is n^5. The n-th row contains 5n-4 entries. Largest coefficients of each row are listed in A077044.
The n-th row is the fifth row of the n-nomial triangle. For example, row 2 (1,5,10,10,5,1) is the fifth row in the binomial triangle.
T(n,k) gives the number of possible ways of randomly selecting k cards from n-1 sets, each with five different playing cards. It is also the number of lattice paths from (0,0) to (5,k) using steps (1,0), (1,1), (1,2), ..., (1,n-1).

Examples

			Triangle starts:
1;
1, 5, 10, 10, 5, 1;
1, 5, 15, 30, 45, 51, 45, 30, 15, 5, 1;
1, 5, 15, 35, 65, 101, 135, 155, 155, 135, 101, 65, 35, 15, 5, 1;
1, 5, 15, 35, 70, 121, 185, 255, 320, 365, 381, 365, 320, 255, 185, 121, 70, 35, 15, 5, 1;
		

Crossrefs

Programs

  • Mathematica
    Table[CoefficientList[Series[((x^n - 1)/(x - 1))^5, {x, 0, 5 n}], x], {n, 10}] // Flatten
  • PARI
    row(n) = Vec(((1 - x^n)/(1 - x))^5); tabf(nn) = for (n=1, nn, print(row(n)));

Formula

T(n,k) = Sum_{i=k-n+1..k} A277949(T(n,i)).
From Juan Pablo Herrera P., Dec 20 2016: (Start)
T(n,k) = A000332(k+4) = (k+4)!/(k!*24) if 0 =< k < n.
T(n,k) = ((k+4)!/k!-5*(k-n+4)!/(k-n)!)/24 if n =< k < 2*n.
T(n,k) = ((k+4)!/k!-5*(k-n+4)!/(k-n)!+10*(k-2*n+4)!/(k-2*n)!)/24 if 2*n =< k < 3*n.
T(n,k) = ((5*n-k-1)!/(5*n-k-5)!-5*(4*n-k-1)!/(4*n-k-5)!)/24 if 3*n-4 =< k < 4*n-4.
T(n,k) = A000332(5*n-k-1) = (5*n-k-1)!/(5*n-k-5)!*24 4*n-4 =< k < 5*n-4. (End)

A277951 Triangle read by rows, in which row n gives coefficients in expansion of ((x^n - 1)/(x - 1))^6.

Original entry on oeis.org

1, 1, 6, 15, 20, 15, 6, 1, 1, 6, 21, 50, 90, 126, 141, 126, 90, 50, 21, 6, 1, 1, 6, 21, 56, 120, 216, 336, 456, 546, 580, 546, 456, 336, 216, 120, 56, 21, 6, 1, 1, 6, 21, 56, 126, 246, 426, 666, 951, 1246, 1506, 1686, 1751, 1686, 1506, 1246, 951, 666, 426, 246, 126, 56, 21, 6, 1
Offset: 1

Views

Author

Juan Pablo Herrera P., Nov 18 2016

Keywords

Comments

Sum of n-th row is n^6. The n-th row contains 6n-5 entries. Largest coefficients of each row are listed in A071816.
The n-th row is the sixth row of the n-nomial triangle. For example, row 2 (1,6,15,20,15,6,1) is the sixth row in the binomial triangle
T(n,k) gives the number of possible ways of randomly selecting k cards from n-1 sets, each with six different playing cards. It is also the number of lattice paths from (0,0) to (6,k) using steps (1,0), (1,1), (1,2), ..., (1,n-1).

Examples

			Triangle starts:
1;
1, 6, 15, 20, 15, 6, 1;
1, 6, 21, 50, 90, 126, 141, 126, 90, 50, 21, 6, 1;
1, 6, 21, 56, 120, 216, 336, 456, 546, 580, 546, 456, 336, 216, 120, 56, 21, 6, 1.
		

Crossrefs

Programs

  • Mathematica
    Table[CoefficientList[Series[((x^n - 1)/(x - 1))^6, {x, 0, 6 n}], x], {n, 10}] // Flatten
  • PARI
    row(n) = Vec(((1 - x^n)/(1 - x))^6);
    tabf(nn) = for (n=1, nn, print(row(n)));

Formula

T(n,k) = Sum_{i=k-n+1..k} A277950(T(n,i))

A266387 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 322560.

Original entry on oeis.org

0, 0, 0, 0, 0, 7, 42, 147, 392, 882, 1764, 3234, 5544, 9009, 14014, 21021, 30576, 43316, 59976, 81396, 108528, 142443, 184338, 235543, 297528, 371910, 460460, 565110, 687960, 831285, 997542, 1189377, 1409632, 1661352, 1947792, 2272424, 2638944, 3051279
Offset: 1

Views

Author

Keywords

Comments

The sequence was discovered by enumerating all orbits of Aut(Z^7) and sorting the orbits as function of the infinity norm of the representative integer lattice points. This sequence is one of the 30 sequences that are obtained by classifying the orbits in a table with the rows being the infinity norm and the columns being the 30 cardinalities (1, 14, 84, 128, 168, 280, 448, 560, 672, 840, 896, 1680, 2240, 2688, 3360, 4480, 5376, 6720, 8960, 13440, 17920, 20160, 26880, 40320, 53760, 80640, 107520, 161280, 322560, 645120) generated by signed permutations of integer lattice points of Z^7.
The continued fraction expansion of this sequence is finite and simplifies to the g.f. 7*x^6/(1-x)^6 (see Mathematica). - Benedict W. J. Irwin, Feb 09 2016

Crossrefs

Other sequences that give the number of orbits of Aut(Z^7) as function of the infinity norm for different cardinalities of these orbits: A000579, A154286, A102860, A002412, A045943, A115067, A008586, A008585, A005843, A001477, A000217.

Programs

  • Mathematica
    Join[{0, 0, 0, 0, 0},Table[Abs[SeriesCoefficient[Series[7/(x+6/(x - 5/2/(x + ContinuedFractionK[If[Mod[k, 2] ==0, (7 + k/2)/(6 + 2 k), ((k + 1)/2 - 5)/(2 (k - 1) +6)], x, {k, 0, 8}]))), {x, Infinity, 101}],2 n + 1]], {n, 0, 50}]] - (* Benedict W. J. Irwin, Feb 09 2016 *)
  • PARI
    concat(vector(5), Vec(7*x^6/(1-x)^6 + O(x^50))) \\ Colin Barker, May 04 2016

Formula

From Colin Barker, Dec 29 2015: (Start)
a(n) = 7*(n-1)*(n-2)*(n-3)*(n-4)*(n-5)/120.
a(n) = 6*a(n-1)-15*a(n-2)+20*a(n-3)-15*a(n-4)+6*a(n-5)-a(n-6) for n>6.
G.f.: 7*x^6 / (1-x)^6.
(End)

A266395 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 161280.

Original entry on oeis.org

0, 0, 0, 0, 15, 75, 225, 525, 1050, 1890, 3150, 4950, 7425, 10725, 15015, 20475, 27300, 35700, 45900, 58140, 72675, 89775, 109725, 132825, 159390, 189750, 224250, 263250, 307125, 356265, 411075, 471975, 539400, 613800, 695640, 785400, 883575, 990675, 1107225
Offset: 1

Views

Author

Keywords

Crossrefs

Number of orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008586, A008585, A005843, A001477, A000217.

Programs

  • PARI
    concat(vector(4), Vec(15*x^5/(1-x)^5 + O(x^50))) \\ Colin Barker, May 05 2016

Formula

From Colin Barker, Dec 29 2015: (Start)
a(n) = 5*(n-1)*(n-2)*(n-3)*(n-4)/8 = 15*A000332(n-1).
a(n) = 5*a(n-1)-10*a(n-2)+10*a(n-3)-5*a(n-4)+a(n-5) for n>5.
G.f.: 15*x^5 / (1-x)^5.
(End)
Showing 1-10 of 13 results. Next