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

A133494 Diagonal of the array of iterated differences of A047848.

Original entry on oeis.org

1, 1, 3, 9, 27, 81, 243, 729, 2187, 6561, 19683, 59049, 177147, 531441, 1594323, 4782969, 14348907, 43046721, 129140163, 387420489, 1162261467, 3486784401, 10460353203, 31381059609, 94143178827, 282429536481, 847288609443, 2541865828329, 7625597484987, 22876792454961, 68630377364883
Offset: 0

Views

Author

Paul Barry, Paul Curtz, Dec 23 2007

Keywords

Comments

a(n) is the number of ways to choose a composition C, and then choose a composition of each part of C. - Geoffrey Critzer, Mar 19 2012
a(n) is the top left entry of the n-th power of the 3 X 3 matrix [1, 1, 1; 1, 1, 1; 1, 1, 1]. - R. J. Mathar, Feb 03 2014
a(n) is the reptend length of 1/3^(n+1) in decimal. - Jianing Song, Nov 14 2018
Also the number of pairs of integer compositions, the first summing to n and the second with sum equal to the length of the first. If an integer composition is regarded as an arrow from sum to length, these are composable pairs, and the obvious composition operation founds a category of integer compositions. For example, we have (2,1,1,4) . (1,2,1) . (1,2) = (2,6), where dots represent the composition operation. The version without empty compositions is A000244. Composable triples are counted by 1 followed by A000302. The unordered version is A022811. - Gus Wiseman, Jul 14 2022

Examples

			From _Gus Wiseman_, Jul 15 2020: (Start)
The a(0) = 1 through a(3) = 9 ways to choose a composition of each part of a composition:
  ()  (1)  (2)      (3)
           (1,1)    (1,2)
           (1),(1)  (2,1)
                    (1,1,1)
                    (1),(2)
                    (2),(1)
                    (1),(1,1)
                    (1,1),(1)
                    (1),(1),(1)
(End)
		

Crossrefs

The strict version is A336139.
Splittings of partitions are A323583.
Multiset partitions of partitions are A001970.
Partitions of each part of a partition are A063834.
Compositions of each part of a partition are A075900.
Strict partitions of each part of a strict partition are A279785.
Compositions of each part of a strict partition are A304961.
Strict compositions of each part of a composition are A307068.
Compositions of each part of a strict composition are A336127.

Programs

Formula

