OFFSET
1,2
PROG
(PARI) sd(n) = digs = digits(n); sum(i=1, #digs, digs[i]);
isok(n) = sd(n) == sd(sigma(n)); \\ Michel Marcus, Dec 26 2013
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Labos Elemer, Jan 16 2002
STATUS
approved