OFFSET
1,1
COMMENTS
Property of the sequence :
{a(n)} = E1 union E2 where E1 = {38, 58, 66, 87, 205, 714} and E2 = {110, 210, 310, 410, 510, 610, 710, 810, 910, 1010, 2010,...}.
E2 = F1 union F2 union... union Fk union... where :
F1 = {110, 210,..., 910} with one zero;
F2 = {1010, 2010,..., 9010} with two zeros;
F3 = {10010, 20010,..., 90010} with three zeros;
......................................................................
Fk = {100…0010, 200…0010,..., 900...0010} with k zeros;
......................................................................
Hence the proposition :
If n is of the form n = a0000...010 with k zeros and a =1,2,..., 9 then sigma(n) = sigma(x)*sigma(y) where x = a0000...01 with k-1 zeros and y = 10.
EXAMPLE
3010 is in the sequence because sigma(3010) = 6336; sigma(301) = 352 and sigma(10) = 18 => 6336 = 352*18.
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):xn:=sigma(n):
if xn = x1*x2
then
printf(`%d, `, n):
else
fi:
od:
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Jun 25 2014
STATUS
approved