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

A000567 Octagonal numbers: n*(3*n-2). Also called star numbers.

Original entry on oeis.org

0, 1, 8, 21, 40, 65, 96, 133, 176, 225, 280, 341, 408, 481, 560, 645, 736, 833, 936, 1045, 1160, 1281, 1408, 1541, 1680, 1825, 1976, 2133, 2296, 2465, 2640, 2821, 3008, 3201, 3400, 3605, 3816, 4033, 4256, 4485, 4720, 4961, 5208, 5461
Offset: 0

Views

Author

Keywords

Comments

From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0,1,....
The spiral begins:
.
85--84--83--82--81--80
/ \
86 56--55--54--53--52 79
/ / \ \
87 57 33--32--31--30 51 78
/ / / \ \ \
88 58 34 16--15--14 29 50 77
/ / / / \ \ \ \
89 59 35 17 5---4 13 28 49 76
/ / / / / \ \ \ \ \
90 60 36 18 6 0 3 12 27 48 75
/ / / / / / / / / / /
91 61 37 19 7 1---2 11 26 47 74
\ \ \ \ \ . / / / /
92 62 38 20 8---9--10 25 46 73
\ \ \ \ . / / /
93 63 39 21--22--23--24 45 72
\ \ \ . / /
94 64 40--41--42--43--44 71
\ \ . /
95 65--66--67--68--69--70
\ .
96
.
(End)
From Lekraj Beedassy, Oct 02 2003: (Start)
Also the number of distinct three-cell blocks that may be removed out of A000217(n+1) square cells arranged in a stepping triangular array of side (n+1). A 5-layer triangular array of square cells, for instance, has vertices outlined thus:
x x
x x x
x x x x
x x x x x
x x x x x x
x x x x x x (End)
First derivative at n of A045991. - Ross La Haye, Oct 23 2004
Starting from n=1, the sequence corresponds to the Wiener index of K_{n,n} (the complete bipartite graph wherein each independent set has n vertices). - Kailasam Viswanathan Iyer, Mar 11 2009
Number of divisors of 24^(n-1) for n > 0 (cf A009968). - J. Lowell, Aug 30 2008
a(n) = A001399(6n-5), number of partitions of 6*n - 5 into parts < 4. For example a(2)=8 and partitions of 6*2 - 5 = 7 into parts < 4 are: [1,1,1,1,1,1,1], [1,1,1,1,1,2],[1,1,1,1,3], [1,1,1,2,2], [1,1,2,3], [1,2,2,2], [1,3,3], [2,2,3]. - Adi Dani, Jun 07 2011
Also, sequence found by reading the line from 0 in the direction 0, 8, ..., and the parallel line from 1 in the direction 1, 21, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. - Omar E. Pol, Sep 10 2011
Partial sums give A002414. - Omar E. Pol, Jan 12 2013
Generate a Pythagorean triple using Euclid's formula with (n, n-1) to give A,B,C. a(n) = B + (A + C)/2. - J. M. Bergot, Jul 13 2013
The number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 773", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
For n >= 1, the continued fraction expansion of sqrt(27*a(n)) is [9n-4; {1, 2n-2, 3, 2n-2, 1, 18n-8}]. For n=1, this collapses to [5; {5, 10}]. - Magus K. Chu, Oct 10 2022
a(n)*a(n+1) + 1 = (3n^2 + n - 1)^2. In general, a(n)*a(n+k) + k^2 = (3n^2 + (3k-2)n - k)^2. - Charlie Marion, May 23 2023

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
  • 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. 1.
  • 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, pages 19-20.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 123.

Crossrefs

Cf. A014641, A014642, A014793, A014794, A001835, A016777, A045944, A093563 ((6, 1) Pascal, column m=2). A016921 (differences).
Cf. A005408 (the odd numbers).

