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.

Previous Showing 51-60 of 382 results. Next

A000435 Normalized total height of all nodes in all rooted trees with n labeled nodes.

Original entry on oeis.org

0, 1, 8, 78, 944, 13800, 237432, 4708144, 105822432, 2660215680, 73983185000, 2255828154624, 74841555118992, 2684366717713408, 103512489775594200, 4270718991667353600, 187728592242564421568, 8759085548690928992256, 432357188322752488126152, 22510748754252398927872000
Offset: 1

Views

Author

Keywords

Comments

This is the sequence that started it all: the first sequence in the database!
The height h(V) of a node V in a rooted tree is its distance from the root. a(n) = Sum_{all nodes V in all n^(n-1) rooted trees on n nodes} h(V)/n.
In the trees which have [0, n-1] = (0, 1, ..., n-1) as their ordered set of nodes, the number of nodes at distance i from node 0 is f(n,i) = (n-1)...(n-i)(i+1)n^(j-1), 0 <= i < n-1, i+j = n-1 (and f(n,n-1) = (n-1)!): (n-1)...(n-i) counts the words coding the paths of length i from any node to 0, n^(j-1) counts the Pruefer codes of the rest, words build by iterated deletion of the greater node of degree 1 ... except the last one, (i+1), necessary pointing at the path. If g(n,i) = (n-1)...(n-i)n^j, i+j = n-1, f(n,i) = g(n,i) - g(n,i+1), g(n,i) = Sum_{k>=i} f(n,k), the sequence is Sum_{i=1..n-1} g(n,i). - Claude Lenormand (claude.lenormand(AT)free.fr), Jan 26 2001
If one randomly selects one ball from an urn containing n different balls, with replacement, until exactly one ball has been selected twice, the probability that this ball was also the second ball to be selected once is a(n)/n^n. See also A001865. - Matthew Vandermast, Jun 15 2004
a(n) is the number of connected endofunctions with no fixed points. - Geoffrey Critzer, Dec 13 2011
a(n) is the number of weakly connected simple digraphs on n labeled nodes where every node has out-degree 1. A digraph where all out-degrees are 1 can be called a functional digraph due to the correspondence with endofunctions. - Andrew Howroyd, Feb 06 2024

Examples

			For n = 3 there are 3^2 = 9 rooted labeled trees on 3 nodes, namely (with o denoting a node, O the root node):
   o
   |
   o     o   o
   |      \ /
   O       O
The first can be labeled in 6 ways and contains nodes at heights 1 and 2 above the root, so contributes 6*(1+2) = 18 to the total; the second can be labeled in 3 ways and contains 2 nodes at height 1 above the root, so contributes 3*2=6 to the total, giving 24 in all. Dividing by 3 we get a(3) = 24/3 = 8.
For n = 4 there are 4^3 = 64 rooted labeled trees on 4 nodes, namely (with o denoting a node, O the root node):
   o
   |
   o     o        o   o
   |     |         \ /
   o     o   o      o     o o o
   |      \ /       |      \|/
   O       O        O       O
  (1)     (2)      (3)     (4)
Tree (1) can be labeled in 24 ways and contains nodes at heights 1, 2, 3 above the root, so contributes 24*(1+2+3) = 144 to the total;
tree (2) can be labeled in 24 ways and contains nodes at heights 1, 1, 2 above the root, so contributes 24*(1+1+2) = 96 to the total;
tree (3) can be labeled in 12 ways and contains nodes at heights 1, 2, 2 above the root, so contributes 12*(1+2+2) = 60 to the total;
tree (4) can be labeled in 4 ways and contains nodes at heights 1, 1, 1 above the root, so contributes 4*(1+1+1) = 12 to the total;
giving 312 in all. Dividing by 4 we get a(4) = 312/4 = 78.
		

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. A001863, A001864, A001854, A002862 (unlabeled version), A234953, A259334.
Column k=1 of A350452.

