OFFSET
1,2
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The terms together with their prime indices begin:
1: {}
2: {1}
4: {1,1}
6: {1,2}
8: {1,1,1}
10: {1,3}
12: {1,1,2}
16: {1,1,1,1}
18: {1,2,2}
20: {1,1,3}
22: {1,5}
24: {1,1,1,2}
30: {1,2,3}
32: {1,1,1,1,1}
34: {1,7}
36: {1,1,2,2}
40: {1,1,1,3}
42: {1,2,4}
44: {1,1,5}
48: {1,1,1,1,2}
MAPLE
pindices:= proc(n) uses NumberTheory; map(pi, PrimeFactors(n)) end proc:
filter:= proc(n) local F;
F:= pindices(n);
andmap(t -> NumberTheory:-Divisors(t) subset F, F)
end proc:
select(filter, [1, seq(i, i=2..200, 2)]); # Robert Israel, Dec 26 2025
MATHEMATICA
Select[Range[100], PrimeNu[#]==Length[Union @@ Divisors/@PrimePi/@First/@If[#==1, {}, FactorInteger[#]]]&]
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 18 2024
STATUS
approved