Programs

  • GAP
    List([0..50], n -> n*(3*n-2)); # G. C. Greubel, Nov 15 2018
    
  • Haskell
    a000567 n = n * (3 * n - 2)  -- Reinhard Zumkeller, Dec 20 2012
    
  • Magma
    [n*(3*n-2) : n in [0..50]]; // Wesley Ivan Hurt, Oct 10 2021
  • Maple
    A000567 := proc(n)
        n*(3*n-2) ;
    end proc:
    seq(A000567(n), n=1..50) ;
  • Mathematica
    Table[n (3 n - 2), {n, 0, 50}] (* Harvey P. Dale, May 06 2012 *)
    Table[PolygonalNumber[RegularPolygon[8], n], {n, 0, 43}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
    PolygonalNumber[8, Range[0, 20]] (* Eric W. Weisstein, Sep 07 2017 *)
    LinearRecurrence[{3, -3, 1}, {1, 8, 21}, {0, 20}] (* Eric W. Weisstein, Sep 07 2017 *)
  • PARI
    a(n)=n*(3*n-2) \\ Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    vector(50, n, n--; n*(3*n-2)) \\ G. C. Greubel, Nov 15 2018
    
  • 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 + 6, y + 6
    A000567 = aList()
    print([next(A000567) for i in range(49)]) # Peter Luschny, Aug 04 2019
    
  • Python
    [n*(3*n-2) for n in range(50)] # Gennady Eremin, Mar 10 2022
    
  • Sage
    [n*(3*n-2) for n in range(50)] # G. C. Greubel, Nov 15 2018
    

Formula

a(n) = n*(3*n-2).
a(n) = (3n-2)*(3n-1)*(3n)/((3n-1) + (3n-2) + (3n)), i.e., (the product of three consecutive numbers)/(their sum). a(1) = 1*2*3/(1+2+3), a(2) = 4*5*6/(4+5+6), etc. - Amarnath Murthy, Aug 29 2002
E.g.f.: exp(x)*(x+3*x^2). - Paul Barry, Jul 23 2003
G.f.: x*(1+5*x)/(1-x)^3. Simon Plouffe in his 1992 dissertation
a(n) = Sum_{k=1..n} (5*n - 4*k). - Paul Barry, Sep 06 2005
a(n) = n + 6*A000217(n-1). - Floor van Lamoen, Oct 14 2005
a(n) = C(n+1,2) + 5*C(n,2).
Starting (1, 8, 21, 40, 65, ...) = binomial transform of [1, 7, 6, 0, 0, 0, ...]. - Gary W. Adamson, Apr 30 2008
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=0, a(1)=1, a(2)=8. - Jaume Oliver Lafont, Dec 02 2008
a(n) = A000578(n) - A007531(n). - Reinhard Zumkeller, Sep 18 2009
a(n) = a(n-1) + 6*n - 5 (with a(0)=0). - Vincenzo Librandi, Nov 20 2010
a(n) = 2*a(n-1) - a(n-2) + 6. - Ant King, Sep 01 2011
a(n) = A000217(n) + 5*A000217(n-1). - Vincenzo Librandi, Nov 20 2010
a(n) = (A185212(n) - 1) / 4. - Reinhard Zumkeller, Dec 20 2012
a(n) = A174709(6n). - Philippe Deléham, Mar 26 2013
a(n) = (2*n-1)^2 - (n-1)^2. - Ivan N. Ianakiev, Apr 10 2013
a(6*a(n) + 16*n + 1) = a(6*a(n) + 16*n) + a(6*n + 1). - Vladimir Shevelev, Jan 24 2014
a(0) = 0, a(n) = Sum_{k=0..n-1} A005408(A051162(n-1,k)), n >= 1. - L. Edson Jeffery, Jul 28 2014
Sum_{n>=1} 1/a(n) = (sqrt(3)*Pi + 9*log(3))/12 = 1.2774090575596367311949534921... . - Vaclav Kotesovec, Apr 27 2016
From Ilya Gutkovskiy, Jul 29 2016: (Start)
Inverse binomial transform of A084857.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(2*sqrt(3)) = A093766. (End)
a(n) = n * A016777(n-1) = A053755(n) - A000290(n+1). - Bruce J. Nicholson, Aug 10 2017
Product_{n>=2} (1 - 1/a(n)) = 3/4. - Amiram Eldar, Jan 21 2021
P(4k+4,n) = ((k+1)*n - k)^2 - (k*n - k)^2 where P(m,n) is the n-th m-gonal number (a generalization of the Apr 10 2013 formula, a(n) = (2*n-1)^2 - (n-1)^2). - Charlie Marion, Oct 07 2021
From Leo Tavares, Oct 31 2021: (Start)
a(n) = A000290(n) + 4*A000217(n-1). See Square Rays illustration.
a(n) = A000290(n) + A046092(n-1)
a(n) = A000384(n) + 2*A000217(n-1). See Twin Rectangular Rays illustration.
a(n) = A000384(n) + A002378(n-1)
a(n) = A003154(n) - A045944(n-1). See Star Rows illustration. (End)

Extensions

Incorrect example removed by Joerg Arndt, Mar 11 2010

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

A275062 Number A(n,k) of permutations p of [n] such that p(i)-i is a multiple of k for all i in [n]; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 6, 1, 1, 1, 1, 2, 24, 1, 1, 1, 1, 1, 4, 120, 1, 1, 1, 1, 1, 2, 12, 720, 1, 1, 1, 1, 1, 1, 4, 36, 5040, 1, 1, 1, 1, 1, 1, 2, 8, 144, 40320, 1, 1, 1, 1, 1, 1, 1, 4, 24, 576, 362880, 1, 1, 1, 1, 1, 1, 1, 2, 8, 72, 2880, 3628800, 1
Offset: 0

Views

Author

Alois P. Heinz, Jul 15 2016

Keywords

Examples

			A(5,0) = A(5,5) = 1: 12345.
A(5,1) = 5! = 120: all permutations of {1,2,3,4,5}.
A(5,2) = 12: 12345, 12543, 14325, 14523, 32145, 32541, 34125, 34521, 52143, 52341, 54123, 54321.
A(5,3) = 4: 12345, 15342, 42315, 45312.
A(5,4) = 2: 12345, 52341.
A(7,4) = 8: 1234567, 1274563, 1634527, 1674523, 5234167, 5274163, 5634127, 5674123.
Square array A(n,k) begins:
  1,       1,     1,   1,   1,  1,  1, 1, 1, 1, 1, ...
  1,       1,     1,   1,   1,  1,  1, 1, 1, 1, 1, ...
  1,       2,     1,   1,   1,  1,  1, 1, 1, 1, 1, ...
  1,       6,     2,   1,   1,  1,  1, 1, 1, 1, 1, ...
  1,      24,     4,   2,   1,  1,  1, 1, 1, 1, 1, ...
  1,     120,    12,   4,   2,  1,  1, 1, 1, 1, 1, ...
  1,     720,    36,   8,   4,  2,  1, 1, 1, 1, 1, ...
  1,    5040,   144,  24,   8,  4,  2, 1, 1, 1, 1, ...
  1,   40320,   576,  72,  16,  8,  4, 2, 1, 1, 1, ...
  1,  362880,  2880, 216,  48, 16,  8, 4, 2, 1, 1, ...
  1, 3628800, 14400, 864, 144, 32, 16, 8, 4, 2, 1, ...
		

Crossrefs

A(k*n,n) for k=1..4 give: A000012, A000079, A000400, A009968.
Cf. A225816.

Programs

  • Maple
    A:= (n, k)-> mul(floor((n+i)/k)!, i=0..k-1):
    seq(seq(A(n, d-n), n=0..d), d=0..14);
  • Mathematica
    A[n_, k_] := Product[Floor[(n+i)/k]!, {i, 0, k-1}];
    Table[A[n, d-n], {d, 0, 14}, {n, 0, d}] // Flatten (* Jean-François Alcover, May 26 2019, from Maple *)

Formula

A(n,k) = Product_{i=0..k-1} floor((n+i)/k)!.
A(k*n,k) = (n!)^k = A225816(k,n).
For k > 0, A(n, k) ~ (2*Pi*n)^((k - 1)/2) * n! / k^(n + k/2). - Vaclav Kotesovec, Oct 02 2018

A009992 Powers of 48: a(n) = 48^n.

Original entry on oeis.org

1, 48, 2304, 110592, 5308416, 254803968, 12230590464, 587068342272, 28179280429056, 1352605460594688, 64925062108545024, 3116402981210161152, 149587343098087735296, 7180192468708211294208, 344649238497994142121984, 16543163447903718821855232
Offset: 0

Views

Author

Keywords

Comments

Same as Pisot sequences E(1, 48), L(1, 48), P(1, 48), T(1, 48). Essentially same as Pisot sequences E(48, 2304), L(48, 2304), P(48, 2304), T(48, 2304). See A008776 for definitions of Pisot sequences.
If X_1, X_2, ..., X_n is a partition of the set {1,2,...,2*n} into blocks of size 2 then, for n>=1, a(n) is equal to the number of functions f : {1,2,..., 2*n}->{1,2,3,4,5,6,7} such that for fixed y_1,y_2,...,y_n in {1,2,3,4,5,6,7} we have f(X_i)<>{y_i}, (i=1,2,...,n). - Milan Janjic, May 24 2007
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 48-colored compositions of n such that no adjacent parts have the same color. - Milan Janjic, Nov 17 2011

Crossrefs

Cf. A001018 (powers of 8), ..., A001029 (powers of 19), A009964 (powers of 20), ..., A009991 (powers of 47), A087752 (powers of 49).
Cf. A000079 (2^n), A000244 (3^n), A000302 (4^n), A000400 (6^n), A001018 (8^n), A001021 (12^n), A001025 (16^n), A009968 (24^n).

Programs

Formula

G.f.: 1/(1-48*x). - Philippe Deléham, Nov 24 2008
a(n) = 48^n; a(n) = 48*a(n-1), a(0)=1. - Vincenzo Librandi, Nov 21 2010
E.g.f.: exp(48*x). - Muniru A Asiru, Nov 21 2018

Extensions

Edited by M. F. Hasler, Apr 19 2015

A322583 a(n) is the number of factorizations of n into factorial numbers larger than one; a(1) = 1.

Original entry on oeis.org

1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Antti Karttunen, Dec 25 2018

Keywords

Comments

Values 1..7 occur for the first time at n = 1, 24, 576, 13824, 69120, 414720, 1658880.
In range 1..69120 differs from A034876 only at positions n = 1, 2, 9, 10 and 16.

Examples

			a(4) = 1 because 4 = 2! * 2!.
a(24) = 2 because 24 = 4! = 3! * 2! * 2!.
a(576) = 3 because 576 = 4! * 4! = 4! * 3! * 2! * 2! = (3!)^2 * (2!)^4.
a(13824) = 4 because 13824 = (4!)^3 = (4!)^2 * 3! * (2!)^2 = 4! * (3!)^2 * (2!)^4 = (3!)^3 * (2!)^6.
a(69120) = 5 because 69120 = 6! * 4! * 2! * 2! = 6! * 3! * 2! * 2! * 2! * 2! = 5! * 3! * 3! * 2! * 2! * 2! * 2! = 5! * 4! * 4! = 5! * 4! * 3! * 2! * 2!. Note that 69120 < 24^4 = 331776.
		

Crossrefs

Cf. A000142, A001013 (positions of nonzero terms), A009968, A034876, A034968.

Programs

  • PARI
    A322583(n, m=n) = if(1==n, 1, my(s=0); for(i=2, oo, my(f=i!); if(f>m, return(s)); if(!(n%f), s += A322583(n/f, f))));
    
  • PARI
    A034968(n) = { my(s=0, b=2, d); while(n, d = (n%b); s += d; n = (n-d)/b; b++); (s); };
    A322583(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&(1==A034968(d)), s += A322583(n/d, d))); (s));

