OFFSET
1,1
COMMENTS
Let x(0)x(1)... x(q-1)x(q) denote the decimal expansion of a number n. The sequence lists the numbers n such that sigma(p) = sigma(s) where p is the prefix x(0)x(1)... x(q-1) and s is the suffix x(1)... x(q-1)x(q).
EXAMPLE
5376 is in the sequence because sigma(537) = sigma(376) = 720.
MAPLE
with(numtheory):
for n from 10 to 10000 do:
x:=convert(n, base, 10):n1:=nops(x):
s1:=sum('x[i]*10^(i-1) ', 'i'=1..n1-1):
s2:=(n-irem(n, 10))/10:
x1:=sigma(s1):x2:=sigma(s2):
if x1 = x2
then
printf(`%d, `, n):
else
fi:
od:
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Jun 27 2014
STATUS
approved