A341089 - OEIS

%I #16 Feb 19 2022 04:54:07

%S 1,2,3,4,5,7,8,9,11,13,16,17,23,31,32,37,61,71,73,79,97,101,107,113,

%T 121,125,128,131,149,151,157,163,167,169,179,181,191,199,311,313,337,

%U 343,347,353,359,361,373,383,389,521,701,709,727,733,739,743,751,757

%N Powers of primes whose reversal in base 10 is also a power of prime.

%C A subsequence of A000961 that contains A007500.

%e 16 = 2^4 and its reversal in base 10, 61, is prime, so 16 and 61 belong to this sequence.

%t Select[Range[800], AllTrue[{#, IntegerReverse[#]}, PrimeNu[#1] < 2 &] &] (* _Amiram Eldar_, Feb 15 2022 *)

%o (PARI) is(n) = omega(n)<=1 && omega(fromdigits(Vecrev(digits(n))))<=1

%Y Cf. A000961, A004086, A007500.

%K nonn,base

%O 1,2

%A _Rémy Sigrist_, Feb 13 2022