Programs

  • Maple
    A000435 := n-> (n-1)!*add (n^k/k!, k=0..n-2);
    seq(simplify((n-1)*GAMMA(n-1,n)*exp(n)),n=1..20); # Vladeta Jovovic, Jul 21 2005
  • Mathematica
    f[n_] := (n - 1)! Sum [n^k/k!, {k, 0, n - 2}]; Array[f, 18] (* Robert G. Wilson v, Aug 10 2010 *)
    nx = 18; Rest[ Range[0, nx]! CoefficientList[ Series[ LambertW[-x] - Log[1 + LambertW[-x]], {x, 0, nx}], x]] (* Robert G. Wilson v, Apr 13 2013 *)
  • PARI
    x='x+O('x^30); concat(0, Vec(serlaplace(lambertw(-x)-log(1+lambertw(-x))))) \\ Altug Alkan, Sep 05 2018
    
  • PARI
    A000435(n)=(n-1)*A001863(n) \\ M. F. Hasler, Dec 10 2018
    
  • Python
    from math import comb
    def A000435(n): return ((sum(comb(n,k)*(n-k)**(n-k)*k**k for k in range(1,(n+1>>1)))<<1) + (0 if n&1 else comb(n,m:=n>>1)*m**n))//n # Chai Wah Wu, Apr 25-26 2023

Formula

a(n) = (n-1)! * Sum_{k=0..n-2} n^k/k!.
a(n) = A001864(n)/n.
E.g.f.: LambertW(-x) - log(1+LambertW(-x)). - Vladeta Jovovic, Apr 10 2001
a(n) = A001865(n) - n^(n-1).
a(n) = A001865(n) - A000169(n). - Geoffrey Critzer, Dec 13 2011
a(n) ~ sqrt(Pi/2)*n^(n-1/2). - Vaclav Kotesovec, Aug 07 2013
a(n)/A001854(n) ~ 1/2 [See Renyi-Szekeres, (4.7)]. Also a(n) = Sum_{k=0..n-1} k*A259334(n,k). - David desJardins, Jan 20 2017
a(n) = (n-1)*A001863(n). - M. F. Hasler, Dec 10 2018

Extensions

Additional references from Valery A. Liskovets
Editorial changes by N. J. A. Sloane, Feb 03 2012
Edited by M. F. Hasler, Dec 10 2018

A007830 a(n) = (n+3)^n.

Original entry on oeis.org

1, 4, 25, 216, 2401, 32768, 531441, 10000000, 214358881, 5159780352, 137858491849, 4049565169664, 129746337890625, 4503599627370496, 168377826559400929, 6746640616477458432, 288441413567621167681, 13107200000000000000000, 630880792396715529789561
Offset: 0

Views

Author

Peter J. Cameron, Mar 15 1996

Keywords

Comments

a(n-2) is the number of trees with n+1 unlabeled vertices and n labeled edges for n > 1. - Christian G. Bower, 12/99 [corrected by Jonathan Vos Post, Sep 22 2012]
a(n) is the number of nonequivalent primitive meromorphic functions with one pole of order n+3 on a Riemann surface of genus 0. - Noam Katz (noamkj(AT)hotmail.com), Mar 30 2001
Pikhurko writes: "Cameron demonstrated that the total number of edge-labeled trees with n >= 2 edges is (n+1)^(n-2) by showing that the number of vertex-labeled trees of size n is n+1 times larger than the number of edge-labeled ones." - Jonathan Vos Post, Sep 22 2012
With offset 1, a(n) is the number of ways to build a rooted labeled forest with some (possibly all or none) of the nodes from {1,2,...,n} and then build another forest with the remaining nodes. - Geoffrey Critzer, May 10 2013

References

  • M. Shapiro, B. Shapiro and A. Vainshtein - Ramified coverings of S^2 with one degenerate branching point and enumeration of edge-ordered graphs, Amer. Math. Soc. Transl., Vol. 180 (1997), pp. 219-227.
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.27.

