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-6 of 6 results.

A022811 Number of terms in n-th derivative of a function composed with itself 3 times.

Original entry on oeis.org

1, 1, 3, 6, 13, 23, 44, 74, 129, 210, 345, 542, 858, 1310, 2004, 2996, 4467, 6540, 9552, 13744, 19711, 27943, 39452, 55172, 76865, 106200, 146173, 199806, 272075, 368247, 496642, 666201, 890602, 1184957, 1571417, 2075058, 2731677, 3582119, 4683595, 6102256
Offset: 0

Views

Author

Winston C. Yang (yang(AT)math.wisc.edu)

Keywords

Comments

This also counts a restricted set of plane partitions of n. Each element of the set which contains the A000041(n) partitions of n can be converted into plane partitions by insertion of line feeds at some or all places of the "pluses." Since the length of rows in plane partitions must be nonincreasing, there are only A000041(L(P)) ways to comply with this rule, where L(P) is the number of terms in that particular partition. Example for n=4: consider all five partitions 4 = 3+1 = 2+2 = 2+1+1 = 1+1+1+1 of four. The associated a(4)=13 plane partitions are 4, 31, 3|1, 22, 2|2, 211, 21|1, 2|1|1, 1111, 111|1, 11|11, 11|1|1 and 1|1|1|1, where the bar represents start of the next row, where a(4) = A000041(L(4)) + A000041(L(3+1)) + A000041(L(2+2)) + A000041(L(2+1+1))+ A000041(L(1+1+1+1)) = A000041(1) + A000041(2) + A000041(2) + A000041(3) + A000041(4). By construction from sorted partitions, all the plane partitions are strictly decreasing along each row and each column. - R. J. Mathar, Aug 12 2008
Also the number of pairs of integer partitions, the first with sum n and the second with sum equal to the length of the first. - Gus Wiseman, Jul 19 2018

Examples

			From _Gus Wiseman_, Jul 19 2018: (Start)
Using the chain rule, we compute the second derivative of f(f(f(x))) to be the following sum of a(2) = 3 terms.
  d^2/dx^2 f(f(f(x))) =
  f'(f(x)) f'(f(f(x))) f''(x) +
  f'(x)^2 f'(f(f(x))) f''(f(x)) +
  f'(x)^2 f'(f(x))^2 f''(f(f(x))).
(End)
		

References

  • W. C. Yang, Derivatives of self-compositions of functions, preprint, 1997.

Crossrefs

Column k=3 of A022818.
First column of A039805.
A row or column of A081718.

