OFFSET
1,1
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]
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
PROG
(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
CROSSREFS
KEYWORD
nonn,base
AUTHOR
N. J. A. Sloane, Jun 22 2009
EXTENSIONS
More terms from Hagen von Eitzen and R. J. Mathar, Jun 23 2009
STATUS
approved