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

A000169 Number of labeled rooted trees with n nodes: n^(n-1).

Original entry on oeis.org

1, 2, 9, 64, 625, 7776, 117649, 2097152, 43046721, 1000000000, 25937424601, 743008370688, 23298085122481, 793714773254144, 29192926025390625, 1152921504606846976, 48661191875666868481, 2185911559738696531968, 104127350297911241532841, 5242880000000000000000000
Offset: 1

Views

Author

Keywords

Comments

Also the number of connected transitive subtree acyclic digraphs on n vertices. - Robert Castelo, Jan 06 2001
For any given integer k, a(n) is also the number of functions from {1,2,...,n} to {1,2,...,n} such that the sum of the function values is k mod n. - Sharon Sela (sharonsela(AT)hotmail.com), Feb 16 2002
The n-th term of a geometric progression with first term 1 and common ratio n: a(1) = 1 -> 1,1,1,1,... a(2) = 2 -> 1,2,... a(3) = 9 -> 1,3,9,... a(4) = 64 -> 1,4,16,64,... - Amarnath Murthy, Mar 25 2004
All rational solutions to the equation x^y = y^x, with x < y, are given by x = A000169(n+1)/A000312(n), y = A000312(n+1)/A007778(n), where n = 1, 2, 3, ... . - Nick Hobson, Nov 30 2006
a(n+1) is also the number of partial functions on n labeled objects. - Franklin T. Adams-Watters, Dec 25 2006
In other words, if A is a finite set of size n-1, then a(n) is the number of binary relations on A that are also functions. Note that a(n) = Sum_{k=0..n-1} binomial(n-1,k)*(n-1)^k = n^(n-1), where binomial(n-1,k) is the number of ways to select a domain D of size k from A and (n-1)^k is the number of functions from D to A. - Dennis P. Walsh, Apr 21 2011
This is the fourth member of a set of which the other members are the symmetric group, full transformation semigroup, and symmetric inverse semigroup. For the first three, see A000142, A000312, A002720. - Peter J. Cameron, Nov 03 2024.
More generally, consider the class of sequences of the form a(n) = (n*c(1)*...*c(i))^(n-1). This sequence has c(1)=1. A052746 has a(n) = (2*n)^(n-1), A052756 has a(n) = (3*n)^(n-1), A052764 has a(n) = (4*n)^(n-1), A052789 has a(n) = (5*n)^(n-1) for n>0. These sequences have a combinatorial structure like simple grammars. - Ctibor O. Zizka, Feb 23 2008
a(n) is equal to the logarithmic transform of the sequence b(n) = n^(n-2) starting at b(2). - Kevin Hu (10thsymphony(AT)gmail.com), Aug 23 2010
Also, number of labeled connected multigraphs of order n without cycles except one loop. See link below to have a picture showing the bijection between rooted trees and multigraphs of this kind. (Note that there are no labels in the picture, but the bijection remains true if we label the nodes.) - Washington Bomfim, Sep 04 2010
a(n) is also the number of functions f:{1,2,...,n} -> {1,2,...,n} such that f(1) = 1.
For a signed version of A000169 arising from the Vandermonde determinant of (1,1/2,...,1/n), see the Mathematica section. - Clark Kimberling, Jan 02 2012
Numerator of (1+1/(n-1))^(n-1) for n>1. - Jean-François Alcover, Jan 14 2013
Right edge of triangle A075513. - Michel Marcus, May 17 2013
a(n+1) is the number of n x n binary matrices with no more than a single one in each row. Partitioning the set of such matrices by the number k of rows with a one, we obtain a(n+1) = Sum_{k=0..n} binomial(n,k)*n^k = (n+1)^n. - Dennis P. Walsh, May 27 2014
Central terms of triangle A051129: a(n) = A051129(2*n-1,n). - Reinhard Zumkeller, Sep 14 2014
a(n) is the row sum of the n-th rows of A248120 and A055302, so it enumerates the monomials in the expansion of [x(1) + x(2) + ... + x(n)]^(n-1). - Tom Copeland, Jul 17 2015
For any given integer k, a(n) is the number of sums x_1 + ... + x_m = k (mod n) such that: x_1, ..., x_m are nonnegative integers less than n, the order of the summands does not matter, and each integer appears fewer than n times as a summand. - Carlo Sanna, Oct 04 2015
a(n) is the number of words of length n-1 over an alphabet of n letters. - Joerg Arndt, Oct 07 2015
a(n) is the number of parking functions whose largest element is n and length is n. For example, a(3) = 9 because there are nine such parking functions, namely (1,2,3), (1,3,2), (2,3,1), (2,1,3), (3,1,2), (3,2,1), (1,1,3), (1,3,1), (3,1,1). - Ran Pan, Nov 15 2015
Consider the following problem: n^2 cells are arranged in a square array. A step can be defined as going from one cell to the one directly above it, to the right of it or under it. A step above cannot be followed by a step below and vice versa. Once the last column of the square array is reached, you can only take steps down. a(n) is the number of possible paths (i.e., sequences of steps) from the cell on the bottom left to the cell on the bottom right. - Nicolas Nagel, Oct 13 2016
The rationals c(n) = a(n+1)/a(n), n >= 1, appear in the proof of G. Pólya's "elementary, but not too elementary, theorem": Sum_{n>=1} (Product_{k=1..n} a_k)^(1/n) < exp(1)*Sum_{n>=1} a_n, for n >= 1, with the sequence {a_k}{k>=1} of nonnegative terms, not all equal to 0. - _Wolfdieter Lang, Mar 16 2018
Coefficients of the generating series for the preLie operadic algebra. Cf. p. 417 of the Loday et al. paper. - Tom Copeland, Jul 08 2018
a(n)/2^(n-1) is the square of the determinant of the n X n matrix M_n with elements m(j,k) = cos(Pi*j*k/n). See Zhi-Wei Sun, Petrov link. - Hugo Pfoertner, Sep 19 2021
a(n) is the determinant of the n X n matrix P_n such that, when indexed [0, n), P(0, j) = 1, P(i <= j) = i, and P(i > j) = i-n. - C.S. Elder, Mar 11 2024

