OFFSET
1,1
COMMENTS
a(21) > 10^11. - Hiroaki Yamanouchi, Sep 28 2015
EXAMPLE
sigma*(10) = 3+4+7 = 14.
sigma(14)-14 = 1+2+7 = 10.
sigma*(32)= 3+5+7+9+13+21 = 58.
sigma(58)-58 = 1+2+29 = 32.
MAPLE
with(numtheory);
P:= proc(n)
local a, b, c, i, ks;
for i from 3 to n do
a:={};
for k from 2 to i-1 do
if abs((i mod k)- k/2) < 1 then
a:=a union {k};
fi;
od;
b:=nops(a); c:=op(a); s:=0;
for k from 1 to b do
s:=s+c[k];
od;
if sigma(s)-s=i then
print(i, s);
fi;
od;
end:
P(10000);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, Jun 29 2011
EXTENSIONS
a(11)-a(20) from Hiroaki Yamanouchi, Sep 28 2015
STATUS
approved