Formula

a(24^n) = a(A009968(n)) = n+1.

A225816 Square array A(n,k) = (k!)^n, n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 6, 4, 1, 1, 1, 24, 36, 8, 1, 1, 1, 120, 576, 216, 16, 1, 1, 1, 720, 14400, 13824, 1296, 32, 1, 1, 1, 5040, 518400, 1728000, 331776, 7776, 64, 1, 1, 1, 40320, 25401600, 373248000, 207360000, 7962624, 46656, 128, 1, 1
Offset: 0

Views

Author

Alois P. Heinz, Jul 29 2013

Keywords

Comments

A(n,k) is the determinant of the k X k matrix M = [Stirling2(n+i,j)] for 1<=i,j<=k. A(2,3) = det([1,3,1; 1,7,6; 1,15,25]) = 36.
A(n,k) is the determinant of the symmetric k X k matrix M = [sigma_n(gcd(i,j))] for 1<=i,j<=k. A(2,3) = det([1,1,1; 1,5,1; 1,1,10]) = 36.
A(n,k) is (-1)^(n*k) times the determinant of the n X n matrix M = [Stirling1(k+i,j)] for 1<=i,j<=n. A(2,3) = (-1)^(2+3) * det([-6,11; 24,-50]) = 36.
A(n,k) is the number of lattice paths from {n}^k to {0}^k using steps that decrement one component by 1 such that for each point (p_1,p_2,...,p_k) we have abs(p_i-p_j) <= 1 for 1<=i,j<=k. A(2,3) = 36:
(1,2,2)-(1,1,2) (0,1,1)-(0,0,1)
/ X \ / X \
(2,2,2)-(2,1,2) (1,2,1)-(1,1,1)-(1,0,1) (0,1,0)-(0,0,0).
\ X / \ X /
(2,2,1) (2,1,1) (1,1,0) (1,0,0)
A(n,k) is the number of set partitions of [k*(n+1)] into k blocks of size n+1 such that the elements of each block are distinct mod n+1. A(2,3) = 36: 123|456|789, 126|345|789, ..., 189|234|567, 189|246|357.

