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 11-20 of 111 results. Next

A188431 The number of n-full sets, F(n).

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 2, 2, 1, 2, 1, 2, 3, 4, 5, 7, 7, 8, 9, 11, 10, 13, 14, 17, 20, 25, 28, 34, 40, 46, 54, 62, 69, 80, 90, 102, 115, 131, 144, 167, 186, 213, 239, 273, 304, 349, 388, 441, 495, 563, 625, 710, 790, 890, 990, 1114, 1232, 1387, 1530, 1713, 1894, 2119, 2330, 2605, 2866, 3192, 3512, 3910, 4289, 4774, 5237, 5809, 6377, 7068, 7739
Offset: 0

Views

Author

Madjid Mirzavaziri, Mar 31 2011

Keywords

Comments

Let A be a set of positive integers. We say that A is n-full if (sum A)=[n] for a positive integer n, where (sum A) is the set of all positive integers which are a sum of distinct elements of A and [n]={1,2,...,n}. Then F(n) denotes the number of n-full sets.
Also the number of distinct and complete partitions of n, by definition, which are counted by A000009 and A126796. - George Beck, Nov 06 2017
An integer partition of n is complete (see also A325781) if every number from 0 to n is the sum of some submultiset of the parts. The Heinz numbers of these partitions are given by A325986. - Gus Wiseman, May 31 2019

Examples

			a(26) = 10, because there are 10 26-full sets: {1,2,4,5,6,8}, {1,2,3,5,7,8}, {1,2,3,5,6,9}, {1,2,3,4,7,9}, {1,2,3,4,6,10}, {1,2,3,4,5,11}, {1,2,4,8,11}, {1,2,4,7,12}, {1,2,4,6,13}, {1,2,3,7,13}.
G.f.: 1 = 1/(1+x) + 1*x/((1+x)*(1+x^2)) + 0*x^2/((1+x)*(1+x^2)*(1+x^3)) + 1*x^3/((1+x)*(1+x^2)*(1+x^3)*(1+x^4)) +...+ a(n)*x^n / Product_{k=1..n+1} (1+x^k) +...
		

Crossrefs

