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

A176710 Mix A001021, 2*A001021.

Original entry on oeis.org

1, 2, 12, 24, 144, 288, 1728, 3456, 20736, 41472, 248832, 497664, 2985984, 5971968, 35831808, 71663616, 429981696, 859963392, 5159780352, 10319560704, 61917364224, 123834728448, 743008370688, 1486016741376, 8916100448256, 17832200896512, 106993205379072, 213986410758144
Offset: 0

Views

Author

Paul Curtz, Apr 24 2010

Keywords

Comments

Divides A091137. See A165313.

Programs

  • Mathematica
    LinearRecurrence[{0, 12}, {1, 2}, 30] (* Paolo Xausa, Jan 29 2024 *)

Formula

a(n)= +12*a(n-2). G.f.: ( -1-2*x ) / ( -1+12*x^2 ). [R. J. Mathar, Apr 26 2010]

Extensions

a(21)-a(27) from Paolo Xausa, Jan 29 2024

A000384 Hexagonal numbers: a(n) = n*(2*n-1).

Original entry on oeis.org

0, 1, 6, 15, 28, 45, 66, 91, 120, 153, 190, 231, 276, 325, 378, 435, 496, 561, 630, 703, 780, 861, 946, 1035, 1128, 1225, 1326, 1431, 1540, 1653, 1770, 1891, 2016, 2145, 2278, 2415, 2556, 2701, 2850, 3003, 3160, 3321, 3486, 3655, 3828, 4005, 4186, 4371, 4560
Offset: 0

Views

Author

Keywords

Comments

