A372121 Row sums of A371783 and A371954 (k-quanimous partitions).
1, 3, 4, 9, 8, 22, 16, 42, 41, 74, 57, 183, 102, 233, 263, 463, 298, 875, 491, 1350, 1172, 1775, 1256, 4273, 2225, 4399, 4584, 8049, 4566, 14913, 6843, 18539, 15831, 22894, 18196, 53323, 21638, 48947, 50281, 94500, 44584, 144976, 63262, 173436, 169361, 202153
Offset: 1
Keywords
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[Sum[Length[Select[IntegerPartitions[n], Select[facs[Times@@Prime/@#], Length[#]==k&&SameQ@@hwt/@#&]!={}&]],{k,Divisors[n]}],{n,1,10}]
-
PARI
T(n, d) = my(v=partitions(n/d), w=List([])); forvec(s=vector(d, i, [1, #v]), listput(w, vecsort(concat(vector(d, i, v[s[i]])))), 1); #Set(w); a(n) = sumdiv(n, d, T(n, d)); \\ Jinyuan Wang, Feb 13 2025
Extensions
More terms from Jinyuan Wang, Feb 13 2025
Comments