Examples

			For n=3, a(3)=9 because there are exactly 9 binary relations on A={1, 2} that are functions, namely: {}, {(1,1)}, {(1,2)}, {(2,1)}, {(2,2)}, {(1,1),(2,1)}, {(1,1),(2,2)}, {(1,2),(2,1)} and {(1,2),(2,2)}. - _Dennis P. Walsh_, Apr 21 2011
G.f. = x + 2*x^2 + 9*x^3 + 64*x^4 + 625*x^5 + 7776*x^6 + 117649*x^7 + ...
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 169.
  • Jonathan L. Gross and Jay Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 524.
  • Hannes Heikinheimo, Heikki Mannila and Jouni K. Seppnen, Finding Trees from Unordered 01 Data, in Knowledge Discovery in Databases: PKDD 2006, Lecture Notes in Computer Science, Volume 4213/2006, Springer-Verlag. - N. J. A. Sloane, Jul 09 2009
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 63.
  • John Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 128.
  • 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).
  • Richard P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see page 25, Prop. 5.3.2, and p. 37, (5.52).

Crossrefs

Programs

  • Haskell
    a000169 n = n ^ (n - 1)  -- Reinhard Zumkeller, Sep 14 2014
    
  • Magma
    [n^(n-1): n in [1..20]]; // Vincenzo Librandi, Jul 17 2015
    
  • Maple
    A000169 := n -> n^(n-1);
    # second program:
    spec := [A, {A=Prod(Z, Set(A))}, labeled]; [seq(combstruct[count](spec, size=n), n=1..20)];
    # third program:
    A000169 := n -> add((-1)^(n+k-1)*pochhammer(n, k)*Stirling2(n-1, k), k = 0..n-1):
    seq(A000169(n), n = 1 .. 23);  # Mélika Tebni, May 07 2023
  • Mathematica
    Table[n^(n - 1), {n, 1, 20}] (* Stefan Steinerberger, Apr 01 2006 *)
    Range[0, 18]! CoefficientList[ Series[ -LambertW[-x], {x, 0, 18}], x] // Rest (* Robert G. Wilson v, updated by Jean-François Alcover, Oct 14 2019 *)
    (* Next, a signed version A000169 from the Vandermonde determinant of (1,1/2,...,1/n) *)
    f[j_] := 1/j; z = 12;
    v[n_] := Product[Product[f[k] - f[j], {j, 1, k - 1}], {k, 2, n}]
    Table[v[n], {n, 1, z}]
    1/%  (* A203421 *)
    Table[v[n]/v[n + 1], {n, 1, z - 1}]  (* A000169 signed *)
    (* Clark Kimberling, Jan 02 2012 *)
    a[n_]:=Det[Table[If[i==0,1,If[i<=j,i,i-n]],{i,0,n-1},{j,0,n-1}]]; Array[a,20] (* Stefano Spezia, Mar 12 2024 *)
  • MuPAD
    n^(n-1) $ n=1..20 /* Zerinvary Lajos, Apr 01 2007 */
    
  • PARI
    a(n) = n^(n-1)
    
  • Python
    def a(n): return n**(n-1)
    print([a(n) for n in range(1, 21)]) # Michael S. Branicky, Sep 19 2021
    
  • Python
    from sympy import Matrix
    def P(n): return [[ (i-n if i > j else i) + (i == 0) for j in range(n) ] for i in range(n)]
    print(*(Matrix(P(n)).det() for n in range(1, 21)), sep=', ') # C.S. Elder, Mar 12 2024

