A009998 -id:A009998 - cp's OEIS frontend

A026898 a(n) = Sum_{k=0..n} (n-k+1)^k.

1, 2, 4, 9, 23, 66, 210, 733, 2781, 11378, 49864, 232769, 1151915, 6018786, 33087206, 190780213, 1150653921, 7241710930, 47454745804, 323154696185, 2282779990495, 16700904488706, 126356632390298, 987303454928973, 7957133905608837, 66071772829247410

Offset: 0

N. J. A. Sloane

nonn

Row sums of A004248, A009998, A009999.
First differences are in A047970.
First differences of A103439.
Antidiagonal sums of array A003992.
a(n-1), for n>=1, is the number of length-n restricted growth strings (RGS) [s(0),s(1),...,s(n-1)] where s(0)=0 and s(k)<=1+max(prefix) for k>=1, that are simultaneously projections as maps f: [n] -> [n] where f(x)<=x and f(f(x))=f(x); see example and the two comments (Arndt, Apr 30 2011 Jan 04 2013) in A000110. - Joerg Arndt, Mar 07 2015
Number of finite sequences s of length n+1 whose discriminator sequence is s itself. Here the discriminator sequence of s is the one where the n-th term (n>=1) is the least positive integer k such that the first n terms are pairwise incongruent, modulo k. - Jeffrey Shallit, May 17 2016
From Gus Wiseman, Jan 08 2019: (Start)
Also the number of set partitions of {1,...,n+1} whose minima form an initial interval of positive integers. For example, the a(3) = 9 set partitions are:
  {{1},{2},{3},{4}}
  {{1},{2},{3,4}}
  {{1},{2,4},{3}}
  {{1,4},{2},{3}}
  {{1},{2,3,4}}
  {{1,3},{2,4}}
  {{1,4},{2,3}}
  {{1,3,4},{2}}
  {{1,2,3,4}}
Missing from this list are:
  {{1},{2,3},{4}}
  {{1,2},{3},{4}}
  {{1,3},{2},{4}}
  {{1,2},{3,4}}
  {{1,2,3},{4}}
  {{1,2,4},{3}}
(End)
a(n) is the number of m-tuples of nonnegative integers less than or equal to n-m (including the "0-tuple"). - Mathew Englander, Apr 11 2021

			G.f.: A(x) = 1 + 2*x + 4*x^2 + 9*x^3 + 23*x^4 + 66*x^5 + 210*x^6 + ...
where we have the identity:
A(x) = 1/(1-x) + x/(1-2*x) + x^2/(1-3*x) + x^3/(1-4*x) + x^4/(1-5*x) + ...
is equal to
A(x) = 1/(1-x) + x/((1-x)^2*(1+x)) + 2!*x^2/((1-x)^3*(1+x)*(1+2*x)) + 3!*x^3/((1-x)^4*(1+x)*(1+2*x)*(1+3*x)) + 4!*x^4/((1-x)^5*(1+x)*(1+2*x)*(1+3*x)*(1+4*x)) + ...
From _Joerg Arndt_, Mar 07 2015: (Start)
The a(5-1) = 23 RGS described in the comment are (dots denote zeros):
01:  [ . . . . . ]
02:  [ . 1 . . . ]
03:  [ . 1 . . 1 ]
04:  [ . 1 . 1 . ]
05:  [ . 1 . 1 1 ]
06:  [ . 1 1 . . ]
07:  [ . 1 1 . 1 ]
08:  [ . 1 1 1 . ]
09:  [ . 1 1 1 1 ]
10:  [ . 1 2 . . ]
11:  [ . 1 2 . 1 ]
12:  [ . 1 2 . 2 ]
13:  [ . 1 2 1 . ]
14:  [ . 1 2 1 1 ]
15:  [ . 1 2 1 2 ]
16:  [ . 1 2 2 . ]
17:  [ . 1 2 2 1 ]
18:  [ . 1 2 2 2 ]
19:  [ . 1 2 3 . ]
20:  [ . 1 2 3 1 ]
21:  [ . 1 2 3 2 ]
22:  [ . 1 2 3 3 ]
23:  [ . 1 2 3 4 ]
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 0..500
Fufa Beyene, Jörgen Backelin, Roberto Mantaci, and Samuel A. Fufa, Set Partitions and Other Bell Number Enumerated Objects, J. Int. Seq., Vol. 26 (2023), Article 23.1.8.
Giulio Cerbai, Pattern-avoiding modified ascent sequences, arXiv:2401.10027 [math.CO], 2024. See p. 12.
Giulio Cerbai and Anders Claesson, Counting fixed-point-free Cayley permutations, arXiv:2507.09304 [math.CO], 2025. See p. 25.
Sajed Haque, Discriminators of Integer Sequences, 2017, See p. 33 Corollary 29.
Mathematics Stack Exchange, Asymptotics of ..., 2011.
Chunyan Yan and Zhicong Lin, Inversion sequences avoiding pairs of patterns, arXiv:1912.03674 [math.CO], 2019. See p. 3.
Robin D. P. Zhou, Pattern avoidance in revised ascent sequences, arXiv:2505.05171 [math.CO], 2025. See p. 4.

Cf. A000110, A000258, A000670, A003101, A008277, A038125, A062810.

Cf. A105795, A179928, A287215, A287216.

a026898 n = sum $ zipWith (^) [n + 1, n .. 1] [0 ..]
-- Reinhard Zumkeller, Sep 14 2014

[(&+[(n-k+1)^k: k in [0..n]]): n in [0..50]]; // Stefano Spezia, Jan 09 2019

a:= n-> add((n+1-j)^j, j=0..n): seq(a(n), n=0..23); # Zerinvary Lajos, Apr 18 2009

Table[Sum[(n-k+1)^k, {k,0,n}], {n, 0, 25}] (* Michael De Vlieger, Apr 01 2015 *)

{a(n)=polcoeff(sum(m=0,n,x^m/(1-(m+1)*x+x*O(x^n))),n)} /* Paul D. Hanna, Sep 13 2011 */

