|
PROG
|
(Python)
from sympy import primefactors, prime
A242092 = [n for n in range(1, 10**7) if primefactors(n) == primefactors(int(str(prime(n))[::-1]))]
(PARI)
rev(n)=r=""; d=digits(n); for(i=1, #d, r=concat(Str(d[i]), r)); eval(r)
for(n=1, 10^7, p=rev(prime(n)); if(omega(n)==omega(p), if(gcd(n, p)==min(n, p), print1(n, ", ")))) \\ Derek Orr, Aug 14 2014
|