A151768 Complement of A071786.

19, 23, 29, 41, 43, 46, 47, 53, 57, 58, 59, 61, 67, 69, 82, 83, 86, 87, 89, 94, 103, 109, 115, 116, 122, 123, 127, 129, 137, 138, 139, 141, 159, 161, 163, 171, 173, 174, 177, 178, 183, 193, 197, 201, 203, 205, 206, 207, 209, 211, 215, 218, 223, 227, 229, 230

Offset: 1

Author

N. J. A. Sloane, Jun 22 2009

Keywords

• nonn
• base

Comments

A number n > 1 is in the sequence if and only if n/A004086(p) is in the sequence for all primes p with A004086(p) dividing n. [Hagen von Eitzen, Jun 23 2009]

A076056 is a subsequence. [Reinhard Zumkeller, Jul 06 2009]

Programs

Maple

read("transforms") ; tdpr := proc(nd) local p,n ; p := [] ; for n from 1 do if ithprime(n) > 10^nd then break; else p := [op(p),digrev(ithprime(n))] ; fi; od: sort(p) ; end: A071786 := proc(L,nd) local tmp,tmp2,j,k,i ; tmp := [] ; for j from 0 do if op(1,L)^j > 10^nd then break; fi; tmp := [op(tmp),op(1,L)^j] ; od: for i from 2 to nops(L) do tmp2 := {} ; for k from 1 to nops(tmp) do for j from 0 do if op(k,tmp)*op(i,L)^j > 10^nd then break; fi; tmp2 := tmp2 union { op(k,tmp)*op(i,L)^j} ; od: od: tmp := convert(tmp2,list) ; od: tmp ; end: maxp10 := 3 : L := tdpr(maxp10) : a151768c := A071786(L,maxp10) : for n from 1 to 10^maxp10 do if not n in a151768c then printf("%d,",n) ; fi; od: # R. J. Mathar, Jun 24 2009

PARI

A=Set([]);for(n=2,300,ok=0;fordiv(n,d,if(!setsearch(A,d)&&isprime(rev(n/d)),ok=1;break));if(!ok,print1(n,",");A=setunion(A,Set([n])))) \\ Hagen von Eitzen, Jun 23 2009

Extensions

More terms from Hagen von Eitzen and R. J. Mathar, Jun 23 2009