OFFSET
1,1
COMMENTS
EXAMPLE
2063 is in the sequence because 1/2 + 1/6 + 1/3 = 1.
MAPLE
with(numtheory):nn:=500000:for m from 1 to nn do:n:=ithprime(m):y:=convert(n, base, 10):n2:=nops(y):s:=0:for i from 1 to n2 do: if y[i]<>0 then s:=s+1/y[i]:else fi:od:if s=1 then printf(`%d, `, n):else fi:od:
MATHEMATICA
Select[Prime[Range[42000]], Total[1/Select[IntegerDigits[#], #!=0&]]==1&] (* Harvey P. Dale, May 31 2019 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Mar 24 2014
STATUS
approved