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 21-30 of 239 results. Next

A305551 Number of partitions of partitions of n where all partitions have the same sum.

Original entry on oeis.org

1, 1, 3, 4, 9, 8, 22, 16, 43, 41, 77, 57, 201, 102, 264, 282, 564, 298, 1175, 491, 1878, 1509, 2611, 1256, 7872, 2421, 7602, 8026, 16304, 4566, 38434, 6843, 48308, 41078, 56582, 28228, 221115, 21638, 146331, 208142, 453017, 44584, 844773, 63262, 1034193, 1296708
Offset: 0

Views

Author

Gus Wiseman, Jun 20 2018

Keywords

Examples

			The a(4) = 9 partitions of partitions where all partitions have the same sum:
(4), (31), (22), (211), (1111),
(2)(2), (2)(11), (11)(11),
(1)(1)(1)(1).
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[PartitionsP[n/k]+k-1,k],{k,Divisors[n]}],{n,60}]
  • PARI
    a(n)={if(n<1, n==0, sumdiv(n, d, binomial(numbpart(n/d) + d - 1, d)))} \\ Andrew Howroyd, Jun 22 2018

Formula

a(n) = Sum_{d|n} binomial(A000041(n/d) + d - 1, d).

A302478 Products of prime numbers of squarefree index.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 20, 22, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 36, 39, 40, 41, 43, 44, 45, 47, 48, 50, 51, 52, 54, 55, 58, 59, 60, 62, 64, 65, 66, 67, 68, 72, 73, 75, 78, 79, 80, 81, 82, 83, 85, 86, 87, 88, 90, 93, 94
Offset: 1

Views

Author

Gus Wiseman, Apr 08 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.

