OFFSET
1,2
LINKS
EXAMPLE
The sequence of terms together with their prime indices begins:
1: {}
667: {9,10}
2021: {14,15}
4331: {18,20}
5767: {21,22}
7081: {21,25}
7663: {22,25}
11021: {27,28}
11639: {27,30}
12091: {28,30}
13837: {26,33}
15049: {26,35}
20413: {33,35}
23213: {34,39}
26123: {36,40}
31553: {34,49}
34933: {42,44}
35657: {42,45}
35723: {34,55}
37909: {39,49}
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[10000], !PrimeQ[#]&&SquareFreeQ[#]&&SameQ@@PrimeOmega/@primeMS[#]&&SameQ@@Total/@primeMS/@primeMS[#]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 30 2019
STATUS
approved