|
|
A289558
|
|
Numbers with two distinct prime factors not divisible by a square larger than 4.
|
|
1
|
|
|
6, 10, 12, 14, 15, 20, 21, 22, 24, 26, 28, 33, 34, 35, 38, 39, 40, 44, 46, 48, 51, 52, 55, 56, 57, 58, 62, 65, 68, 69, 74, 76, 77, 80, 82, 85, 86, 87, 88, 91, 92, 93, 94, 95, 96, 104, 106, 111, 112, 115, 116, 118, 119, 122, 123, 124, 129, 133, 134, 136, 141, 142, 143, 145, 146, 148, 152
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
EXAMPLE
|
10 is in this sequence because A001221(10) = 2 and 10 is not divisible by a square greater than 4.
|
|
MATHEMATICA
|
With[{nn = 10^4, r = 2}, Take[#, 67] &@ Union@ Flatten@ Table[Function[p, {p*2^Range[Log2[nn/p]], p DeleteCases[Prime@ Range[r, PrimePi[nn/p]], q_ /; q == p]}]@ Prime@ n, {n, r, PrimePi[nn/Prime[r]]}]] (* Michael De Vlieger, Jul 10 2017 *)
|
|
PROG
|
(PARI) list(lim)=my(v=List()); forprime(p=5, lim\3, forprime(q=3, min(lim\p, p-2), listput(v, p*q))); for(e=1, logint(lim\3, 2), forprime(p=3, lim>>e, listput(v, p<<e))); Set(v) \\ Charles R Greathouse IV, Aug 09 2017
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|