Examples

			Square array A(n,k) begins:
  1, 1,  1,    1,       1,           1, ...
  1, 1,  2,    6,      24,         120, ...
  1, 1,  4,   36,     576,       14400, ...
  1, 1,  8,  216,   13824,     1728000, ...
  1, 1, 16, 1296,  331776,   207360000, ...
  1, 1, 32, 7776, 7962624, 24883200000, ...
		

Crossrefs

Columns k=0+1, 2-4 give: A000012, A000079, A000400, A009968.
Rows n=0-4 give: A000012, A000142, A001044, A000442, A134375.
Main diagonal gives: A036740.

Programs

  • Maple
    A:= (n, k)-> k!^n:
    seq(seq(A(n,d-n), n=0..d), d=0..12);

Formula

A(n,k) = (k!)^n.
A(n,k) = k^n * A(n,k-1) for k>0, A(n,0) = 1.
A(n,k) = k! * A(n-1,k) for n>0, A(0,k) = 1.
G.f. of column k: 1/(1-k!*x).

A090215 A generalization of triangles A071951 (Legendre-Stirling) and A089504.

Original entry on oeis.org

1, 24, 1, 576, 144, 1, 13824, 17856, 504, 1, 331776, 2156544, 199296, 1344, 1, 7962624, 259117056, 73903104, 1328256, 3024, 1, 191102976, 31102009344, 26864234496, 1189638144, 6408576, 6048, 1, 4586471424, 3732432224256, 9702226427904, 1026160275456, 11956045824, 24697728, 11088, 1
Offset: 1