Formula

The e.g.f. T(x) = Sum_{n>=1} n^(n-1)*x^n/n! satisfies T(x) = x*exp(T(x)), so T(x) is the functional inverse (series reversion) of x*exp(-x).
Also T(x) = -LambertW(-x) where W(x) is the principal branch of Lambert's function.
T(x) is sometimes called Euler's tree function.
a(n) = A000312(n-1)*A128434(n,1)/A128433(n,1). - Reinhard Zumkeller, Mar 03 2007
E.g.f.: LambertW(x)=x*G(0); G(k) = 1 - x*((2*k+2)^(2*k))/(((2*k+1)^(2*k)) - x*((2*k+1)^(2*k))*((2*k+3)^(2*k+1))/(x*((2*k+3)^(2*k+1)) - ((2*k+2)^(2*k+1))/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Dec 30 2011
a(n) = Sum_{i=1..n} binomial(n-1,i-1)*i^(i-2)*(n-i)^(n-i). - Dmitry Kruchinin, Oct 28 2013
Limit_{n->oo} a(n)/A000312(n-1) = e. - Daniel Suteu, Jul 23 2016
From Amiram Eldar, Nov 20 2020: (Start)
Sum_{n>=1} 1/a(n) = A098686.
Sum_{n>=1} (-1)^(n+1)/a(n) = A262974. (End)
a(n) = Sum_{k=0..n-1} (-1)^(n+k-1)*Pochhammer(n, k)*Stirling2(n-1, k). - Mélika Tebni, May 07 2023
In terms of Eulerian numbers A340556(n,k) of the second order Sum_{m>=1} m^(m+n) z^m/m! = 1/(1-T(z))^(2n+1) * Sum_{k=0..n} A2(n,k) T(z)^k. - Marko Riedel, Jan 10 2024

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

Original entry on oeis.org

0, 1, 3, 8, 22, 65, 209, 732, 2780, 11377, 49863, 232768, 1151914, 6018785, 33087205, 190780212, 1150653920, 7241710929, 47454745803, 323154696184, 2282779990494, 16700904488705, 126356632390297, 987303454928972, 7957133905608836, 66071772829247409
Offset: 0

Views

Author

Keywords

Comments

For n > 0: a(n) = sum of row n of triangles A051129 and A247358. - Reinhard Zumkeller, Sep 14 2014
a(n-1) is the number of set partitions of [n] into two or more blocks such that all absolute differences between least elements of consecutive blocks are 1. a(3) = 8: 134|2, 13|24, 14|23, 1|234, 14|2|3, 1|24|3, 1|2|34, 1|2|3|4. - Alois P. Heinz, May 22 2017
Min_{n >= 1} a(n+1)/a(n) = 8/3. This is the answer to the 4th problem proposed during the first day of the final round of the 16th Austrian Mathematical Olympiad in 1985 (see link IMO Compendium). - Bernard Schott, Jan 07 2019
If n rings of different internal diameter can fit close together on a tapering column, a(n) is the number of different arrangements of at least one ring. For example, if the rings increasing in size are 1, 2 and 3, then a(3) = 8 corresponding to the possible arrangements from the point on the column of smallest diameter (1XX), (X2X), (XX3), (12X), (32X), (1X3), (X23) and (123), where X denotes a space on the column. - Ian Duff, Jun 23 2025

Examples

			For n = 3 we get a(3) = 3^1 + 2^2 + 1^3 = 8. For n = 4 we get a(4) = 4^1 + 3^2 + 2^3 + 1^4 = 22.
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

First differences are in A047970.

Programs

  • Haskell
    a003101 n = sum $ zipWith (^) [0 ..] [n + 1, n .. 1]
    -- Reinhard Zumkeller, Sep 14 2014
    
  • Magma
    [n eq 0 select 0 else (&+[(n-j+1)^j: j in [1..n]]): n in [0..50]]; // G. C. Greubel, Oct 26 2022
    
  • Maple
    A003101 := n->add((n-k+1)^k, k=1..n);
    a:= n-> add((n-j+1)^j, j=1..n): seq(a(n), n=0..30); # Zerinvary Lajos, Jun 07 2008
  • Mathematica
    Table[Sum[(n-k+1)^k,{k,n}],{n,0,25}] (* Harvey P. Dale, Aug 14 2011 *)
  • PARI
    a(n)=sum(k=1,n,(n-k+1)^k) \\ Charles R Greathouse IV, Oct 31 2011
    
  • SageMath
    def A003101(n): return sum( (n-k+1)^k for k in range(1,n+1))
    [A003101(n) for n in range(50)] # G. C. Greubel, Oct 26 2022

Formula

a(n) = A026898(n) - 1.
G.f.: G(0)/x-1/(1-x)/x 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) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 26 2013
G.f.: Sum_{k>=1} x^k/(1 - (k + 1)*x). - Ilya Gutkovskiy, Oct 09 2018
a(n) = n^1 + (n-1)^2 + (n-2)^3 + ... + 3^(n-2) + 2^(n-1) + 1^n. - Bernard Schott, Jan 07 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 + w(n))) * w(n)^(n+2 - w(n)), where w(n) = (n+1)/LambertW(exp(1)*(n+1)). - Vaclav Kotesovec, Jun 25 2021, after user "leonbloy", see Mathematics Stack Exchange link.

