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

A066186 Sum of all parts of all partitions of n.

Original entry on oeis.org

0, 1, 4, 9, 20, 35, 66, 105, 176, 270, 420, 616, 924, 1313, 1890, 2640, 3696, 5049, 6930, 9310, 12540, 16632, 22044, 28865, 37800, 48950, 63336, 81270, 104104, 132385, 168120, 212102, 267168, 334719, 418540, 520905, 647172, 800569, 988570, 1216215, 1493520
Offset: 0

Views

Author

Wouter Meeussen, Dec 15 2001

Keywords

Comments

Sum of the zeroth moments of all partitions of n.
Also the number of one-element transitions from the integer partitions of n to the partitions of n-1 for labeled parts with the assumption that any part z is composed of labeled elements of amount 1, i.e., z = 1_1 + 1_2 + ... + 1_z. Then one can take from z a single element in z different ways. E.g., for n=3 to n=2 we have A066186(3) = 9 and [111] --> [11], [111] --> [11], [111] --> [11], [12] --> [111], [12] --> [111], [12] --> [2], [3] --> 2, [3] --> 2, [3] --> 2. For the unlabeled case, one can take a single element from z in only one way. Then the number of one-element transitions from the integer partitions of n to the partitions of n-1 is given by A000070. E.g., A000070(3) = 4 and for the transition from n=3 to n=2 one has [111] --> [11], [12] --> [11], [12] --> [2], [3] --> [2]. - Thomas Wieder, May 20 2004
Also sum of all parts of all regions of n (Cf. A206437). - Omar E. Pol, Jan 13 2013
From Omar E. Pol, Jan 19 2021: (Start)
Apart from initial zero this is also as follows:
Convolution of A000203 and A000041.
Convolution of A024916 and A002865.
For n >= 1, a(n) is also the number of cells in a symmetric polycube in which the terraces are the symmetric representation of sigma(k), for k = n..1, (cf. A237593) starting from the base and located at the levels A000041(0)..A000041(n-1) respectively. The polycube looks like a symmetric tower (cf. A221529). A dissection is a three-dimensional spiral whose top view is described in A239660. The growth of the volume of the polycube represents each convolution mentioned above. (End)
From Omar E. Pol, Feb 04 2021: (Start)
a(n) is also the sum of all divisors of all positive integers in a sequence with n blocks where the m-th block consists of A000041(n-m) copies of m, with 1 <= m <= n. The mentioned divisors are also all parts of all partitions of n.
Apart from initial zero this is also the convolution of A340793 and A000070. (End)

Examples

			a(3)=9 because the partitions of 3 are: 3, 2+1 and 1+1+1; and (3) + (2+1) + (1+1+1) = 9.
a(4)=20 because A000041(4)=5 and 4*5=20.
		

Crossrefs

Cf. A000041, A093694, A000070, A132825, A001787 (same for ordered partitions), A277029, A000203, A221529, A237593, A239660.
First differences give A138879. - Omar E. Pol, Aug 16 2013

Programs

  • Haskell
    a066186 = sum . concat . ps 1 where
       ps _ 0 = [[]]
       ps i j = [t:ts | t <- [i..j], ts <- ps t (j - t)]
    -- Reinhard Zumkeller, Jul 13 2013
    
  • Maple
    with(combinat): a:= n-> n*numbpart(n): seq(a(n), n=0..50); # Zerinvary Lajos, Apr 25 2007
  • Mathematica
    PartitionsP[ Range[0, 60] ] * Range[0, 60]
  • PARI
    a(n)=numbpart(n)*n \\ Charles R Greathouse IV, Mar 10 2012
    
  • Python
    from sympy import npartitions
    def A066186(n): return n*npartitions(n) # Chai Wah Wu, Oct 22 2023
  • Sage
    [n*Partitions(n).cardinality() for n in range(41)] # Peter Luschny, Jul 29 2014
    

Formula

