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 81-90 of 579 results. Next

A087897 Number of partitions of n into odd parts greater than 1.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 5, 5, 6, 8, 8, 10, 12, 13, 15, 18, 20, 23, 27, 30, 34, 40, 44, 50, 58, 64, 73, 83, 92, 104, 118, 131, 147, 166, 184, 206, 232, 256, 286, 320, 354, 394, 439, 485, 538, 598, 660, 730, 809, 891, 984, 1088, 1196, 1318, 1454, 1596, 1756
Offset: 0

Views

Author

N. J. A. Sloane, Dec 04 2003

Keywords

Comments

Also number of partitions of n into distinct parts which are not powers of 2.
Also number of partitions of n into distinct parts such that the two largest parts differ by 1.
Also number of partitions of n such that the largest part occurs an odd number of times that is at least 3 and every other part occurs an even number of times. Example: a(10) = 2 because we have [2,2,2,1,1,1,1] and [2,2,2,2,2]. - Emeric Deutsch, Mar 30 2006
Also difference between number of partitions of 1+n into distinct parts and number of partitions of n into distinct parts. - Philippe LALLOUET, May 08 2007
In the Berndt reference replace {a -> -x, q -> x} in equation (3.1) to get f(x). G.f. is 1 - x * (1 - f(x)).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Also number of symmetric unimodal compositions of n+3 where the maximal part appears three times. - Joerg Arndt, Jun 11 2013
Let c(n) = number of palindromic partitions of n whose greatest part has multiplicity 3; then c(n) = a(n-3) for n>=3. - Clark Kimberling, Mar 05 2014
From Gus Wiseman, Aug 22 2021: (Start)
Also the number of integer partitions of n - 1 whose parts cover an interval of positive integers starting with 2. These partitions are ranked by A339886. For example, the a(6) = 1 through a(16) = 5 partitions are:
32 222 322 332 432 3322 3332 4332 4432 5432 43332
2222 3222 22222 4322 33222 33322 33332 44322
32222 222222 43222 43322 333222
322222 332222 432222
2222222 3222222
(End)

Examples

			1 + x^3 + x^5 + x^6 + x^7 + x^8 + 2*x^9 + 2*x^10 + 2*x^11 + 3*x^12 + 3*x^13 + ...
q + q^73 + q^121 + q^145 + q^169 + q^193 + 2*q^217 + 2*q^241 + 2*q^265 + ...
a(10)=2 because we have [7,3] and [5,5].
From _Joerg Arndt_, Jun 11 2013: (Start)
There are a(22)=13 symmetric unimodal compositions of 22+3=25 where the maximal part appears three times:
01:  [ 1 1 1 1 1 1 1 1 3 3 3 1 1 1 1 1 1 1 1 ]
02:  [ 1 1 1 1 1 1 2 3 3 3 2 1 1 1 1 1 1 ]
03:  [ 1 1 1 1 1 5 5 5 1 1 1 1 1 ]
04:  [ 1 1 1 1 2 2 3 3 3 2 2 1 1 1 1 ]
05:  [ 1 1 1 2 5 5 5 2 1 1 1 ]
06:  [ 1 1 2 2 2 3 3 3 2 2 2 1 1 ]
07:  [ 1 1 3 5 5 5 3 1 1 ]
08:  [ 1 1 7 7 7 1 1 ]
09:  [ 1 2 2 5 5 5 2 2 1 ]
10:  [ 1 4 5 5 5 4 1 ]
11:  [ 2 2 2 2 3 3 3 2 2 2 2 ]
12:  [ 2 3 5 5 5 3 2 ]
13:  [ 2 7 7 7 2 ]
(End)
From _Gus Wiseman_, Feb 16 2021: (Start)
The a(7) = 1 through a(19) = 8 partitions are the following (A..J = 10..19). The Heinz numbers of these partitions are given by A341449.
  7  53  9    55  B    75    D    77    F      97    H      99      J
         333  73  533  93    553  95    555    B5    755    B7      775
                       3333  733  B3    753    D3    773    D5      955
                                  5333  933    5533  953    F3      973
                                        33333  7333  B33    5553    B53
                                                     53333  7533    D33
                                                            9333    55333
                                                            333333  73333
(End)
		

References

  • J. W. L. Glaisher, Identities, Messenger of Mathematics, 5 (1876), pp. 111-112. see Eq. I

Crossrefs

The ordered version is A000931.
Partitions with no ones are counted by A002865, ranked by A005408.
The even version is A035363, ranked by A066207.
The version for factorizations is A340101.
Partitions whose only even part is the smallest are counted by A341447.
The Heinz numbers of these partitions are given by A341449.
A000009 counts partitions into odd parts, ranked by A066208.
A025147 counts strict partitions with no 1's.
A025148 counts strict partitions with no 1's or 2's.
A026804 counts partitions whose smallest part is odd, ranked by A340932.
A027187 counts partitions with even length/maximum, ranks A028260/A244990.
A027193 counts partitions with odd length/maximum, ranks A026424/A244991.
A058695 counts partitions of odd numbers, ranked by A300063.
A058696 counts partitions of even numbers, ranked by A300061.
A340385 counts partitions with odd length and maximum, ranked by A340386.