A089072 Triangle read by rows: T(n,k) = k^n, n >= 1, 1 <= k <= n.

Original entry on oeis.org

1, 1, 4, 1, 8, 27, 1, 16, 81, 256, 1, 32, 243, 1024, 3125, 1, 64, 729, 4096, 15625, 46656, 1, 128, 2187, 16384, 78125, 279936, 823543, 1, 256, 6561, 65536, 390625, 1679616, 5764801, 16777216, 1, 512, 19683, 262144, 1953125, 10077696, 40353607, 134217728, 387420489
Offset: 1

Views

Author

Alford Arnold, Dec 04 2003

Keywords

Comments

T(n, k) = number of mappings from an n-element set into a k-element set. - Clark Kimberling, Nov 26 2004
Let S be the semigroup of (full) transformations on [n]. Let a be in S with rank(a) = k. Then T(n,k) = |a S|, the number of elements in the right principal ideal generated by a. - Geoffrey Critzer, Dec 30 2021
From Manfred Boergens, Jun 23 2024: (Start)
In the following two comments the restriction k<=n can be lifted, allowing all k>=1.
T(n,k) is the number of n X k binary matrices with row sums = 1.
T(n,k) is the number of coverings of [n] by tuples (A_1,...,A_k) in P([n])^k with disjoint A_j, with P(.) denoting the power set.
For nonempty A_j see A019538.
For tuples with "disjoint" dropped see A092477.
For tuples with nonempty A_j and with "disjoint" dropped see A218695. (End)

