OFFSET
1,4
COMMENTS
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
FORMULA
a(2n) = A056239(2n) - 1 for n > 0.
EXAMPLE
The prime indices of 33 are {2,5}, with combinations
2 = 2
4 = 2+2
5 = 5
6 = 2+2+2
7 = 5+2
Hence a(33) = 5.
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
combs[n_, y_]:=With[{s=Table[{k, i}, {k, y}, {i, 0, Floor[n/k]}]}, Select[Tuples[s], Total[Times@@@#]==n&]];
Table[Length[Select[Range[Total[prix[n]]], combs[#, prix[n]]!={}&]], {n, 30}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Oct 19 2023
STATUS
approved