Binomial transform of A078008. - Paul Curtz, Aug 04 2008
From R. J. Mathar, Nov 11 2008: (Start)
G.f.: (1 - 2*x)/(1 - 3*x).
a(n) = A000244(n-1), n > 0. (End)
From Philippe Deléham, Nov 13 2008: (Start)
a(n) = Sum_{k=0..n} A112467(n,k)*2^k.
a(n) = Sum_{k=0..n} A071919(n,k)*2^k. (End)
Let A(x) be the g.f. Then B(x) = x*A(x) satisfies B(x/(1-x)) = x/(1 - 2*B(x)). - Vladimir Kruchinin, Dec 05 2011
G.f.: 1/(1 - (Sum_{k>=1} (x/(1 - x))^k)). - Joerg Arndt, Sep 30 2012
For n > 0, a(n) = 2*(Sum_{k=0..n-1} a(k)) - 1 = 3^(n-1). - J. Conrad, Oct 29 2015
G.f.: 1 + x/(1 + x)*(1 + 4*x/(1 + 4*x)*(1 + 7*x/(1 + 7*x)*(1 + 10*x/(1 + 10*x)*(1 + .... - Peter Bala, May 27 2017
Invert transform of A011782(n) = 2^(n-1). Second invert transform of A000012. - Gus Wiseman, Jul 19 2020
a(n) = ceiling(3^(n-1)). - Alois P. Heinz, Jul 26 2020
From Elmo R. Oliveira, Mar 31 2025: (Start)
E.g.f.: (2 + exp(3*x))/3.
a(n) = 3*a(n-1) for n > 1. (End)

Extensions

Definition clarified by R. J. Mathar, Nov 11 2008

A075900 Expansion of g.f.: Product_{n>0} 1/(1 - 2^(n-1)*x^n).

Original entry on oeis.org

1, 1, 3, 7, 19, 43, 115, 259, 659, 1523, 3731, 8531, 20883, 47379, 113043, 259219, 609683, 1385363, 3245459, 7344531, 17028499, 38579603, 88585619, 199845267, 457864595, 1028904339, 2339763603, 5256820115, 11896157587, 26626389395
Offset: 0

Views

Author

N. J. A. Sloane, Oct 15 2002

Keywords

Comments

Number of compositions of partitions of n. a(3) = 7: 3, 21, 12, 111, 2|1, 11|1, 1|1|1. - Alois P. Heinz, Sep 16 2019
Also the number of ways to split an integer composition of n into consecutive subsequences with weakly decreasing (or increasing) sums. - Gus Wiseman, Jul 13 2020
This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = 1, g(n) = 2^(n-1). - Seiichi Manyama, Aug 22 2020

Examples

			From _Gus Wiseman_, Jul 13 2020: (Start)
The a(0) = 1 through a(4) = 19 splittings:
  ()  (1)  (2)      (3)          (4)
           (1,1)    (1,2)        (1,3)
           (1),(1)  (2,1)        (2,2)
                    (1,1,1)      (3,1)
                    (2),(1)      (1,1,2)
                    (1,1),(1)    (1,2,1)
                    (1),(1),(1)  (2,1,1)
                                 (2),(2)
                                 (3),(1)
                                 (1,1,1,1)
                                 (1,1),(2)
                                 (1,2),(1)
                                 (2),(1,1)
                                 (2,1),(1)
                                 (1,1),(1,1)
                                 (1,1,1),(1)
                                 (2),(1),(1)
                                 (1,1),(1),(1)
                                 (1),(1),(1),(1)
(End)
		

Crossrefs

Row sums of A327549.
The strict case is A304961.
Partitions of partitions are A001970.
Splittings with equal sums are A074854.
Splittings of compositions are A133494.
Splittings of partitions are A323583.
Splittings with distinct sums are A336127.
Starting with a reversed partition gives A316245.
Starting with a partition instead of composition gives A336136.

Programs

  • Magma
    m:=80;
    R:=PowerSeriesRing(Integers(), m);
    Coefficients(R!( 1/(&*[1-2^(j-1)*x^j: j in [1..m+2]]) )); // G. C. Greubel, Jan 25 2024
    
  • Maple
    oo := 101; t1 := mul(1/(1-x^n/2),n=1..oo): t2 := series(t1,x,oo-1): t3 := seriestolist(t2): A075900 := n->2^n*t3[n+1];
    with(combinat); A075900 := proc(n) local i,t1,t2,t3; t1 := partition(n); t2 := 0; for i from 1 to nops(t1) do t3 := t1[i]; t2 := t2+2^(n-nops(t3)); od: t2; end;
  • Mathematica
    b[n_]:= b[n]= Sum[d*2^(n - n/d), {d, Divisors[n]}];
    a[0]= 1; a[n_]:= a[n]= 1/n*Sum[b[k]*a[n-k], {k,n}];
    Table[a[n], {n,0,30}] (* Jean-François Alcover, Mar 20 2014, after Vladeta Jovovic, fixed by Vaclav Kotesovec, Mar 08 2018 *)
  • Maxima
    s(m,n):=if nVladimir Kruchinin, Sep 06 2014 */
    
  • PARI
    {a(n)=polcoeff(prod(k=1,n,1/(1-2^(k-1)*x^k+x*O(x^n))),n)} \\ Paul D. Hanna, Jan 13 2013
    
  • PARI
    {a(n)=polcoeff(exp(sum(k=1,n+1,x^k/(k*(1-2^k*x^k)+x*O(x^n)))),n)} \\ Paul D. Hanna, Jan 13 2013
    
  • SageMath
    m=80;
    def A075900_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( 1/product(1-2^(j-1)*x^j for j in range(1,m+1)) ).list()
    A075900_list(m) # G. C. Greubel, Jan 25 2024

Formula

a(n) = Sum_{ partitions n = c_1 + ... + c_k } 2^(n-k). If p(n, m) = number of partitions of n into m parts, a(n) = Sum_{m=1..n} p(n, m)*2^(n-m).
G.f.: Sum_{n>=0} (a(n)/2^n)*x^n = Product_{n>0} 1/(1-x^n/2). - Vladeta Jovovic, Feb 11 2003
a(n) = 1/n*Sum_{k=1..n} A080267(k)*a(n-k). - Vladeta Jovovic, Feb 11 2003
G.f.: exp( Sum_{n>=1} x^n / (n*(1 - 2^n*x^n)) ). - Paul D. Hanna, Jan 13 2013
a(n) = s(1,n), a(0)=1, where s(m,n) = Sum_{k=m..n/2} 2^(k-1)*s(k, n-k) + 2^(n-1), s(n,n) = 2^(n-1), s(m,n)=0, m>. - Vladimir Kruchinin, Sep 06 2014
a(n) ~ 2^(n-2) * (Pi^2 - 6*log(2)^2)^(1/4) * exp(sqrt((Pi^2 - 6*log(2)^2)*n/3)) / (3^(1/4) * sqrt(Pi) * n^(3/4)). - Vaclav Kotesovec, Mar 09 2018

Extensions

More terms from Vladeta Jovovic, Feb 11 2003

A098407 Number of different hierarchical orderings that can be formed from n unlabeled elements with no repetition of subhierarchies.

Original entry on oeis.org

1, 1, 2, 6, 13, 33, 78, 186, 436, 1028, 2394, 5566, 12877, 29689, 68198, 156194, 356599, 811959, 1843956, 4177436, 9442166, 21295934, 47932572, 107677140, 241443980, 540441068, 1207689636, 2694452060, 6002389882, 13351958546, 29659179804, 65794744420, 145768641091
Offset: 0

Views

Author

Thomas Wieder, Sep 07 2004; corrected Sep 09 2004

Keywords

Comments

a(n) is the number of finite sets of compositions with total sum n. The case of constant sums is A358904, cf. A074854. The case of distinct sums is A304961, ordered A336127. The ordered version (sequences of distinct compositions) is A358907. - Gus Wiseman, Dec 12 2022

Examples

			Let a pair of parentheses () indicate a subhierarchy and let square brackets [] denote a set of subhierarchies, that is, a hierarchy (also called a society). Let the ranks be ordered from left to right and separated by a colon; e.g., (2:3) is a subhierarchy with three elements ("individuals") on top and two elements on the bottom rank.
Then the hierarchical ordering for n = 4 is composed of the following sets: [(1:1),(2)]; [(1),(3)]; [(1),(1:1:1)]; [(1),(2:1)]; [(1),(1:2)]; [(4)]; [(2:2)]; [(1:3)]; [(3:1)]; [(1:1:2)]; [(1:2:1)]; [(2:1:1)]; [(1:1:1:1)]; thus a(4) = 13.
For example, the following hierarchy is not allowed: [(1),(1),(1),(1)] because of the repetition of (1).
		

Crossrefs

A034691 counts multisets of compositions, ordered A133494.
A261049 counts sets of partitions, ordered A358906.

Programs

  • Maple
    main := proc(n::integer) local a, ListOfPartitions, NumberOfPartitions, APartition, APart, ASet, MultipliticityOfAPart, ndxprttn, ndxprt, Term, Produkt; with(combinat): with(ListTools): a := 0; ListOfPartitions := partition(n); NumberOfPartitions := nops(ListOfPartitions); for ndxprttn from 1 to NumberOfPartitions do APartition := ListOfPartitions[ndxprttn]; ASet := convert(APartition,set); Produkt := 1; for ndxprt from 1 to nops(ASet) do APart := op(ndxprt,ASet); MultipliticityOfAPart := Occurrences(APart, APartition); Term := 2^(APart-1); Term := binomial(Term,MultipliticityOfAPart); Produkt := Produkt * Term; # End of do-loop *** ndxprt ***. end do; a := a + Produkt; # End of do-loop *** ndxprttn ***. end do; print("n, a(n):",n,a); end proc;
    PartitionList := proc (n, k) # Authors: # Herbert S. Wilf and Joanna Nordlicht, # Source: # Lecture Notes "East Side West Side,..." # University of Pennsylvania, USA, 2002. # Available from http://www.cis.upenn.edu/~wilf/lecnotes.html # Berechnet die Partitionen von n mit k Summanden. local East, West; if n < 1 or k < 1 or n < k then RETURN([]) elif n = 1 then RETURN([[1]]) else if n < 2 or k < 2 or n < k then West := [] else West := map(proc (x) options operator, arrow; [op(x), 1] end proc, PartitionList(n-1, k-1)) end if; if k <= n-k then East := map(proc(y) options operator, arrow; map(proc (x) options operator, arrow; x+1 end proc, y) end proc, PartitionList(n-k, k)) else East := [] end if; RETURN([op(West), op(East)]) end if end proc;
    # second Maple program:
    series(exp(add((-1)^(j-1)/j*z^j/(1-2*z^j), j=1..40)), z, 40); # Cf. A102866; Vladeta Jovovic, Feb 19 2008
    # alternative Maple program:
    b:= proc(n, i) option remember; `if`(n=0 or i=1, `if`(n>1, 0, 1),
          add(b(n-i*j, i-1)*binomial(2^(i-1), j), j=0..n/i))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..32);  # Alois P. Heinz, May 22 2018
  • Mathematica
    terms = 32; CoefficientList[Product[(1 + x^k)^(2^(k-1)), {k, 1, terms+1}] + O[x]^(terms+1), x] // Rest (* Jean-François Alcover, Nov 10 2017, after Vladeta Jovovic *)
    nmax = 40; CoefficientList[Series[Exp[Sum[-(-1)^k*x^k/(k*(1 - 2 x^k)), {k, 1, nmax}]], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jun 08 2018 *)

Formula

a(n) = Sum_{ partitions n = s_1 + ... + s_n } Product_{ Set{s_i} } C(2^(s_i - 1), m(s_i)), where the sum runs over all partitions of n, the product runs over the set of parts of a given partition, s_i is the i-th part in the set of parts, C(k, l) denotes the binomial coefficient and m(s_i) is the multiplicity of part s_i in the given partition.
G.f.: Product_{k>=1} (1+x^k)^(2^(k-1)). - Vladeta Jovovic, Feb 19 2008
a(n) ~ 2^n * exp(sqrt(2*n) - 1/4 + c) / (sqrt(2*Pi) * 2^(3/4) * n^(3/4)), where c = Sum_{k>=2} -(-1)^k / (k*(2^k-2)) = -0.207530918644117743551169251314627032059... - Vaclav Kotesovec, Jun 08 2018
Weigh transform of A011782. - Alois P. Heinz, Jun 25 2018

Extensions

More terms from Alois P. Heinz, Apr 21 2012
a(0)=1 prepended by Alois P. Heinz, May 22 2018

A074854 a(n) = Sum_{d|n} (2^(n-d)).

Original entry on oeis.org

1, 3, 5, 13, 17, 57, 65, 209, 321, 801, 1025, 3905, 4097, 12417, 21505, 53505, 65537, 233985, 262145, 885761, 1327105, 3147777, 4194305, 16060417, 17825793, 50339841, 84148225, 220217345, 268435457, 990937089, 1073741825, 3506503681
Offset: 1

Views

Author

Miklos Kristof, Sep 11 2002

Keywords

Comments

A034729 = Sum_{d|n} (2^(d-1)).
If p is a prime, then a(p) = A034729(p) = 2^(p-1)+1.
From Gus Wiseman, Jul 14 2020: (Start)
Number of ways to tile a rectangle of size n using horizontal strips. Also the number of ways to choose a composition of each part of a constant partition of n. The a(0) = 1 through a(5) = 17 splittings are:
() (1) (2) (3) (4) (5)
(1,1) (1,2) (1,3) (1,4)
(1),(1) (2,1) (2,2) (2,3)
(1,1,1) (3,1) (3,2)
(1),(1),(1) (1,1,2) (4,1)
(1,2,1) (1,1,3)
(2,1,1) (1,2,2)
(2),(2) (1,3,1)
(1,1,1,1) (2,1,2)
(1,1),(2) (2,2,1)
(2),(1,1) (3,1,1)
(1,1),(1,1) (1,1,1,2)
(1),(1),(1),(1) (1,1,2,1)
(1,2,1,1)
(2,1,1,1)
(1,1,1,1,1)
(1),(1),(1),(1),(1)
(End)

Examples

			Divisors of 6 = 1,2,3,6 and 6-1 = 5, 6-2 = 4, 6-3 = 3, 6-6 = 0. a(6) = 2^5 + 2^4 + 2^3 + 2^0 = 32 + 16 + 8 + 1 = 57.
G.f. = x + 3*x^2 + 5*x^3 + 13*x^4 + 17*x^5 + 57*x^6 + 65*x^7 + ...
a(14) = 1 + 2^7 + 2^12 + 2^13 = 12417. - _Gus Wiseman_, Jun 20 2018
		

Crossrefs

Cf. A080267.
Cf. A051731.
The version looking at lengths instead of sums is A101509.
The strictly increasing (or strictly decreasing) version is A304961.
Starting with a partition gives A317715.
Starting with a strict partition gives A318683.
Requiring distinct instead of equal sums gives A336127.
Starting with a strict composition gives A336130.
Partitions of partitions are A001970.
Splittings of compositions are A133494.
Splittings of partitions are A323583.

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, Sum[ 2^(n - d), {d, Divisors[n]}]] (* Michael Somos, Mar 28 2013 *)
  • PARI
    a(n)=if(n<1,0,2^n*polcoeff(sum(k=1,n,2/(2-x^k),x*O(x^n)),n))
    
  • PARI
    a(n) = sumdiv(n,d, 2^(n-d) ); /* Joerg Arndt, Mar 28 2013 */

Formula

G.f.: 2^n times coefficient of x^n in Sum_{k>=1} x^k/(2-x^k). - Benoit Cloitre, Apr 21 2003; corrected by Joerg Arndt, Mar 28 2013
G.f.: Sum_{k>0} 2^(k-1)*x^k/(1-2^(k-1)*x^k). - Vladeta Jovovic, Jun 24 2003
G.f.: Sum_{n>=1} a*z^n/(1-a*z^n) (generalized Lambert series) where z=2*x and a=1/2. - Joerg Arndt, Jan 30 2011
Triangle A051731 mod 2 converted to decimal. - Philippe Deléham, Oct 04 2003
G.f.: Sum_{k>0} 1 / (2 / (2*x)^k - 1). - Michael Somos, Mar 28 2013

Extensions

a(14) corrected from 9407 to 12417 by Gus Wiseman, Jun 20 2018

A336135 Number of ways to split an integer partition of n into contiguous subsequences with strictly decreasing sums.

Original entry on oeis.org

1, 1, 2, 5, 8, 16, 29, 50, 79, 135, 213, 337, 522, 796, 1191, 1791, 2603, 3799, 5506, 7873, 11154, 15768, 21986, 30565, 42218, 57917, 78968, 107399, 144932, 194889, 261061, 347773, 461249, 610059, 802778, 1053173, 1377325, 1793985, 2329009, 3015922, 3891142
Offset: 0

Views

Author

Gus Wiseman, Jul 11 2020

Keywords

Examples

			The a(1) = 1 through a(5) = 16 splittings:
  (1)  (2)    (3)        (4)          (5)
       (1,1)  (2,1)      (2,2)        (3,2)
              (1,1,1)    (3,1)        (4,1)
              (2),(1)    (2,1,1)      (2,2,1)
              (1,1),(1)  (3),(1)      (3,1,1)
                         (1,1,1,1)    (3),(2)
                         (2,1),(1)    (4),(1)
                         (1,1,1),(1)  (2,1,1,1)
                                      (2,2),(1)
                                      (3),(1,1)
                                      (3,1),(1)
                                      (1,1,1,1,1)
                                      (2,1),(1,1)
                                      (2,1,1),(1)
                                      (1,1,1),(1,1)
                                      (1,1,1,1),(1)
		

Crossrefs

The version with equal sums is A317715.
The version with strictly increasing sums is A336134.
The version with weakly increasing sums is A336136.
The version with weakly decreasing sums is A316245.
The version with different sums is A336131.
Starting with a composition gives A304961.
Starting with a strict partition gives A318684.
Partitions of partitions are A001970.
Partitions of compositions are A075900.
Compositions of compositions are A133494.
Compositions of partitions are A323583.

Programs

  • Mathematica
    splits[dom_]:=Append[Join@@Table[Prepend[#,Take[dom,i]]&/@splits[Drop[dom,i]],{i,Length[dom]-1}],{dom}];
    Table[Sum[Length[Select[splits[ctn],Greater@@Total/@#&]],{ctn,IntegerPartitions[n]}],{n,0,10}]
  • PARI
    a(n)={my(recurse(r,m,s,t,f)=if(m==0, r==0, if(f, self()(r,min(m,t-1),t-1,0,0)) + self()(r,m-1,s,t,0) + if(t+m<=s, self()(r-m,min(m,r-m),s,t+m,1)))); recurse(n,n,n,0)} \\ Andrew Howroyd, Jan 18 2024

Extensions

a(21) onwards from Andrew Howroyd, Jan 18 2024

A336342 Number of ways to choose a partition of each part of a strict composition of n.

Original entry on oeis.org

1, 1, 2, 7, 11, 29, 81, 155, 312, 708, 1950, 3384, 7729, 14929, 32407, 81708, 151429, 305899, 623713, 1234736, 2463743, 6208978, 10732222, 22487671, 43000345, 86573952, 160595426, 324990308, 744946690, 1336552491, 2629260284, 5050032692, 9681365777
Offset: 0

Views

Author

Gus Wiseman, Jul 18 2020

Keywords

Comments

A strict composition of n is a finite sequence of distinct positive integers summing to n.
Is there a simple generating function?

Examples

			The a(1) = 1 through a(4) = 11 ways:
  (1)  (2)    (3)        (4)
       (1,1)  (2,1)      (2,2)
              (1,1,1)    (3,1)
              (1),(2)    (1),(3)
              (2),(1)    (2,1,1)
              (1),(1,1)  (3),(1)
              (1,1),(1)  (1,1,1,1)
                         (1),(2,1)
                         (2,1),(1)
                         (1),(1,1,1)
                         (1,1,1),(1)
		

Crossrefs

Multiset partitions of partitions are A001970.
Strict compositions are counted by A032020, A072574, and A072575.
Splittings of partitions are A323583.
Splittings of partitions with distinct sums are A336131.
Partitions:
- Partitions of each part of a partition are A063834.
- Compositions of each part of a partition are A075900.
- Strict partitions of each part of a partition are A270995.
- Strict compositions of each part of a partition are A336141.
Strict partitions:
- Partitions of each part of a strict partition are A271619.
- Compositions of each part of a strict partition are A304961.
- Strict partitions of each part of a strict partition are A279785.
- Strict compositions of each part of a strict partition are A336142.
Compositions:
- Partitions of each part of a composition are A055887.
- Compositions of each part of a composition are A133494.
- Strict partitions of each part of a composition are A304969.
- Strict compositions of each part of a composition are A307068.
Strict compositions:
- Partitions of each part of a strict composition are A336342.
- Compositions of each part of a strict composition are A336127.
- Strict partitions of each part of a strict composition are A336343.
- Strict compositions of each part of a strict composition are A336139.

Programs

  • Mathematica
    Table[Length[Join@@Table[Tuples[IntegerPartitions/@ctn],{ctn,Join@@Permutations/@Select[IntegerPartitions[n],UnsameQ@@#&]}]],{n,0,10}]
  • PARI
    seq(n)={[subst(serlaplace(p),y,1) | p<-Vec(prod(k=1, n, 1 + y*x^k*numbpart(k) + O(x*x^n)))]} \\ Andrew Howroyd, Apr 16 2021

Formula

G.f.: Sum_{k>=0} k! * [y^k](Product_{j>=1} 1 + y*x^j*A000041(j)). - Andrew Howroyd, Apr 16 2021

A344062 Expansion of Product_{k>=1} (1 + 3^(k-1)*x^k).

Original entry on oeis.org

1, 1, 3, 12, 36, 135, 432, 1539, 4860, 17496, 55404, 192456, 623295, 2125764, 6849684, 23442453, 75110328, 252965916, 822670668, 2735858268, 8838926712, 29501352792, 95090206689, 314068876416, 1018141045092, 3342663979092, 10798571289897, 35481518064576
Offset: 0

Views

Author

Ilya Gutkovskiy, May 08 2021

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 27; CoefficientList[Series[Product[(1 + 3^(k - 1) x^k), {k, 1, nmax}], {x, 0, nmax}], x]
    Table[Sum[Length[Select[IntegerPartitions[n, {k}], UnsameQ @@ # &]] 3^(n - k), {k, 0, Floor[(Sqrt[8 n + 1] - 1)/2]}], {n, 0, 27}]
  • PARI
    seq(n)={Vec(prod(k=1, n, 1 + 3^(k-1)*x^k + O(x*x^n)))} \\ Andrew Howroyd, May 08 2021

Formula

a(n) = Sum_{k=0..A003056(n)} q(n,k) * 3^(n-k), where q(n,k) is the number of partitions of n into k distinct parts.
a(n) ~ (-polylog(2, -1/3))^(1/4) * 3^n * exp(2*sqrt(-polylog(2, -1/3)*n)) / (4*sqrt(Pi/3)*n^(3/4)). - Vaclav Kotesovec, May 09 2021

A336134 Number of ways to split an integer partition of n into contiguous subsequences with strictly increasing sums.

Original entry on oeis.org

1, 1, 2, 4, 6, 11, 17, 27, 37, 62, 82, 125, 168, 246, 320, 462, 585, 839, 1078, 1466, 1830, 2528, 3136, 4188, 5210, 6907, 8498, 11177, 13570, 17668, 21614, 27580, 33339, 42817, 51469, 65083, 78457, 98409, 117602, 147106, 174663, 217400, 259318, 319076, 377707
Offset: 0

Views

Author

Gus Wiseman, Jul 11 2020

Keywords

Examples

			The a(1) = 1 through a(6) = 17 splits:
  (1)  (2)    (3)        (4)          (5)            (6)
       (1,1)  (2,1)      (2,2)        (3,2)          (3,3)
              (1,1,1)    (3,1)        (4,1)          (4,2)
              (1),(1,1)  (2,1,1)      (2,2,1)        (5,1)
                         (1,1,1,1)    (3,1,1)        (2,2,2)
                         (1),(1,1,1)  (2,1,1,1)      (3,2,1)
                                      (2),(2,1)      (4,1,1)
                                      (1,1,1,1,1)    (2,2,1,1)
                                      (2),(1,1,1)    (2),(2,2)
                                      (1),(1,1,1,1)  (3,1,1,1)
                                      (1,1),(1,1,1)  (2,1,1,1,1)
                                                     (2),(2,1,1)
                                                     (1,1,1,1,1,1)
                                                     (2),(1,1,1,1)
                                                     (1),(1,1,1,1,1)
                                                     (1,1),(1,1,1,1)
                                                     (1),(1,1),(1,1,1)
		

Crossrefs

The version with equal sums is A317715.
The version with strictly decreasing sums is A336135.
The version with weakly decreasing sums is A316245.
The version with different sums is A336131.
Starting with a composition gives A304961.
Starting with a strict partition gives A336133.
Partitions of partitions are A001970.
Partitions of compositions are A075900.
Compositions of compositions are A133494.
Compositions of partitions are A323583.

Programs

  • Mathematica
    splits[dom_]:=Append[Join@@Table[Prepend[#,Take[dom,i]]&/@splits[Drop[dom,i]],{i,Length[dom]-1}],{dom}];
    Table[Sum[Length[Select[splits[ctn],Less@@Total/@#&]],{ctn,IntegerPartitions[n]}],{n,0,10}]
  • PARI
    a(n)={my(recurse(r,m,s,t,f)=if(m==0, r==0, if(f && r > t && t >= s, self()(r,m,t+1,0,0)) + self()(r,m-1,s,t,0) + self()(r-m,min(m,r-m), s,t+m,1))); recurse(n,n,0,0,0)} \\ Andrew Howroyd, Jan 18 2024

Extensions

a(21) onwards from Andrew Howroyd, Jan 18 2024

A336139 Number of ways to choose a strict composition of each part of a strict composition of n.

Original entry on oeis.org

1, 1, 1, 5, 9, 17, 45, 81, 181, 397, 965, 1729, 3673, 7313, 15401, 34065, 68617, 135069, 266701, 556969, 1061921, 2434385, 4436157, 9120869, 17811665, 35651301, 68949549, 136796317, 283612973, 537616261, 1039994921, 2081261717, 3980842425, 7723253181, 15027216049
Offset: 0

Views

Author

Gus Wiseman, Jul 16 2020

Keywords

Comments

A strict composition of n is a finite sequence of distinct positive integers summing to n.

Examples

			The a(1) = 1 through a(5) = 17 splittings:
  (1)  (2)  (3)      (4)        (5)
            (1,2)    (1,3)      (1,4)
            (2,1)    (3,1)      (2,3)
            (1),(2)  (1),(3)    (3,2)
            (2),(1)  (3),(1)    (4,1)
                     (1),(1,2)  (1),(4)
                     (1),(2,1)  (2),(3)
                     (1,2),(1)  (3),(2)
                     (2,1),(1)  (4),(1)
                                (1),(1,3)
                                (1,2),(2)
                                (1),(3,1)
                                (1,3),(1)
                                (2),(1,2)
                                (2,1),(2)
                                (2),(2,1)
                                (3,1),(1)
		

Crossrefs

The version for partitions is A063834.
Row sums of A072574.
The version for non-strict compositions is A133494.
The version for strict partitions is A279785.
Multiset partitions of partitions are A001970.
Strict compositions are A032020.
Taking a composition of each part of a partition: A075900.
Taking a composition of each part of a strict partition: A304961.
Taking a strict composition of each part of a composition: A307068.
Splittings of partitions are A323583.
Compositions of parts of strict compositions are A336127.
Set partitions of strict compositions are A336140.

Programs

  • Mathematica
    strs[n_]:=Join@@Permutations/@Select[IntegerPartitions[n],UnsameQ@@#&];
    Table[Length[Join@@Table[Tuples[strs/@ctn],{ctn,strs[n]}]],{n,0,15}]

A307068 Expansion of 1/(1 - Sum_{k>=1} k!*x^(k*(k+1)/2) / Product_{j=1..k} (1 - x^j)).

Original entry on oeis.org

1, 1, 2, 6, 14, 34, 88, 216, 532, 1322, 3290, 8142, 20192, 50080, 124144, 307878, 763474, 1893038, 4694060, 11639580, 28861736, 71567206, 177460750, 440037738, 1091134276, 2705618900, 6708953156, 16635775698, 41250705518, 102286806130, 253634237896, 628921097352, 1559496588628
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 22 2019

Keywords

Comments

Invert transform of A032020.
Number of ways to choose a strict composition of each part of a composition of n. - Gus Wiseman, Jul 18 2020
The Invert transform T(a) of a sequence a is given by T(a)n = Sum_c Product_i a(c_i), where the sum is over all compositions c of n. - _Gus Wiseman, Aug 01 2020

Examples

			From _Gus Wiseman_, Jul 18 2020: (Start)
The a(1) = 1 through a(4) = 14 ways to choose a strict composition of each part of a composition:
    (1)  (2)      (3)          (4)
         (1),(1)  (1,2)        (1,3)
                  (2,1)        (3,1)
                  (1),(2)      (1),(3)
                  (2),(1)      (2),(2)
                  (1),(1),(1)  (3),(1)
                               (1),(1,2)
                               (1),(2,1)
                               (1,2),(1)
                               (2,1),(1)
                               (1),(1),(2)
                               (1),(2),(1)
                               (2),(1),(1)
                               (1),(1),(1),(1)
(End)
		

Crossrefs

The version for partitions is A270995.
Starting with a strict composition gives A336139.
Strict compositions are counted by A032020.
Partitions of each part of a partition are A063834.
Compositions of each part of a partition are A075900.
Compositions of each part of a composition are A133494.
Strict partitions of each part of a strict partition are A279785.
Compositions of each part of a strict partition are A304961.
Strict partitions of each part of a composition are A304969.
Compositions of each part of a strict composition are A336127.
Set partitions of strict compositions are A336140.
Strict compositions of each part of a partition are A336141.

Programs

  • Magma
    m:=80;
    R:=PowerSeriesRing(Integers(), m);
    Coefficients(R!( 1/(1 - (&+[Factorial(k)*x^Binomial(k+1,2)/(&*[ 1-x^j: j in [1..k]]): k in [1..m+2]]) ) )); // G. C. Greubel, Jan 25 2024
    
  • Maple
    T:= proc(n, k) option remember; `if`(k<0 or n<0, 0,
          `if`(k=0, `if`(n=0, 1, 0), T(n-k, k) +k*T(n-k, k-1)))
        end:
    g:= proc(n) option remember; add(T(n, k), k=0..floor((sqrt(8*n+1)-1)/2)) end:
    a:= proc(n) option remember; `if`(n<1, 1,
          add(a(n-i)*g(i), i=1..n))
        end:
    seq(a(n), n=0..32);  # Alois P. Heinz, Dec 16 2022
  • Mathematica
    nmax = 32; CoefficientList[Series[1/(1 - Sum[k!*x^(k*(k+1)/2)/Product[ (1-x^j), {j,k}], {k,nmax}]), {x, 0, nmax}], x]
  • SageMath
    m=80;
    def p(x, j): return product(1-x^k for k in range(1,j+1))
    def f(x): return 1/(1 - sum(factorial(j)*x^binomial(j+1,2)/p(x,j) for j in range(1, m+3)) )
    def A307068_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( f(x) ).list()
    A307068_list(m) # G. C. Greubel, Jan 25 2024

Formula

a(0) = 1; a(n) = Sum_{k=1..n} A032020(k)*a(n-k).
Showing 1-10 of 36 results. Next