a(n) = n * A000041(n). - Omar E. Pol, Oct 10 2011
G.f.: x * (d/dx) Product_{k>=1} 1/(1-x^k), i.e., derivative of g.f. for A000041. - Jon Perry, Mar 17 2004 (adjusted to match the offset by Geoffrey Critzer, Nov 29 2014)
Equals A132825 * [1, 2, 3, ...]. - Gary W. Adamson, Sep 02 2007
a(n) = A066967(n) + A066966(n). - Omar E. Pol, Mar 10 2012
a(n) = A207381(n) + A207382(n). - Omar E. Pol, Mar 13 2012
a(n) = A006128(n) + A196087(n). - Omar E. Pol, Apr 22 2012
a(n) = A220909(n)/2. - Omar E. Pol, Jan 13 2013
a(n) = Sum_{k=1..n} A000203(k)*A000041(n-k), n >= 1. - Omar E. Pol, Jan 20 2013
a(n) = Sum_{k=1..n} k*A036043(n,n-k+1). - L. Edson Jeffery, Aug 03 2013
a(n) = Sum_{k=1..n} A024916(k)*A002865(n-k), n >= 1. - Omar E. Pol, Jul 13 2014
a(n) ~ exp(Pi*sqrt(2*n/3))/(4*sqrt(3)) * (1 - (sqrt(3/2)/Pi + Pi/(24*sqrt(6))) / sqrt(n)). - Vaclav Kotesovec, Oct 24 2016
a(n) = Sum_{k=1..n} A340793(k)*A000070(n-k), n >= 1. - Omar E. Pol, Feb 04 2021

Extensions

a(0) added by Franklin T. Adams-Watters, Jul 28 2014

A066897 Total number of odd parts in all partitions of n.

Original entry on oeis.org

1, 2, 5, 8, 15, 24, 39, 58, 90, 130, 190, 268, 379, 522, 722, 974, 1317, 1754, 2330, 3058, 4010, 5200, 6731, 8642, 11068, 14076, 17864, 22528, 28347, 35490, 44320, 55100, 68355, 84450, 104111, 127898, 156779, 191574, 233625, 284070, 344745, 417292, 504151
Offset: 1

Views

Author

Naohiro Nomoto, Jan 24 2002

Keywords

Comments

Also sum of all odd-indexed parts minus the sum of all even-indexed parts of all partitions of n (Cf. A206563). - Omar E. Pol, Feb 12 2012
Column 1 of A206563. - Omar E. Pol, Feb 15 2012
Suppose that p=[p(1),p(2),p(3),...] is a partition of n with parts in nonincreasing order. Let f(p) = p(1) - p(2) + p(3) - ... be the alternating sum of parts of p and let F(n) = sum of alternating sums of all partitions of n. Conjecture: F(n) = A066897(n) for n >= 1. - Clark Kimberling, May 17 2019
From Omar E. Pol, Apr 02 2023: (Start)
Convolution of A000041 and A001227.
Convolution of A002865 and A060831.
a(n) is also the total number of odd divisors of all positive integers in a sequence with n blocks where the m-th block consists of A000041(n-m) copies of m, with 1 <= m <= n. The mentioned odd divisors are also all odd parts of all partitions of n. (End)
a(n) is odd iff n is a term of A067567 (proof: n*p(n) = the sum of the parts in all the partitions of n == the number of odd parts in all partitions of n (mod 2). Hence the number of odd parts in all partitions of n is odd iff n*p(n) is odd, equivalently, iff both n and p(n) are odd). - Peter Bala, Jan 11 2025

Examples

			a(4) = 8 because in the partitions of 4, namely [4],[3,1],[2,2],[2,1,1],[1,1,1,1], we have a total of 0+2+0+2+4=8 odd parts.
		

Crossrefs

