OFFSET
1,1
COMMENTS
LINKS
Paul Weisenhorn, Table of n, a(n) for n = 1..70
EXAMPLE
R(3)=111 multiplied by a(3)=6 yields z=666=2*3*3*37 = A006753(34): 6+6+6 = 2+3+3+3+7 = 18.
MAPLE
# digits-sum of primfactors of z=dsp(z)
for n from 2 to 70 do f(n):=1: test:=false:
while (f(n) < 420000) and (test=false) do
f(n):=f(n)+1: z:=f(n)*r(n): ds(z):=0:
dsp(z):=dsp(r(n))+dsp(f(n)):
while (z>0) do z:=iquo(z, 10, 'm'): ds(z):=ds(z)+m: end do:
if(ds(z)=dsp(z)) then test:=true: print(n, f(n)): end if:
end do: end do:
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paul Weisenhorn, Apr 16 2010, Apr 23 2010
EXTENSIONS
Keyword:base added by R. J. Mathar, Apr 24 2010
STATUS
approved