Number of edges in the join of two complete graphs, each of order n, K_n * K_n. - Roberto E. Martinez II, Jan 07 2002
The power series expansion of the entropy function H(x) = (1+x)log(1+x) + (1-x)log(1-x) has 1/a_i as the coefficient of x^(2i) (the odd terms being zero). - Tommaso Toffoli (tt(AT)bu.edu), May 06 2002
Partial sums of A016813 (4n+1). Also with offset = 0, a(n) = (2n+1)(n+1) = A005408 * A000027 = 2n^2 + 3n + 1, i.e., a(0) = 1. - Jeremy Gardiner, Sep 29 2002
Sequence also gives the greatest semiperimeter of primitive Pythagorean triangles having inradius n-1. Such a triangle has consecutive longer sides, with short leg 2n-1, hypotenuse a(n) - (n-1) = A001844(n), and area (n-1)*a(n) = 6*A000330(n-1). - Lekraj Beedassy, Apr 23 2003
Number of divisors of 12^(n-1), i.e., A000005(A001021(n-1)). - Henry Bottomley, Oct 22 2001
More generally, if p1 and p2 are two arbitrarily chosen distinct primes then a(n) is the number of divisors of (p1^2*p2)^(n-1) or equivalently of any member of A054753^(n-1). - Ant King, Aug 29 2011
Number of standard tableaux of shape (2n-1,1,1) (n>=1). - Emeric Deutsch, May 30 2004
It is well known that for n>0, A014105(n) [0,3,10,21,...] is the first of 2n+1 consecutive integers such that the sum of the squares of the first n+1 such integers is equal to the sum of the squares of the last n; e.g., 10^2 + 11^2 + 12^2 = 13^2 + 14^2.
Less well known is that for n>1, a(n) [0,1,6,15,28,...] is the first of 2n consecutive integers such that sum of the squares of the first n such integers is equal to the sum of the squares of the last n-1 plus n^2; e.g., 15^2 + 16^2 + 17^2 = 19^2 + 20^2 + 3^2. - Charlie Marion, Dec 16 2006
a(n) is also a perfect number A000396 when n is an even superperfect number A061652. - Omar E. Pol, Sep 05 2008
Sequence found by reading the line from 0, in the direction 0, 6, ... and the line from 1, in the direction 1, 15, ..., in the square spiral whose vertices are the generalized hexagonal numbers A000217. - Omar E. Pol, Jan 09 2009
For n>=1, 1/a(n) = Sum_{k=0..2*n-1} ((-1)^(k+1)*binomial(2*n-1,k)*binomial(2*n-1+k,k)*H(k)/(k+1)) with H(k) harmonic number of order k.
The number of possible distinct colorings of any 2 colors chosen from n colors of a square divided into quadrants. - Paul Cleary, Dec 21 2010
Central terms of the triangle in A051173. - Reinhard Zumkeller, Apr 23 2011
For n>0, a(n-1) is the number of triples (w,x,y) with all terms in {0,...,n} and max(|w-x|,|x-y|) = |w-y|. - Clark Kimberling, Jun 12 2012
a(n) is the number of positions of one domino in an even pyramidal board with base 2n. - César Eliud Lozada, Sep 26 2012
Partial sums give A002412. - Omar E. Pol, Jan 12 2013
Let a triangle have T(0,0) = 0 and T(r,c) = |r^2 - c^2|. The sum of the differences of the terms in row(n) and row(n-1) is a(n). - J. M. Bergot, Jun 17 2013
With T_(i+1,i)=a(i+1) and all other elements of the lower triangular matrix T zero, T is the infinitesimal generator for A176230, analogous to A132440 for the Pascal matrix. - Tom Copeland, Dec 11 2013
a(n) is the number of length 2n binary sequences that have exactly two 1's. a(2) = 6 because we have: {0,0,1,1}, {0,1,0,1}, {0,1,1,0}, {1,0,0,1}, {1,0,1,0}, {1,1,0,0}. The ordinary generating function with interpolated zeros is: (x^2 + 3*x^4)/(1-x^2)^3. - Geoffrey Critzer, Jan 02 2014
For n > 0, a(n) is the largest integer k such that k^2 + n^2 is a multiple of k + n. More generally, for m > 0 and n > 0, the largest integer k such that k^(2*m) + n^(2*m) is a multiple of k + n is given by k = 2*n^(2*m) - n. - Derek Orr, Sep 04 2014
Binomial transform of (0, 1, 4, 0, 0, 0, ...) and second partial sum of (0, 1, 4, 4, 4, ...). - Gary W. Adamson, Oct 05 2015
a(n) also gives the dimension of the simple Lie algebras D_n, for n >= 4. - Wolfdieter Lang, Oct 21 2015
For n > 0, a(n) equals the number of compositions of n+11 into n parts avoiding parts 2, 3, 4. - Milan Janjic, Jan 07 2016
Also the number of minimum dominating sets and maximal irredundant sets in the n-cocktail party graph. - Eric W. Weisstein, Jun 29 and Aug 17 2017
As Beedassy's formula shows, this Hexagonal number sequence is the odd bisection of the Triangle number sequence. Both of these sequences are figurative number sequences. For A000384, a(n) can be found by multiplying its triangle number by its hexagonal number. For example let's use the number 153. 153 is said to be the 17th triangle number but is also said to be the 9th hexagonal number. Triangle(17) Hexagonal(9). 17*9=153. Because the Hexagonal number sequence is a subset of the Triangle number sequence, the Hexagonal number sequence will always have both a triangle number and a hexagonal number. n* (2*n-1) because (2*n-1) renders the triangle number. - Bruce J. Nicholson, Nov 05 2017
Also numbers k with the property that in the symmetric representation of sigma(k) the smallest Dyck path has a central valley and the largest Dyck path has a central peak, n >= 1. Thus all hexagonal numbers > 0 have middle divisors. (Cf. A237593.) - Omar E. Pol, Aug 28 2018
k^a(n-1) mod n = 1 for prime n and k=2..n-1. - Joseph M. Shunia, Feb 10 2019
Consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z: a(n+1) gives the semiperimeter of related triangles; A005408, A046092 and A001844 give the X, Y and Z values. - Ralf Steiner, Feb 25 2020
See A002939(n) = 2*a(n) for the corresponding perimeters. - M. F. Hasler, Mar 09 2020
It appears that these are the numbers k with the property that the smallest subpart in the symmetric representation of sigma(k) is 1. - Omar E. Pol, Aug 28 2021
The above conjecture is true. See A280851 for a proof. - Hartmut F. W. Hoft, Feb 02 2022
The n-th hexagonal number equals the sum of the n consecutive integers with the same parity starting at n; for example, 1, 2+4, 3+5+7, 4+6+8+10, etc. In general, the n-th 2k-gonal number is the sum of the n consecutive integers with the same parity starting at (k-2)*n - (k-3). When k = 1 and 2, this result generates the positive integers, A000027, and the squares, A000290, respectively. - Charlie Marion, Mar 02 2022
Conjecture: For n>0, min{k such that there exist subsets A,B of {0,1,2,...,a(n)} such that |A|=|B|=k and A+B={0,1,2,...,2*a(n)}} = 2*n. - Michael Chu, Mar 09 2022

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, pp. 77-78. (In the integral formula on p. 77 a left bracket is missing for the cosine argument.)
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 38.
  • 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. 2.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 53-54, 129-130, 132.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 21.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See pp. 122-123.