Programs

  • Haskell
    a066897 = p 0 1 where
       p o _             0 = o
       p o k m | m < k     = 0
               | otherwise = p (o + mod k 2) k (m - k) + p o (k + 1) m
    -- Reinhard Zumkeller, Mar 09 2012
    
  • Haskell
    a066897 = length . filter odd . concat . ps 1 where
       ps _ 0 = [[]]
       ps i j = [t:ts | t <- [i..j], ts <- ps t (j - t)]
    -- Reinhard Zumkeller, Jul 13 2013
  • Maple
    g:=sum(x^(2*j-1)/(1-x^(2*j-1)),j=1..70)/product(1-x^j,j=1..70): gser:=series(g,x=0,45): seq(coeff(gser,x^n),n=1..44);
    # Emeric Deutsch, Mar 13 2006
    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]+ (i mod 2)*g[1]]
          fi
        end:
    a:= n-> b(n, n)[2]:
    seq(a(n), n=1..50);
    # Alois P. Heinz, Mar 22 2012
  • Mathematica
    f[n_, i_] := Count[Flatten[IntegerPartitions[n]], i]
    o[n_] := Sum[f[n, i], {i, 1, n, 2}]
    e[n_] := Sum[f[n, i], {i, 2, n, 2}]
    Table[o[n], {n, 1, 45}]  (* A066897 *)
    Table[e[n], {n, 1, 45}]  (* A066898 *)
    %% - %                   (* A209423 *)
    (* Clark Kimberling, Mar 08 2012 *)
    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]] + Mod[i, 2]*g[[1]]}] ]; a[n_] := b[n, n][[2]]; Table[a[n], {n, 1, 50}] (* Jean-François Alcover, Sep 26 2015, after Alois P. Heinz *)

Formula

a(n) = Sum_{k=1..n} b(k)*numbpart(n-k), where b(k)=A001227(k)=number of odd divisors of k and numbpart() is A000041. - Vladeta Jovovic, Jan 26 2002
a(n) = Sum_{k=0..n} k*A103919(n,k). - Emeric Deutsch, Mar 13 2006
G.f.: Sum_{j>=1}(x^(2j-1)/(1-x^(2j-1)))/Product_{j>=1}(1-x^j). - Emeric Deutsch, Mar 13 2006
a(n) = A066898(n) + A209423(n) = A006128(n) - A066898(n). [Reinhard Zumkeller, Mar 09 2012]
a(n) = A207381(n) - A207382(n). - Omar E. Pol, Mar 11 2012
a(n) = (A006128(n) + A209423(n))/2. - Vaclav Kotesovec, May 25 2018
a(n) ~ exp(Pi*sqrt(2*n/3)) * (2*gamma + log(24*n/Pi^2)) / (8*Pi*sqrt(2*n)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, May 25 2018

Extensions

More terms from Vladeta Jovovic, Jan 26 2002

A207381 Total sum of the odd-indexed parts of all partitions of n.

Original entry on oeis.org

1, 3, 7, 14, 25, 45, 72, 117, 180, 275, 403, 596, 846, 1206, 1681, 2335, 3183, 4342, 5820, 7799, 10321, 13622, 17798, 23221, 30009, 38706, 49567, 63316, 80366, 101805, 128211, 161134, 201537, 251495, 312508, 387535, 478674, 590072, 724920, 888795, 1086324
Offset: 1

Views

Author

Omar E. Pol, Feb 17 2012

Keywords

Comments

For more information see A206563.

Examples

			For n = 5, write the partitions of 5 and below write the sums of their odd-indexed parts:
.    5
.    3+2
.    4+1
.    2+2+1
.    3+1+1
.    2+1+1+1
.    1+1+1+1+1
.  ------------
.   20 + 4 + 1 = 25
The total sum of the odd-indexed parts is 25 so a(5) = 25.
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; local g, h;
          if n=0 then [1, 0$2]
        elif i<1 then [0$3]
        else g:= b(n, i-1); h:= `if`(i>n, [0$3], b(n-i, i));
             [g[1]+h[1], g[2]+h[3], g[3]+h[2]+i*h[1]]
          fi
        end:
    a:= n-> b(n,n)[3]:
    seq(a(n), n=1..50); # Alois P. Heinz, Mar 12 2012
  • Mathematica
    b[n_, i_] := b[n, i] = Module[{g, h}, If[n == 0 , {1, 0, 0}, If[i < 1, {0, 0, 0},  g = b[n, i - 1]; h = If[i > n, {0, 0, 0}, b[n - i, i]]; {g[[1]] + h[[1]], g[[2]] + h[[3]], g[[3]] + h[[2]] + i*h[[1]]}]]]; a[n_] := b[n, n][[3]]; Table [a[n], {n, 1, 50}] (* Jean-François Alcover, Dec 09 2016 after Alois P. Heinz *)

Formula

a(n) = A066186(n) - A207382(n) = A066897(n) + A207382(n).

Extensions

More terms from Alois P. Heinz, Mar 12 2012
Showing 1-3 of 3 results.