OFFSET
1,1
COMMENTS
FORMULA
a(n) = 2 * A301988(n).
EXAMPLE
The sequence of terms together with their prime indices begins:
18: {1,2,2}
60: {1,1,2,3}
168: {1,1,1,2,4}
216: {1,1,1,2,2,2}
400: {1,1,1,1,3,3}
528: {1,1,1,1,2,5}
1248: {1,1,1,1,1,2,6}
2240: {1,1,1,1,1,1,3,4}
2880: {1,1,1,1,1,1,2,2,3}
3264: {1,1,1,1,1,1,2,7}
7296: {1,1,1,1,1,1,1,2,8}
14080: {1,1,1,1,1,1,1,1,3,5}
17664: {1,1,1,1,1,1,1,1,2,9}
25088: {1,1,1,1,1,1,1,1,1,4,4}
32256: {1,1,1,1,1,1,1,1,1,2,2,4}
41472: {1,1,1,1,1,1,1,1,1,2,2,2,2}
44544: {1,1,1,1,1,1,1,1,1,2,10}
66560: {1,1,1,1,1,1,1,1,1,1,3,6}
95232: {1,1,1,1,1,1,1,1,1,1,2,11}
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[10000], And[PrimeOmega[#]>2, Times@@primeMS[#]==Total[primeMS[#]]-1]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 25 2019
EXTENSIONS
More terms from Jinyuan Wang, Jun 27 2020
STATUS
approved