Crossrefs

a(n)= A093561(n+1, 2), (4, 1)-Pascal column.
a(n) = A100345(n, n-1) for n>0.
Cf. A002939 (twice a(n): sums of Pythagorean triples (X, Y, Z=Y+1)).
Cf. A280851.

Programs

  • Haskell
    a000384 n = n * (2 * n - 1)
    a000384_list = scanl (+) 0 a016813_list
    -- Reinhard Zumkeller, Dec 16 2012
    
  • Maple
    A000384:=n->n*(2*n-1); seq(A000384(k), k=0..100); # Wesley Ivan Hurt, Sep 27 2013
  • Mathematica
    Table[n*(2 n - 1), {n, 0, 100}] (* Wesley Ivan Hurt, Sep 27 2013 *)
    LinearRecurrence[{3, -3, 1}, {0, 1, 6}, 50] (* Harvey P. Dale, Sep 10 2015 *)
    Join[{0}, Accumulate[Range[1, 312, 4]]] (* Harvey P. Dale, Mar 26 2016 *)
    (* For Mathematica 10.4+ *) Table[PolygonalNumber[RegularPolygon[6], n], {n, 0, 48}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
    PolygonalNumber[6, Range[0, 20]] (* Eric W. Weisstein, Aug 17 2017 *)
    CoefficientList[Series[x*(1 + 3*x)/(1 - x)^3 , {x, 0, 100}], x] (* Stefano Spezia, Sep 02 2018 *)
  • PARI
    a(n)=n*(2*n-1)
    
  • PARI
    a(n) = binomial(2*n,2) \\ Altug Alkan, Oct 06 2015
    
  • 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 + 4, y + 4
    A000384 = aList()
    print([next(A000384) for i in range(49)]) # Peter Luschny, Aug 04 2019

Formula

a(n) = Sum_{k=1..n} tan^2((k - 1/2)*Pi/(2n)). - Ignacio Larrosa Cañestro, Apr 17 2001
E.g.f.: exp(x)*(x+2x^2). - Paul Barry, Jun 09 2003
G.f.: x*(1+3*x)/(1-x)^3. - Simon Plouffe in his 1992 dissertation, dropping the initial zero
a(n) = A000217(2*n-1) = A014105(-n).
a(n) = 4*A000217(n-1) + n. - Lekraj Beedassy, Jun 03 2004
a(n) = right term of M^n * [1,0,0], where M = the 3 X 3 matrix [1,0,0; 1,1,0; 1,4,1]. Example: a(5) = 45 since M^5 *[1,0,0] = [1,5,45]. - Gary W. Adamson, Dec 24 2006
Row sums of triangle A131914. - Gary W. Adamson, Jul 27 2007
Row sums of n-th row, triangle A134234 starting (1, 6, 15, 28, ...). - Gary W. Adamson, Oct 14 2007
Starting with offset 1, = binomial transform of [1, 5, 4, 0, 0, 0, ...]. Also, A004736 * [1, 4, 4, 4, ...]. - Gary W. Adamson, Oct 25 2007
a(n)^2 + (a(n)+1)^2 + ... + (a(n)+n-1)^2 = (a(n)+n+1)^2 + ... + (a(n)+2n-1)^2 + n^2; e.g., 6^2 + 7^2 = 9^2 + 2^2; 28^2 + 29^2 + 30^2 + 31^2 = 33^2 + 34^2 + 35^2 + 4^2. - Charlie Marion, Nov 10 2007
a(n) = binomial(n+1,2) + 3*binomial(n,2).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=0, a(1)=1, a(2)=6. - Jaume Oliver Lafont, Dec 02 2008
a(n) = T(n) + 3*T(n-1), where T(n) is the n-th triangular number. - Vincenzo Librandi, Nov 10 2010
a(n) = a(n-1) + 4*n - 3 (with a(0)=0). - Vincenzo Librandi, Nov 20 2010
a(n) = A007606(A000290(n)). - Reinhard Zumkeller, Feb 12 2011
a(n) = 2*a(n-1) - a(n-2) + 4. - Ant King, Aug 26 2011
a(n+1) = A045896(2*n). - Reinhard Zumkeller, Dec 12 2011
a(2^n) = 2^(2n+1) - 2^n. - Ivan N. Ianakiev, Apr 13 2013
a(n) = binomial(2*n,2). - Gary Detlefs, Jul 28 2013
a(n+1) = A128918(2*n+1). - Reinhard Zumkeller, Oct 13 2013
a(4*a(n)+7*n+1) = a(4*a(n)+7*n) + a(4*n+1). - Vladimir Shevelev, Jan 24 2014
Sum_{n>=1} 1/a(n) = 2*log(2) = 1.38629436111989...= A016627. - Vaclav Kotesovec, Apr 27 2016
Sum_{n>=1} (-1)^n/a(n) = log(2) - Pi/2. - Vaclav Kotesovec, Apr 20 2018
a(n+1) = trinomial(2*n+1, 2) = trinomial(2*n+1, 4*n), for n >= 0, with the trinomial irregular triangle A027907. a(n+1) = (n+1)*(2*n+1) = (1/Pi)*Integral_{x=0..2} (1/sqrt(4 - x^2))*(x^2 - 1)^(2*n+1)*R(4*n-2, x) with the R polynomial coefficients given in A127672. [Comtet, p. 77, the integral formula for q=3, n -> 2*n+1, k = 2, rewritten with x = 2*cos(phi)]. - Wolfdieter Lang, Apr 19 2018
Sum_{n>=1} 1/(a(n))^2 = 2*Pi^2/3-8*log(2) = 1.0345588... = 10*A182448 - A257872. - R. J. Mathar, Sep 12 2019
a(n) = (A005408(n-1) + A046092(n-1) + A001844(n-1))/2. - Ralf Steiner, Feb 27 2020
Product_{n>=2} (1 - 1/a(n)) = 2/3. - Amiram Eldar, Jan 21 2021
a(n) = floor(Sum_{k=(n-1)^2..n^2} sqrt(k)), for n >= 1. - Amrit Awasthi, Jun 13 2021
a(n+1) = A084265(2*n), n>=0. - Hartmut F. W. Hoft, Feb 02 2022
a(n) = A000290(n) + A002378(n-1). - Charles Kusniec, Sep 11 2022

Extensions

Partially edited by Joerg Arndt, Mar 11 2010

A000351 Powers of 5: a(n) = 5^n.

Original entry on oeis.org

1, 5, 25, 125, 625, 3125, 15625, 78125, 390625, 1953125, 9765625, 48828125, 244140625, 1220703125, 6103515625, 30517578125, 152587890625, 762939453125, 3814697265625, 19073486328125, 95367431640625, 476837158203125, 2384185791015625, 11920928955078125
Offset: 0

Views

Author

Keywords

Comments

Same as Pisot sequences E(1, 5), L(1, 5), P(1, 5), T(1, 5). Essentially same as Pisot sequences E(5, 25), L(5, 25), P(5, 25), T(5, 25). See A008776 for definitions of Pisot sequences.
a(n) has leading digit 1 if and only if n = A067497 - 1. - Lekraj Beedassy, Jul 09 2002
With interpolated zeros 0, 1, 0, 5, 0, 25, ... (g.f.: x/(1 - 5*x^2)) second inverse binomial transform of Fibonacci(3n)/Fibonacci(3) (A001076). Binomial transform is A085449. - Paul Barry, Mar 14 2004
Sums of rows of the triangles in A013620 and A038220. - Reinhard Zumkeller, May 14 2006
Sum of coefficients of expansion of (1 + x + x^2 + x^3 + x^4)^n. a(n) is number of compositions of natural numbers into n parts less than 5. a(2) = 25 there are 25 compositions of natural numbers into 2 parts less than 5. - Adi Dani, Jun 22 2011
The compositions of n in which each natural number is colored by one of p different colors are called p-colored compositions of n. For n >= 1, a(n) equals the number of 5-colored compositions of n such that no adjacent parts have the same color. - Milan Janjic, Nov 17 2011
Numbers n such that sigma(5n) = 5n + sigma(n). In fact we have this theorem: p is a prime if and only if all solutions of the equation sigma(p*x) = p*x + sigma(x) are powers of p. - Jahangeer Kholdi, Nov 23 2013
From Doug Bell, Jun 22 2015: (Start)
Empirical observation: Where n is an odd multiple of 3, let x = (a(n) + 1)/9 and let y be the decimal expansion of x/a(n); then y*(x+1)/x + 1 = y rotated to the left.
Example:
a(3) = 125;
x = (125 + 1)/9 = 14;
y = 112, which is the decimal expansion of 14/125 = 0.112;
112*(14 + 1)/14 + 1 = 121 = 112 rotated to the left.
(End)
a(n) is the number of n-digit integers that contain only odd digits (A014261). - Bernard Schott, Nov 12 2022
Number of pyramids in the Sierpinski fractal square-based pyramid at the n-th step, while A279511 gives the corresponding number of vertices (see IREM link with drawings). - Bernard Schott, Nov 29 2022

References

  • 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. A009969 (even bisection), A013710 (odd bisection), A005054 (first differences), A003463 (partial sums).
Sierpinski fractal square-based pyramid: A020858 (Hausdorff dimension), A279511 (number of vertices), this sequence (number of pyramids).

Programs

Formula

a(n) = 5^n.
a(0) = 1; a(n) = 5*a(n-1) for n > 0.
G.f.: 1/(1 - 5*x).
E.g.f.: exp(5*x).
a(n) = A006495(n)^2 + A006496(n)^2.
a(n) = A159991(n) / A001021(n). - Reinhard Zumkeller, May 02 2009
From Bernard Schott, Nov 12 2022: (Start)
Sum_{n>=0} 1/a(n) = 5/4.
Sum_{n>=0} (-1)^n/a(n) = 5/6. (End)
a(n) = Sum_{k=0..n} C(2*n+1,n-k)*A000045(2*k+1). - Vladimir Kruchinin, Jan 14 2025

A001022 Powers of 13.

Original entry on oeis.org

1, 13, 169, 2197, 28561, 371293, 4826809, 62748517, 815730721, 10604499373, 137858491849, 1792160394037, 23298085122481, 302875106592253, 3937376385699289, 51185893014090757, 665416609183179841, 8650415919381337933, 112455406951957393129, 1461920290375446110677, 19004963774880799438801
Offset: 0

Views

Author

Keywords

Comments

Same as Pisot sequences E(1, 13), L(1, 13), P(1, 13), T(1, 13). Essentially same as Pisot sequences E(13, 169), L(13, 169), P(13, 169), T(13, 169). See A008776 for definitions of Pisot sequences.
The compositions of n in which each natural number is colored by one of p different colors are called p-colored compositions of n. For n>=1, a(n) equals the number of 13-colored compositions of n such that no adjacent parts have the same color. - Milan Janjic, Nov 17 2011
Numbers n such that sigma(13*n) = 13*n+sigma(n). - Jahangeer Kholdi, Nov 23 2013

Examples

			For the fifth formula: a(7) = 1*1 + 12*7 + 144*21 + 1728*35 + 20736*35 + 248832*21 + 2985984*7 + 35831808*1 = 62748517. - _Bruno Berselli_, Aug 06 2013
		

References

  • 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. A001021.

Programs

Formula

G.f.: 1/(1-13*x).
E.g.f.: exp(13*x).
a(n) = 13^n. - Vincenzo Librandi, Nov 21 2010
a(n) = 13*a(n-1) n > 0, a(0)=1. - Vincenzo Librandi, Nov 21 2010
a(n) = Sum_{k=0..n} A001021(k)*binomial(n,k). It is well known that Sum_{k=0..n} (h-1)^k*binomial(n,k) = h^n. - Bruno Berselli, Aug 06 2013

Extensions

More terms from James Sellers, Sep 19 2000

A016125 Expansion of 1/((1-x)*(1-12*x)).

Original entry on oeis.org

1, 13, 157, 1885, 22621, 271453, 3257437, 39089245, 469070941, 5628851293, 67546215517, 810554586205, 9726655034461, 116719860413533, 1400638324962397, 16807659899548765, 201691918794585181
Offset: 0

Views

Author

Keywords

Comments

Let A be the Hessenberg matrix of the order n, defined by: A[1,j]=1, A[i,i]:=12, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)=det(A). - Milan Janjic, Feb 21 2010
Let A be the Hessenberg matrix of the order n, defined by: A[1,j]=1, A[i,i]:=13, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=2, a(n-2)=(-1)^n*charpoly(A,1). - Milan Janjic, Feb 21 2010
Numbers that are repunits in duodecimal representation. - Reinhard Zumkeller, Dec 12 2012
a(n) is the total number of holes in a certain box fractal (start with 12 boxes, 1 hole) after n iterations. See illustration in links. - Kival Ngaokrajang, Jan 28 2015