Examples

			Triangle begins:
  1;
  1,  4;
  1,  8,  27;
  1, 16,  81,  256;
  1, 32, 243, 1024,  3125;
  1, 64, 729, 4096, 15625, 46656;
  ...
		

Crossrefs

Related to triangle of Eulerian numbers A008292.

Programs

  • Haskell
    a089072 = flip (^)
    a089072_row n = map (a089072 n) [1..n]
    a089072_tabl = map a089072_row [1..]  -- Reinhard Zumkeller, Mar 18 2013
    
  • Magma
    [k^n: k in [1..n], n in [1..12]]; // G. C. Greubel, Nov 01 2022
    
  • Mathematica
    Column[Table[k^n, {n, 8}, {k, n}], Center] (* Alonso del Arte, Nov 14 2011 *)
  • SageMath
    flatten([[k^n for k in range(1,n+1)] for n in range(1,12)]) # G. C. Greubel, Nov 01 2022

Formula

Sum_{k=1..n} T(n, k) = A031971(n).
T(n, n) = A000312(n).
T(2*n, n) = A062206(n).
a(n) = (n + T*(1-T)/2)^T, where T = round(sqrt(2*n),0). - Gerald Hillier, Apr 12 2015
T(n,k) = A051129(n,k). - R. J. Mathar, Dec 10 2015
T(n,k) = Sum_{i=0..k} Stirling2(n,i)*binomial(k,i)*i!. - Geoffrey Critzer, Dec 30 2021
From G. C. Greubel, Nov 01 2022: (Start)
T(n, n-1) = A007778(n-1), n >= 2.
T(n, n-2) = A008788(n-2), n >= 3.
T(2*n+1, n) = A085526(n).
T(2*n-1, n) = A085524(n).
T(2*n-1, n-1) = A085526(n-1), n >= 2.
T(3*n, n) = A083282(n).
Sum_{k=1..n} (-1)^k * T(n, k) = (-1)^n * A120485(n).
Sum_{k=1..floor(n/2)} T(n-k, k) = A226065(n).
Sum_{k=1..floor(n/2)} T(n, k) = A352981(n).
Sum_{k=1..floor(n/3)} T(n, k) = A352982(n). (End)

Extensions

More terms and better definition from Herman Jamke (hermanjamke(AT)fastmail.fm), Jul 10 2004
Offset corrected by Reinhard Zumkeller, Mar 18 2013

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

Views

Author

Keywords

Comments

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

Examples

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

Crossrefs

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

Programs

  • Mathematica
    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 *)
  • PARI
    T(x, y) = x^y \\ Charles R Greathouse IV, Feb 07 2017
    
  • SageMath
    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

Formula

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)

Extensions

New name by Peter Luschny, Apr 16 2024.

A051128 Table T(n,k) = n^k read by upwards antidiagonals (n >= 1, k >= 1).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Sum of antidiagonals is A003101(n) for n>0. - Alford Arnold, Jan 14 2007

Examples

			Table begins
1,    1,    1,    1,    1, ...
2,    4,    8,   16,   32, ...
3,    9,   27,   81,  243, ...
4,   16,   64,  256, 1024, ...
		

Crossrefs

Programs

Formula

a(n) = A004736(n)^A002260(n) or ((t*t+3*t+4)/2-n)^(n-(t*(t+1)/2)), where t=floor((-1+sqrt(8*n-7))/2). - Boris Putievskiy, Dec 14 2012

Extensions

More terms from James Sellers, Dec 11 1999

A003320 a(n) = max_{k=0..n} k^(n-k).

Original entry on oeis.org

1, 1, 1, 2, 4, 9, 27, 81, 256, 1024, 4096, 16384, 78125, 390625, 1953125, 10077696, 60466176, 362797056, 2176782336, 13841287201, 96889010407, 678223072849, 4747561509943, 35184372088832, 281474976710656, 2251799813685248
Offset: 0

Views

Author

Keywords

Comments

For n > 0, a(n+1) = largest term of row n in triangles A051129 and A247358. - Reinhard Zumkeller, Sep 14 2014

