(VBA)
sub calcul()
sheets("Result").select
range("A1").select
for i=1 to 10^13
sum=0
for k=1 to len(i)
sum=sum+(len(i)-k+1)^mid(i, k, 1)
next
if i=sum then
activecell.value=i
activesheet.offset(1, 0).select
end if
next
end sub
(PARI) isok(n) = my(d=digits(n)); sum(k=1, #d, (#d-k+1)^d[k]) == n; \\
Michel Marcus, Feb 24 2017