Examples

			For n=5, a(5) = 1*6 + 11*15 + 121*20 + 1331*15 + 14641*6 + 161051*1 = 271453. - _Bruno Berselli_, Nov 11 2015
		

Crossrefs

Programs

Formula

a(n) = (12^(n+1) - 1)/11.
a(n) = 12*a(n-1)+1 for n>0, a(0)=1. - Vincenzo Librandi, Nov 19 2010
a(n) = Sum_{i=0...n} 11^i*binomial(n+1,n-i). - Bruno Berselli, Nov 11 2015
E.g.f.: exp(x)*(12*exp(11*x) - 1)/11. - Stefano Spezia, Mar 11 2023

A091030 Partial sums of powers of 13 (A001022).

Original entry on oeis.org

1, 14, 183, 2380, 30941, 402234, 5229043, 67977560, 883708281, 11488207654, 149346699503, 1941507093540, 25239592216021, 328114698808274, 4265491084507563, 55451384098598320, 720867993281778161
Offset: 1

Views

Author

Wolfdieter Lang, Jan 23 2004

Keywords

Comments

13^a(n) is highest power of 13 dividing (13^n)!.
For analogs with primes 2, 3, 5, 7 and 11 see A000225, A003462, A003463, A023000 and A016123 respectively.
Let A be the Hessenberg matrix of the order n, defined by: A[1,j]=1,A[i,i]:=13, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n)=det(A). - Milan Janjic, Feb 21 2010
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=14, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n)=(-1)^(n)*charpoly(A,1). - Milan Janjic, Feb 21 2010