Programs

  • Haskell
    import Data.MemoCombinators (memo2, integral, Memo)
    a188431 n = a188431_list !! (n-1)
    a188431_list = map
       (\x -> sum [fMemo x i | i <- [a188429 x .. a188430 x]]) [1..] where
       fMemo = memo2 integral integral f
       f _ 1 = 1
       f m i = sum [fMemo (m - i) j |
                    j <- [a188429 (m - i) .. min (a188430 (m - i)) (i - 1)]]
    -- Reinhard Zumkeller, Aug 06 2015
  • Maple
    sums:= proc(s) local i, m;
              m:= max(s[]);
             `if`(m<1, {}, {m, seq([i, i+m][], i=sums(s minus {m}))})
           end:
    a:= proc(n) local b;
          b:= proc(i,s) local si;
                if i=1 then `if`(sums(s)={$1..n}, 1, 0)
              else si:= s union {i};
                   b(i-1, s)+ `if`(max(sums(si)[])>n, 0, b(i-1, si))
                fi
              end; b(n, {1})
        end:
    seq(a(n), n=1..40);  # Alois P. Heinz, Apr 03 2011
    # second Maple program:
    b:= proc(n, i) option remember; `if`(i*(i+1)/2n or i>n-i+1, 0, b(n-i, i-1))))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..80);  # Alois P. Heinz, May 20 2017
  • Mathematica
    Sums[s_] := Sums[s] = With[{m = Max[s]}, If[m < 1, {}, Union @ Flatten @ Join[{m}, Table[{i, i + m}, {i, Sums[s ~Complement~ {m}]}]]]];
    a[n_] := Module[{b}, b[i_, s_] := b[i, s] = Module[{si}, If[i == 1, If[Sums[s] == Range[n], 1, 0], si = s ~Union~ {i}; b[i-1, s] + If[Max[ Sums[si]] > n, 0, b[i - 1, si]]]]; b[n, {1}]];
    Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 1, 80}] (* Jean-François Alcover, Apr 12 2017, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&Union[Total/@Union[Subsets[#]]]==Range[0,n]&]],{n,30}] (* Gus Wiseman, May 31 2019 *)
  • PARI
    /* As coefficients in g.f. */
    {a(n)=local(A=[1]); for(i=1, n+1, A=concat(A,0); A[#A]=polcoeff(1 - sum(m=1,#A,A[m]*x^m/prod(k=1, m, 1+x^k +x*O(x^#A) )), #A) ); A[n+1]}
    for(n=0, 50, print1(a(n),", ")) /* Paul D. Hanna, Mar 06 2012 */
    

Formula

F(n) = Sum_(i=L(n) .. U(n), F(n,i)), where F(n,i) = Sum_(j=L(n-i) .. min(U(n-i),i-1), F(n-i,j) ) and L(n), U(n) are defined in A188429 and A188430, respectively.
G.f.: 1 = Sum_{n>=0} a(n)*x^n / Product_{k=1..n+1} (1+x^k), with a(0)=1. - Paul D. Hanna, Mar 08 2012
a(n) ~ c * exp(Pi*sqrt(n/3)) / n^(3/4), where c = 0.03316508... - Vaclav Kotesovec, Oct 21 2019

Extensions

More terms from Alois P. Heinz, Apr 03 2011
a(0)=1 prepended by Alois P. Heinz, May 20 2017

A316313 Number of integer partitions of n such that every distinct submultiset has a different average.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 5, 6, 6, 9, 11, 10, 15, 17, 18, 22, 28, 26, 34, 37, 44, 50, 61, 53, 77, 82, 93, 89, 120, 120, 149, 138, 174, 180, 222, 193, 257, 262, 305, 281, 367, 359, 424, 398, 487, 507, 590, 526, 662, 666, 782, 729, 894, 892, 995, 987, 1154, 1188, 1370
Offset: 1

Views

Author

Gus Wiseman, Jun 29 2018

Keywords

Comments

Note that such a partition is necessarily strict.

Examples

			The a(8) = 6 integer partitions are (8), (71), (62), (53), (521), (431).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@Mean/@Union[Subsets[#]]&]],{n,20}]

A326020 Number of complete subsets of {1..n}.

Original entry on oeis.org

1, 2, 3, 4, 6, 9, 15, 27, 50, 95, 185, 365, 724, 1441, 2873, 5735, 11458, 22902, 45789, 91561, 183102, 366180, 732331, 1464626, 2929209, 5858367, 11716674, 23433277, 46866473, 93732852, 187465596, 374931067, 749861989, 1499723808, 2999447418
Offset: 0

Views

Author

Gus Wiseman, Jun 04 2019

Keywords

Comments

A set of positive integers summing to n is complete if every nonnegative integer up to n is the sum of some subset.

Examples

			The a(0) = 1 through a(6) = 15 subsets:
  {}  {}   {}     {}       {}         {}           {}
      {1}  {1}    {1}      {1}        {1}          {1}
           {1,2}  {1,2}    {1,2}      {1,2}        {1,2}
                  {1,2,3}  {1,2,3}    {1,2,3}      {1,2,3}
                           {1,2,4}    {1,2,4}      {1,2,4}
                           {1,2,3,4}  {1,2,3,4}    {1,2,3,4}
                                      {1,2,3,5}    {1,2,3,5}
                                      {1,2,4,5}    {1,2,3,6}
                                      {1,2,3,4,5}  {1,2,4,5}
                                                   {1,2,4,6}
                                                   {1,2,3,4,5}
                                                   {1,2,3,4,6}
                                                   {1,2,3,5,6}
                                                   {1,2,4,5,6}
                                                   {1,2,3,4,5,6}
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],Union[Plus@@@Subsets[#]]==Range[0,Total[#]]&]],{n,0,10}]

Extensions

a(17)-a(34) from Charlie Neder, Jun 05 2019

A237258 Number of strict partitions of 2n that include a partition of n.

Original entry on oeis.org

1, 0, 0, 1, 1, 3, 4, 7, 9, 16, 21, 32, 43, 63, 84, 122, 158, 220, 293, 393, 511, 685, 881, 1156, 1485, 1925, 2445, 3147, 3952, 5019, 6323, 7924, 9862, 12336, 15259, 18900, 23294, 28646, 35091, 42985, 52341, 63694, 77336, 93588, 112973, 136367, 163874, 196638
Offset: 0

Views

Author

Clark Kimberling, Feb 05 2014

Keywords

Comments

A strict partition is a partition into distinct parts.

Examples

			a(5) counts these partitions of 10: [5,4,1], [5,3,2], [4,3,2,1].
		

Crossrefs

The non-strict version is A002219, ranked by A357976.
These partitions are ranked by A357854.
A000712 counts distinct submultisets of partitions, strict A032302.
A304792 counts subset-sums of partitions, positive A276024, strict A284640.

Programs

  • Mathematica
    z = 24; Table[theTotals = Map[{#, Map[Total, Subsets[#]]} &,  Select[IntegerPartitions[2 nn], # == DeleteDuplicates[#] &]]; Length[Map[#[[1]] &, Select[theTotals, Length[Position[#[[2]], nn]] >= 1 &]]], {nn, z}] (* Peter J. C. Moses, Feb 04 2014 *)

Formula

a(n) = A237194(2n,n).

Extensions

a(31)-a(47) from Alois P. Heinz, Feb 07 2014

A333224 Number of distinct positive consecutive subsequence-sums of the k-th composition in standard order.

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 3, 3, 1, 3, 2, 4, 3, 4, 4, 4, 1, 3, 3, 5, 3, 5, 4, 5, 3, 4, 5, 5, 5, 5, 5, 5, 1, 3, 3, 5, 2, 5, 5, 6, 3, 6, 3, 6, 5, 6, 5, 6, 3, 4, 6, 6, 5, 6, 6, 6, 5, 6, 6, 6, 6, 6, 6, 6, 1, 3, 3, 5, 3, 6, 6, 7, 3, 5, 5, 7, 4, 6, 6, 7, 3, 6, 4, 7, 5, 7, 6
Offset: 0

Views

Author

Gus Wiseman, Mar 18 2020

Keywords

Comments

The k-th composition in standard order (row k of A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again.

Examples

			The composition (4,3,1,2) has positive subsequence-sums 1, 2, 3, 4, 6, 7, 8, 10, so a(550) = 8.
		

Crossrefs

Dominated by A124770.
Compositions where every subinterval has a different sum are counted by A169942 and A325677 and ranked by A333222. The case of partitions is counted by A325768 and ranked by A325779.
Positive subset-sums of partitions are counted by A276024 and A299701.
Knapsack partitions are counted by A108917 and A325592 and ranked by A299702.
Strict knapsack partitions are counted by A275972 and ranked by A059519 and A301899.
Knapsack compositions are counted by A325676 and A325687 and ranked by A333223. The case of partitions is counted by A325769 and ranked by A325778, for which the number of distinct consecutive subsequences is given by A325770.
Allowing empty subsequences gives A333257.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Table[Length[Union[ReplaceList[stc[n],{_,s__,_}:>Plus[s]]]],{n,0,100}]

Formula

a(n) = A333257(n) - 1.

A333257 Number of distinct consecutive subsequence-sums of the k-th composition in standard order.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 4, 4, 2, 4, 3, 5, 4, 5, 5, 5, 2, 4, 4, 6, 4, 6, 5, 6, 4, 5, 6, 6, 6, 6, 6, 6, 2, 4, 4, 6, 3, 6, 6, 7, 4, 7, 4, 7, 6, 7, 6, 7, 4, 5, 7, 7, 6, 7, 7, 7, 6, 7, 7, 7, 7, 7, 7, 7, 2, 4, 4, 6, 4, 7, 7, 8, 4, 6, 6, 8, 5, 7, 7, 8, 4, 7, 5, 8, 6, 8, 7
Offset: 0

Views

Author

Gus Wiseman, Mar 20 2020

Keywords

Comments

A composition of n is a finite sequence of positive integers summing to n. The k-th composition in standard order (row k of A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again.

Examples

			The ninth composition in standard order is (3,1), which has consecutive subsequences (), (1), (3), (3,1), with sums 0, 1, 3, 4, so a(9) = 4.
		

Crossrefs

Dominated by A124771.
Compositions where every subinterval has a different sum are counted by A169942 and A325677 and ranked by A333222, while the case of partitions is counted by A325768 and ranked by A325779.
Positive subset-sums of partitions are counted by A276024 and A299701.
Knapsack partitions are counted by A108917 and ranked by A299702.
Knapsack compositions are counted by A325676 and A325687 and ranked by A333223.
The version for Heinz numbers of partitions is A325770.
Not allowing empty subsequences gives A333224.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Table[Length[Union[ReplaceList[stc[n],{_,s___,_}:>Plus[s]]]],{n,0,100}]

Formula

a(n) = A333224(n) + 1.

A304793 Number of distinct positive subset-sums of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 3, 2, 3, 1, 4, 1, 3, 3, 4, 1, 5, 1, 5, 3, 3, 1, 5, 2, 3, 3, 5, 1, 6, 1, 5, 3, 3, 3, 6, 1, 3, 3, 6, 1, 7, 1, 5, 5, 3, 1, 6, 2, 5, 3, 5, 1, 7, 3, 7, 3, 3, 1, 7, 1, 3, 4, 6, 3, 7, 1, 5, 3, 6, 1, 7, 1, 3, 5, 5, 3, 7, 1, 7, 4, 3, 1, 8, 3, 3, 3, 7, 1, 8, 3, 5, 3, 3, 3, 7, 1, 5, 5, 8, 1, 7, 1, 7, 7
Offset: 1

Views

Author

Gus Wiseman, May 18 2018

Keywords

Comments

A positive integer n is a positive subset-sum of an integer partition y if there exists a submultiset of y with sum n. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
a(n) <= A000005(n).
One less than the number of distinct values obtained when A056239 is applied to all divisors of n. - Antti Karttunen, Jul 01 2018

Examples

			The positive subset-sums of (4,3,1) are {1, 3, 4, 5, 7, 8} so a(70) = 6.
The positive subset-sums of (5,1,1,1) are {1, 2, 3, 5, 6, 7, 8} so a(88) = 7.
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Union[Total/@Rest[Subsets[Join@@Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]]],{n,100}]
  • PARI
    up_to = 65537;
    A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
    v056239 = vector(up_to,n,A056239(n));
    A304793(n) = { my(m=Map(),s,k=0); fordiv(n,d,if(!mapisdefined(m,s = v056239[d]), mapput(m,s,s); k++)); (k-1); }; \\ Antti Karttunen, Jul 01 2018

Extensions

More terms from Antti Karttunen, Jul 01 2018

A299729 Heinz numbers of non-knapsack partitions.

Original entry on oeis.org

12, 24, 30, 36, 40, 48, 60, 63, 70, 72, 80, 84, 90, 96, 108, 112, 120, 126, 132, 140, 144, 150, 154, 156, 160, 165, 168, 180, 189, 192, 198, 200, 204, 210, 216, 220, 224, 228, 240, 252, 264, 270, 273, 276, 280, 286, 288, 300, 308, 312, 315, 320, 324, 325
Offset: 1

Views

Author

Gus Wiseman, Feb 17 2018

Keywords

Comments

An integer partition is non-knapsack if there exist two different submultisets with the same sum. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			12 is the Heinz number of (2,1,1) which is not knapsack because 2 = 1 + 1.
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],!UnsameQ@@Plus@@@Union[Rest@Subsets[primeMS[#]]]&]

A321452 Number of integer partitions of n that can be partitioned into two or more blocks with equal sums.

Original entry on oeis.org

0, 0, 1, 1, 3, 1, 7, 1, 14, 10, 26, 1, 55, 1, 90, 68, 167, 1, 292, 1, 482, 345, 761, 1, 1291, 266, 1949, 1518, 3091, 1, 4793, 1, 7177, 5612, 10566, 2623, 16007, 1, 22912, 18992, 33619, 1, 48529, 1, 68758, 59187, 96571, 1, 137489, 11418, 189979, 167502, 264299
Offset: 0

Views

Author

Gus Wiseman, Nov 10 2018

Keywords

Comments

a(n) = 1 if and only if n is prime. - Chai Wah Wu, Nov 12 2018

Examples

			The a(2) = 1 through a(9) = 10 partitions:
  (11)  (111)  (22)    (11111)  (33)      (1111111)  (44)        (333)
               (211)            (222)                (422)       (3321)
               (1111)           (321)                (431)       (32211)
                                (2211)               (2222)      (33111)
                                (3111)               (3221)      (222111)
                                (21111)              (3311)      (321111)
                                (111111)             (4211)      (2211111)
                                                     (22211)     (3111111)
                                                     (32111)     (21111111)
                                                     (41111)     (111111111)
                                                     (221111)
                                                     (311111)
                                                     (2111111)
                                                     (11111111)
The partition (32111) can be partitioned as ((13)(112)), and the blocks both sum to 4, so (32111) is counted under a(8).
		

Crossrefs

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[IntegerPartitions[n],Length[Select[facs[Times@@Prime/@#],SameQ@@hwt/@#&]]>1&]],{n,10}]

Formula

a(n) = A000041(n) - A321451(n).

Extensions

a(26)-a(52) from Alois P. Heinz, Nov 11 2018

A321451 Number of integer partitions of n that cannot be partitioned into two or more blocks with equal sums.

Original entry on oeis.org

1, 1, 1, 2, 2, 6, 4, 14, 8, 20, 16, 55, 22, 100, 45, 108, 64, 296, 93, 489, 145, 447, 241, 1254, 284, 1692, 487, 1492, 627, 4564, 811, 6841, 1172, 4531, 1744, 12260, 1970, 21636, 3103, 12193, 3719, 44582, 4645, 63260, 6417, 29947, 8987, 124753, 9784, 162107, 14247
Offset: 0

Views

Author

Gus Wiseman, Nov 10 2018

Keywords

Examples

			The a(1) = 1 through a(9) = 20 partitions:
  (1)  (2)  (3)   (4)   (5)     (6)    (7)       (8)     (9)
            (21)  (31)  (32)    (42)   (43)      (53)    (54)
                        (41)    (51)   (52)      (62)    (63)
                        (221)   (411)  (61)      (71)    (72)
                        (311)          (322)     (332)   (81)
                        (2111)         (331)     (521)   (432)
                                       (421)     (611)   (441)
                                       (511)     (5111)  (522)
                                       (2221)            (531)
                                       (3211)            (621)
                                       (4111)            (711)
                                       (22111)           (3222)
                                       (31111)           (4221)
                                       (211111)          (4311)
                                                         (5211)
                                                         (6111)
                                                         (22221)
                                                         (42111)
                                                         (51111)
                                                         (411111)
A complete list of all multiset partitions of the partition (2111) into two or more blocks is: ((1)(112)), ((2)(111)), ((11)(12)), ((1)(1)(12)), ((1)(2)(11)), ((1)(1)(1)(2)). None of these has equal block-sums, so (2111) is counted toward a(5).
On the other hand, the partition (321) can be partitioned as ((12)(3)), which has two or more blocks and equal block-sums, so (321) is not counted toward a(6).
		

Crossrefs

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[IntegerPartitions[n],Length[Select[facs[Times@@Prime/@#],SameQ@@hwt/@#&]]==1&]],{n,10}]

Formula

a(n) = A000041(n) - A321452(n).

Extensions

a(33)-a(50) from Alois P. Heinz, Nov 11 2018
Previous Showing 11-20 of 111 results. Next