OFFSET
1,1
COMMENTS
Numbers k such that the bitwise OR of k and all prime factors of k is not equal to k. - Chai Wah Wu, Dec 18 2022
MAPLE
isA102554 := proc(n)
local p;
for p in numtheory[factorset](n) do
if p <> ANDnos(p, n) then
return true
end if;
end do:
false ;
end proc:
for n from 1 to 500 do
if isA102554(n) then
printf("%d, ", n) ;
end if;
end do: # R. J. Mathar, Jan 20 2023
PROG
(Python)
from itertools import count, islice
from functools import reduce
from operator import ior
from sympy import primefactors
def A102554_gen(startvalue=2): # generator of terms >= startvalue
return filter(lambda n:n|reduce(ior, primefactors(n))!=n, count(max(startvalue, 2)))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Reinhard Zumkeller, Jan 14 2005
STATUS
approved