A063834
Twice partitioned numbers: the number of ways a number can be partitioned into not necessarily different parts and each part is again so partitioned.
Original entry on oeis.org
1, 1, 3, 6, 15, 28, 66, 122, 266, 503, 1027, 1913, 3874, 7099, 13799, 25501, 48508, 88295, 165942, 299649, 554545, 997281, 1817984, 3245430, 5875438, 10410768, 18635587, 32885735, 58399350, 102381103, 180634057, 314957425, 551857780, 958031826, 1667918758
Offset: 0
G.f. = 1 + x + 3*x^2 + 6*x^3 + 15*x^4 + 28*x^5 + 66*x^6 + 122*x^7 + 266*x^8 + ...
If n=6, a possible first partitioning is (3+3), resulting in the following second partitionings: ((3),(3)), ((3),(2+1)), ((3),(1+1+1)), ((2+1),(3)), ((2+1),(2+1)), ((2+1),(1+1+1)), ((1+1+1),(3)), ((1+1+1),(2+1)), ((1+1+1),(1+1+1)).
A001970 counts multiset partitions of integer partitions.
-
with(combinat):
b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
b(n, i-1)+`if`(i>n, 0, numbpart(i)*b(n-i, i)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..50); # Alois P. Heinz, Nov 26 2015
-
Table[Plus @@ Apply[Times, IntegerPartitions[i] /. i_Integer :> PartitionsP[i], 2], {i, 36}]
(* second program: *)
b[n_, i_] := b[n, i] = If[n==0 || i==1, 1, b[n, i-1] + If[i > n, 0, PartitionsP[i]*b[n-i, i]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Jan 20 2016, after Alois P. Heinz *)
-
{a(n) = if( n<0, 0, polcoeff( 1 / prod(k=1, n, 1 - numbpart(k) * x^k, 1 + x * O(x^n)), n))}; /* Michael Somos, Dec 19 2016 */
A381454
Number of multisets that can be obtained by choosing a strict integer partition of each prime index of n and taking the multiset union.
Original entry on oeis.org
1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 3, 1, 4, 2, 2, 1, 5, 1, 6, 2, 2, 3, 8, 1, 3, 4, 1, 2, 10, 2, 12, 1, 3, 5, 4, 1, 15, 6, 4, 2, 18, 2, 22, 3, 2, 8, 27, 1, 3, 3, 5, 4, 32, 1, 6, 2, 6, 10, 38, 2, 46, 12, 2, 1, 8, 3, 54, 5, 8, 4, 64, 1, 76, 15, 3, 6, 6, 4, 89, 2, 1
Offset: 1
The a(25) = 3 multisets are: {3,3}, {1,2,3}, {1,1,2,2}.
Multiset partitions of prime indices:
- For strict multiset partitions with distinct sums (
A321469) see
A381637.
A003963 gives product of prime indices.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
A358914 counts twice-partitions into distinct strict partitions.
Cf.
A000720,
A001222,
A002846,
A005117,
A066328,
A213242,
A213385,
A213427,
A293511,
A299200,
A299201,
A299202,
A300385,
A317142.
-
prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Union[Sort/@Join@@@Tuples[Select[IntegerPartitions[#],UnsameQ@@#&]&/@prix[n]]]],{n,100}]
A381633
Number of ways to partition the prime indices of n into sets with distinct sums.
Original entry on oeis.org
1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 1, 1, 2, 2, 0, 1, 1, 1, 1, 2, 2, 1, 0, 0, 2, 0, 1, 1, 4, 1, 0, 2, 2, 2, 1, 1, 2, 2, 0, 1, 5, 1, 1, 1, 2, 1, 0, 0, 1, 2, 1, 1, 0, 2, 0, 2, 2, 1, 3, 1, 2, 1, 0, 2, 5, 1, 1, 2, 4, 1, 0, 1, 2, 1, 1, 2, 5, 1, 0, 0, 2, 1, 4, 2, 2, 2
Offset: 1
The A050320(60) = 6 ways to partition {1,1,2,3} into sets are:
{{1},{1,2,3}}
{{1,2},{1,3}}
{{1},{1},{2,3}}
{{1},{2},{1,3}}
{{1},{3},{1,2}}
{{1},{1},{2},{3}}
Of these, only the following have distinct block-sums:
{{1},{1,2,3}}
{{1,2},{1,3}}
{{1},{2},{1,3}}
So a(60) = 3.
Taking block-sums (and sorting) gives
A381634.
A003963 gives product of prime indices.
A265947 counts refinement-ordered pairs of integer partitions.
-
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
sfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[sfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
Table[Length[Select[sfacs[n],UnsameQ@@hwt/@#&]],{n,100}]
A293243
Numbers that cannot be written as a product of distinct squarefree numbers.
Original entry on oeis.org
4, 8, 9, 16, 24, 25, 27, 32, 40, 48, 49, 54, 56, 64, 72, 80, 81, 88, 96, 104, 108, 112, 121, 125, 128, 135, 136, 144, 152, 160, 162, 169, 176, 184, 189, 192, 200, 208, 216, 224, 232, 240, 243, 248, 250, 256, 272, 288, 289, 296, 297, 304, 320, 324, 328, 336
Offset: 1
120 is not in the sequence because 120 = 2*6*10. 3600 is not in the sequence because 3600 = 2*6*10*30.
Multiset partitions of this type (set of sets) are counted by
A050342.
Twice-partitions of this type (set of sets) are counted by
A279785, see also
A358914.
Partitions of this type are counted by
A382078.
A050320 counts factorizations into squarefree numbers.
A050345 counts factorizations partitioned into into distinct sets.
-
N:= 1000: # to get all terms <= N
A:= Vector(N):
A[1]:= 1:
for n from 2 to N do
if numtheory:-issqrfree(n) then
S:= [$1..N/n]; T:= n*S; A[T]:= A[T]+A[S]
fi;
od:
select(t -> A[t]=0, [$1..N]); # Robert Israel, Oct 10 2017
-
nn=500;
sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
Select[Range[nn],Length[sqfacs[#]]===0&]
A381806
Numbers that cannot be written as a product of squarefree numbers with distinct sums of prime indices.
Original entry on oeis.org
4, 8, 9, 16, 24, 25, 27, 32, 40, 48, 49, 54, 56, 64, 72, 80, 81, 88, 96, 104, 108, 112, 121, 125, 128, 135, 136, 144, 152, 160, 162, 169, 176, 184, 189, 192, 200, 208, 216, 224, 232, 240, 243, 248, 250, 256, 272, 288, 289, 296, 297, 304, 320, 324, 328, 336
Offset: 1
There are 4 factorizations of 18000 into squarefree numbers:
(2*2*3*5*10*30)
(2*2*5*6*10*15)
(2*2*10*15*30)
(2*5*6*10*30)
but none of these has all distinct sums of prime indices, so 18000 is in the sequence.
Strongly normal multisets of this type are counted by
A292444.
For more on set multipartitions with distinct sums see
A279785,
A381718.
A003963 gives product of prime indices.
-
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
sqfics[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfics[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]]
Select[Range[nn],Length[Select[sqfics[#],UnsameQ@@hwt/@#&]]==0&]
A381992
Number of integer partitions of n that can be partitioned into sets with distinct sums.
Original entry on oeis.org
1, 1, 1, 2, 3, 5, 6, 9, 13, 17, 25, 33, 44, 59, 77, 100, 134, 170, 217, 282, 360, 449, 571, 719, 899, 1122, 1391, 1727, 2136, 2616, 3209, 3947, 4800, 5845, 7094, 8602, 10408, 12533, 15062, 18107, 21686, 25956, 30967, 36936, 43897, 52132, 61850, 73157, 86466, 101992, 120195
Offset: 0
There are 6 ways to partition (3,2,2,1) into sets:
{{2},{1,2,3}}
{{1,2},{2,3}}
{{1},{2},{2,3}}
{{2},{2},{1,3}}
{{2},{3},{1,2}}
{{1},{2},{2},{3}}
Of these, 3 have distinct block sums:
{{2},{1,2,3}}
{{1,2},{2,3}}
{{1},{2},{2,3}}
so (3,2,2,1) is counted under a(8).
The a(1) = 1 through a(8) = 13 partitions:
(1) (2) (3) (4) (5) (6) (7) (8)
(2,1) (3,1) (3,2) (4,2) (4,3) (5,3)
(2,1,1) (4,1) (5,1) (5,2) (6,2)
(2,2,1) (3,2,1) (6,1) (7,1)
(3,1,1) (4,1,1) (3,2,2) (3,3,2)
(2,2,1,1) (3,3,1) (4,2,2)
(4,2,1) (4,3,1)
(5,1,1) (5,2,1)
(3,2,1,1) (6,1,1)
(3,2,2,1)
(3,3,1,1)
(4,2,1,1)
(3,2,1,1,1)
Twice-partitions of this type are counted by
A279785.
Multiset partitions of this type are counted by
A381633, zeros of
A381634.
Normal multiset partitions of this type are counted by
A381718, see
A116539.
These partitions are ranked by
A382075.
For distinct blocks instead of sums we have
A382077, complement
A382078.
For a unique choice we have
A382079.
A050320 counts multiset partitions of prime indices into sets.
A050326 counts multiset partitions of prime indices into distinct sets.
A265947 counts refinement-ordered pairs of integer partitions.
A382201 lists MM-numbers of sets with distinct sums.
-
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
Table[Length[Select[IntegerPartitions[n],Length[Select[mps[#], And@@UnsameQ@@@#&&UnsameQ@@Total/@#&]]>0&]],{n,0,10}]
A381990
Number of integer partitions of n that cannot be partitioned into a set (or multiset) of sets with distinct sums.
Original entry on oeis.org
0, 0, 1, 1, 2, 2, 5, 6, 9, 13, 17, 23, 33, 42, 58, 76, 97, 127, 168, 208, 267, 343, 431, 536, 676, 836, 1045, 1283, 1582, 1949, 2395, 2895, 3549, 4298, 5216, 6281, 7569, 9104, 10953, 13078, 15652, 18627, 22207, 26325, 31278, 37002, 43708, 51597, 60807, 71533, 84031
Offset: 0
The partition y = (3,3,3,2,2,1,1,1,1) has only one multiset partition into a set of sets, namely {{1},{3},{1,2},{1,3},{1,2,3}}, but this does not have distinct sums, so y is counted under a(17).
The a(2) = 1 through a(8) = 9 partitions:
(11) (111) (22) (2111) (33) (2221) (44)
(1111) (11111) (222) (4111) (2222)
(3111) (22111) (5111)
(21111) (31111) (22211)
(111111) (211111) (41111)
(1111111) (221111)
(311111)
(2111111)
(11111111)
Twice-partitions of this type are counted by
A279785.
Normal multiset partitions of this type are counted by
A381718, see
A116539.
MM-numbers of these multiset partitions (strict blocks with distinct sum) are
A382201.
A050320 counts multiset partitions of prime indices into sets.
A050326 counts multiset partitions of prime indices into distinct sets.
A265947 counts refinement-ordered pairs of integer partitions.
-
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
Table[Length[Select[IntegerPartitions[n],Length[Select[mps[#],And@@UnsameQ@@@#&&UnsameQ@@Total/@#&]]==0&]],{n,0,10}]
A382077
Number of integer partitions of n that can be partitioned into a set of sets.
Original entry on oeis.org
1, 1, 1, 2, 3, 5, 6, 9, 13, 17, 25, 33, 44, 59, 77, 100, 134, 171, 217, 283, 361, 449, 574, 721, 900, 1126, 1397, 1731, 2143, 2632, 3223, 3961, 4825, 5874, 7131, 8646, 10452, 12604, 15155, 18216, 21826, 26108, 31169, 37156, 44202, 52492, 62233, 73676, 87089, 102756, 121074
Offset: 0
For y = (3,2,2,2,1,1,1), we have the multiset partition {{1},{2},{1,2},{1,2,3}}, so y is counted under a(12).
The a(1) = 1 through a(8) = 13 partitions:
(1) (2) (3) (4) (5) (6) (7) (8)
(2,1) (3,1) (3,2) (4,2) (4,3) (5,3)
(2,1,1) (4,1) (5,1) (5,2) (6,2)
(2,2,1) (3,2,1) (6,1) (7,1)
(3,1,1) (4,1,1) (3,2,2) (3,3,2)
(2,2,1,1) (3,3,1) (4,2,2)
(4,2,1) (4,3,1)
(5,1,1) (5,2,1)
(3,2,1,1) (6,1,1)
(3,2,2,1)
(3,3,1,1)
(4,2,1,1)
(3,2,1,1,1)
Factorizations of this type are counted by
A050345.
Normal multiset partitions of this type are counted by
A116539.
The MM-numbers of these multiset partitions are
A302494.
Twice-partitions of this type are counted by
A358914.
For distinct block-sums instead of blocks we have
A381992, ranked by
A382075.
For normal multisets instead of integer partitions we have
A382214, complement
A292432.
A050320 counts multiset partitions of prime indices into sets.
A050326 counts multiset partitions of prime indices into distinct sets.
A265947 counts refinement-ordered pairs of integer partitions.
-
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]& /@ sps[Range[Length[set]]]];
Table[Length[Select[IntegerPartitions[n], Length[Select[mps[#],UnsameQ@@#&&And@@UnsameQ@@@#&]]>0&]],{n,0,9}]
A381634
Number of multisets that can be obtained by taking the sum of each block of a set multipartition (multiset of sets) of the prime indices of n with distinct block-sums.
Original entry on oeis.org
1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 1, 1, 2, 2, 0, 1, 1, 1, 1, 2, 2, 1, 0, 0, 2, 0, 1, 1, 4, 1, 0, 2, 2, 2, 1, 1, 2, 2, 0, 1, 5, 1, 1, 1, 2, 1, 0, 0, 1, 2, 1, 1, 0, 2, 0, 2, 2, 1, 3, 1, 2, 1, 0, 2, 5, 1, 1, 2, 4, 1, 0, 1, 2, 1, 1, 2, 5, 1, 0, 0, 2, 1, 4, 2, 2, 2
Offset: 1
The prime indices of 120 are {1,1,2,3}, with 3 ways:
{{1},{1,2,3}}
{{1,2},{1,3}}
{{1},{2},{1,3}}
with block-sums: {1,6}, {3,4}, {1,2,4}, so a(120) = 3.
The prime indices of 210 are {1,2,3,4}, with 12 ways:
{{1,2,3,4}}
{{1},{2,3,4}}
{{2},{1,3,4}}
{{3},{1,2,4}}
{{4},{1,2,3}}
{{1,2},{3,4}}
{{1,3},{2,4}}
{{1},{2},{3,4}}
{{1},{3},{2,4}}
{{1},{4},{2,3}}
{{2},{3},{1,4}}
{{1},{2},{3},{4}}
with block-sums: {10}, {1,9}, {2,8}, {3,7}, {4,6}, {3,7}, {4,6}, {1,2,7}, {1,3,6}, {1,4,5}, {2,3,5}, {1,2,3,4}, of which 10 are distinct, so a(210) = 10.
A003963 gives product of prime indices.
A265947 counts refinement-ordered pairs of integer partitions.
Cf.
A000720,
A001222,
A002846,
A005117,
A116540,
A213242,
A213385,
A213427,
A299202,
A300385,
A317142,
A317143,
A318360.
-
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
sfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[sfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
Table[Length[Union[Sort[hwt/@#]&/@Select[sfacs[n],UnsameQ@@hwt/@#&]]],{n,100}]
A382078
Number of integer partitions of n that cannot be partitioned into a set of sets.
Original entry on oeis.org
0, 0, 1, 1, 2, 2, 5, 6, 9, 13, 17, 23, 33, 42, 58, 76, 97, 126, 168, 207, 266, 343, 428, 534, 675, 832, 1039, 1279, 1575, 1933, 2381, 2881, 3524, 4269, 5179, 6237, 7525, 9033, 10860, 12969, 15512, 18475, 22005, 26105, 30973, 36642, 43325, 51078, 60184, 70769, 83152
Offset: 0
The partition y = (2,2,1,1,1) can be partitioned into sets in the following ways:
{{1},{1,2},{1,2}}
{{1},{1},{2},{1,2}}
{{1},{1},{1},{2},{2}}
But none of these is itself a set, so y is counted under a(7).
The a(2) = 1 through a(8) = 9 partitions:
(11) (111) (22) (2111) (33) (2221) (44)
(1111) (11111) (222) (4111) (2222)
(3111) (22111) (5111)
(21111) (31111) (22211)
(111111) (211111) (41111)
(1111111) (221111)
(311111)
(2111111)
(11111111)
The MM-numbers of these multiset partitions (set of sets) are
A302494.
Twice-partitions of this type are counted by
A358914.
A050320 counts multiset partitions of prime indices into sets.
A050326 counts multiset partitions into distinct sets, complement
A050345.
A265947 counts refinement-ordered pairs of integer partitions.
-
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
Table[Length[Select[IntegerPartitions[n],Length[Select[mps[#],UnsameQ@@#&&And@@UnsameQ@@@#&]]==0&]],{n,0,9}]
Showing 1-10 of 35 results.
Comments