A368099
Triangle read by rows where T(n,k) is the number of non-isomorphic k-element sets of finite nonempty multisets with cardinalities summing to n, or strict multiset partitions of weight n and length k.
Original entry on oeis.org
1, 0, 1, 0, 2, 1, 0, 3, 4, 1, 0, 5, 12, 5, 1, 0, 7, 28, 22, 5, 1, 0, 11, 66, 83, 31, 5, 1, 0, 15, 134, 252, 147, 34, 5, 1, 0, 22, 280, 726, 620, 203, 35, 5, 1, 0, 30, 536, 1946, 2283, 1069, 235, 35, 5, 1, 0, 42, 1043, 4982, 7890, 5019, 1469, 248, 35, 5, 1
Offset: 0
Triangle begins:
1
0 1
0 2 1
0 3 4 1
0 5 12 5 1
0 7 28 22 5 1
0 11 66 83 31 5 1
0 15 134 252 147 34 5 1
0 22 280 726 620 203 35 5 1
0 30 536 1946 2283 1069 235 35 5 1
0 42 1043 4982 7890 5019 1469 248 35 5 1
...
Row n = 4 counts the following representatives:
. {{1,1,1,1}} {{1},{1,1,1}} {{1},{2},{1,1}} {{1},{2},{3},{4}}
{{1,1,1,2}} {{1},{1,1,2}} {{1},{2},{1,2}}
{{1,1,2,2}} {{1},{1,2,2}} {{1},{2},{1,3}}
{{1,1,2,3}} {{1},{1,2,3}} {{1},{2},{3,3}}
{{1,2,3,4}} {{1},{2,2,2}} {{1},{2},{3,4}}
{{1},{2,2,3}}
{{1},{2,3,4}}
{{1,1},{1,2}}
{{1,1},{2,2}}
{{1,1},{2,3}}
{{1,2},{1,3}}
{{1,2},{3,4}}
Counting connected components instead of edges gives
A321194.
For set multipartitions we have
A334550.
Cf.
A255903,
A296122,
A302545,
A306005,
A317532,
A317775,
A317794,
A317795,
A319560,
A368094,
A368095.
-
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mpm[n_]:=Join@@Table[Union[Sort[Sort /@ (#/.x_Integer:>s[[x]])]&/@sps[Range[n]]],{s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
Table[Length[Union[brute /@ Select[mpm[n],UnsameQ@@#&&Length[#]==k&]]], {n,0,5},{k,0,n}]
-
EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
K(q, t, k)={EulerT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}
G(n)={my(s=0); forpart(q=n, my(p=sum(t=1, n, y^t*subst(x*Ser(K(q, t, n\t))/t, x, x^t))); s+=permcount(q)*exp(p-subst(subst(p, x, x^2), y, y^2))); s/n!}
T(n)={[Vecrev(p) | p <- Vec(G(n))]}
{ my(A=T(10)); for(n=1, #A, print(A[n])) } \\ Andrew Howroyd, Jan 11 2024
A383310
Number of ways to choose a strict multiset partition of a factorization of n into factors > 1.
Original entry on oeis.org
1, 1, 1, 2, 1, 3, 1, 5, 2, 3, 1, 8, 1, 3, 3, 9, 1, 8, 1, 8, 3, 3, 1, 20, 2, 3, 5, 8, 1, 12, 1, 19, 3, 3, 3, 24, 1, 3, 3, 20, 1, 12, 1, 8, 8, 3, 1, 46, 2, 8, 3, 8, 1, 20, 3, 20, 3, 3, 1, 38, 1, 3, 8, 37, 3, 12, 1, 8, 3, 12, 1, 67, 1, 3, 8, 8, 3, 12, 1, 46, 9, 3
Offset: 1
The a(36) = 24 choices:
{{2,2,3,3}} {{2},{2,3,3}} {{2},{3},{2,3}}
{{2,2,9}} {{3},{2,2,3}} {{2},{3},{6}}
{{2,3,6}} {{2,2},{3,3}}
{{2,18}} {{2},{2,9}}
{{3,3,4}} {{9},{2,2}}
{{3,12}} {{2},{3,6}}
{{4,9}} {{3},{2,6}}
{{6,6}} {{6},{2,3}}
{{36}} {{2},{18}}
{{3},{3,4}}
{{4},{3,3}}
{{3},{12}}
{{4},{9}}
The case of a unique choice (positions of 1) is
A008578.
This is the strict case of
A050336.
For distinct strict blocks we have
A050345.
For strict blocks that are not necessarily distinct we have
A296119.
Twice-partitions of this type are counted by
A296122.
A050320 counts factorizations into squarefree numbers, distinct
A050326.
Cf.
A000009,
A005117,
A045782,
A050342,
A279785,
A293243,
A293511,
A302494,
A316439,
A358914,
A381992,
A382201.
-
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
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[Sum[Length[Select[mps[y],UnsameQ@@#&]],{y,facs[n]}],{n,30}]
A387116
Number of ways to choose a constant sequence of integer partitions, one of each prime index of n.
Original entry on oeis.org
1, 1, 2, 1, 3, 0, 5, 1, 2, 0, 7, 0, 11, 0, 0, 1, 15, 0, 22, 0, 0, 0, 30, 0, 3, 0, 2, 0, 42, 0, 56, 1, 0, 0, 0, 0, 77, 0, 0, 0, 101, 0, 135, 0, 0, 0, 176, 0, 5, 0, 0, 0, 231, 0, 0, 0, 0, 0, 297, 0, 385, 0, 0, 1, 0, 0, 490, 0, 0, 0, 627, 0, 792, 0, 0, 0, 0, 0
Offset: 1
The a(49) = 5 choices:
((4),(4))
((3,1),(3,1))
((2,2),(2,2))
((2,1,1),(2,1,1))
((1,1,1,1),(1,1,1,1))
Twice-partitions of this type are counted by
A047968, see also
A296122.
For initial intervals instead of partitions we have
A055396, see also
A387111.
For disjoint instead of constant we have
A383706.
For distinct instead of constant we have
A387110.
For strict partitions instead of partitions we have
A387117.
A003963 multiplies together prime indices.
A120383 lists numbers divisible by all of their prime indices.
A289509 lists numbers with relatively prime prime indices.
Cf.
A008284,
A063834,
A276078,
A335433,
A335448,
A355529,
A355741,
A387115,
A387120,
A387133,
A387136.
Comments