Crossrefs

Programs

  • Magma
    [(n+3)^n: n in [0..20]]; // G. C. Greubel, Mar 06 2020
    
  • Maple
    A007830:=n->(n+3)^n; seq(A007830(n), n=0..20);
    T := -LambertW(-x): ser := series(exp(3*T)/(1-T), x, 20):
    seq(n!*coeff(ser, x, n), n = 0..18); # Peter Luschny, Jan 20 2023
  • Mathematica
    Table[(n+3)^n, {n, 0, 18}]
  • PARI
    a(n)=(n+3)^n \\ Charles R Greathouse IV, Feb 06 2017
    
  • Sage
    [(n+3)^n for n in (0..20)] # G. C. Greubel, Mar 06 2020

Formula

E.g.f. for b(n) = a(n-3): T(x) - (3/4)*T^2(x) + (1/6)*T^3(x), where T(x) is Euler's tree function (see A000169). - Len Smiley, Nov 17 2001
E.g.f.: -LambertW(-x)^3/(x^3 * (1+LambertW(-x))). - Vladeta Jovovic, Nov 07 2003
With offset 1: E.g.f.: exp(T(x))^2/2 where T(x) is the e.g.f. for A000169. - Geoffrey Critzer, May 10 2013
E.g.f.: (1/2)*d/dx (LambertW(-x)/(-x))^2. - Wolfdieter Lang, Oct 25 2022

Extensions

More terms from Wesley Ivan Hurt, May 05 2014

A036679 a(n) = n^n - n!.

Original entry on oeis.org

0, 0, 2, 21, 232, 3005, 45936, 818503, 16736896, 387057609, 9996371200, 285271753811, 8915621446656, 302868879571453, 11111919647266816, 437892582706491375, 18446723150919663616, 827239906198908668177, 39346401672922831847424, 1978419534015213180291979
Offset: 0

Views

Author

Keywords

Comments

a(n) = |non-injective functions [n]->[n]| = |non-surjective functions [n]->[n]|.
Fit a polynomial f of degree n-1 to the first n n-th powers of nonnegative integers. Then a(n) = f(n). - Franklin T. Adams-Watters, Dec 28 2006
n^n > n! for n >= 3. [Mitrinovic]

References

  • D. S. Mitrinovic, Analytic Inequalities, Springer-Verlag, 1970; p. 193, 3.1.22.

Crossrefs

Cf. A126130, diagonal of A101030.

Programs

Formula

E.g.f.: 1/(1-T(x))-1/(1-x) where T(x) is the e.g.f. for A000169. - Geoffrey Critzer, Dec 10 2012

A065440 a(n) = (n-1)^n.

Original entry on oeis.org

1, 0, 1, 8, 81, 1024, 15625, 279936, 5764801, 134217728, 3486784401, 100000000000, 3138428376721, 106993205379072, 3937376385699289, 155568095557812224, 6568408355712890625, 295147905179352825856, 14063084452067724991009, 708235345355337676357632
Offset: 0

Views

Author

Len Smiley, Nov 17 2001

Keywords

Comments

a(n) is the number of functions from {1,2,...,n} into {1,2,...,n} that have no fixed points.
The probability that a random function from {1,2,...,n} into {1,2,...,n} has no fixed point is equal to a(n)/n^n; it tends to 1/e when n tends to infinity. - Robert FERREOL, Mar 29 2017

Crossrefs

Essentially the same as A007778 - note T(x) = -W(-x).
Column k=0 of A055134.
Row sums of A350452.

Programs

Formula

