cp's OEIS Frontend

A046349 Composite numbers with only palindromic prime factors.

%I A046349 #22 May 15 2022 16:24:34
%S A046349 4,6,8,9,10,12,14,15,16,18,20,21,22,24,25,27,28,30,32,33,35,36,40,42,
%T A046349 44,45,48,49,50,54,55,56,60,63,64,66,70,72,75,77,80,81,84,88,90,96,98,
%U A046349 99,100,105,108,110,112,120,121,125,126,128,132,135,140,144,147,150
%N A046349 Composite numbers with only palindromic prime factors.
%H A046349 Robert Israel, <a href="/A046349/b046349.txt">Table of n, a(n) for n = 1..3020</a>
%F A046349 A033620 INTERSECT A002808. - _R. J. Mathar_, Sep 09 2015
%p A046349 isA046349 := proc(n)
%p A046349     simplify(isA033620(n) and not isprime(n)) ;
%p A046349 end proc:
%p A046349 for n from 2 to 300 do
%p A046349     if isA046349(n) then
%p A046349         printf("%d,",n) ;
%p A046349     end if;
%p A046349 end do: # _R. J. Mathar_, Sep 09 2015
%t A046349 palQ[n_]:=Reverse[x=IntegerDigits[n]]==x; Select[Range[4,150],!PrimeQ[#]&&And@@palQ/@First/@FactorInteger[#]&] (* _Jayanta Basu_, Jun 05 2013 *)
%t A046349 Select[Range[200],CompositeQ[#]&&AllTrue[FactorInteger[#][[All,1]],PalindromeQ]&] (* _Harvey P. Dale_, May 15 2022 *)
%o A046349 (Python)
%o A046349 from sympy import isprime, primefactors
%o A046349 def pal(n): s = str(n); return s == s[::-1]
%o A046349 def ok(n): return not isprime(n) and all(pal(f) for f in primefactors(n))
%o A046349 print(list(filter(ok, range(4, 151)))) # _Michael S. Branicky_, Apr 06 2021
%Y A046349 Cf. A002385, A033620, A046350, A046351.
%K A046349 nonn,base,easy
%O A046349 1,1
%A A046349 _Patrick De Geest_, Jun 15 1998