Programs

  • Haskell
    a087897 = p [3,5..] where
       p [] _ = 0
       p _  0 = 1
       p ks'@(k:ks) m | m < k     = 0
                      | otherwise = p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Aug 12 2011
    
  • Maple
    To get 128 terms: t4 := mul((1+x^(2^n)),n=0..7); t5 := mul((1+x^k),k=1..128): t6 := series(t5/t4,x,100); t7 := seriestolist(t6);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, 1,
          `if`(i<3, 0, b(n, i-2)+`if`(i>n, 0, b(n-i, i))))
        end:
    a:= n-> b(n, n-1+irem(n, 2)):
    seq(a(n), n=0..80);  # Alois P. Heinz, Jun 11 2013
  • Mathematica
    max = 65; f[x_] := Product[ 1/(1 - x^(2k+1)), {k, 1, max}]; CoefficientList[ Series[f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 16 2011, after Emeric Deutsch *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<3, 0, b[n, i-2]+If[i>n, 0, b[n-i, i]]] ]; a[n_] := b[n, n-1+Mod[n, 2]]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Apr 01 2015, after Alois P. Heinz *)
    Flatten[{1, Table[PartitionsQ[n+1] - PartitionsQ[n], {n, 0, 80}]}] (* Vaclav Kotesovec, Dec 01 2015 *)
    Table[Length[Select[IntegerPartitions[n],FreeQ[#,1]&&OddQ[Times@@#]&]],{n,0,30}] (* Gus Wiseman, Feb 16 2021 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - x) * eta(x^2 + A) / eta(x + A), n))} /* Michael Somos, Nov 13 2011 */
    
  • Python
    from functools import lru_cache
    @lru_cache(maxsize=None)
    def A087897_T(n,k):
        if n==0: return 1
        if k<3 or n<0: return 0
        return A087897_T(n,k-2)+A087897_T(n-k,k)
    def A087897(n): return A087897_T(n,n-(n&1^1)) # Chai Wah Wu, Sep 23 2023, after Alois P. Heinz

Formula

Expansion of q^(-1/24) * (1 - q) * eta(q^2) / eta(q) in powers of q.
Expansion of (1 - x) / chi(-x) in powers of x where chi() is a Ramanujan theta function.
G.f.: 1 + x^3 + x^5*(1 + x) + x^7*(1 + x)*(1 + x^2) + x^9*(1 + x)*(1 + x^2)*(1 + x^3) + ... [Glaisher 1876]. - Michael Somos, Jun 20 2012
G.f.: Product_{k >= 1} 1/(1-x^(2*k+1)).
G.f.: Product_{k >= 1, k not a power of 2} (1+x^k).
G.f.: Sum_{k >= 1} x^(3*k)/Product_{j = 1..k} (1 - x^(2*j)). - Emeric Deutsch, Mar 30 2006
a(n) ~ exp(Pi*sqrt(n/3)) * Pi / (8 * 3^(3/4) * n^(5/4)) * (1 - (15*sqrt(3)/(8*Pi) + 11*Pi/(48*sqrt(3)))/sqrt(n) + (169*Pi^2/13824 + 385/384 + 315/(128*Pi^2))/n). - Vaclav Kotesovec, Aug 30 2015, extended Nov 04 2016
G.f.: 1/(1 - x^3) * Sum_{n >= 0} x^(5*n)/Product_{k = 1..n} (1 - x^(2*k)) = 1/((1 - x^3)*(1 - x^5)) * Sum_{n >= 0} x^(7*n)/Product_{k = 1..n} (1 - x^(2*k)) = ..., extending Deutsch's result dated Mar 30 2006. - Peter Bala, Jan 15 2021
G.f.: Sum_{n >= 0} x^(n*(2*n+1))/Product_{k = 2..2*n+1} (1 - x^k). (Set z = x^3 and q = x^2 in Mc Laughlin et al., Section 1.3, Entry 7.) - Peter Bala, Feb 02 2021
a(2*n+1) = Sum{j>=1} A008284(n+1-j,2*j - 1) and a(2*n) = Sum{j>=1} A008284(n-j, 2*j). - Gregory L. Simay, Sep 22 2023

A359889 Numbers that are 1 or whose prime indices have the same mean as median.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 64, 65, 67, 69, 71, 73, 74, 77, 79, 81, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94
Offset: 1

Views

Author

Gus Wiseman, Jan 22 2023

Keywords

Comments

First differs from A236510 in having 252 (prime indices {1,1,2,2,4}).
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).

Examples

			The prime indices of 900 are {1,1,2,2,3,3}, with mean 2 and median 2, so 900 is in the sequence.
		

Crossrefs

These partitions are counted by A240219, strict A359897.
The LHS (mean of prime indices) is A326567/A326568.
The complement is A359890, counted by A359894.
The odd-length case is A359891, complement A359892, counted by A359895.
The RHS (median of prime indices) is A360005/2.
A058398 counts partitions by mean, see also A008284, A327482.
A088529/A088530 gives mean of prime signature A124010.
A112798 lists prime indices, length A001222, sum A056239.
A316413 lists numbers whose prime indices have integer mean.
A359893 and A359901 count partitions by median, odd-length A359902.
A359908 lists numbers whose prime indices have integer median.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],#==1||Mean[prix[#]]==Median[prix[#]]&]

Formula

Numbers n such that A326567(n)/A326568(n) = A360005(n)/2.

A362615 Irregular triangle read by rows where T(n,k) is the number of integer partitions of n with k co-modes.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 2, 1, 0, 4, 1, 0, 5, 2, 0, 7, 3, 1, 0, 10, 4, 1, 0, 13, 7, 2, 0, 16, 11, 3, 0, 23, 14, 4, 1, 0, 30, 19, 6, 1, 0, 35, 29, 11, 2, 0, 50, 34, 14, 3, 0, 61, 46, 23, 5, 0, 73, 69, 27, 6, 1, 0, 95, 81, 44, 10, 1, 0, 123, 105, 53, 14, 2
Offset: 0

Views

Author

Gus Wiseman, May 04 2023

Keywords

Comments

We define a co-mode in a multiset to be an element that appears at most as many times as each of the others. For example, the co-modes of {a,a,b,b,b,c,c} are {a,c}.

Examples

			Triangle begins:
   1
   0   1
   0   2
   0   2   1
   0   4   1
   0   5   2
   0   7   3   1
   0  10   4   1
   0  13   7   2
   0  16  11   3
   0  23  14   4   1
   0  30  19   6   1
   0  35  29  11   2
   0  50  34  14   3
   0  61  46  23   5
   0  73  69  27   6   1
   0  95  81  44  10   1
Row n = 8 counts the following partitions:
  (8)         (53)     (431)
  (44)        (62)     (521)
  (332)       (71)
  (422)       (3221)
  (611)       (3311)
  (2222)      (4211)
  (5111)      (32111)
  (22211)
  (41111)
  (221111)
  (311111)
  (2111111)
  (11111111)
		

Crossrefs

Row sums are A000041.
Row lengths are A002024.
Removing columns 0 and 1 and taking sums gives A362609, ranks A362606.
Column k = 1 is A362610, ranks A359178.
This statistic (co-mode count) is ranked by A362613.
For mode instead of co-mode we have A362614, ranked by A362611.
A008284 counts partitions by length.
A096144 counts partitions by number of minima, A026794 by maxima.
A238342 counts compositions by number of minima, A238341 by maxima.
A275870 counts collapsible partitions.

Programs

  • Mathematica
    comsi[ms_]:=Select[Union[ms],Count[ms,#]<=Min@@Length/@Split[ms]&];
    Table[Length[Select[IntegerPartitions[n],Length[comsi[#]]==k&]],{n,0,15},{k,0,Floor[(Sqrt[1+8n]-1)/2]}]

Formula

Sum_{k=0..A003056(n)} k * T(n,k) = A372632(n). - Alois P. Heinz, May 07 2024

A362610 Number of integer partitions of n having a unique part of least multiplicity.

Original entry on oeis.org

0, 1, 2, 2, 4, 5, 7, 10, 13, 16, 23, 30, 35, 50, 61, 73, 95, 123, 139, 187, 216, 269, 328, 411, 461, 594, 688, 836, 980, 1211, 1357, 1703, 1936, 2330, 2697, 3253, 3649, 4468, 5057, 6005, 6841, 8182, 9149, 10976, 12341, 14508, 16447, 19380, 21611, 25553, 28628
Offset: 0

Views

Author

Gus Wiseman, Apr 30 2023

Keywords

Comments

Alternatively, these are partitions with a part appearing fewer times than each of the others.

Examples

			The partition (3,3,2,2,2,1,1,1) has least multiplicity 2, and only one part of multiplicity 2 (namely 3), so is counted under a(15).
The a(1) = 1 through a(8) = 13 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (221)    (33)      (322)      (44)
                    (211)   (311)    (222)     (331)      (332)
                    (1111)  (2111)   (411)     (511)      (422)
                            (11111)  (3111)    (2221)     (611)
                                     (21111)   (4111)     (2222)
                                     (111111)  (22111)    (5111)
                                               (31111)    (22211)
                                               (211111)   (41111)
                                               (1111111)  (221111)
                                                          (311111)
                                                          (2111111)
                                                          (11111111)
		

Crossrefs

For parts instead of multiplicities we have A002865, ranks A247180.
For median instead of co-mode we have A238478, complement A238479.
These partitions have ranks A359178.
For mode complement of co-mode we have A362607, ranks A362605.
For mode instead of co-mode we have A362608, ranks A356862.
The complement is counted by A362609, ranks A362606.
A000041 counts integer partitions.
A275870 counts collapsible partitions.
A359893 counts partitions by median.
A362611 counts modes in prime factorization, co-modes A362613.
A362614 counts partitions by number of modes, co-modes A362615.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Count[Length/@Split[#],Min@@Length/@Split[#]]==1&]],{n,0,30}]
  • PARI
    seq(n) = my(A=O(x*x^n)); Vec(sum(m=2, n+1, sum(j=1, n, x^(j*(m-1))/(1 + if(j*m<=n, x^(j*m)/(1-x^j) )) + A)*prod(j=1, n\m, 1 + x^(j*m)/(1 - x^j) + A)), -(n+1)) \\ Andrew Howroyd, May 04 2023

Formula

G.f.: Sum_{m>=2} (Sum_{j>=1} x^(j*(m-1))/(1 + x^(j*m)/(1 - x^j))) * (Product_{j>=1} (1 + x^(j*m)/(1 - x^j))). - Andrew Howroyd, May 04 2023

A209816 Number of partitions of 2n in which every part is

Original entry on oeis.org

1, 3, 7, 15, 30, 58, 105, 186, 318, 530, 863, 1380, 2164, 3345, 5096, 7665, 11395, 16765, 24418, 35251, 50460, 71669, 101050, 141510, 196888, 272293, 374423, 512081, 696760, 943442, 1271527, 1706159, 2279700, 3033772, 4021695, 5311627, 6990367, 9168321
Offset: 1

Views

Author

Clark Kimberling, Mar 13 2012

Keywords

Comments

Also, the number of partitions of 3n in which n is the maximal part.
Also, the number of partitions of 3n into n parts. - Seiichi Manyama, May 07 2018
Also the number of multigraphical partitions of 2n, i.e., integer partitions that comprise the multiset of vertex-degrees of some multigraph. - Gus Wiseman, Oct 24 2018
Also number of partitions of 2n with at most n parts. Conjugate partitions map one to one to partitions of 2*n with each part <= n. - Wolfdieter Lang, May 21 2019

Examples

			The 7 partitions of 6 with parts <4 are as follows:
3+3, 3+2+1, 3+1+1+1
2+2+2, 2+2+1+1, 2+1+1+1+1
1+1+1+1+1+1.
Matching partitions of 2 into rationals as described:
1 + 1
1 + 3/3 + 1/3
1 + 1/3 + 1/3 + 1/3
2/3 + 2/3 + 2/3
2/3 + 2/3 + 1/3 + 1/3
2/3 + 1/3 + 1/3 + 1/3 + 1/3
1/3 + 1/3 + 1/3 + 1/3 + 1/3 + 1/3.
From _Seiichi Manyama_, May 07 2018: (Start)
n | Partitions of 3n into n parts
--+-------------------------------------------------
1 | 3;
2 | 5+1, 4+2, 3+3;
3 | 7+1+1, 6+2+1, 5+3+1, 5+2+2, 4+4+1, 4+3+2, 3+3+3; (End)
From _Gus Wiseman_, Oct 24 2018: (Start)
The a(1) = 1 through a(4) = 15 partitions:
  (11)  (22)    (33)      (44)
        (211)   (222)     (332)
        (1111)  (321)     (422)
                (2211)    (431)
                (3111)    (2222)
                (21111)   (3221)
                (111111)  (3311)
                          (4211)
                          (22211)
                          (32111)
                          (41111)
                          (221111)
                          (311111)
                          (2111111)
                          (11111111)
(End)
		

Crossrefs

Programs

  • Haskell
    a209816 n = p [1..n] (2*n) where
       p _          0 = 1
       p []         _ = 0
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Nov 14 2013
  • Maple
    b:= proc(n, i) option remember;
          `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i))))
        end:
    a:= n-> b(2*n, n):
    seq(a(n), n=1..50);  # Alois P. Heinz, Jul 09 2012
  • Mathematica
    f[n_] := Length[Select[IntegerPartitions[2 n], First[#] <= n &]]; Table[f[n], {n, 1, 30}] (* A209816 *)
    Table[SeriesCoefficient[Product[1/(1-x^k),{k,1,n}],{x,0,2*n}],{n,1,20}] (* Vaclav Kotesovec, May 25 2015 *)
    Table[Length@IntegerPartitions[3n, {n}], {n, 25}] (* Vladimir Reshetnikov, Jul 24 2016 *)
    b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]; a[n_] := b[2*n, n]; Table[a[n], {n, 1, 50}] (* Jean-François Alcover, Aug 29 2016, after Alois P. Heinz *)

Formula

a(n) = A000041(2*n)-A000070(n-1). - Matthew Vandermast, Jul 16 2012
a(n) = Sum_{k=1..n} A008284(2*n, k) = A000041(2*n) - A000070(n-1), for n >= 1. - Wolfdieter Lang, May 21 2019

Extensions

More terms from Alois P. Heinz, Jul 09 2012

A026810 Number of partitions of n in which the greatest part is 4.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 2, 3, 5, 6, 9, 11, 15, 18, 23, 27, 34, 39, 47, 54, 64, 72, 84, 94, 108, 120, 136, 150, 169, 185, 206, 225, 249, 270, 297, 321, 351, 378, 411, 441, 478, 511, 551, 588, 632, 672, 720, 764, 816, 864, 920, 972, 1033, 1089, 1154, 1215, 1285, 1350
Offset: 0

Views

Author

Keywords

Comments

Also number of partitions of n into exactly 4 parts.
Also the number of weighted cubic graphs on 4 nodes (=the tetrahedron) with weight n. - R. J. Mathar, Nov 03 2018
From Gus Wiseman, Jun 27 2021: (Start)
Also the number of strict integer partitions of 2n with alternating sum 4, or (by conjugation) partitions of 2n covering an initial interval of positive integers with exactly 4 odd parts. The strict partitions with alternating sum 4 are:
(4) (5,1) (6,2) (7,3) (8,4) (9,5) (10,6)
(5,2,1) (5,3,2) (5,4,3) (6,5,3) (7,6,3)
(6,3,1) (6,4,2) (7,5,2) (8,6,2)
(7,4,1) (8,5,1) (9,6,1)
(6,3,2,1) (6,4,3,1) (6,5,4,1)
(7,4,2,1) (7,4,3,2)
(7,5,3,1)
(8,5,2,1)
(6,4,3,2,1)
(End)

Examples

			From _Gus Wiseman_, Jun 27 2021: (Start)
The a(4) = 1 through a(10) = 9 partitions of length 4:
  (1111)  (2111)  (2211)  (2221)  (2222)  (3222)  (3322)
                  (3111)  (3211)  (3221)  (3321)  (3331)
                          (4111)  (3311)  (4221)  (4222)
                                  (4211)  (4311)  (4321)
                                  (5111)  (5211)  (4411)
                                          (6111)  (5221)
                                                  (5311)
                                                  (6211)
                                                  (7111)
(End)
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 275.
  • D. E. Knuth, The Art of Computer Programming, vol. 4,fascicle 3, Generating All Combinations and Partitions, Section 7.2.1.4., p. 56, exercise 31.

Crossrefs

Cf. A001400, A026811, A026812, A026813, A026814, A026815, A026816, A069905 (3 positive parts), A002621 (partial sums), A005044 (first differences).
A non-strict version is A000710 or A088218.
This is column k = 2 of A152146.
A reverse version is A343941.

Programs

  • Magma
    [Round((n^3+3*n^2-9*n*(n mod 2))/144): n in [0..60]]; // Vincenzo Librandi, Oct 14 2015
  • Maple
    A049347 := proc(n)
            op(1+(n mod 3),[1,-1,0]) ;
    end proc:
    A056594 := proc(n)
            op(1+(n mod 4),[1,0,-1,0]) ;
    end proc:
    A026810 := proc(n)
            1/288*(n+1)*(2*n^2+4*n-13+9*(-1)^n) ;
            %-A049347(n)/9 ;
            %+A056594(n)/8 ;
    end proc: # R. J. Mathar, Jul 03 2012
  • Mathematica
    Table[Count[IntegerPartitions[n], {4, _}], {n, 0, 60}]
    LinearRecurrence[{1, 1, 0, 0, -2, 0, 0, 1, 1, -1}, {0, 0, 0, 0, 1, 1, 2, 3, 5, 6}, 60] (* Vincenzo Librandi, Oct 14 2015 *)
    Table[Length[IntegerPartitions[n, {4}]], {n, 0, 60}] (* Eric Rowland, Mar 02 2017 *)
    CoefficientList[Series[x^4/Product[1 - x^k, {k, 1, 4}], {x, 0, 60}], x] (* Robert A. Russell, May 13 2018 *)
  • PARI
    for(n=0, 60, print(n, " ", round((n^3 + 3*n^2 -9*n*(n % 2))/144))); \\ Washington Bomfim, Jul 03 2012
    
  • PARI
    x='x+O('x^60); concat([0, 0, 0, 0], Vec(x^4/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)))) \\ Altug Alkan, Oct 14 2015
    
  • PARI
    vector(60, n, n--; (n+1)*(2*n^2+4*n-13+9*(-1)^n)/288 + real(I^n)/8 - ((n+2)%3-1)/9) \\ Altug Alkan, Oct 26 2015
    
  • PARI
    print1(0,", "); for(n=1,60,j=0;forpart(v=n,j++,,[4,4]); print1(j,", ")) \\ Hugo Pfoertner, Oct 01 2018
    

Formula

G.f.: x^4/((1-x)*(1-x^2)*(1-x^3)*(1-x^4)) = x^4/((1-x)^4*(1+x)^2*(1+x+x^2)*(1+x^2)).
a(n+4) = A001400(n). - Michael Somos, Apr 07 2012
a(n) = round( (n^3 + 3*n^2 -9*n*(n mod 2))/144 ). - Washington Bomfim, Jan 06 2021 and Jul 03 2012
a(n) = (n+1)*(2*n^2+4*n-13+9*(-1)^n)/288 -A049347(n)/9 +A056594(n)/8. - R. J. Mathar, Jul 03 2012
From Gregory L. Simay, Oct 13 2015: (Start)
a(n) = (n^3 + 3*n^2 - 9*n)/144 + a(m) - (m^3 + 3*m^2 - 9*m)/144 if n = 12k + m and m is odd. For example, a(23) = a(12*1 + 11) = (23^3 + 3*23^2 - 9*23)/144 + a(11) - (11^3 + 3*11^2 - 9*11)/144 = 94.
a(n) = (n^3 + 3*n^2)/144 + a(m) - (m^3 + 3*m^2)/144 if n = 12k + m and m is even. For example, a(22) = a(12*1 + 10) = (22^3 + 3*22^2)/144 + a(10) - (10^3 + 3*10^2)/144 = 84. (End)
a(n) = A008284(n,4). - Robert A. Russell, May 13 2018
From Gregory L. Simay, Jul 28 2019: (Start)
a(2n+1) = a(2n) + a(n+1) - a(n-3) and
a(2n) = a(2n-1) + a(n+2) - a(n-2). (End)

A364916 Array read by antidiagonals downwards where A(n,k) is the number of ways to write n as a nonnegative linear combination of the parts of a strict integer partition of k.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 17 2023

Keywords

Comments

A way of writing n as a (nonnegative) linear combination of a finite sequence y is any sequence of pairs (k_i,y_i) such that k_i >= 0 and Sum k_i*y_i = n. For example, the pairs ((3,1),(1,1),(1,1),(0,2)) are a way of writing 5 as a linear combination of (1,1,1,2), namely 5 = 3*1 + 1*1 + 1*1 + 0*2. Of course, there are A000041(n) ways to write n as a linear combination of (1..n).
As a triangle, also the number of ways to write n as a *positive* linear combination of the parts of a strict integer partition of k.

Examples

			Array begins:
  1  1  1  2  2  3  4   5   6   8   10   12  15   18   22   27
  0  1  0  1  1  1  2   2   3   3   5    5   7    8    10   12
  0  1  1  2  1  2  4   4   5   6   9    10  13   15   19   23
  0  1  0  3  2  2  4   4   6   7   11   11  15   17   22   27
  0  1  1  3  3  3  7   7   8   10  16   17  23   27   33   42
  0  1  0  3  2  4  7   6   9   9   17   17  23   26   33   43
  0  1  1  5  3  4  12  10  13  16  26   27  36   42   52   68
  0  1  0  4  3  3  10  11  13  13  27   25  35   40   51   67
  0  1  1  5  4  5  15  13  19  20  36   37  51   58   72   97
  0  1  0  6  4  5  14  13  18  23  42   39  54   61   78   105
  0  1  1  6  4  6  20  17  23  25  54   50  69   80   98   138
  0  1  0  6  4  5  19  16  23  24  54   55  71   80   103  144
  0  1  1  8  6  7  27  23  30  35  72   70  103  113  139  199
  0  1  0  7  5  6  24  21  29  31  75   68  95   115  139  201
  0  1  1  8  5  7  31  27  36  39  90   86  122  137  178  255
  0  1  0  9  6  8  31  27  38  42  100  93  129  148  187  289
Triangle begins:
   1
   1  0
   1  1  0
   2  0  1  0
   2  1  1  1  0
   3  1  2  0  1  0
   4  1  1  3  1  1  0
   5  2  2  2  3  0  1  0
   6  2  4  2  3  3  1  1  0
   8  3  4  4  3  2  5  0  1  0
  10  3  5  4  7  4  3  4  1  1  0
  12  5  6  6  7  7  4  3  5  0  1  0
  15  5  9  7  8  6 12  3  4  6  1  1  0
  18  7 10 11 10  9 10 10  5  4  6  0  1  0
  22  8 13 11 16  9 13 11 15  5  4  6  1  1  0
  27 10 15 15 17 17 16 13 13 14  6  4  8  0  1  0
		

Crossrefs

Same as A116861 with offset 0 and rows reversed, non-strict version A364912.
Row n = 0 is A000009.
Row n = 1 is A096765.
Row n = 2 is A365005.
Column k = 0 is A000007.
Column k = 1 is A000012.
Column k = 2 is A000035.
Column k = 3 is A137719.
The main diagonal is A364910.
Left half has row sums A365002.
For not just strict partitions we have A365004, diagonal A364907.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A066328 adds up distinct prime indices.
A364350 counts combination-free strict partitions, complement A364839.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,0,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]];
    t[n_,k_]:=Length[Join@@Table[combs[n,ptn],{ptn,Select[IntegerPartitions[k],UnsameQ@@#&]}]];
    Table[t[k,n-k],{n,0,15},{k,0,n}]

A325242 Irregular triangle read by rows with zeros removed where T(n,k) is the number of integer partitions of n with k distinct multiplicities, n > 0.

Original entry on oeis.org

1, 2, 3, 4, 1, 4, 3, 8, 3, 6, 9, 10, 12, 11, 19, 15, 26, 1, 13, 39, 4, 25, 47, 5, 19, 70, 12, 29, 89, 17, 33, 115, 28, 42, 148, 41, 39, 189, 69, 62, 235, 88, 55, 294, 141, 81, 362, 183, 1, 84, 450, 253, 5, 103, 558, 333, 8, 105, 669, 464, 17, 153, 817, 576, 29
Offset: 1

Views

Author

Gus Wiseman, Apr 15 2019

Keywords

Comments

For example, the partition (32111) has multiplicities {1,1,3}, of which 2 are distinct, so is counted under T(8,2).

Examples

			Triangle begins:
   1
   2
   3
   4   1
   4   3
   8   3
   6   9
  10  12
  11  19
  15  26   1
  13  39   4
  25  47   5
  19  70  12
  29  89  17
  33 115  28
  42 148  41
  39 189  69
  62 235  88
  55 294 141
  81 362 183   1
Row n = 8 counts the following partitions:
  (8)         (332)
  (44)        (422)
  (53)        (611)
  (62)        (3221)
  (71)        (4211)
  (431)       (5111)
  (521)       (22211)
  (2222)      (32111)
  (3311)      (41111)
  (11111111)  (221111)
              (311111)
              (2111111)
		

Crossrefs

Row lengths are A056556. Row sums are A000041. Column k = 1 is A047966. Column k = 2 is A325243.

Programs

  • Mathematica
    DeleteCases[Table[Length[Select[IntegerPartitions[n],Length[Union[Length/@Split[#]]]==k&]],{n,20},{k,n}],0,2]

A004250 Number of partitions of n into 3 or more parts.

Original entry on oeis.org

0, 0, 1, 2, 4, 7, 11, 17, 25, 36, 50, 70, 94, 127, 168, 222, 288, 375, 480, 616, 781, 990, 1243, 1562, 1945, 2422, 2996, 3703, 4550, 5588, 6826, 8332, 10126, 12292, 14865, 17958, 21618, 25995, 31165, 37317, 44562
Offset: 1

Views

Author

Keywords

Comments

Number of (n+1)-vertex spider graphs: trees with n+1 vertices and exactly 1 vertex of degree at least 3 (i.e. branching vertex). There is a trivial bijection with the objects described in the definition. - Emeric Deutsch, Feb 22 2014
Also the number of graphical partitions of 2n into n parts. - Gus Wiseman, Jan 08 2021

Examples

			a(6)=7 because there are three partitions of n=6 with i=3 parts: [4, 1, 1], [3, 2, 1], [2, 2, 2] and two partitions with i=4 parts: [3, 1, 1, 1], [2, 2, 1, 1] and one partition with i=5 parts: [2, 1, 1, 1, 1] and one partition with i=6 parts: [1, 1, 1, 1, 1, 1].
From _Gus Wiseman_, Jan 18 2021: (Start)
The a(3) = 1 through a(7) = 11 graphical partitions of 2n into n parts:
  (222)  (2222)  (22222)  (222222)  (2222222)
         (3221)  (32221)  (322221)  (3222221)
                 (33211)  (332211)  (3322211)
                 (42211)  (333111)  (3332111)
                          (422211)  (4222211)
                          (432111)  (4322111)
                          (522111)  (4331111)
                                    (4421111)
                                    (5222111)
                                    (5321111)
                                    (6221111)
(End)
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • P. R. Stein, On the number of graphical partitions, pp. 671-684 of Proc. 9th S-E Conf. Combinatorics, Graph Theory, Computing, Congr. Numer. 21 (1978).

Crossrefs

Rightmost column of A259873.
Central column of A339659.
A000041 counts partitions of 2n into n parts, ranked by A340387.
A000569 counts graphical partitions, ranked by A320922.
A008284 counts partitions by sum and length.
A027187 counts partitions of even length.
A309356 ranks simple covering graphs.
The following count vertex-degree partitions and give their Heinz numbers:
- A209816 counts multigraphical partitions (A320924).
- A320921 counts connected graphical partitions (A320923).
- A339617 counts non-graphical partitions of 2n (A339618).
- A339656 counts loop-graphical partitions (A339658).
Partial sums of A117995.

Programs

  • Maple
    with(combinat);
    for i from 1 to 15 do pik(i,3) od;
    pik:= proc(n::integer, k::integer)
    # Thomas Wieder, Jan 30 2007
    local i, Liste, Result;
    if k > n or n < 0 or k < 1 then
    return fail
    end if;
    Result := 0;
    for i from k to n do
    Liste:= PartitionList(n,i);
    #print(Liste);
    Result := Result + nops(Liste);
    end do;
    return Result;
    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 at: http://www.cis.upenn.edu/~wilf/lecnotes.html
    # Calculates the partition of n into k parts.
    # E.g. PartitionList(5,2) --> [[4, 1], [3, 2]].
    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;
    #  Thomas Wieder, Feb 01 2007
    ZL :=[S, {S = Set(Cycle(Z),3 <= card)}, unlabelled]: seq(combstruct[count](ZL, size=n), n=1..41); # Zerinvary Lajos, Mar 25 2008
    B:=[S,{S = Set(Sequence(Z,1 <= card),card >=3)},unlabelled]: seq(combstruct[count](B, size=n), n=1..41); # Zerinvary Lajos, Mar 21 2009
  • Mathematica
    Length /@ Table[Select[Partitions[n], Length[#] > 2 &], {n, 20}] (* Eric W. Weisstein, May 16 2007 *)
    Table[Count[Length /@ Partitions[n], ?(# > 2 &)], {n, 20}] (* _Eric W. Weisstein, May 16 2017 *)
    Table[PartitionsP[n] - Floor[n/2] - 1, {n, 20}] (* Eric W. Weisstein, May 16 2017 *)
    Length /@ Table[IntegerPartitions[n, {3, n}], {n, 20}] (* Eric W. Weisstein, May 16 2017 *)
  • PARI
    a(n) = numbpart(n) - (n+2)\2; /* Joerg Arndt, Apr 03 2013 */

Formula

G.f.: Sum_{n>=0} (q^n / Product_{k=1..n+3} (1 - q^k)). - N. J. A. Sloane
a(n) = A000041(n) - floor((n+2)/2) = A000041(n)-A004526(n+2) = A058984(n)-1. - Vladeta Jovovic, Jun 18 2003
Let P(n,i) denote the number of partitions of n into i parts. Then a(n) = Sum_{i=3..n} P(n,i). - Thomas Wieder, Feb 01 2007
a(n) = A259873(n,n). - Gus Wiseman, Jan 08 2021

Extensions

Definition corrected by Thomas Wieder, Feb 01 2007 and by Eric W. Weisstein, May 16 2007

A220377 Number of partitions of n into three distinct and mutually relatively prime parts.

Original entry on oeis.org

1, 0, 2, 1, 3, 1, 6, 1, 7, 3, 7, 3, 14, 3, 15, 6, 14, 6, 25, 6, 22, 10, 25, 9, 42, 8, 34, 15, 37, 15, 53, 13, 48, 22, 53, 17, 78, 17, 65, 30, 63, 24, 99, 24, 88, 35, 84, 30, 126, 34, 103, 45, 103, 38, 166, 35, 124, 57, 128, 51, 184, 44, 150, 67, 172, 52, 218
Offset: 6

Views

Author

Carl Najafi, Dec 13 2012

Keywords

Comments

The Heinz numbers of these partitions are the intersection of A005117 (strict), A014612 (triples), and A302696 (coprime). - Gus Wiseman, Oct 14 2020

Examples

			For n=10 we have three such partitions: 1+2+7, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 14 2020: (Start)
The a(6) = 1 through a(20) = 15 triples (empty column indicated by dot, A..H = 10..17):
321  .  431  531  532  731  543  751  743  753  754  971  765  B53  875
        521       541       651       752  951  853  B51  873  B71  974
                  721       732       761  B31  871  D31  954  D51  A73
                            741       851       952       972       A91
                            831       941       B32       981       B54
                            921       A31       B41       A71       B72
                                      B21       D21       B43       B81
                                                          B52       C71
                                                          B61       D43
                                                          C51       D52
                                                          D32       D61
                                                          D41       E51
                                                          E31       F41
                                                          F21       G31
                                                                    H21
(End)
		

Crossrefs

A023022 is the 2-part version.
A101271 is the relative prime instead of pairwise coprime version.
A220377*6 is the ordered version.
A305713 counts these partitions of any length, with Heinz numbers A302797.
A307719 is the non-strict version.
A337461 is the non-strict ordered version.
A337563 is the case with no 1's.
A337605 is the pairwise non-coprime instead of pairwise coprime version.
A001399(n-6) counts strict 3-part partitions, with Heinz numbers A007304.
A008284 counts partitions by sum and length, with strict case A008289.
A318717 counts pairwise non-coprime strict partitions.
A326675 ranks pairwise coprime sets.
A327516 counts pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.

Programs

  • Mathematica
    Table[Length@Select[ IntegerPartitions[ n, {3}], #[[1]] != #[[2]] != #[[3]] && GCD[#[[1]], #[[2]]] == 1 && GCD[#[[1]], #[[3]]] == 1 && GCD[#[[2]], #[[3]]] == 1 &], {n, 6, 100}]
    Table[Count[IntegerPartitions[n,{3}],?(CoprimeQ@@#&&Length[ Union[#]] == 3&)],{n,6,100}] (* _Harvey P. Dale, May 22 2020 *)
  • PARI
    a(n)=my(P=partitions(n));sum(i=1,#P,#P[i]==3&&P[i][1]Charles R Greathouse IV, Dec 14 2012

Formula

a(n > 2) = A307719(n) - 1. - Gus Wiseman, Oct 15 2020
Previous Showing 81-90 of 579 results. Next