OFFSET
1,1
COMMENTS
a(32) > 10^10. a(33) = 1167555543. a(34) > 10^10. - Donovan Johnson, May 08 2010
a(35) > 3*10^11. - Oscar Volpatti, Aug 06 2020
MAPLE
A048853 := proc(n::integer) local resul, ddigs, d, c, tmp ; resul := 0 ; ddigs := convert(n, base, 10) ; for d from 1 to nops(ddigs) do for c from 0 to 9 do if c = 0 and d = nops(ddigs) then continue ; else if c <> op(d, ddigs) then tmp := [op(1..d-1, ddigs), c, op(d+1..nops(ddigs), ddigs)] ; tst := sum(op(i, tmp)*10^(i-1), i=1..nops(tmp)) ; if isprime(tst) then resul := resul+1 ; fi ; fi ; fi ; od : od ; RETURN(resul) ; end: A050662 := proc(n::integer) local i; for i from 1 to 10000000 do if A048853(i) = n then RETURN(i) ; fi ; od ; RETURN(-1) ; end : for n from 1 to 90 do print(A050662(n)) ; od ; # R. J. Mathar, Apr 25 2006
CROSSREFS
KEYWORD
nonn,base,more
AUTHOR
Patrick De Geest, Jul 15 1999
EXTENSIONS
a(25)-a(31) and a(33) from Donovan Johnson, May 08 2010
a(32) and a(34) from Oscar Volpatti, Aug 06 2020
STATUS
approved