|
|
A320836
|
|
a(n) = Sum (-1)^k where the sum is over all strict multiset partitions of a multiset whose multiplicities are the prime indices of n and k is the number of parts, or strict factorizations of A181821(n).
|
|
3
|
|
|
1, -1, -1, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, -1, -1, 0, -2, -1, 0, -2, 0, -2, -1, -1, -1, -4, -1, -1, -1, -3, 0, -3, 0, -2, -4, -1, -1, -6, -2, -3, -2, -2, 0, -6, -2, -4, -1, -1, 0, -5, 0, -1, -3, -9, -2, -3, 0, -2, -1, -3, 0, -7, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,27
|
|
COMMENTS
|
This multiset (row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.
|
|
LINKS
|
|
|
FORMULA
|
|
|
MAPLE
|
with(numtheory):
b:= proc(n, k) option remember; `if`(n>k, 0, -1)+`if`(isprime(n), 0,
-add(`if`(d>k, 0, b(n/d, d-1)), d=divisors(n) minus {1, n}))
end:
a:= n-> `if`(n=1, 1, b(((l-> mul(ithprime(i)^l[i], i=1..nops(l)))(
sort(map(i-> pi(i[1])$i[2], ifactors(n)[2]), `>`)))$2)):
|
|
MATHEMATICA
|
nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]], {#1}]&, If[n==1, {}, Flatten[Cases[FactorInteger[n]//Reverse, {p_, k_}:>Table[PrimePi[p], {k}]]]]];
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[(-1)^Length[m], {m, Select[mps[nrmptn[n]], UnsameQ@@#&]}], {n, 30}]
|
|
CROSSREFS
|
Cf. A001055, A001222, A007716, A045778, A114592, A162247, A181821, A305936, A316441, A318284, A319237, A319238, A320835.
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|