Examples

			For n=6, a(6) = 1*6 + 12*15 + 144*20 + 1728*15 + 20736*6 + 248832*1 = 402234. - _Bruno Berselli_, Nov 12 2015
		

Crossrefs

Programs

Formula

G.f.: x/((1-13*x)*(1-x)) = (1/(1-13*x) - 1/(1-x))/12.
a(n) = Sum_{k=0..n-1} 13^k = (13^n-1)/12.
a(n) = 13*a(n-1)+1 for n>1, a(1)=1. - Vincenzo Librandi, Feb 05 2011
a(n) = Sum_{k=0...n-1} 12^k*binomial(n,n-1-k). - Bruno Berselli, Nov 12 2015
E.g.f.: exp(x)*(exp(12*x) - 1)/12. - Stefano Spezia, Mar 11 2023

A181818 Products of superprimorials (A006939).

Original entry on oeis.org

1, 2, 4, 8, 12, 16, 24, 32, 48, 64, 96, 128, 144, 192, 256, 288, 360, 384, 512, 576, 720, 768, 1024, 1152, 1440, 1536, 1728, 2048, 2304, 2880, 3072, 3456, 4096, 4320, 4608, 5760, 6144, 6912, 8192, 8640, 9216, 11520, 12288, 13824, 16384, 17280, 18432, 20736, 23040, 24576, 27648, 32768
Offset: 1

