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

A074141 Sum of products of parts increased by 1 in all partitions of n.

Original entry on oeis.org

1, 2, 7, 18, 50, 118, 301, 684, 1621, 3620, 8193, 17846, 39359, 84198, 181313, 383208, 811546, 1695062, 3546634, 7341288, 15207022, 31261006, 64255264, 131317012, 268336125, 545858260, 1110092387, 2250057282, 4558875555, 9213251118, 18613373708, 37529713890
Offset: 0

Views

Author

Amarnath Murthy, Aug 28 2002

Keywords

Comments

Replace each term in A036035 by the number of its divisors as in A074139; sequence gives sum of terms in the n-th row.
This is the sum of the number of submultisets of the multisets with n elements; a part of a partition is a frequency of such an element. - George Beck, Nov 01 2011

Examples

			The partitions of 4 are 4, 3+1, 2+2, 2+1+1, 1+1+1+1, the corresponding products when parts are increased by 1 are 5,8,9,12,16 and their sum is a(4) = 50.
		

Crossrefs

Row sums of A074139 and of A079025 and of A079308 and of A238963.
Column k=2 of A261718.
Cf. A267008.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i=1,
          2^n, b(n, i-1) +(1+i)*b(n-i, min(n-i, i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..50); # Alois P. Heinz, Sep 07 2014
  • Mathematica
    Table[Plus @@ Times @@@ (IntegerPartitions[n] + 1), {n, 0, 28}] (* T. D. Noe, Nov 01 2011 *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, (1+i) * b[n-i, i]]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Oct 08 2015, after Alois P. Heinz *)
  • Maxima
    S(n,m):=if n=0 then 1 else if nVladimir Kruchinin, Sep 07 2014 */

Formula

G.f.: 1/Product_{m>0} (1-(m+1)*x^m).
a(n) = 1/n*Sum_{k=1..n} b(k)*a(n-k), where b(k) = Sum_{d divides k} d*(d+1)^(k/d).
a(n) = S(n,1), where S(n,m) = sum(k=m..n/2, (k+1)*S(n-k,k))+(n+1), S(n,n)=n+1, S(0,m)=1, S(n,m)=0 for nVladimir Kruchinin, Sep 07 2014
a(n) ~ c * 2^n, where c = Product_{k>=2} 1/(1-(k+1)/2^k) = 18.56314656361011472747535423226928404842588594722907068201... = A256155. - Vaclav Kotesovec, Sep 11 2014, updated May 10 2021

Extensions

More terms from Alford Arnold, Sep 17 2002
More terms, better description and formulas from Vladeta Jovovic, Vladimir Baltic, Nov 28 2002

A293366 Number of partitions of n where each part i is marked with a word of length i over a binary alphabet whose letters appear in alphabetical order and both letters occur at least once in the partition.

Original entry on oeis.org

3, 12, 40, 104, 279, 654, 1577, 3560, 8109, 17734, 39205, 83996, 181043, 382856, 811084, 1694468, 3545864, 7340308, 15205768, 31259422, 64253260, 131314502, 268332975, 545854344, 1110087515, 2250051262, 4558868119, 9213241988, 18613362500, 37529700206
Offset: 2

Views

Author

Alois P. Heinz, Oct 07 2017

Keywords

Crossrefs

Column k=2 of A261719.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1, k)+`if`(i>n, 0, b(n-i, i, k)*binomial(i+k-1, k-1))))
        end:
    a:= n-> (k-> add(b(n$2, k-i)*(-1)^i*binomial(k, i), i=0..k))(2):
    seq(a(n), n=2..35);
  • Mathematica
    b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i < 1, 0, b[n, i - 1, k] + If[i > n, 0, b[n - i, i, k] Binomial[i + k - 1, k - 1]]]];
    a[n_] := With[{k = 2}, Sum[b[n, n, k - i] (-1)^i Binomial[k, i], {i, 0, k}]];
    a /@ Range[2, 35] (* Jean-François Alcover, Dec 08 2020, after Alois P. Heinz *)

Formula

a(n) ~ c * 2^n, where c = A256155 = 18.563146563610114727475354232269284... - Vaclav Kotesovec, Oct 11 2017
Showing 1-2 of 2 results.