|
%I #11 Jan 20 2023 08:23:31
%S 4,6,8,9,10,12,14,16,18,20,21,22,24,25,26,28,30,32,33,34,35,36,38,39,
%T 40,42,44,45,46,48,49,50,52,54,55,56,57,58,60,62,64,65,66,68,69,70,72,
%U 74,75,76,77,78,80,81,82,84,86,87,88,90,91,92,93,94,96,98,99,100,102,104
%N Numbers k such that p <> (k AND p) for at least one prime-factor p.
%C Numbers k such that A102550(k) < A001221(k).
%C 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
%p isA102554 := proc(n)
%p local p;
%p for p in numtheory[factorset](n) do
%p if p <> ANDnos(p,n) then
%p return true
%p end if;
%p end do:
%p false ;
%p end proc:
%p for n from 1 to 500 do
%p if isA102554(n) then
%p printf("%d,",n) ;
%p end if;
%p end do: # _R. J. Mathar_, Jan 20 2023
%o (Python)
%o from itertools import count, islice
%o from functools import reduce
%o from operator import ior
%o from sympy import primefactors
%o def A102554_gen(startvalue=2): # generator of terms >= startvalue
%o return filter(lambda n:n|reduce(ior,primefactors(n))!=n,count(max(startvalue,2)))
%o A102554_list = list(islice(A102554_gen(),20)) # _Chai Wah Wu_, Dec 18 2022
%Y Cf. A001221, A102550, A102553, A007088, A004676.
%K nonn,base
%O 1,1
%A _Reinhard Zumkeller_, Jan 14 2005
|