Views

Author

Matthew Vandermast, Nov 30 2010

Keywords

Comments

Sorted list of positive integers with a factorization Product p(i)^e(i) such that (e(1) - e(2)) >= (e(2) - e(3)) >= ... >= (e(k-1) - e(k)) >= e(k), with k = A001221(n), and p(k) = A006530(n) = A000040(k), i.e., the prime factors p(1) .. p(k) must be consecutive primes from 2 onward. - Comment clarified by Antti Karttunen, Apr 28 2022
Subsequence of A025487. A025487(n) belongs to this sequence iff A181815(n) is a member of A025487.
If prime signatures are considered as partitions, these are the members of A025487 whose prime signature is conjugate to the prime signature of a member of A182863. - Matthew Vandermast, May 20 2012

Examples

			2, 12, and 360 are all superprimorials (i.e., members of A006939). Therefore, 2*2*12*360 = 17280 is included in the sequence.
From _Gus Wiseman_, Aug 12 2020 (Start):
The sequence of factorizations (which are unique) begins:
    1 = empty product
    2 = 2
    4 = 2*2
    8 = 2*2*2
   12 = 12
   16 = 2*2*2*2
   24 = 2*12
   32 = 2*2*2*2*2
   48 = 2*2*12
   64 = 2*2*2*2*2*2
   96 = 2*2*2*12
  128 = 2*2*2*2*2*2*2
  144 = 12*12
  192 = 2*2*2*2*12
  256 = 2*2*2*2*2*2*2*2