a(n) = A007778(n-1).
E.g.f.: x/(T(x)*(1-T(x))) (where T(x) is Euler's tree function, the E.g.f. for n^(n-1)) (see A000169).
a(n) = Sum_{k=0..n} (-1)^k*binomial(n,k)*n^(n-k). - Robert FERREOL, Mar 28 2017
a(n) = Sum_{k=0..n} (-1)^k*binomial(n+2,k+2)*(k+1)*(2*k+n+3)^n. - Vladimir Kruchinin, Aug 13 2025

A048802 Number of labeled rooted trees of nonempty sets with n points. (Each node is a set of 1 or more points.)

Original entry on oeis.org

1, 3, 16, 133, 1521, 22184, 393681, 8233803, 198342718, 5408091155, 164658043397, 5537255169582, 203840528337291, 8153112960102283, 352079321494938344, 16325961781591781401, 809073412162081974237, 42674870241038732398720, 2386963662244981472850709
Offset: 1

Views

Author

Christian G. Bower, Mar 15 1999

Keywords

Examples

			G.f. = x + 3*x^2 + 16*x^3 + 133*x^4 + 1521*x^5 + 22184*x^6 + 393681*x^7 + ...
		

Crossrefs

Programs

  • Mathematica
    nn=20;t=Sum[n^(n-1)x^n/n!,{n,1,nn}];Range[0,nn]!CoefficientList[ ComposeSeries[ Series[t,{x,0,nn}],Series[Exp[x]-1 ,{x,0,nn}]],x]  (* Geoffrey Critzer, Sep 16 2012 *)
  • PARI
    {a(n) = sum( k=1, n, stirling(n, k, 2) * k^(k - 1))}; /* Michael Somos, Jun 09 2012 */
    
  • PARI
    {a(n) = n! * polcoeff( serreverse( log(1 + x*exp(-x +x*O(x^n))) ),n)}
    for(n=1,30,print1(a(n),", ")) \\ Paul D. Hanna, Jan 24 2016

Formula

E.g.f.: B(exp(x)-1) where B is e.g.f. of A000169.
E.g.f.: Series_Reversion( log(1 + x*exp(-x)) ). - Paul D. Hanna, Jan 24 2016
a(n) = Sum_{k=1..n} Stirling2(n, k)*k^(k-1). - Vladeta Jovovic, Sep 17 2003
Stirling transform of A000169. - Michael Somos, Jun 09 2012
a(n) ~ sqrt(1+exp(1)) * n^(n-1) / (exp(n) * (log(1+exp(-1)))^(n-1/2)). - Vaclav Kotesovec, Feb 17 2014

A251592 Triangle of coefficients of polynomials P(n,t) related to the Mittag-Leffler function, where P(n,t) = Product_{k=0..n-2} n*t-k.

Original entry on oeis.org

1, 0, 2, 0, -3, 9, 0, 8, -48, 64, 0, -30, 275, -750, 625, 0, 144, -1800, 7560, -12960, 7776, 0, -840, 13426, -77175, 204085, -252105, 117649, 0, 5760, -112896, 831488, -3010560, 5734400, -5505024, 2097152, 0, -45360, 1058508, -9573228
Offset: 1

Views

Author

Jean-François Alcover, Dec 05 2014

Keywords

Comments

Second column (unsigned) 2, 3, 8, 30, 144, ... is A001048.
Diagonal 1, 2, 9, 64, 625, 7776, ... is A000169.

Examples

			Triangle begins :
  1;
  0,   2;
  0,  -3,     9;
  0,   8,   -48,   64;
  0, -30,   275, -750,    625;
  0, 144, -1800, 7560, -12960, 7776;
  ...
		

References

  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, Reading, MA, 2nd ed. 1998

Crossrefs

Cf. A000169, A001048, A156136, A000108 (B_2(x)), A001764 (B_3(x)), A002293 (B_4(x)), A002294 (B_5(x)), A002295 (B_6(x)), A002296 (B_7(x)), A007556 (B_8(x)), A062994 (B_9(x)), A059968 (B_10(x)), A230388 (B_11(x)), A139526, A260687.

Programs

  • Mathematica
    P[n_, t_] := Product[n*t - k, {k, 0, n-2}]; row[n_] := CoefficientList[P[n, t], t]; Table[row[n], {n, 1, 10}] // Flatten

Formula

P(n,t) = (n-1)!*binomial(n*t, n-1).
From Peter Bala, Nov 15 2015: (Start)
E.g.f. (with constant term 1): B_t(x) = Sum_{n >= 0} 1/(n*t + 1)*binomial(n*t + 1,n)*x^n = 1 + x + 2*t*x^2/2! + 3*t(3*t - 1)*x^3/3! + 4*t*(4*t - 1)*(4*t - 2)*x^4/4! + ... is the generalized binomial series of Lambert. See Graham et al., Section 5.4 and Section 7.5.
In the notation of the Bala link, B_t(x) = I^t(1 + x) where I^t is a fractional inversion operator. B_(1+t)(x) is the e.g.f. for A260687.
B_t(x) = 1 + x*B_t(x)^t.
For complex r, B_t(x)^r = Sum_{n >= 0} r/(n*t + r)*binomial(n*t + r,n)*x^n.
log (B_t(x)) = Sum_{n >= 1} 1/(n*t)*binomial(n*t,n)*x^n.
B_2(x) is the o.g.f. for the Catalan numbers A000108. B_t(x) for t = 3,4,5,... gives the o.g.f. for various Fuss-Catalan sequences. See the cross references. (End)

A007925 a(n) = n^(n+1) - (n+1)^n.

Original entry on oeis.org

-1, -1, -1, 17, 399, 7849, 162287, 3667649, 91171007, 2486784401, 74062575399, 2395420006033, 83695120256591, 3143661612445145, 126375169532421599, 5415486851106043649, 246486713303685957375, 11877172892329028459041, 604107995057426434824791
Offset: 0

Views

Author

Dennis S. Kluk (mathemagician(AT)ameritech.net)

Keywords

Comments

From Mathew Englander, Jul 07 2020: (Start)
All a(n) are odd and for n even, a(n) == 3 (mod 4); for n odd and n != 1, a(n) == 1 (mod 4).
The correspondence between n and a(n) when considered mod 6 is as follows: for n == 0, 1, 2, or 3, a(n) == 5; for n == 4, a(n) == 3; for n == 5, a(n) == 1.
For all n, a(n)+1 is a multiple of n^2.
For n odd and n >= 3, a(n)-1 is a multiple of (n+1)^2.
For n even and n >= 0, a(n)+1 is a multiple of (n+1)^2.
For proofs of the above, see the Englander link. (End)

Examples

			a(2) = 1^2 - 2^1 = -1,
a(4) = 3^4 - 4^3 = 17.
		

References

  • G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.

Crossrefs

Programs

Formula

Asymptotic expression for a(n) is a(n) ~ n^n * (n - e). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 19 2001
From Mathew Englander, Jul 07 2020: (Start)
a(n) = A111454(n+4) - 1.
a(n) = A055651(n, n+1).
a(n) = A220417(n+1, n) for n >= 1.
a(n) = A007778(n) - A000169(n+1).
(End)
E.g.f.: LambertW(-x)/((1+LambertW(-x))*x)-LambertW(-x)/(1+LambertW(-x))^3. - Alois P. Heinz, Jul 04 2022

A081721 Number of bracelets of n beads in up to n colors.

Original entry on oeis.org

1, 3, 10, 55, 377, 4291, 60028, 1058058, 21552969, 500280022, 12969598086, 371514016094, 11649073935505, 396857785692525, 14596464294191704, 576460770691256356, 24330595997127372497, 1092955780817066765469, 52063675152021153895330, 2621440000054016000176044
Offset: 1

Views

Author

N. J. A. Sloane, based on information supplied by Gary W. Adamson, Apr 05 2003

Keywords

Comments

T(n,n), T given in A081720.
From Olivier Gérard, Aug 01 2016: (Start)
Number of classes of functions of [n] to [n] under rotation and reversal.
.
Classes can be of size between 1 and 2n
depending on divisibility properties of n.
.
n 1 2 3 4 5 n 2n
----------------------------------------
1 1
2 2 1
3 3 0 6 1
4 4 6 0 30 15
5 5 0 0 120 252
6 6 15 30 725 3515
7 7 0 0 2394 57627
.
(End)

Crossrefs

Cf. A000312 All endofunctions
Cf. A000169 Classes under translation mod n
Cf. A001700 Classes under sort
Cf. A056665 Classes under rotation
Cf. A168658 Classes under complement to n+1
Cf. A130293 Classes under translation and rotation
Cf. A275549 Classes under reversal
Cf. A275550 Classes under reversal and complement
Cf. A275551 Classes under translation and reversal
Cf. A275552 Classes under translation and complement
Cf. A275553 Classes under translation, complement and reversal
Cf. A275554 Classes under translation, rotation and complement
Cf. A275555 Classes under translation, rotation and reversal
Cf. A275556 Classes under translation, rotation, complement and reversal
Cf. A275557 Classes under rotation and complement
Cf. A275558 Classes under rotation, complement and reversal
Row sums of partition array A213941.
Main diagonal of A321791.

Programs

  • Mathematica
    Table[CycleIndex[DihedralGroup[n],s]/.Table[s[i]->n,{i,1,n}],{n,1,20}] (* Geoffrey Critzer, Jun 18 2013 *)
    t[n_, k_] := (For[t1 = 0; d = 1, d <= n, d++, If[Mod[n, d] == 0, t1 = t1 + EulerPhi[d]*k^(n/d)]]; If[EvenQ[n], (t1 + (n/2)*(1 + k)*k^(n/2))/(2*n), (t1 + n*k^((n + 1)/2))/(2*n)]); a[n_] := t[n, n]; Array[a, 20] (* Jean-François Alcover, Nov 02 2017, after Maple code for A081720 *)

Formula

a(n) ~ n^(n-1) / 2. - Vaclav Kotesovec, Mar 18 2017

Extensions

Name changed by Olivier Gérard, Aug 05 2016
Name revised by Álvar Ibeas, Apr 20 2018

A277473 E.g.f.: -exp(x)*LambertW(-x).

Original entry on oeis.org

0, 1, 4, 18, 116, 1060, 12702, 187810, 3296120, 66897288, 1540762010, 39693752494, 1130866726596, 35300006582620, 1198036854980630, 43921652697277170, 1729775120233353968, 72831210167041246480, 3264674481128340280242, 155220967397580333229270
Offset: 0

Views

Author

Vaclav Kotesovec, Oct 17 2016

Keywords

Crossrefs

Partial sums of A038051.

Programs

  • Mathematica
    CoefficientList[Series[-Exp[x]*LambertW[-x], {x, 0, 20}], x] * Range[0, 20]!
    Table[Sum[Binomial[n, k]*k^(k-1), {k, 1, n}], {n, 0, 20}]
  • PARI
    x='x+O('x^50); concat([0], Vec(serlaplace(-exp(x)*lambertw(-x)))) \\ G. C. Greubel, Jun 11 2017

Formula

a(n) = Sum_{k=1..n} binomial(n,k) * k^(k-1).
a(n) ~ exp(exp(-1)) * n^(n-1).

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

Views

Author

Gus Wiseman, Apr 28 2021

Keywords

Comments

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

Examples

			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
		

Crossrefs

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

Programs

  • Mathematica
    Table[DivisorSigma[0,k^(n-k)],{n,10},{k,n}]
  • PARI
    A(n, k) = numdiv(n^k); \\ Seiichi Manyama, May 15 2021

Formula

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
Previous Showing 51-60 of 382 results. Next