OFFSET
1,2
COMMENTS
Also Heinz numbers of integer partitions whose parts have (1) squarefree product and (2) binary indices covering an initial interval.
A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
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.
EXAMPLE
The terms together with their binary indices of prime indices begin:
1: {}
2: {{1}}
4: {{1},{1}}
5: {{1,2}}
6: {{1},{2}}
8: {{1},{1},{1}}
10: {{1},{1,2}}
12: {{1},{1},{2}}
15: {{2},{1,2}}
16: {{1},{1},{1},{1}}
17: {{1,2,3}}
20: {{1},{1},{1,2}}
24: {{1},{1},{1},{2}}
26: {{1},{2,3}}
30: {{1},{2},{1,2}}
32: {{1},{1},{1},{1},{1}}
33: {{2},{1,3}}
34: {{1},{1,2,3}}
40: {{1},{1},{1},{1,2}}
47: {{1,2,3,4}}
48: {{1},{1},{1},{1},{2}}
51: {{2},{1,2,3}}
MATHEMATICA
normQ[m_]:=Or[m=={}, Union[m]==Range[Max[m]]];
bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n, 2]], 1];
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[1000], SquareFreeQ[Times@@prix[#]]&&normQ[Join@@bpe/@prix[#]]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 31 2024
STATUS
approved