login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A325223
Sum of the prime indices of n minus the greater of the number of prime factors of n counted with multiplicity and the largest prime index of n.
3
0, 0, 0, 0, 0, 1, 0, 0, 2, 1, 0, 1, 0, 1, 2, 0, 0, 2, 0, 2, 2, 1, 0, 1, 3, 1, 3, 2, 0, 3, 0, 0, 2, 1, 3, 2, 0, 1, 2, 2, 0, 3, 0, 2, 4, 1, 0, 1, 4, 4, 2, 2, 0, 3, 3, 3, 2, 1, 0, 3, 0, 1, 4, 0, 3, 3, 0, 2, 2, 4, 0, 2, 0, 1, 5, 2, 4, 3, 0, 2, 4, 1, 0, 4, 3, 1, 2
OFFSET
1,9
COMMENTS
A prime index of n is a number m such that prime(m) divides n.
Also the number of squares in the Young diagram of the integer partition with Heinz number n after the first row or the first column, whichever is larger, is removed. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
FORMULA
a(n) = A056239(n) - max(A001222(n), A061395(n)) = A056239(n) - A263297(n).
EXAMPLE
88 has 4 prime indices {1,1,1,5} with sum 8 and maximum 5, so a(88) = 8 - max(4,5) = 3.
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[Total[primeMS[n]]-Max[Length[primeMS[n]], Max[primeMS[n]]], {n, 100}]
CROSSREFS
Positions of 0's are A174090. Positions of 1's are A325231.
Sequence in context: A307776 A341027 A050317 * A141095 A175599 A179759
KEYWORD
nonn
AUTHOR
Gus Wiseman, Apr 12 2019
STATUS
approved