Programs

  • Maple
    A022811 := proc(n) local a,P,p,lp ; a := 0 ; P := combinat[partition](n) ; for p in P do lp := nops(p) ; a := a+combinat[numbpart](lp) ; od: RETURN(a) ; end: for n from 1 do print(n,A022811(n)) ; od: # R. J. Mathar, Aug 12 2008
  • Mathematica
    a[n_] := Total[PartitionsP[Length[#]]& /@ IntegerPartitions[n]];
    Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 0, 80}] (* Jean-François Alcover, Apr 28 2017 *)
    Table[Length[1+D[f[f[f[x]]],{x,n}]]-1,{n,10}] (* Gus Wiseman, Jul 19 2018 *)

Formula

If a(n,m) = number of terms in m-derivative of a function composed with itself n times, p(n,k) = number of partitions of n into k parts, then a(n,m) = Sum_{i=0..m} p(m,i)*a(n-1,i).
G.f.: Sum_{k>=0} p(k) * x^k / Product_{j=1..k} (1 - x^j), where p(k) = number of partitions of k. - Ilya Gutkovskiy, Jan 28 2020

Extensions

Typo corrected by Neven Juric, Mar 25 2013

A131407 Repeated set partitions or nested set partitions. Possible coalitions among n persons.

Original entry on oeis.org

1, 1, 2, 11, 95, 1307, 27035, 788279, 30812087, 1554832679, 98387784047, 7628836816295, 711320467520855, 78520062277781087, 10127079289703949695, 1508987827451079129599, 257250406707409951420079, 49750955749787132205813743, 10833471589449269308161546191
Offset: 0

Views

Author

Thomas Wieder, Jul 09 2007, Jul 20 2007

Keywords

Comments

Consider a set N={1,2,3,...n}. We can apply the operation S~(N) on N which gives us the set partitions S~(N)=SP(N) of N. Let denote SP_i(N) such a set partition, then SP(N)={SP_1(N), SP_2(N)...,SP_B(n)} (There are B(n) set partitions of N with B(n) as the Bell number). Observe that in each SP(N) we have SP(1)={{1,2,3,...,n}} and SP(B(n))={{1},{2},{3},...,{n}} and their magnitudes are |SP(1)|=1 and |SP(B(n)|=n.
Now we perform an iteration on the set partitions SP_i(N). We set partition each SP_i(N), thus we perform S~(SP_i(N), but we exclude SP(1)={{1,2,3,...,n}} and SP(B(n))={{1},{2},{3},...,{n}} from this repetition. Otherwise an infinite recursion arises. Thus if 1 < |SP_i(N)|=m < n, then we apply S~ on SP_i(N) again and get S~(SP_i(N))= SP(SP_i(N))={SP_1(SP_i(N)),...,SP_B(m)(SP_i(N))}. We repeat this partition operation S~ on every set partition we encounter. Let denote U_k a subset of SP_i(X) were X is a set. X may be any of the subsequent set partitions. Since |U_k| < X (under the condition above on m) the repeated application of S~ will end in set partitions SP(X) with |SP(X)| = 1.
Let us consider the example N={1,2,3}. The S~(N) gives us {{1,2,3}}, {{1,2},{3}}, {{1,3},{2}}, {{2,3},{13}} and {{1},{2},{1}}. We exclude {{1,2,3}} and {{1},{2},{1}} from further partitioning. From {{1,2},{3}} we get {{{1,2},{3}}} and {{{1,2}},{{3}}}. Consider the last two partitions. They correspond to N'={1',2'} and are thus {{1',2'}} and {{1'},{2'}}. Since |{{1',2'}}|=1 and |{{1'},{2'}}|=2 these last two set partitions cannot be partitioned any further according to our condition above. In total we get {{1,2,3}}, {{1,2},{3}}, {{1,3},{2}}, {{2,3},{1}}, {{{1,2},{3}}}, {{{1,2}},{{3}}}, {{{1,3},{2}}}, {{{1,3}},{{2}}}, {{{2,3},{1}}}, {{{2,3}},{{1}}}, {{1},{2},{3}} and we have a(3)=11.
A possible application are the number of coalitions among the set N={1,2,...,n} of n persons. These persons will split into parties = subsets U_k of N. Then coalitions will form among these parties, thus we encounter sets of subsets. It is even possible that coalitions form coalitions in turn. We thus define a coalition structure as a set of repeated set partitions. For example if n=6 we could have {{1,2},{3}},{{4,5,6}}, the parties {1,2} and {3} form the coalition {{1,2},{3}}. Since {{456}}={4,5,6} one might not want to consider a single set as a coalition, but formally it is possible to do so. However, if in the example all three parties are patriotic, they may stand together in questions of national interest and the coalition structure would be {{{1,2},{3}},{{4,5,6}}}.
However, in my opinion, the usual definition of a coalition as a partition of a set falls too short.
See also A005121 = Ultradissimilarity relations on an n-set. The paper "On the Asymptotic Analysis of a Class of Linear Recurrences" (by Thomas Prellberg) outlines how to find an asymptotic formula for A005121. Perhaps this method is applicable to the present sequence as well, but one needs to have the generating function as starting point.

Examples

			a(3)=11 because we have
  {{1,2,3}},
  {{1,2},{3}},
  {{1,3},{2}},
  {{2,3},{1}},
  {{{1,2},{3}}},
  {{{1,2}},{{3}}},
  {{{1,3},{2}}},
  {{{1,3}},{{2}}},
  {{{2,3},{1}}},
  {{{2,3}},{{1}}},
  {{1},{2},{3}}.
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, p. 319 and 556.

Crossrefs

Programs

  • Maple
    rctlnn := proc(n::nonnegint) # Thanks to Joe Riel, who suggested the use of # "procname" instead of "rctlnn" within the program.
    local j; option remember; if n = 0 then 1; else bell(n)+add(stirling2(n,j)*procname(j), j=2..n-1); end if; end proc:
    # second Maple program:
    a:= proc(n) option remember; uses combinat;
          bell(n) + add(stirling2(n, i)*a(i), i=2..n-1)
        end:
    seq(a(n), n=0..20);  # Alois P. Heinz, Apr 05 2012
  • Mathematica
    a[n_] := a[n] = If[n<2, 1, BellB[n] + Sum[StirlingS2[n, i]*a[i], {i, 2, n-1}]]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jul 15 2015, after Alois P. Heinz *)

Formula

Recurrence: a(n) = Bell(n) + Sum_{i=2..n-1} S2(n,i)*a(i). E.g.: a(n=4) = Bell(4) + S2(4,2) a(2) + S2(4,3) a(3) = 15+2+7*2+6*11 = 95. "closed" formula: a(n=4) = Bell(n=4) + Sum_{i1=2..(n=4)-1} Bell(i1) + S2(n,i1)*Sum_{i2=2..i1-1} Bell(i2) + S2(i1,i2)*Sum_{i3=2..i2-1} Bell(i3) + S2(i2,i3)*Sum_{i4=2..i3-1} Stirling2(i3,i4).
a(n) ~ 3 * L * (n!)^2 / (n^(1+log(2)/3) * (2*log(2))^n), where L = Lengyel's constant A086053 = 1.0986858055... . - Vaclav Kotesovec, Sep 04 2014

Extensions

a(0)=1 prepended by Alois P. Heinz, Sep 02 2020

A246828 Decimal expansion of a constant related to A055887.

Original entry on oeis.org

2, 6, 9, 8, 3, 2, 9, 1, 0, 6, 4, 7, 4, 2, 1, 1, 2, 3, 1, 2, 6, 3, 9, 9, 8, 6, 6, 6, 1, 8, 8, 3, 7, 6, 3, 3, 0, 7, 1, 3, 4, 6, 5, 1, 2, 5, 9, 1, 3, 9, 8, 6, 3, 5, 6, 7, 6, 9, 0, 1, 2, 3, 1, 1, 7, 8, 1, 9, 8, 6, 5, 9, 3, 6, 6, 9, 5, 0, 5, 5, 9, 4, 5, 1, 3, 6, 6, 4, 7, 6, 6, 5, 2, 0, 2, 2, 0, 3, 5, 5, 8, 0, 0, 7, 7
Offset: 1

Views

Author

Vaclav Kotesovec, Sep 04 2014

Keywords

Examples

			2.698329106474211231263998666188376330713465125913986356769...
		

Crossrefs

Programs

  • Mathematica
    RealDigits[1/x /. FindRoot[QPochhammer[x] == 1/2, {x, 1/2}, WorkingPrecision -> 120]][[1]] (* Vaclav Kotesovec, May 23 2018 *)

Formula

Equals lim n -> infinity A055887(n)^(1/n).
Equals lim n -> infinity A095975(n)^(1/n).
Equals lim n -> infinity A141799(n)^(1/n).
Equals lim n -> infinity A131408(n)^(1/n).
Root of the equation QPochhammer[x] = 1/2. - Vaclav Kotesovec, May 23 2018

A141799 Number of repeated integer partitions of n.

Original entry on oeis.org

1, 3, 8, 25, 66, 192, 511, 1418, 3812, 10383, 27958, 75758, 204215, 551821, 1488561, 4018722, 10842422, 29262357, 78955472, 213063551, 574905487, 1551325859, 4185959285, 11295211039, 30478118079, 82240300045, 221911189754, 598790247900, 1615732588962
Offset: 1

Views

Author

Thomas Wieder, Jul 05 2008

Keywords

Comments

An integer n can be partitioned into P(n) partitions P([n],i) where i=1,...,P(n) counts the partitions. The partition P([n],i) consists of T(n,i) integer parts t(i,j) with j=1,...,T(n,i). Now we perform on each t(i,j) an integer partition again and arrive at new partitions. Their parts can be partitioned again and so forth. We count such repeated partitions of n. One convention is necessary to avoid an infinite loop: The trivial partition P([n],1)=[n] will not be partitioned again but just counted once (and therefore we also have a(1)=1).

Examples

			For the integers 1, 2, 3 and 4 we have
[1] -> 1,
thus a(1)=1.
[2] -> 1,
[1,1] => [1] ->, [1] -> 1.
thus a(2)=3.
[3] -> 1,
[1,2] => [1] -> 1, [2] -> 3,
[1,1,1] => [1] -> 1, [1] -> 1, [1] -> 1,
thus a(3)=8.
[4] -> 1,
[1,3] => [1] -> 1, [3] -> 8,
[2,2] => [2] -> 3, [2] -> 3,
[1,1,2] => [1] -> 1, [1] -> 1, [2] -> 3,
[1,1,1,1] => [1] -> 1, [1] -> 1, [1] -> 1, [1] -> 1,
thus a(4)=25.
		

Crossrefs

Programs

  • Maple
    A141799 := proc(n) option remember ; local a,P,i,p ; if n =1 then 1; else a := 0 ; for P in combinat[partition](n) do if nops(P) > 1 then for i in P do a := a+procname(i) ; od: else a := a+1 ; fi; od: RETURN(a) ; fi ; end: for n from 1 to 40 do printf("%d,",A141799(n)) ; od: # R. J. Mathar, Aug 25 2008
    # second Maple program
    a:= proc(n) option remember;
          1+ `if`(n>1, b(n, n-1)[2], 0)
        end:
    b:= proc(n, i) option remember; local f, g;
          if n=0 or i=1 then [1, n]
        else f:= b(n, i-1); g:= `if`(i>n, [0, 0], b(n-i, i));
             [f[1]+g[1], f[2]+g[2] +g[1]*a(i)]
          fi
        end:
    seq(a(n), n=1..40); # Alois P. Heinz, Apr 05 2012
  • Mathematica
    a[n_] := a[n] = 1 + If[n>1, b[n, n-1][[2]], 0]; b[n_, i_] := b[n, i] = Module[{f, g}, If[n == 0 || i == 1, {1, n}, f = b[n, i-1]; g = If[i>n, {0, 0}, b[n-i, i]]; {f[[1]] + g[[1]], f[[2]] + g[[2]] + g[[1]]*a[i]}]]; Table[a[n], {n, 1, 40}] (* Jean-François Alcover, Oct 29 2015, after Alois P. Heinz *)

Formula

Let sum_{i=1}^P(n) denote the sum over all integer partitions P([n],i) of n. Let sum_{j=1}^T(i,j) denote the sum over all parts of the i-th integer partition. Then we have the recursive formula 1 if t(i,j)=n a(n) = sum_{i=1}^P(n) sum_{j=1}^T(i,j) { a(t(i,j)) else. E.g. a(4)=25 because [4] contributes 1, [1,3] contributes a(1)+a(3)=1+8=9, [2,2] contributes a(2)+a(2)=3+3=6, [1,1,2] contributes a(1)+a(1)+a(2)=1+1+3=5, [1,1,1,1] contributes a(1)+a(1)+a(1)+a(1)=1+1+1+1=4 which gives in total 25.
a(n) ~ c * d^n, where d = 2.69832910647421123126399866... (see A246828), c = 0.5088820425072641934222229579416714164592334575899644931509447692360546... . - Vaclav Kotesovec, Sep 04 2014

Extensions

Extended by R. J. Mathar, Aug 25 2008

A143141 Total number of all repeated partitions of the integer n and its parts down to parts equal to 1. Essentially first differences of A055887.

Original entry on oeis.org

1, 2, 5, 14, 37, 101, 271, 733, 1976, 5334, 14390, 38833, 104779, 282734, 762903, 2058571, 5554692, 14988400, 40443620, 109130216, 294469216, 794574883, 2144024501, 5785283758, 15610599502, 42122535067, 113660462337, 306693333868, 827559549428, 2233028019698
Offset: 1

Views

Author

Thomas Wieder, Jul 27 2008

Keywords

Comments

Start from the A000041(n) integer partitions P(n,i,s) of the integer n at stage s=1.
The index i=1,...,A000041(n) denotes the different partitions.
We call the index s the partition stage and increase it by one as we sub-partition the partitions of a previous stage.
Each P(n,i,s) is a set P(n,i,s)={t(n,1,j,s)),...,t(P,i,j,s),...,t(P,i,J,s)} of parts t(P,i,j,s) of S.
The index j is attached to the parts of a partition P(n,i,s). 1<=j<=n since there are at most n parts.
Now apply the set partition process on every P(n,i,s=1).
That is, each t(n,i,j,s=1) is subjected to a further partitioning.
We get partitions P(t'(n,i,j,1),i',j',2)={t'(t(n,i,j,1),i',1,2),...,t'(t(n,i,j,1), i',j',2),...,t'(t(n,i,j,1),i',J',2)} of the second partition stage.
We repeat this partitioning process on each part t'(i,j',2) until we arrive at parts equal to 1 which cannot be partitioned any further.
We may speak of the full decomposition F of n into parts.
The sequence counts the total number of partitions of all stages of the full decomposition of n.
Note that n is its own partition, e.g. P(n=3,i=1,s=1)={3} is an integer partition of n=3.
We do not apply the repeated partitioning on the partition P(n,i,s)={n} (otherwise an infinite loop would arise).
For n=1 and n=2 there is no second partition stage: s stays at s=1.
The corresponding labeled counterpart is sequence A143140.

Examples

			n=1:
[[1]]
n=2:
[[2], [1, 1]]
n=3:
[[3], [2, 1], [1, 1, 1]], [[2], [1, 1]]
n=4 in more detail:
[4], [3, 1], [2, 2], [2, 1, 1], [1, 1, 1, 1]], <- stage s=1, partition of 4
[[3], [2, 1], [1, 1, 1]], <- stage s=2 partitioning the first 3 of the 2nd partition
[[2], [1, 1]], <- stage s=2 partitioning the first 2 of the 3rd partition
[[2], [1, 1]], <- stage s=2 partitioning the second 2 of the 3rd partition
[[2], [1, 1]] <- stage s=2 partitioning the first 2 of the 4th partition
a(4) = 14 = 5 (from s=1)+9 (from s=2).
		

Crossrefs

Programs

  • Maple
    A055887 := proc(n) option remember ; if n = 0 then 1; else add(combinat[numbpart](k)*procname(n-k),k=1..n) ; fi; end: A143141 := proc(n) if n = 1 then 1; else A055887(n)-A055887(n-1) ; fi; end: seq(A143141(n),n=1..20) ;
  • Mathematica
    b[n_] := b[n] = Sum[PartitionsP[k]*b[n-k], {k, 1, n}]; b[0]=1; A055887 = Table[b[n], {n, 0, 30}]; Join[{1}, Rest[Differences[A055887]]] (* Jean-François Alcover, Feb 05 2017 *)

Formula

a(n) = A055887(n) - A055887(n-1), n>1.

Extensions

Edited and extended by R. J. Mathar, Aug 25 2008

A131965 a(n) = 1 + Sum_{i=2..n-1} n*a(i).

Original entry on oeis.org

1, 1, 1, 4, 21, 131, 943, 7701, 70409, 712891, 7921011, 95844233, 1254688141, 17670191319, 266412115271, 4281623281141, 73073037331473, 1319881736799731, 25155393101359579, 504505383866156001, 10621165976129600021, 234196709773657680463, 5397676549069062730671
Offset: 0

Views

Author

Thomas Wieder, Aug 02 2007

Keywords

Comments

a(n) = 1 + Sum_{i=2..n-1} 1*a(i) = 2^n; a(n) = 1 + Sum_{i=2..n-1} 2*a(i) = 3^n; etc. It seems that a(n+1)/(n*a(n)) -> 1 for n -> oo. [Comment corrected by Emeric Deutsch, Aug 10 2007]
Let M(n) denote the n X n matrix with ones along the subdiagonal, ones everywhere above the main diagonal, the integers 4, 5, etc., along the main diagonal, and zeros everywhere else. Then a(n) equals the permanent of M(n-2) for n >= 3. - John M. Campbell, Apr 20 2021

Examples

			a(4)=21 because 1 + 4*1 + 4*4 = 21.
		

Crossrefs

Programs

  • Magma
    m:=25; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( (-8*(1+x) + 2*(3-x)*Exp(x) + (4+3*x^2-x^3))/(2*(1-x)^3) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Mar 09 2019
    
  • Maple
    rctlnn := proc(n::nonnegint) local j; option remember; if n = 0 then 0; else 1+add(n*procname(j), j=2..n-1); end if; end proc:
    a[1] := 1; for n from 2 to 18 do a[n] := 1+sum(n*a[i], i = 2 .. n-1) end do: seq(a[n], n = 1 .. 18); # Emeric Deutsch, Aug 10 2007
    # third Maple program:
    a:= proc(n) option remember;
          1+add(n*a(i), i=2..n-1)
        end:
    seq(a(n), n=0..30);  # Alois P. Heinz, Sep 03 2020
  • Mathematica
    a[1] = a[2] = 1; a[n_] := a[n] = (n^2*a[n-1]-1)/(n-1); Array[a, 30] (* Jean-François Alcover, Feb 08 2017 *)
  • Sage
    m = 25; T = taylor((-8*(1+x) + 2*(3-x)*exp(x) + (4+3*x^2-x^3))/(2*(1-x)^3), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Mar 09 2019

Formula

a(n) = 1 + Sum_{i=2..n-1} n*a(i).
E.g.f.: 1/2 * (x + (2*exp(x)-5)/(x-1)^2 -5/(x-1)).
Asymptotic expansion: a(n)/n! = (5/2 + e)*n^2 + O(n).
a(n) = (n+1)*a(n-1) + a(n-2) + ... + a(2), e.g., a(5) = 6*21 + 4 + 1 = 131.

Extensions

More terms from Emeric Deutsch, Aug 10 2007
a(0)=1 prepended and edited by Alois P. Heinz, Sep 03 2020
Showing 1-6 of 6 results.