OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
4 is not in the sequence since 4 is a prime power, and all divisors d | k of prime power k = p^e are also prime powers.
6 is not in the sequence since 6 is squarefree, and all divisors d | k of squarefree k are also squarefree.
12 is not in the sequence since 12 is in A366825, and there is only 1 divisor in A126706, which is 12 itself.
24 is in the sequence since the intersection of A126706 and row 24 of A027750, indicated by bracketed numbers, is {1, 2, 3, 4, 6, [12, 24]}, etc.
24: {12, 24}
36: {12, 18, 36}
40: {20, 40}
48: {12, 24, 48}
54: {18, 54}
56: {28, 56}
60: {12, 20, 60}
72: {12, 18, 24, 36, 72}
80: {20, 40, 80}
84: {12, 28, 84}
88: {44, 88}
90: {18, 45, 90}
MATHEMATICA
Select[Range[300], Function[k, DivisorSum[k, 1 &, Nor[PrimePowerQ[#], SquareFreeQ[#]] &] > 1]]
PROG
(PARI) list(lim)=my(v=List()); forfactored(k=24, lim\1, my(e=k[2][, 2]); if(#e>1 && vecmax(e)>1 && (#e>2 || vecsum(e)>3), listput(v, k[1]))); Vec(v) \\ Charles R Greathouse IV, Oct 01 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Sep 28 2024
STATUS
approved