{INTEGRATE(n,F)=local(G=F);for(i=1,n,G=intformal(G));G}
{a(n)=local(A=1+x);A=sum(k=0,n,INTEGRATE(k,exp((k+1)*x+x*O(x^n))));n!*polcoeff(A,n)} \\ Paul D. Hanna, Dec 28 2013
for(n=0,30,print1(a(n),", "))

{a(n)=polcoeff( sum(m=0, n, m!*x^m/(1-x +x*O(x^n))^(m+1)/prod(k=1, m, 1+k*x +x*O(x^n))), n)}  /* From o.g.f. (Paul D. Hanna, Jul 20 2014) */
for(n=0, 25, print1(a(n), ", "))

[sum((n-j+1)^j for j in (0..n)) for n in (0..30)] # G. C. Greubel, Jun 15 2021

a(n) = A003101(n) + 1.
G.f.: Sum_{n>=0} x^n/(1 - (n+1)*x). - Paul D. Hanna, Sep 13 2011
G.f.: G(0) where G(k) = 1 + x*(2*k*x-1)/((2*k*x+x-1) - x*(2*k*x+x-1)^2/(x*(2*k*x+x-1) + (2*k*x+2*x-1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 26 2013
E.g.f.: Sum_{n>=0} Integral^n exp((n+1)*x) dx^n, where Integral^n F(x) dx^n is the n-th integration of F(x) with no constant of integration. - Paul D. Hanna, Dec 28 2013
O.g.f.: Sum_{n>=0} n! * x^n/(1-x)^(n+1) / Product_{k=1..n} (1 + k*x). - Paul D. Hanna, Jul 20 2014
a(n) = A101494(n+1,0). - Vladimir Kruchinin, Apr 01 2015
a(n-1) = Sum_{k = 1..n} k^(n-k). - Gus Wiseman, Jan 08 2019
log(a(n)) ~ (1 - 1/LambertW(exp(1)*n)) * n * log(1 + n/LambertW(exp(1)*n)). - Vaclav Kotesovec, Jun 15 2021
a(n) ~ sqrt(2*Pi/(n+1 + (n+1)/w(n))) * ((n+1)/w(n))^(n+2 - (n+1)/w(n)), where w(n) = LambertW(exp(1)*(n+1)). - Vaclav Kotesovec, Jun 25 2021, after user "leonbloy", see Mathematics Stack Exchange link.

a(23)-a(25) from Paul D. Hanna, Dec 28 2013

A047969 Square array of nexus numbers a(n,k) = (n+1)^(k+1) - n^(k+1) (n >= 0, k >= 0) read by upwards antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 7, 1, 1, 7, 19, 15, 1, 1, 9, 37, 65, 31, 1, 1, 11, 61, 175, 211, 63, 1, 1, 13, 91, 369, 781, 665, 127, 1, 1, 15, 127, 671, 2101, 3367, 2059, 255, 1, 1, 17, 169, 1105, 4651, 11529, 14197, 6305, 511, 1, 1, 19, 217, 1695, 9031

Offset: 0

N. J. A. Sloane

If each row started with an initial 0 (i.e., a(n,k) = (n+1)^k - n^k) then each row would be the binomial transform of the preceding row. - Henry Bottomley, May 31 2001
a(n-1, k-1) is the number of ordered k-tuples of positive integers such that the largest of these integers is n. - Alford Arnold, Sep 07 2005
From Alford Arnold, Jul 21 2006: (Start)
The sequences in A047969 can also be calculated using the Eulerian Array (A008292) and Pascal's Triangle (A007318) as illustrated below: (cf. A101095).
  1       1       1       1       1       1
  1       1       1       1       1       1
  -----------------------------------------
  1       2       3       4       5       6
          1       2       3       4       5
  1       3       5       7       9      11
  -----------------------------------------
  1       3       6      10      15      21
          4      12      24      40      60
                  1       3       6      10
  1       7      19      37      61      91
  -----------------------------------------
  1       4      10      20      35      56
         11      44     110     220     385
                 11      44     110     220
                          1       4      10
  1      15      65     175     369     671
  -----------------------------------------  (End)
From Peter Bala, Oct 26 2008: (Start)
The above remarks of Alford Arnold may be summarized by saying that (the transpose of) this array is the Hilbert transform of the triangle of Eulerian numbers A008292 (see A145905 for the definition of the Hilbert transform). In this context, A008292 is best viewed as the array of h-vectors of permutohedra of type A. See A108553 for the Hilbert transform of the array of h-vectors of type D permutohedra. Compare this array with A009998.
The polynomials n^k - (n-1)^k, k = 1,2,3,..., which give the nonzero entries in the columns of this array, satisfy a Riemann hypothesis: their zeros lie on the vertical line Re s = 1/2 in the complex plane. See A019538 for the connection between the polynomials n^k - (n-1)^k and the Stirling polynomials of the simplicial complexes dual to the type A permutohedra.
(End)
Empirical: (n+1)^(k+1) - n^(k+1) is the number of first differences of length k+1 arrays of numbers in 0..n, k > 0. - R. H. Hardin, Jun 30 2013
a(n-1, k-1) is the number of bargraphs of width k and height n. Examples: a(1,2) = 7 because we have [1,1,2], [1,2,1], [2,1,1], [1,2,2], [2,1,2], [2,2,1], and [2,2,2]; a(2,1) = 5 because we have [1,3], [2,3], [3,1], [3,2], and [3,3] (bargraphs are given as compositions). This comment is equivalent to A. Arnold's Sep 2005 comment. - Emeric Deutsch, Jan 30 2017

			Array a begins:
  [n\k][0  1   2    3    4   5  6  ...
  [0]   1  1   1    1    1   1  1  ...
  [1]   1  3   7   15   31  63  ...
  [2]   1  5  19   65  211  ...
  [3]   1  7  37  175  ...
  ...
Triangle T begins:
  n\m   0   1    2     3     4      5      6      7      8     9  10 ...
  0:    1
  1:    1   1
  2:    1   3    1
  3:    1   5    7     1
  4:    1   7   19    15     1
  5:    1   9   37    65    31      1
  6:    1  11   61   175   211     63      1
  7:    1  13   91   369   781    665    127      1
  8:    1  15  127   671  2101   3367   2059    255      1
  9:    1  17  169  1105  4651  11529  14197   6305    511     1
  10:   1  19  217  1695  9031  31031  61741  58975  19171  1023   1
  ...  - _Wolfdieter Lang_, May 07 2021

J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 54.

T. D. Noe, Rows n = 0..100 of triangle, flattened
A. Blecher, C. Brennan, A. Knopfmacher and Helmut Prodinger, The height and width of bargraphs, Discrete Applied Math. 180, (2015), 36-44.
Eric Weisstein's World of Mathematics, Nexus Number

Cf. A047970.
Cf. A009998, A108553 (Hilbert transform of array of h-vectors of type D permutohedra), A145904, A145905.
Row n sequences of array a: A000012, A000225(k+1), A001047(k+1), A005061(k+1), A005060(k+1), A005062(k+1), A016169(k+1), A016177(k+1), A016185(k+1), A016189(k+1), A016195(k+1), A016197(k+1).
Column k sequences of array a: (nexus numbers): A000012, A005408, A003215, A005917(n+1), A022521, A022522, A022523, A022524, A022525, A022526, A022527, A022528.
Cf. A343237 (row reversed triangle).

Flatten[Table[n = d - e; k = e; (n + 1)^(k + 1) - n^(k + 1), {d, 0, 100}, {e, 0, d}]] (* T. D. Noe, Feb 22 2012 *)

T(n,m):=if m=0 then 1 else sum(k!*(-1)^(m+k)*stirling2(m,k)*binomial(n+k-1,n),k,0,m); /* Vladimir Kruchinin, Jan 28 2018 */

From Vladimir Kruchinin: (Start)
O.g.f. of e.g.f of rows of array: ((1-x)*exp(y))/(1-x*exp(y))^2.
T(n,m) = Sum_{k=0..m} k!*(-1)^(m+k)*Stirling2(m,k)*C(n+k-1,n), T(n,0)=1.(End)
From Wolfdieter Lang, May 07 2021: (Start)
T(n,m) = a(n-m,m) = (n-m+1)^(m+1) - (n-m)^(m+1), n >= 0, m = 0, 1,..., n.
O.g.f. column k of the array: polylog(-(k+1), x)*(1-x)/x. See the Peter Bala comment above, and the Eulerian triangle A008292 formula by Vladeta Jovovic, Sep 02 2002.
E.g.f. of e.g.f. of row of the array: exp(y)*(1 + x*(exp(y) - 1))*exp(x*exp(y)).
O.g.f. of triangle's exponential row polynomials R(n, y) = Sum_{m=0} T(n, m)*(y^m)/m!: G(x, y) = exp(x*y)*(1 - x)/(1 - x*exp(x*y))^2. (End)

A088956 Triangle, read by rows, of coefficients of the hyperbinomial transform.

Original entry on oeis.org

1, 1, 1, 3, 2, 1, 16, 9, 3, 1, 125, 64, 18, 4, 1, 1296, 625, 160, 30, 5, 1, 16807, 7776, 1875, 320, 45, 6, 1, 262144, 117649, 27216, 4375, 560, 63, 7, 1, 4782969, 2097152, 470596, 72576, 8750, 896, 84, 8, 1, 100000000, 43046721, 9437184, 1411788, 163296, 15750

Offset: 0

Paul D. Hanna, Oct 26 2003

The hyperbinomial transform of a sequence {b} is defined to be the sequence {d} given by d(n) = Sum_{k=0..n} T(n,k)*b(k), where T(n,k) = (n-k+1)^(n-k-1)*C(n,k).
Given a table in which the n-th row is the n-th binomial transform of the first row, then the hyperbinomial transform of any diagonal results in the next lower diagonal in the table.
The simplest example of a table of iterated binomial transforms is A009998, with a main diagonal of {1,2,9,64,625,...}; and the hyperbinomial transform of this diagonal gives the next lower diagonal, {1,3,16,125,1296,...}, since 1=(1)*1, 3=(1)*1+(1)*2, 16=(3)*1+(2)*2+(1)*9, 125=(16)*1+(9)*2+(3)*9+(1)*64, etc.
Another example: the hyperbinomial transform maps A065440 into A055541, since HYPERBINOMIAL([1,1,1,8,81,1024,15625]) = [1,2,6,36,320,3750,54432] where e.g.f.: A065440(x)+x = x-x/( LambertW(-x)*(1+LambertW(-x)) ), e.g.f.: A055541(x) = x-x*LambertW(-x).
The m-th iteration of the hyperbinomial transform is given by the triangle of coefficients defined by T_m(n,k) = m*(n-k+m)^(n-k-1)*binomial(n,k).
Example: PARI code for T_m: {a=[1,1,1,8,81,1024,15625]; m=1; b=vector(length(a)); for(n=0,length(a)-1, b[n+1]=sum(k=0,n, m*(n-k+m)^(n-k-1)*binomial(n,k)*a[k+1]); print1(b[n+1],","))} RETURNS b=[1,2,6,36,320,3750,54432].
The INVERSE hyperbinomial transform is thus given by m=-1: {a=[1,2,6,36,320,3750,54432]; m=-1; b=vector(length(a)); for(n=0,length(a)-1, b[n+1]=sum(k=0,n, m*(n-k+m)^(n-k-1)*binomial(n,k)*a[k+1]); print1(b[n+1],","))} RETURNS b=[1,1,1,8,81,1024,15625].
Simply stated, the HYPERBINOMIAL transform is to -LambertW(-x)/x as the BINOMIAL transform is to exp(x).
Let A[n] be the set of all forests of labeled rooted trees on n nodes. Build a superset B[n] of A[n] by designating "some" (possibly all or none) of the isolated nodes in each forest. T(n,k) is the number of elements in B[n] with exactly k designated nodes. See A219034. - Geoffrey Critzer, Nov 10 2012

			Rows begin:
       {1},
       {1,      1},
       {3,      2,     1},
      {16,      9,     3,    1},
     {125,     64,    18,    4,   1},
    {1296,    625,   160,   30,   5,  1},
   {16807,   7776,  1875,  320,  45,  6, 1},
  {262144, 117649, 27216, 4375, 560, 63, 7, 1}, ...

T. D. Noe, Rows n=0..50 of triangle, flattened
T. Copeland, Composition, Conjugation, and the Umbral Calculus-Part I, 2021.
G. Helms, Pascalmatrix tetrated
E. W. Weisstein, Abel Polynomial, From MathWorld - A Wolfram Web Resource.
Index entries for triangles and arrays related to Pascal's triangle

Cf. A088957 (row sums), A000272 (first column), A009998, A105599, A132440, A215534 (matrix inverse), A215652.
Cf. A227325 (central terms).
Cf. A007318, A059297, A137542.

a088956 n k =  a095890 (n + 1) (k + 1) * a007318' n k `div` (n - k + 1)
a088956_row n = map (a088956 n) [0..n]
a088956_tabl = map a088956_row [0..]
-- Reinhard Zumkeller, Jul 07 2013

nn=8; t=Sum[n^(n-1)x^n/n!, {n,1,nn}]; Range[0,nn]! CoefficientList[Series[Exp[t+y x] ,{x,0,nn}], {x,y}] //Grid (* Geoffrey Critzer, Nov 10 2012 *)

T(n, k) = (n-k+1)^(n-k-1)*C(n, k).
E.g.f.: -LambertW(-x)*exp(x*y)/x. - Vladeta Jovovic, Oct 27 2003
From Peter Bala, Sep 11 2012: (Start)
Let T(x) = Sum_{n >= 0} n^(n-1)*x^n/n! denote the tree function of A000169. The e.g.f. is (T(x)/x)*exp(t*x) = exp(T(x))*exp(t*x) = 1 + (1 + t)*x + (3 + 2*t + t^2)*x^2/2! + .... Hence the triangle is the exponential Riordan array [T(x)/x,x] belonging to the exponential Appell group.
The matrix power (A088956)^r has the e.g.f. exp(r*T(x))*exp(t*x) with triangle entries given by r*(n-k+r)^(n-k-1)*binomial(n,k) for n and k >= 0. See A215534 for the case r = -1.
Let A(n,x) = x*(x+n)^(n-1) be an Abel polynomial. The present triangle is the triangle of connection constants expressing A(n,x+1) as a linear combination of the basis polynomials A(k,x), 0 <= k <= n. For example, A(4,x+1) = 125*A(0,x) + 64*A(1,x) + 18*A(2,x) + 4*A(3,x) + A(4,x) gives row 4 as [125,64,18,4,1].
Let S be the array with the sequence [1,2,3,...] on the main subdiagonal and zeros elsewhere. S is the infinitesimal generator for Pascal's triangle (see A132440). Then the infinitesimal generator for this triangle is S*A088956; that is, A088956 = Exp(S*A088956), where Exp is the matrix exponential.
With T(x) the tree function as above, define E(x) = T(x)/x. Then A088956 = E(S) = Sum_{n>=0} (n+1)^(n-1)*S^n/n!.
For commuting lower unit triangular matrices A and B, we define A raised to the matrix power B, denoted A^^B, to be the matrix Exp(B*log(A)), where the matrix logarithm Log(A) is defined as Sum_{n >= 1} (-1)^(n+1)*(A-1)^n/n. Let P denote Pascal's triangle A007318. Then the present triangle, call it X, solves the matrix equation P^^X = X . See A215652 for the solution to X^^P = P. Furthermore, if we denote the inverse of X by Y then X^^Y = P. As an infinite tower of matrix powers, A088956 = P^^(P^^(P^^(...))).
A088956 augmented with the sequence (x,x,x,...) on the first superdiagonal is the production matrix for the row polynomials of A105599.
(End)
T(n,k) = A095890(n+1,k+1) * A007318(n,k) / (n-k+1), 0 <= k <= n. - Reinhard Zumkeller, Jul 07 2013
Sum_{k = 0..n} T(n,n-k)*(x - k - 1)^(n-k) = x^n. Setting x = n + 1 gives Sum_{k = 0..n} T(n,k)*k^k = (n + 1)^n. - Peter Bala, Feb 17 2017
As lower triangular matrices, this entry, T, equals unsigned A137542 * A007318 * signed A059297. The Pascal matrix is sandwiched between a pair of inverse matrices, so this entry is conjugate to the Pascal matrix, allowing convergent analytic expressions of T, say f(T), to be computed as f(A007318) sandwiched between the inverse pair. - Tom Copeland, Dec 06 2021

A062319 Number of divisors of n^n, or of A000312(n).

Original entry on oeis.org

1, 1, 3, 4, 9, 6, 49, 8, 25, 19, 121, 12, 325, 14, 225, 256, 65, 18, 703, 20, 861, 484, 529, 24, 1825, 51, 729, 82, 1653, 30, 29791, 32, 161, 1156, 1225, 1296, 5329, 38, 1521, 1600, 4961, 42, 79507, 44, 4005, 4186, 2209, 48, 9457, 99, 5151, 2704, 5565, 54

Offset: 0

Jason Earls, Jul 05 2001

From Gus Wiseman, May 02 2021: (Start)
Conjecture: The number of divisors of n^n equals the number of pairwise coprime ordered n-tuples of divisors of n. Confirmed up to n = 30. For example, the a(1) = 1 through a(5) = 6 tuples are:
  (1)  (1,1)  (1,1,1)  (1,1,1,1)  (1,1,1,1,1)
       (1,2)  (1,1,3)  (1,1,1,2)  (1,1,1,1,5)
       (2,1)  (1,3,1)  (1,1,1,4)  (1,1,1,5,1)
              (3,1,1)  (1,1,2,1)  (1,1,5,1,1)
                       (1,1,4,1)  (1,5,1,1,1)
                       (1,2,1,1)  (5,1,1,1,1)
                       (1,4,1,1)
                       (2,1,1,1)
                       (4,1,1,1)
The unordered case (pairwise coprime n-multisets of divisors of n) is counted by A343654.
(End)

			From _Gus Wiseman_, May 02 2021: (Start)
The a(1) = 1 through a(5) = 6 divisors:
  1  1  1   1    1
     2  3   2    5
     4  9   4    25
        27  8    125
            16   625
            32   3125
            64
            128
            256
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Harry J. Smith)

Cf. A046798, A077592, A035116. - Enrique Pérez Herrero, Nov 09 2010
Number of divisors of A000312(n).
Taking Omega instead of sigma gives A066959.
Positions of squares are A173339.
Diagonal n = k of the array A343656.
A000005 counts divisors.
A059481 counts k-multisets of elements of {1..n}.
A334997 counts length-k strict chains of divisors of n.
A343658 counts k-multisets of divisors.
Pairwise coprimality:
- A018892 counts coprime pairs of divisors.
- A084422 counts pairwise coprime subsets of {1..n}.
- A100565 counts pairwise coprime triples of divisors.
- A225520 counts pairwise coprime sets of divisors.
- A343652 counts maximal pairwise coprime sets of divisors.
- A343653 counts pairwise coprime non-singleton sets of divisors > 1.
- A343654 counts pairwise coprime sets of divisors > 1.
Cf. A000169, A000272, A002064, A002109, A009998, A048691, A143773, A146291, A176029, A327527, A343657.

[NumberOfDivisors(n^n): n in  [0..60]]; // Vincenzo Librandi, Nov 09 2014

A062319[n_IntegerQ]:=DivisorSigma[0,n^n]; (* Enrique Pérez Herrero, Nov 09 2010 *)
Join[{1},DivisorSigma[0,#^#]&/@Range[60]] (* Harvey P. Dale, Jun 06 2024 *)

je=[]; for(n=0,200,je=concat(je,numdiv(n^n))); je

{ for (n=0, 1000, write("b062319.txt", n, " ", numdiv(n^n)); ) } \\ Harry J. Smith, Aug 04 2009

a(n)=local(fm);fm=factor(n);prod(k=1,matsize(fm)[1],fm[k,2]*n+1) \\ Franklin T. Adams-Watters, May 03 2011

a(n) = if(n==0, 1, sumdiv(n, d, n^omega(d))); \\ Seiichi Manyama, May 12 2021

from math import prod
from sympy import factorint
def A062319(n): return prod(n*d+1 for d in factorint(n).values()) # Chai Wah Wu, Jun 03 2021

a(n) = A000005(A000312(n)). - Enrique Pérez Herrero, Nov 09 2010
a(2^n) = A002064(n). - Gus Wiseman, May 02 2021
a(prime(n)) = prime(n) + 1. - Gus Wiseman, May 02 2021
a(n) = Product_{i=1..s} (1 + n * m_i) where (m_1,...,m_s) is the sequence of prime multiplicities (prime signature) of n. - Gus Wiseman, May 02 2021
a(n) = Sum_{d|n} n^omega(d) for n > 0. - Seiichi Manyama May 12 2021

A072597 Expansion of 1/(exp(-x) - x) as exponential generating function.

Original entry on oeis.org

1, 2, 7, 37, 261, 2301, 24343, 300455, 4238153, 67255273, 1185860331, 23000296155, 486655768525, 11155073325917, 275364320099807, 7282929854486431, 205462851526617489, 6158705454187353297, 195465061563672788947, 6548320737474275229347, 230922973019493881984021

Offset: 0

Michael Somos, Jun 23 2002

Polynomials from A140749/A141412 are linked to Stirling1 (see A048594, A129841, A140749). See also P. Flajolet, X. Gourdon, B. Salvy in, available on Internet, RR-1857.pdf (preprint of unavailable Gazette des Mathematiciens 55, 1993, pp. 67-78; for graph 2 see also X. Gourdon RR-1852.pdf, pp. 64-65). What is the corresponding graph for A152650/A152656 = simplified A009998/A119502 linked, via A152818, to a(n), then Stirling2? - Paul Curtz, Dec 16 2008
Denominators in rational approximations of Lambert W(1). See Ramanujan, Notebooks, volume 2, page 22: "2. If e^{-x} = x, shew that the convergents to x are 1/2, 4/7, 21/37, 148/261, &c." Numerators in A006153. - Michael Somos, Jan 21 2019
Call an element g in a semigroup a group element if g^j = g for some j > 1. Then a(n) is the number of group elements in the semigroup of partial transformations of an n-set. Hence a(n) = Sum_{k=0..n} A154372(n,k)*k!. - Geoffrey Critzer, Nov 27 2021

			G.f. = 1 + 2*x + 7*x^2 + 37*x^3 + 261*x^4 + 2301*x^5 + 24343*x^6 + ...

O. Ganyushkin and V Mazorchuk, Classical Finite Transformation Semigroups, Springer, 2009, page 70.
S. Ramanujan, Notebooks, Tata Institute of Fundamental Research, Bombay 1957 Vol. 2, see page 22.

Seiichi Manyama, Table of n, a(n) for n = 0..411
W. S. Gray and M. Thitsa, System Interconnections and Combinatorial Integer Sequences, in: System Theory (SSST), 2013 45th Southeastern Symposium on, Date of Conference: 11-11 March 2013.
G. Jiraskova and J. Shallit, The state complexity of star-complement-star, arXiv preprint arXiv:1203.5353 [cs.FL], 2012. - From _N. J. A. Sloane_, Sep 21 2012

Cf. A000110, A006153, A030178, A089148.
Cf. A048993, A052820.
Cf. A006153, A154372.

CoefficientList[Series[1/(Exp[-x]-x), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Jun 26 2013 *)
a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ 1 / (Exp[-x] - x), {x, 0, n}]]; (* Michael Somos, Jan 21 2019 *)
a[ n_] := If[ n < 0, 0, n! Sum[ (n - k + 1)^k / k!, {k, 0, n}]]; (* Michael Somos, Jan 21 2019 *)

{a(n) = if( n<0, 0, n! * polcoeff( 1 / (exp(-x + x * O(x^n)) - x), n))};

{a(n) = if( n<0, 0, n! * sum(k=0, n, (n-k+1)^k / k!))}; /* Michael Somos, Jan 21 2019 */

E.g.f.: 1 / (exp(-x) - x).
a(n) = n!*Sum_{k=0..n} (n-k+1)^k/k!. - Vladeta Jovovic, Aug 31 2003
a(n) = Sum_{k=0..n} (-1)^(n-k)*Stirling2(n, k)*A052820(k). - Vladeta Jovovic, Apr 12 2004
Recurrence: a(n+1) = 1 + Sum_{j=1..n} binomial(n, j)*a(j)*j. - Jon Perry, Apr 25 2005
E.g.f.: 1/(Q(0) - x) where Q(k) = 1 - x/(2*k+1 - x*(2*k+1)/(x - (2*k+2)/Q(k+1) )); (continued fraction ). - Sergei N. Gladkovskii, Apr 04 2013
a(n) ~ n!/((1+c)*c^(n+1)), where c = A030178 = LambertW(1) = 0.5671432904... - Vaclav Kotesovec, Jun 26 2013
O.g.f.: Sum_{k>=0} k!*x^k/(1 - (k + 1)*x)^(k+1). - Ilya Gutkovskiy, Oct 09 2018
a(n) = A006153(n+1)/(n+1). - Seiichi Manyama, Nov 05 2024

A004248 Array read by ascending antidiagonals: A(n, k) = k^n.

Original entry on oeis.org

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

Offset: 0

David W. Wilson

This array transforms into A371761 using the Akiyama-Tanigawa algorithm for powers applied to the rows. - Peter Luschny, Apr 16 2024

This array transforms into A344499 using the Akiyama-Tanigawa algorithm for powers applied to the columns. - Peter Luschny, Apr 27 2024

			Seen as an array that is read by ascending antidiagonals:
[0] 1, 1,   1,    1,     1,     1,      1,      1,       1, ...
[1] 0, 1,   2,    3,     4,     5,      6,      7,       8, ...
[2] 0, 1,   4,    9,    16,    25,     36,     49,      64, ...
[3] 0, 1,   8,   27,    64,   125,    216,    343,     512, ...
[4] 0, 1,  16,   81,   256,   625,   1296,   2401,    4096, ...
[5] 0, 1,  32,  243,  1024,  3125,   7776,  16807,   32768, ...
[6] 0, 1,  64,  729,  4096, 15625,  46656, 117649,  262144, ...
[7] 0, 1, 128, 2187, 16384, 78125, 279936, 823543, 2097152, ...

T. D. Noe, Rows n=0..50 of triangle, flattened
Michelle Rudolph-Lilith, On the Product Representation of Number Sequences, with Application to the Fibonacci Family, arXiv preprint arXiv:1508.07894 [math.NT], 2015.

For other versions see A051129 and A009998.
Row sums are A026898, diagonal sums are A104872. [Paul Barry, Mar 28 2005]
Cf. A003992, A048723, A371761, A344499.

T[x_, y_] := If[y == 0, 1, (x - y)^y];
Table[T[x, y], {x, 0, 11}, {y, x, 0, -1}] // Flatten (* Jean-François Alcover, Dec 15 2017 *)

T(x, y) = x^y \\ Charles R Greathouse IV, Feb 07 2017

def Arow(n, len): return [k**n for k in range(len)]
for n in range(8): print([n], Arow(n, 9))  # Peter Luschny, Apr 16 2024

Table of x^y, where (x,y) = (0,0), (0,1), (1,0), (0,2), (1,1), (2,0), ...
As a number triangle, columns have g.f. x^k/(1 - kx). - Paul Barry, Mar 28 2005
From Paul Barry, Jul 13 2005: (Start)
T(n, k) = if(k <= n, k^(n - k), 0).
T(n, k) = Sum_{j=0..floor((n-k)/2)} (-1)^j*C(n-k, j)*C(n-k-j, n-k)*k^(n-k-2j).
(End)

New name by Peter Luschny, Apr 16 2024.

A343656 Array read by antidiagonals where A(n,k) is the number of divisors of n^k.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 5, 2, 1, 1, 6, 5, 7, 3, 4, 1, 1, 7, 6, 9, 4, 9, 2, 1, 1, 8, 7, 11, 5, 16, 3, 4, 1, 1, 9, 8, 13, 6, 25, 4, 7, 3, 1, 1, 10, 9, 15, 7, 36, 5, 10, 5, 4, 1, 1, 11, 10, 17, 8, 49, 6, 13, 7, 9, 2, 1, 1, 12, 11, 19, 9, 64, 7, 16, 9, 16, 3, 6, 1

Offset: 1

Gus Wiseman, Apr 28 2021

First differs from A343658 at A(4,2) = 5, A343658(4,2) = 6.

As a triangle, T(n,k) = number of divisors of k^(n-k).

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7
  n=1:  1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8
  n=3:  1   2   3   4   5   6   7   8
  n=4:  1   3   5   7   9  11  13  15
  n=5:  1   2   3   4   5   6   7   8
  n=6:  1   4   9  16  25  36  49  64
  n=7:  1   2   3   4   5   6   7   8
  n=8:  1   4   7  10  13  16  19  22
  n=9:  1   3   5   7   9  11  13  15
Triangle begins:
  1
  1  1
  1  2  1
  1  3  2  1
  1  4  3  3  1
  1  5  4  5  2  1
  1  6  5  7  3  4  1
  1  7  6  9  4  9  2  1
  1  8  7 11  5 16  3  4  1
  1  9  8 13  6 25  4  7  3  1
  1 10  9 15  7 36  5 10  5  4  1
  1 11 10 17  8 49  6 13  7  9  2  1
  1 12 11 19  9 64  7 16  9 16  3  6  1
  1 13 12 21 10 81  8 19 11 25  4 15  2  1
For example, row n = 8 counts the following divisors:
  1  64  243  256  125  36  7  1
     32  81   128  25   18  1
     16  27   64   5    12
     8   9    32   1    9
     4   3    16        6
     2   1    8         4
     1        4         3
              2         2
              1         1

Seiichi Manyama, Antidiagonals n = 1..140, flattened

Columns k=1..9 of the array give A000005, A048691, A048785, A344327, A344328, A344329, A343526, A344335, A344336.
Row n = 6 of the array is A000290.
Diagonal n = k of the array is A062319.
Array antidiagonal sums (row sums of the triangle) are A343657.
Dominated by A343658.
A000312 = n^n.
A007318 counts k-sets of elements of {1..n}.
A009998(n,k) = n^k (as an array, offset 1).
A059481 counts k-multisets of elements of {1..n}.
Cf. A000169, A000272, A002064, A002109, A066959, A143773, A146291, A176029, A251683, A282935, A326358, A327527, A334996, A343652.

Table[DivisorSigma[0,k^(n-k)],{n,10},{k,n}]

A(n, k) = numdiv(n^k); \\ Seiichi Manyama, May 15 2021

A(n,k) = A000005(A009998(n,k)), where A009998(n,k) = n^k is the interpretation as an array.

A(n,k) = Sum_{d|n} k^omega(d). - Seiichi Manyama, May 15 2021

A152818 Array read by antidiagonals: A(n,k) = (k+1)^n*(n+k)!/n!.

Original entry on oeis.org

1, 1, 1, 1, 4, 2, 1, 12, 18, 6, 1, 32, 108, 96, 24, 1, 80, 540, 960, 600, 120, 1, 192, 2430, 7680, 9000, 4320, 720, 1, 448, 10206, 53760, 105000, 90720, 35280, 5040, 1, 1024, 40824, 344064, 1050000, 1451520, 987840, 322560, 40320

Offset: 0

Paul Curtz, Dec 13 2008

A009998/A119502 gives triangle of unreduced coefficients of polynomials defined by A152650/A152656. a(n) gives numerators with denominators n! for each row.

Row 0 is A000142. Row 1 is formed from positive members of A001563. Row 2 is A055533. Column 0 is A000012. Column 1 is formed from positive members of A001787. Column 2 is A006043. Column 3 is A006044. - Omar E. Pol, Jan 06 2009

			From _Omar E. Pol_, Jan 06 2009: (Start)
Array begins:
  1,    1,      2,        6,         24,          120, ...
  1,    4,     18,       96,        600,         4320, ...
  1,   12,    108,      960,       9000,        90720, ...
  1,   32,    540,     7680,     105000,      1451520, ...
  1,   80,   2430,    53760,    1050000,     19595520, ...
  1,  192,  10206,   344064,    9450000,    235146240, ...
  1,  448,  40824,  2064384,   78750000,   2586608640, ...
  1, 1024, 157464, 11796480,  618750000,  26605117440, ...
  1, 2304, 590490, 64880640, 4640625000, 259399895040, ... (End)
Antidiagonal triangle:
  1;
  1,   1;
  1,   4,     2;
  1,  12,    18,     6;
  1,  32,   108,    96,     24;
  1,  80,   540,   960,    600,   120;
  1, 192,  2430,  7680,   9000,  4320,   720;
  1, 448, 10206, 53760, 105000, 90720, 35280, 5040;

G. C. Greubel, Antidiagonals n = 0..50, flattened
F. A. Haight, Overflow at a traffic light, Biometrika, 46 (1959), 420-424. See page 422.
F. A. Haight, Overflow at a traffic light, Biometrika, 46 (1959), 420-424. (Annotated scanned copy)
F. A. Haight, Letter to N. J. A. Sloane, n.d.

Cf. A000012, A000142, A001563, A001787, A006043, A006044, A009998.
Cf. A055533, A072597 (antidiagonal sums), A089148, A119502, A152650.
Cf. A152656, A154372.

A152818:= func< n,k | (k+1)^(n-k)*Factorial(k)*Binomial(n,k) >;
[A152818(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Apr 10 2023

len= 45; m= 1 + Ceiling[Sqrt[len]]; Sort[Flatten[#, 1] &[MapIndexed[ {(2 +#2[[1]]^2 +(#2[[2]] -1)*#2[[2]] +#2[[1]]*(2*#2[[2]] -3))/ 2, #1}&, Table[(k+1)^n*(n+k)!/n!, {n,0,m}, {k,0,m}], {2}]]][[All, 2]][[1 ;; len]] (* From Jean-François Alcover, May 27 2011 *)
T[n_, k_]:= (k+1)^(n-k)*k!*Binomial[n, k];
Table[T[n,k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Apr 10 2023 *)

A(n,k) = (k+1)^n*(n+k)!/n! \\ Charles R Greathouse IV, Sep 10 2016

def A152818_row(n):
    R. = ZZ[]
    P = add((n-k+1)^k*x^(n-k+1)*factorial(n)/factorial(k) for k in (0..n))
    return P.coefficients()
for n in (0..12): print(A152818_row(n))  # Peter Luschny, May 03 2013

E.g.f. for array as a triangle: exp(x)/(1-t*x*exp(x)) = 1+(1+t)*x+(1+4*t+2*t^2)*x^2/2! + (1+12*t+18*t^2+6*t^3)*x^3/3! + .... E.g.f. is int {z = 0..inf} exp(-z)*F(x,t*z), (x and t chosen sufficiently small for the integral to converge), where F(x,t) = exp(x*(1+t*exp(x))) is the e.g.f. for A154372. - Peter Bala, Oct 09 2011
From Peter Bala, Oct 09 2011: (Start)
From the e.g.f., the row polynomials R(n,t) satisfy the recursion R(n,t) = 1 + t*sum {k = 0..n-1} n!/(k!*(n-k-1)!)*R(n-k-1,t). The polynomials 1/n!*R(n,x) are the polynomials P(n,x) of A152650.
Sum_{k=0..n} T(n, k) = A072597(n) (antidiagonal sums). (End)
From G. C. Greubel, Apr 10 2023: (Start)
T(n, k) = (k+1)^(n-k) * k! * binomial(n, k) (antidiagonal triangle).
Sum_{k=0..n} (-1)^k*T(n, k) = A089148(n). (End)

Better definition, extended and edited by Omar E. Pol and N. J. A. Sloane, Jan 05 2009

A343658 Array read by antidiagonals where A(n,k) is the number of ways to choose a multiset of k divisors of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 6, 2, 1, 1, 6, 5, 10, 3, 4, 1, 1, 7, 6, 15, 4, 10, 2, 1, 1, 8, 7, 21, 5, 20, 3, 4, 1, 1, 9, 8, 28, 6, 35, 4, 10, 3, 1, 1, 10, 9, 36, 7, 56, 5, 20, 6, 4, 1, 1, 11, 10, 45, 8, 84, 6, 35, 10, 10, 2, 1

Offset: 1

Gus Wiseman, Apr 29 2021

First differs from A343656 at A(4,2) = 6, A343656(4,2) = 5.

As a triangle, T(n,k) = number of ways to choose a multiset of n - k divisors of k.

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8
  n=1:  1   1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8   9
  n=3:  1   2   3   4   5   6   7   8   9
  n=4:  1   3   6  10  15  21  28  36  45
  n=5:  1   2   3   4   5   6   7   8   9
  n=6:  1   4  10  20  35  56  84 120 165
  n=7:  1   2   3   4   5   6   7   8   9
  n=8:  1   4  10  20  35  56  84 120 165
  n=9:  1   3   6  10  15  21  28  36  45
Triangle begins:
   1
   1   1
   1   2   1
   1   3   2   1
   1   4   3   3   1
   1   5   4   6   2   1
   1   6   5  10   3   4   1
   1   7   6  15   4  10   2   1
   1   8   7  21   5  20   3   4   1
   1   9   8  28   6  35   4  10   3   1
   1  10   9  36   7  56   5  20   6   4   1
   1  11  10  45   8  84   6  35  10  10   2   1
For example, row n = 6 counts the following multisets:
  {1,1,1,1,1}  {1,1,1,1}  {1,1,1}  {1,1}  {1}  {}
               {1,1,1,2}  {1,1,3}  {1,2}  {5}
               {1,1,2,2}  {1,3,3}  {1,4}
               {1,2,2,2}  {3,3,3}  {2,2}
               {2,2,2,2}           {2,4}
                                   {4,4}
Note that for n = 6, k = 4 in the triangle, the two multisets {1,4} and {2,2} represent the same divisor 4, so they are only counted once under A343656(4,2) = 5.

Andrew Howroyd, Table of n, a(n) for n = 1..1275 (first 50 antidiagonals)

Row k = 1 of the array is A000005.
Column n = 4 of the array is A000217.
Column n = 6 of the array is A000292.
Row k = 2 of the array is A184389.
The distinct products of these multisets are counted by A343656.
Antidiagonal sums of the array (or row sums of the triangle) are A343661.
A000312 = n^n.
A009998(n,k) = n^k (as an array, offset 1).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.
Cf. A000169, A062319, A067824, A143773, A146291, A176029, A285572, A326077, A327527, A334996, A343652, A343657.

multchoo[n_,k_]:=Binomial[n+k-1,k];
Table[multchoo[DivisorSigma[0,k],n-k],{n,10},{k,n}]

A(n,k) = binomial(numdiv(n) + k - 1, k)
{ for(n=1, 9, for(k=0, 8, print1(A(n,k), ", ")); print ) } \\ Andrew Howroyd, Jan 11 2024

A(n,k) = ((A000005(n), k)) = A007318(A000005(n) + k - 1, k).

T(n,k) = ((A000005(k), n - k)) = A007318(A000005(k) + n - k - 1, n - k).

A067402 Triangle with columns built from certain power sequences.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 12, 5, 1, 1, 48, 45, 7, 1, 1, 192, 405, 112, 9, 1, 1, 768, 3645, 1792, 225, 11, 1, 1, 3072, 32805, 28672, 5625, 396, 13, 1, 1, 12288, 295245, 458752, 140625, 14256, 637, 15, 1, 1, 49152, 2657205, 7340032, 3515625, 513216, 31213, 960, 17, 1

Offset: 0

Wolfdieter Lang, Jan 25 2002

			Triangle starts:
      m=0   1   2   3   4
  n=0:  1;
  n=1:  1,  1;
  n=2:  1,  3,  1;
  n=3:  1, 12,  5,  1;
  n=4:  1, 48, 45,  7,  1;
  ...

Indranil Ghosh, Rows 0..125, flattened

Cf. A009998 (triangle built from powers of (m+1)).

Columns m=0..8: A000012, A002001, A067403-A067409.

A[n_,m_]:=If[n==m,1,(2m+1)(m+1)^(2(n-m-1))]; Flatten[Table[A[n,m],{n,0,9},{m,0,n}]] (* Stefano Spezia, Sep 30 2022 *)

a(n, m) = 1 if n=m; a(n, m) = (2*m+1)*(m+1)^(2*(n-m-1)) if n>m>=0.

G.f. for column m: (x^m)*(1-x*m^2)/(1-x*(m+1)^2).

More terms from Larry Reeves (larryr(AT)acm.org), Jun 11 2002

cp's OEIS Frontend

A026898 a(n) = Sum_{k=0..n} (n-k+1)^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Sage

Formula

Extensions

A047969 Square array of nexus numbers a(n,k) = (n+1)^(k+1) - n^(k+1) (n >= 0, k >= 0) read by upwards antidiagonals.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Maxima

Formula

A088956 Triangle, read by rows, of coefficients of the hyperbinomial transform.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A062319 Number of divisors of n^n, or of A000312(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

PARI

PARI

Python

Formula

A072597 Expansion of 1/(exp(-x) - x) as exponential generating function.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A004248 Array read by ascending antidiagonals: A(n, k) = k^n.

Views

Author

Keywords