Views

Author

Wolfdieter Lang, Dec 01 2003

Keywords

Comments

This triangle underlies the array entry A090214 ((4,4)-generalized Stirling2).

Examples

			[1]; [24,1]; [576,144,1]; [13824,17856,504,1]; ...
		

Crossrefs

Cf. A071951 (Legendre-Stirling, (2, 2) case), A089504 ((3, 3)-case).
The column sequences (without leading zeros) are A009968 (powers of 24), etc.

Programs

  • Mathematica
    max = 10; f[m_] := 1/Product[1-FactorialPower[r+3, 4]*x, {r, 1, m}]; col[m_] := CoefficientList[f[m] + O[x]^(max-m+1), x]; a[n_, m_] := col[m][[n-m+1]]; Table[a[n, m], {n, 1, max}, {m, 1, n}] // Flatten (* Jean-François Alcover, Sep 01 2016 *)

Formula

G.f. for m-th column sequence (without leading zeros and m>=1) is 1/product(1-fallfac(r+3, 4)*x, r=1..m) with fallfac(n, k) := A008279(n, k) (falling factorials).
a(n, m) = sum(A089515(m, p)*fallfac(p, 4)^(n-m), p=1..m)/D(m) if n>=m>=1 else 0; with D(m) := A089516(m).