(End)
		

Crossrefs

A181817 rearranged in numerical order. Also includes all members of A000079, A001021, A006939, A009968, A009992, A066120, A166475, A167448, A181813, A181814, A181816, A182763.
Subsequence of A025487, A055932, A087980, A130091, A181824.
A001013 is the version for factorials.
A336426 is the complement.
A336496 is the version for superfactorials.
A001055 counts factorizations.
A006939 lists superprimorials or Chernoff numbers.
A317829 counts factorizations of superprimorials.
Cf. A022915, A076954, A304686, A325368, A336419, A336420, A336421, A353518 (characteristic function).

Programs

  • Mathematica
    Select[Range[100],PrimePi[First/@If[#==1,{}, FactorInteger[#]]]==Range[ PrimeNu[#]]&&LessEqual@@Differences[ Append[Last/@FactorInteger[#],0]]&] (* Gus Wiseman, Aug 12 2020 *)
  • PARI
    firstdiffs0forward(vec) = { my(v=vector(#vec)); for(n=1,#v,v[n] = vec[n]-if(#v==n,0,vec[1+n])); (v); };
    A353518(n) = if(1==n,1,my(f=factor(n), len=#f~); if(primepi(f[len,1])!=len, return(0), my(diffs=firstdiffs0forward(f[,2])); for(i=1,#diffs-1,if(diffs[i+1]>diffs[i],return(0))); (1)));
    isA181818(n) = A353518(n); \\ Antti Karttunen, Apr 28 2022

A159991 Powers of 60: a(n) = 60^n.

Original entry on oeis.org

1, 60, 3600, 216000, 12960000, 777600000, 46656000000, 2799360000000, 167961600000000, 10077696000000000, 604661760000000000, 36279705600000000000, 2176782336000000000000, 130606940160000000000000, 7836416409600000000000000, 470184984576000000000000000
Offset: 0

Views

Author

Reinhard Zumkeller, May 01 2009

Keywords

Examples

			G.f. = 1 + 60*x + 3600*x^2 + 216000*x^3 + 12960000*x^4 + 77600000*x^5 + ... - _Michael Somos_, Jan 01 2019
		

Crossrefs

Programs

Formula

a(n) = A000400(n)*A011557(n) = A000351(n)*A001021(n) = A000302(n)*A001024(n) = A000244(n)*A009964(n). (Corrected by Robert B Fowler, Jan 25 2023)
From Muniru A Asiru, Nov 21 2018: (Start)
a(n) = 60^n.
a(n) = 60*a(n-1) for n > 0, a(0) = 1.
G.f.: 1/(1-60*x).
E.g.f: exp(60*x). (End)
a(n) = 1/a(-n) for all n in Z. - Michael Somos, Jan 01 2019

A003992 Square array read by upwards antidiagonals: T(n,k) = n^k for n >= 0, k >= 0.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 4, 1, 0, 1, 4, 9, 8, 1, 0, 1, 5, 16, 27, 16, 1, 0, 1, 6, 25, 64, 81, 32, 1, 0, 1, 7, 36, 125, 256, 243, 64, 1, 0, 1, 8, 49, 216, 625, 1024, 729, 128, 1, 0, 1, 9, 64, 343, 1296, 3125, 4096, 2187, 256, 1, 0, 1, 10, 81, 512, 2401, 7776, 15625, 16384, 6561, 512, 1, 0
Offset: 0

Views

Author

Keywords

Comments

If the array is transposed, T(n,k) is the number of oriented rows of n colors using up to k different colors. The formula would be T(n,k) = [n==0] + [n>0]*k^n. The generating function for column k would be 1/(1-k*x). For T(3,2)=8, the rows are AAA, AAB, ABA, ABB, BAA, BAB, BBA, and BBB. - Robert A. Russell, Nov 08 2018
T(n,k) is the number of multichains of length n from {} to [k] in the Boolean lattice B_k. - Geoffrey Critzer, Apr 03 2020

Examples

			Rows begin:
[1, 0,  0,   0,    0,     0,      0,      0, ...],
[1, 1,  1,   1,    1,     1,      1,      1, ...],
[1, 2,  4,   8,   16,    32,     64,    128, ...],
[1, 3,  9,  27,   81,   243,    729,   2187, ...],
[1, 4, 16,  64,  256,  1024,   4096,  16384, ...],
[1, 5, 25, 125,  625,  3125,  15625,  78125, ...],
[1, 6, 36, 216, 1296,  7776,  46656, 279936, ...],
[1, 7, 49, 343, 2401, 16807, 117649, 823543, ...], ...
		

Crossrefs

Main diagonal is A000312. Other diagonals include A000169, A007778, A000272, A008788. Antidiagonal sums are in A026898.
Cf. A099555.
Transpose is A004248. See A051128, A095884, A009999 for other versions.
Cf. A277504 (unoriented), A293500 (chiral).

Programs

  • Magma
    [[(n-k)^k: k in [0..n]]: n in [0..10]]; // G. C. Greubel, Nov 08 2018
  • Mathematica
    Table[If[k == 0, 1, (n - k)^k], {n, 0, 11}, {k, 0, n}]//Flatten
  • PARI
    T(n,k) = (n-k)^k \\ Charles R Greathouse IV, Feb 07 2017
    

Formula

E.g.f.: Sum T(n,k)*x^n*y^k/k! = 1/(1-x*exp(y)). - Paul D. Hanna, Oct 22 2004
E.g.f.: Sum T(n,k)*x^n/n!*y^k/k! = e^(x*e^y). - Franklin T. Adams-Watters, Jun 23 2006

Extensions

More terms from David W. Wilson
Edited by Paul D. Hanna, Oct 22 2004

A178248 a(n) = 12^n + 1.

Original entry on oeis.org

2, 13, 145, 1729, 20737, 248833, 2985985, 35831809, 429981697, 5159780353, 61917364225, 743008370689, 8916100448257, 106993205379073, 1283918464548865, 15407021574586369, 184884258895036417, 2218611106740436993
Offset: 0

Views

Author

Stuart Clary, Dec 20 2010

Keywords

Comments

Prime factors of a(n) are in the Cunningham Project.

Examples

			a(3) = 12^3 + 1 = 1729.
		

Crossrefs

Programs

Formula

O.g.f.: (2-13*x)/(1-13*x+12*x^2) = (2-13*x)/((1-x)(1-12*x)).
E.g.f.: exp(x) + exp(12*x). - Stefano Spezia, Mar 20 2023
From Elmo R. Oliveira, Dec 15 2023: (Start)
a(n) = 12*a(n-1) - 11 for n>0.
a(n) = 13*a(n-1) - 12*a(n-2) for n>1.
a(n) = A001021(n)+1 = A024140(n)+2.
a(n) = (11*A016125(n) + 13)/12. (End)
Showing 1-10 of 58 results. Next