Examples

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of set multisystems.
01:  {}
02:  {{}}
03:  {{1}}
04:  {{},{}}
05:  {{2}}
06:  {{},{1}}
08:  {{},{},{}}
09:  {{1},{1}}
10:  {{},{2}}
11:  {{3}}
12:  {{},{},{1}}
13:  {{1,2}}
15:  {{1},{2}}
16:  {{},{},{},{}}
17:  {{4}}
18:  {{},{1},{1}}
20:  {{},{},{2}}
22:  {{},{3}}
24:  {{},{},{},{1}}
25:  {{2},{2}}
26:  {{},{1,2}}
27:  {{1},{1},{1}}
29:  {{1,3}}
30:  {{},{1},{2}}
31:  {{5}}
32:  {{},{},{},{},{}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,And@@SquareFreeQ/@PrimePi/@FactorInteger[#][[All,1]]]&]
  • PARI
    ok(n)={!#select(p->!issquarefree(primepi(p)), factor(n)[,1])} \\ Andrew Howroyd, Aug 26 2018

A279790 Number of twice-partitions of type (Q,P,Q) and weight n.

Original entry on oeis.org

1, 1, 3, 3, 5, 11, 12, 18, 24, 49, 53, 82, 102, 149, 236, 297, 392, 540, 702, 924, 1423, 1723, 2318, 3016, 3969, 5037, 6647, 9285, 11448, 15048, 19143, 24695, 31288, 40075, 50014, 68583, 83056, 107252, 133796, 171565
Offset: 1

Views

Author

Gus Wiseman, Dec 18 2016

Keywords

Comments

Also number of ways to choose a sum-preserving permutation of a set partition of a strict partition of n.

Examples

			The a(6)=11 twice-partitions are:
((6)),
((5)(1)), ((51)),
((4)(2)), ((42)),
((3)(2)(1)), ((3)(21)), ((21)(3)), ((31)(2)), ((32)(1)), ((321)).
		

Crossrefs

Programs

  • Mathematica
    nn=20;
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Total[Sum[Times@@Factorial/@Length/@Split[Sort[Total/@ptn]],{ptn,sps[Sort[#]]}]&/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,nn}]

A299202 Moebius function of the multiorder of integer partitions indexed by their Heinz numbers.

Original entry on oeis.org

0, 1, 1, -1, 1, -1, 1, 0, -1, -1, 1, 2, 1, -1, -1, -1, 1, 1, 1, 1, -1, -1, 1, -1, -1, -1, 0, 1, 1, 3, 1, 0, -1, -1, -1, -1, 1, -1, -1, -1, 1, 2, 1, 1, 1, -1, 1, 0, -1, 1, -1, 1, 1, -1, -1, -1, -1, -1, 1, -3, 1, -1, 2, 0, -1, 2, 1, 1, -1, 3, 1, 2, 1, -1, 1, 1, -1, 2, 1, 1, -1, -1, 1, -5, -1, -1, -1, -1, 1, -4
Offset: 1

Views

Author

Gus Wiseman, Feb 05 2018

Keywords

Comments

By convention, mu() = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Heinz number of (2,1,1) is 12, so mu(2,1,1) = a(12) = 2.
		

Crossrefs

Programs

  • Mathematica
    nn=120;
    ptns=Table[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]],{n,nn}];
    tris=Join@@Map[Tuples[IntegerPartitions/@#]&,ptns];
    mu[y_]:=mu[y]=If[Length[y]===1,1,-Sum[Times@@mu/@t,{t,Select[tris,And[Length[#]>1,Sort[Join@@#,Greater]===y]&]}]];
    mu/@ptns

Formula

mu(y) = Sum_{g(t)=y} (-1)^d(t), where the sum is over all enriched p-trees (A289501, A299203) whose multiset of leaves is the integer partition y, and d(t) is the number of non-leaf nodes in t.

A300063 Heinz numbers of integer partitions of odd numbers.

Original entry on oeis.org

2, 5, 6, 8, 11, 14, 15, 17, 18, 20, 23, 24, 26, 31, 32, 33, 35, 38, 41, 42, 44, 45, 47, 50, 51, 54, 56, 58, 59, 60, 65, 67, 68, 69, 72, 73, 74, 77, 78, 80, 83, 86, 92, 93, 95, 96, 97, 98, 99, 103, 104, 105, 106, 109, 110, 114, 119, 122, 123, 124, 125, 126, 127
Offset: 1

Views

Author

Gus Wiseman, Feb 23 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			15 is the Heinz number of (3,2), which has odd weight, so 15 belongs to the sequence.
Sequence of odd-weight partitions begins: (1) (3) (2,1) (1,1,1) (5) (4,1) (3,2) (7) (2,2,1) (3,1,1) (9) (2,1,1,1) (6,1).
		

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; local k; for k from 1+
         `if`(n=1, 0, a(n-1)) while add(numtheory[pi]
          (i[1])*i[2], i=ifactors(k)[2])::even do od; k
        end:
    seq(a(n), n=1..100);  # Alois P. Heinz, May 22 2018
  • Mathematica
    Select[Range[200],OddQ[Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]]&]

A296188 Number of normal semistandard Young tableaux whose shape is the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 2, 1, 4, 4, 8, 1, 6, 12, 16, 6, 32, 32, 28, 1, 64, 16, 128, 24, 96, 80, 256, 8, 44, 192, 22, 80, 512, 96, 1024, 1, 288, 448, 224, 30, 2048, 1024, 800, 40, 4096, 400, 8192, 240, 168, 2304, 16384, 10, 360, 204, 2112, 672, 32768, 68, 832, 160, 5376, 5120
Offset: 1

Views

Author

Gus Wiseman, Feb 14 2018

Keywords

Comments

A tableau is normal if its entries span an initial interval of positive integers. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The a(9) = 6 tableaux:
1 3   1 2   1 2   1 2   1 1   1 1
2 4   3 4   3 3   2 3   2 3   2 2
		

References

  • Richard P. Stanley, Enumerative Combinatorics Volume 2, Cambridge University Press, 1999, Chapter 7.10.

Crossrefs

Programs

  • Mathematica
    conj[y_List]:=If[Length[y]===0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    conj[n_Integer]:=Times@@Prime/@conj[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]];
    ssyt[n_]:=If[n===1,1,Sum[ssyt[n/q*Times@@Cases[FactorInteger[q],{p_,k_}:>If[p===2,1,NextPrime[p,-1]^k]]],{q,Rest[Divisors[n]]}]];
    Table[ssyt[conj[n]],{n,50}]

Formula

Let b(n) = Sum_{d|n, d>1} b(n * d' / d) where if d = Product_i prime(s_i)^m(i) then d' = Product_i prime(s_i - 1)^m(i) and prime(0) = 1. Then a(n) = b(conj(n)) where conj = A122111.

A304969 Expansion of 1/(1 - Sum_{k>=1} q(k)*x^k), where q(k) = number of partitions of k into distinct parts (A000009).

Original entry on oeis.org

1, 1, 2, 5, 11, 25, 57, 129, 292, 662, 1500, 3398, 7699, 17443, 39519, 89536, 202855, 459593, 1041267, 2359122, 5344889, 12109524, 27435660, 62158961, 140828999, 319065932, 722884274, 1637785870, 3710611298, 8406859805, 19046805534, 43152950024, 97768473163
Offset: 0

Views

Author

Ilya Gutkovskiy, May 22 2018

Keywords

Comments

Invert transform of A000009.
From Gus Wiseman, Jul 31 2022: (Start)
Also the number of ways to choose a multiset partition into distinct constant multisets of a multiset of length n that covers an initial interval of positive integers. This interpretation involves only multisets, not sequences. For example, the a(1) = 1 through a(4) = 11 multiset partitions are:
{{1}} {{1,1}} {{1,1,1}} {{1,1,1,1}}
{{1},{2}} {{1},{1,1}} {{1},{1,1,1}}
{{1},{2,2}} {{1,1},{2,2}}
{{2},{1,1}} {{1},{2,2,2}}
{{1},{2},{3}} {{2},{1,1,1}}
{{1},{2},{1,1}}
{{1},{2},{2,2}}
{{1},{2},{3,3}}
{{1},{3},{2,2}}
{{2},{3},{1,1}}
{{1},{2},{3},{4}}
The non-strict version is A055887.
The strongly normal non-strict version is A063834.
The strongly normal version is A270995.
(End)

Examples

			From _Gus Wiseman_, Jul 31 2022: (Start)
a(n) is the number of ways to choose a strict integer partition of each part of an integer composition of n. The a(1) = 1 through a(4) = 11 choices are:
  ((1))  ((2))     ((3))        ((4))
         ((1)(1))  ((21))       ((31))
                   ((1)(2))     ((1)(3))
                   ((2)(1))     ((2)(2))
                   ((1)(1)(1))  ((3)(1))
                                ((1)(21))
                                ((21)(1))
                                ((1)(1)(2))
                                ((1)(2)(1))
                                ((2)(1)(1))
                                ((1)(1)(1)(1))
(End)
		

Crossrefs

Row sums of A308680.
The unordered version is A089259, non-strict A001970 (row-sums of A061260).
For partitions instead of compositions we have A270995, non-strict A063834.
A000041 counts integer partitions, strict A000009.
A072233 counts partitions by sum and length.
Cf. A279784.

Programs

  • Maple
    b:= proc(n) option remember; `if`(n=0, 1, add(b(n-j)*add(
         `if`(d::odd, d, 0), d=numtheory[divisors](j)), j=1..n)/n)
        end:
    a:= proc(n) option remember; `if`(n=0, 1,
          add(b(j)*a(n-j), j=1..n))
        end:
    seq(a(n), n=0..40);  # Alois P. Heinz, May 22 2018
  • Mathematica
    nmax = 32; CoefficientList[Series[1/(1 - Sum[PartitionsQ[k] x^k, {k, 1, nmax}]), {x, 0, nmax}], x]
    nmax = 32; CoefficientList[Series[1/(2 - Product[1 + x^k, {k, 1, nmax}]), {x, 0, nmax}], x]
    nmax = 32; CoefficientList[Series[1/(2 - 1/QPochhammer[x, x^2]), {x, 0, nmax}], x]
    a[0] = 1; a[n_] := a[n] = Sum[PartitionsQ[k] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 32}]

Formula

G.f.: 1/(1 - Sum_{k>=1} A000009(k)*x^k).
G.f.: 1/(2 - Product_{k>=1} (1 + x^k)).
G.f.: 1/(2 - Product_{k>=1} 1/(1 - x^(2*k-1))).
G.f.: 1/(2 - exp(Sum_{k>=1} (-1)^(k+1)*x^k/(k*(1 - x^k)))).
a(n) ~ c / r^n, where r = 0.441378990861652015438479635503868737167721352874... is the root of the equation QPochhammer[-1, r] = 4 and c = 0.4208931614610039677452560636348863586180784719323982664940444607322... - Vaclav Kotesovec, May 23 2018

A066739 Number of representations of n as a sum of products of positive integers. 1 is not allowed as a factor, unless it is the only factor. Representations which differ only in the order of terms or factors are considered equivalent.

Original entry on oeis.org

1, 1, 2, 3, 6, 8, 14, 19, 32, 44, 67, 91, 139, 186, 269, 362, 518, 687, 960, 1267, 1747, 2294, 3106, 4052, 5449, 7063, 9365, 12092, 15914, 20422, 26639, 34029, 44091, 56076, 72110, 91306, 116808, 147272, 187224, 235201, 297594, 372390, 468844, 584644, 732942
Offset: 0

Views

Author

Naohiro Nomoto, Jan 16 2002

Keywords

Examples

			For n=5, 5 = 4+1 = 2*2+1 = 3+2 = 3+1+1 = 2+2+1 = 2+1+1+1 = 1+1+1+1+1, so a(5) = 8.
For n=8, 8 = 4*2 = 2*2*2 = ... = 4+4 = 2*2+4 = 2*2+2*2 = ...; note that there are 3 ways to factor the terms of 4+4. In general, if a partition contains a number k exactly r times, then the number of ways to factor the k's is the binomial coefficient C(A001055(k)+r-1,r).
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n, k) option remember;
          `if`(n>k, 0, 1) +`if`(isprime(n), 0,
          add(`if`(d>k, 0, b(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    a:= proc(n) option remember;
          `if`(n=0, 1, add(add(d*b(d, d), d=divisors(j)) *a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..60); # Alois P. Heinz, Apr 22 2012
  • Mathematica
    p[ n_, 1 ] := If[ n==1, 1, 0 ]; p[ 1, k_ ] := 1; p[ n_, k_ ] := p[ n, k ]=p[ n, k-1 ]+If[ Mod[ n, k ]==0, p[ n/k, k ], 0 ]; A001055[ n_ ] := p[ n, n ]; a[ n_, 1 ] := 1; a[ 0, k_ ] := 1; a[ n_, k_ ] := If[ k>n, a[ n, n ], a[ n, k ]=a[ n, k-1 ]+Sum[ Binomial[ A001055[ k ]+r-1, r ]a[ n-k*r, k-1 ], {r, 1, Floor[ n/k ]} ] ]; a[ n_ ] := a[ n, n ]; (* p[ n, k ]=number of factorizations of n with factors <= k. a[ n, k ]=number of representations of n as a sum of products of positive integers, with summands <= k *)
    b[n_, k_] := b[n, k] = If[n>k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d>k, 0, b[n/d, d]], {d, Divisors[n] ~Complement~ {1, n}}]]; a[0] = 1; a[n_] := a[n] = If[n == 0, 1, Sum[DivisorSum[j, #*b[#, #]&]*a[n-j], {j, 1, n}]/n]; Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Nov 10 2015, after Alois P. Heinz *)
    facs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[facs[n/d],Min@@#1>=d&],{d,Rest[Divisors[n]]}]];
    Table[Length[Union[Sort/@Join@@Table[Tuples[facs/@ptn],{ptn,IntegerPartitions[n]}]]],{n,50}] (* Gus Wiseman, Sep 05 2018 *)
  • Python
    from sympy.core.cache import cacheit
    from sympy import divisors, isprime
    @cacheit
    def b(n, k): return (0 if n>k else 1) + (0 if isprime(n) else sum([0 if d>k else b(n//d, d) for d in divisors(n)[1:-1]]))
    @cacheit
    def a(n): return 1 if n==0 else sum(sum(d*b(d, d) for d in divisors(j))*a(n - j)  for j in range(1, n + 1))//n
    print([a(n) for n in range(61)]) # Indranil Ghosh, Aug 19 2017, after Maple code

Formula

a(n) = Sum_{pi} Product_{m=1..n} binomial(k(m)+A001055(m)-1, k(m)), where pi runs through all partitions k(1) + 2 * k( 2) + ... + n * k(n) = n. a(n)=1/n*Sum_{m=1..n} a(n-m)*b(m), n > 0, a(0)=1, b(m)=Sum_{d|m} d*A001055(d). Euler transform of A001055(n): Product_{m=1..infinity} (1-x^m)^(-A001055(m)). - Vladeta Jovovic, Jan 21 2002

Extensions

Edited by Dean Hickerson, Jan 19 2002

A279375 Number of set partitions of strict integer partitions of n that have distinct block-sums.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 12, 16, 24, 39, 49, 70, 94, 127, 202, 247, 340, 450, 606, 772, 1169, 1407, 1920, 2454, 3267, 4089, 5469, 7293, 9222, 11884, 15291, 19417, 24890, 31469, 39662, 52619, 64764, 82502, 103576, 131169, 162726, 206015, 254233, 318464, 406262, 499210, 620593, 773673, 957073, 1181593
Offset: 0

Views

Author

Gus Wiseman, Dec 11 2016

Keywords

Comments

Also twice partitioned numbers where all partitions are strict. Also triangles of weight n in the multisystem of strict partitions. Strict partitions are an example of a multisystem that is neither transitive nor partitive nor contractible but is decomposable; see link for details.

Examples

			The a(6)=9 set partitions of strict integer partitions of 6 are: ((6)), ((51)), ((5)(1)), ((42)), ((4)(2)), ((321)), ((32)(1)), ((31)(2)), ((3)(2)(1)). The set partition ((3)(21)) is not counted because its blocks do not have distinct sums.
		

Crossrefs

Programs

  • Mathematica
    nn=20;sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Total[Length[Select[sps[Sort[#]],UnsameQ@@Total/@#&]]&/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,nn}]

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
Previous Showing 21-30 of 239 results. Next