Examples

			a(5) = max(5^0, 4^1, 3^2, 2^3, 1^4, 0^5) = max(1,4,9,8,1,0) = 9.
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • I. Tomescu, Introducere in Combinatorica. Editura Tehnica, Bucharest, 1972, p. 231.

Crossrefs

Programs

  • Haskell
    a003320 n = maximum $ zipWith (^) [0 .. n] [n, n-1 ..]
    -- Reinhard Zumkeller, Jun 24 2013
    
  • Mathematica
    Join[{1},Max[#]&/@Table[k^(n-k),{n,25},{k,n}]] (* Harvey P. Dale, Jun 20 2011 *)
  • PARI
    a(n) = vecmax(vector(n+1, k, (k-1)^(n-k+1))); \\ Michel Marcus, Jun 13 2017

Formula

a(n) = A056155(n-1)^(n - A056155(n-1)), for n >= 2. - Ridouane Oudra, Dec 09 2020

Extensions

Easdown reference from Michail Kats (KatsMM(AT)info.sgu.ru)
More terms from James Sellers, Aug 21 2000

A265583 Array T(n,k) = k*(k-1)^(n-1) read by ascending antidiagonals; k,n >= 1.

Original entry on oeis.org

1, 0, 2, 0, 2, 3, 0, 2, 6, 4, 0, 2, 12, 12, 5, 0, 2, 24, 36, 20, 6, 0, 2, 48, 108, 80, 30, 7, 0, 2, 96, 324, 320, 150, 42, 8, 0, 2, 192, 972, 1280, 750, 252, 56, 9, 0, 2, 384, 2916, 5120, 3750, 1512, 392, 72, 10, 0, 2, 768, 8748, 20480, 18750, 9072, 2744, 576, 90, 11
Offset: 1

Views

Author

R. J. Mathar, Dec 10 2015

Keywords

Comments

T(n,k) is the number of n-letter words in a k-letter alphabet with no adjacent letters the same. The factor k represents the number of choices of the first letter, and the n-1 times repeated factor k-1 represents the choices of the next n-1 letters avoiding their predecessor.
The antidiagonal sums are s(d) = 1, 2, 5, 12, 31, 88, 275, 942, 3513, 14158, 61241, 282632, .. for d = n+k >= 2.

Examples

			      1       2       3       4       5       6       7
      0       2       6      12      20      30      42
      0       2      12      36      80     150     252
      0       2      24     108     320     750    1512
      0       2      48     324    1280    3750    9072
      0       2      96     972    5120   18750   54432
      0       2     192    2916   20480   93750  326592
T(3,3)=12 counts aba, abc, aca, acb, bab, bac, bca, bcb, cab, cac, cba, cbc. Words like aab or cbb are not counted.
		

Crossrefs

Cf. A007283 (column 3), A003946 (column 4), A003947 (column 5), A002378 (row 2), A011379 (row 3), A179824 (row 4), A055897 (diagonal), A265584.

Programs

  • GAP
    T:= function(n,k)
        if (n=1 and k=1) then return 1;
        else return k*(k-1)^(n-k-1);
        fi;
      end;
    Flat(List([2..15], n-> List([1..n-1], k-> T(n,k) ))); # G. C. Greubel, Aug 10 2019
  • Magma
    T:= func< n,k | (n eq 1 and k eq 1) select 1 else k*(k-1)^(n-k-1) >;
    [T(n,k): k in [1..n-1], n in [2..15]]; // G. C. Greubel, Aug 10 2019
    
  • Maple
    A265583 := proc(n,k)
        k*(k-1)^(n-1) ;
    end proc:
    seq(seq( A265583(d-k,k),k=1..d-1),d=2..13) ;
  • Mathematica
    T[1,1] = 1; T[n_, k_] := If[k==1, 0, k*(k-1)^(n-1)]; Table[T[n-k,k], {n,2,12}, {k,1,n-1}] // Flatten (* Amiram Eldar, Dec 13 2018 *)
  • PARI
    T(n,k) = if(n==k==1, 1, k*(k-1)^(n-k-1) );
    for(n=2,15, for(k=1,n-1, print1(T(n,k), ", "))) \\ G. C. Greubel, Aug 10 2019
    
  • Sage
    def T(n, k):
        if (n==k==1): return 1
        else: return k*(k-1)^(n-k-1)
    [[T(n, k) for k in (1..n-1)] for n in (2..15)] # G. C. Greubel, Aug 10 2019
    

Formula

T(n,k) = k*A051129(n-1,k-1) = k*A003992(k-1,n-1).
G.f. for column k: k*x/(1-(k-1)*x). - R. J. Mathar, Dec 12 2015
G.f. for array: y/(y-1) - (1+1/x)*y*LerchPhi(y,1,-1/x). - Robert Israel, Dec 13 2018

A292741 Number A(n,k) of partitions of n with k sorts of part 1; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 1, 2, 2, 1, 1, 3, 5, 3, 2, 1, 4, 10, 11, 5, 2, 1, 5, 17, 31, 24, 7, 4, 1, 6, 26, 69, 95, 50, 11, 4, 1, 7, 37, 131, 278, 287, 104, 15, 7, 1, 8, 50, 223, 657, 1114, 865, 212, 22, 8, 1, 9, 65, 351, 1340, 3287, 4460, 2599, 431, 30, 12, 1, 10, 82, 521, 2459, 8042, 16439, 17844, 7804, 870, 42, 14
Offset: 0

Views

Author

Alois P. Heinz, Sep 22 2017

Keywords

Examples

			A(1,3) = 3: 1a, 1b, 1c.
A(2,3) = 10: 2, 1a1a, 1a1b, 1a1c, 1b1a, 1b1b, 1b1c, 1c1a, 1c1b, 1c1c.
A(3,2) = 11: 3, 21a, 21b, 1a1a1a, 1a1a1b, 1a1b1a, 1a1b1b, 1b1a1a, 1b1a1b, 1b1b1a, 1b1b1b.
Square array A(n,k) begins:
  1,  1,   1,    1,     1,     1,      1,      1, ...
  0,  1,   2,    3,     4,     5,      6,      7, ...
  1,  2,   5,   10,    17,    26,     37,     50, ...
  1,  3,  11,   31,    69,   131,    223,    351, ...
  2,  5,  24,   95,   278,   657,   1340,   2459, ...
  2,  7,  50,  287,  1114,  3287,   8042,  17215, ...
  4, 11, 104,  865,  4460, 16439,  48256, 120509, ...
  4, 15, 212, 2599, 17844, 82199, 289540, 843567, ...
		

Crossrefs

Columns k=0-2 give: A002865, A000041, A090764.
Rows n=0-2 give: A000012, A001477, A002522, A071568.
Main diagonal gives A292462.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0 or i<2, k^n,
          add(b(n-i*j, i-1, k), j=0..iquo(n, i)))
        end:
    A:= (n, k)-> b(n$2, k):
    seq(seq(A(n, d-n), n=0..d), d=0..14);
  • Mathematica
    b[0, , ] = 1; b[n_, i_, k_] := b[n, i, k] = If[i < 2, k^n, Sum[b[n - i*j, i - 1, k], {j, 0, Quotient[n, i]}]];
    A[n_, k_] := b[n, n, k];
    Table[A[n, d - n], {d, 0, 14}, {n, 0, d}] // Flatten (* Jean-François Alcover, May 19 2018, translated from Maple *)

Formula

G.f. of column k: 1/(1-k*x) * 1/Product_{j>=2} (1-x^j).
A(n,k) = Sum_{j=0..n} A002865(j) * k^(n-j).

A332101 Least m such that m^n <= Sum_{k

Original entry on oeis.org

2, 3, 5, 6, 8, 9, 11, 12, 14, 15, 16, 18, 19, 21, 22, 24, 25, 27, 28, 29, 31, 32, 34, 35, 37, 38, 40, 41, 42, 44, 45, 47, 48, 50, 51, 52, 54, 55, 57, 58, 60, 61, 63, 64, 65, 67, 68, 70, 71, 73, 74, 76, 77, 78, 80, 81, 83, 84, 86, 87, 89, 90, 91, 93, 94, 96, 97
Offset: 0

Views

Author

M. F. Hasler, Apr 14 2020

Keywords

Comments

In a list (1^n, 2^n, 3^n, ...) (rows of table A051128 or A051129), a(n) is the index of the first term less than or equal to the sum of all earlier terms, cf. example.
Obviously a lower bound for any s solution to s^n = Sum_{x in S} x^n, S subset of {1, ..., s-1}, cf. A030052.

Examples

			For n = 0, m^0 > Sum_{0 < k < m} k^0 = 0 for m = 0, 1 (empty sums), but 2^0 = Sum_{0 < k < 2} k^0 = 1, so a(0) = 2.
For n = 1, 1^1 > Sum_{0 < k < 1} k^1 = 0 (empty sum) and 2^1 > Sum_{0 < k < 2} k^1 = 1, but 3^1 <= Sum_{0 < k < 3} k^1 = 1 + 2, so a(1) = 3.
To find a(n) one can add up terms in row n of the table k^n until the sum equals or exceeds the next term, whose column number k is then a(n):
  n |k: 1  2   3   4    5    6          Comment
  --+---------------------------------------------------------------
  1 |  1   2   3                  1 < 2 but 1 + 2 >= 3, so a(1) = 3.
  2 |  1   4   9  16   25         1 + 4 + 9 + 16 > 25, and a(2) = 5.
  3 |  1   8  27  64  125  216    1 + 8 + 27 + 64 + 125 > 216: a(3) = 6.
		

Crossrefs

Cf. A078607, A332097 (maximum of E(s), cf comments), A030052 (least k such that k^n = sum of distinct n-th powers), A332065 (all k such that k^n is a sum of distinct n-th powers).

Programs

  • Mathematica
    Table[Block[{m = 1, s = 0}, While[m^n > s, s = s + m^n; m++]; m], {n, 0, 66}] (* Michael De Vlieger, Apr 30 2020 *)
  • PARI
    apply( A332101(n,s)=for(m=1,oo, s
    				

Formula

a(n) = round(n / log(2)) + 2. (Conjectured; verified up to 10^4, in particular for 3525/log(2) = 5085.500019... and 7844/log(2) ~ 11316.49990...)
a(n) = A078607(n) + 2 for almost all n > 1. (n = 777451915729368 might be an exception to this equality or the above one.) - M. F. Hasler, May 08 2020

A247358 Triangle read by rows: n-th row contains powers b^e with b + e = n + 1 in natural order.

Original entry on oeis.org

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

Views

Author

Reinhard Zumkeller, Sep 14 2014

Keywords

Comments

Sorted rows of triangle A051129.

Examples

			.  1  |  1                            |  1^1
.  2  |  1 2                          |  1^2 2^1
.  3  |  1 3  4                       |  1^3 3^1 2^2
.  4  |  1 4  8   9                   |  1^4 4^1 2^3 3^2
.  5  |  1 5 16  16  27               |  1^5 5^1 2^4 4^2 3^3
.  6  |  1 6 25  32  64  81           |  1^6 6^1 5^2 2^5 4^3 3^4
.  7  |  1 7 36  64 125 243 256       |  1^7 7^1 6^2 2^6 5^3 3^5 4^4
.  8  |  1 8 49 128 216 625 729 1024  |  1^8 8^1 7^2 2^7 6^3 5^4 3^6 4^5 .
		

Crossrefs

Cf. A051129, A003101 (row sums), A247363 (central terms), A003320 (row maxima).

Programs

  • Haskell
    import Data.List (sort)
    a247358 n k = a247358_tabl !! (n-1) !! (k-1)
    a247358_row n = a247358_tabl !! (n-1)
    a247358_tabl = map sort a051129_tabl
    
  • Mathematica
    Table[Table[k^(n-k+1), {k, 1, n}] // Sort, {n, 1, 11}] // Flatten (* Jean-François Alcover, Nov 18 2019 *)
  • PARI
    row(n) = vecsort(vector(n, k, k^(n-k+1))); \\ Michel Marcus, Jan 24 2022
  • Python
    from itertools import chain
    A247358_list = list(chain.from_iterable(sorted((b+1)**(n-b) for b in range(n)) for n in range(1,8))) # Chai Wah Wu, Sep 14 2014
    
Showing 1-10 of 17 results. Next