Extensions

More terms coming from a-file added by Michel Marcus, Feb 08 2023

A178795 Expansion of the polynomial (x^15-1)*(x^12-1)*(x^10-1)*(x^9-1)*(x^7-1)*(x^6-1)*(x^4-1)*(x-1) in increasing powers of x.

Original entry on oeis.org

1, -1, 0, 0, -1, 1, -1, 0, 1, -1, 1, 1, -2, 3, -1, -1, 3, -3, 1, 1, -4, 3, -1, -3, 4, -4, 1, 3, -5, 4, 0, -3, 6, -3, 0, 4, -5, 3, 1, -4, 4, -3, -1, 3, -4, 1, 1, -3, 3, -1, -1, 3, -2, 1, 1, -1, 1, 0, -1, 1, -1, 0, 0, -1, 1
Offset: 0

Views

Author

N. J. A. Sloane, Dec 26 2010

Keywords

Comments

q^120*(q^30-1)*(q^24-1)*(q^20-1)*(q^18-1)*(q^14-1)*(q^12-1)*(q^8-1)*(q^2-1) is the order of the simple group E_8(q), if q is a prime power.
If f(x) is the x-polynomial and g(q) the q-polynomial, then g(q) = q^120*f(q^2). - Jean-François Alcover, Aug 25 2022

Examples

			With p=2 one gets the order of E_8(2): 337804753143634806261388190614085595079991692242467651576160959909068800000. - _Jean-François Alcover_, Aug 25 2022
		

References

  • R. L. Griess, Jr., Twelve Sporadic Groups, Springer, 1998; see p. 169.

Crossrefs

Programs

  • PARI
    Vec((x^15-1)*(x^12-1)*(x^10-1)*(x^9-1)*(x^7-1)*(x^6-1)*(x^4-1)*(x-1)) \\ Michel Marcus, Aug 25 2022

A067423 Eighth column of triangle A067417.

Original entry on oeis.org

1, 10, 240, 5760, 138240, 3317760, 79626240, 1911029760, 45864714240, 1100753141760, 26418075402240, 634033809653760, 15216811431690240, 365203474360565760, 8764883384653578240, 210357201231685877760
Offset: 0

Views

Author

Wolfdieter Lang, Jan 25 2002

Keywords

Crossrefs

Cf. A067422 (seventh column), A067424 (ninth column), A009968 (powers of 24).

Programs

Formula

a(n) = A067417(n+7, 7).
a(n) = 10*(3*8)^(n-1), n >= 1, a(0)=1.
G.f.: (1-14*x)/(1-24*x).

A089507 Second column of triangle A089504 and second column of array A078741 divided by 18.

Original entry on oeis.org

1, 30, 756, 18360, 441936, 10614240, 254788416, 6115201920, 146766525696, 3522406694400, 84537821131776, 2028908069959680, 48693795855814656, 1168651113600245760, 28047626804770062336, 673143043784666480640
Offset: 0

Views

Author

Wolfdieter Lang, Dec 01 2003

Keywords

Comments

Convolution of A000400 (powers of 6) with A009968 (powers of 24).

Crossrefs

Programs

  • Magma
    [6^n*(4^(n+1)-1)/3: n in [0..15]]; // Vincenzo Librandi, Oct 18 2017
  • Mathematica
    CoefficientList[Series[1/((1-6x)(1-24x)),{x,0,20}],x] (* or *) LinearRecurrence[{30,-144},{1,30},20] (* Harvey P. Dale, Sep 25 2017 *)

Formula

G.f.: 1/((1-3*2*1*x)*(1-4*3*2*x)).
a(n) = A089504(n+2, 2), n>=0.
a(n) = (4*(4*3*2)^n - (3*2*1)^n)/3 = (2^n)*(2^(2*(n+1))-1)*3^(n-1).
a(n) = 6^n*(4^(n+1)-1)/3. - Vincenzo Librandi, Oct 18 2017
Showing 1-10 of 24 results. Next