OFFSET
1,1
LINKS
Lars Blomberg, Table of n, a(n) for n = 1..47
EXAMPLE
phi(22) = 10: 2 + 2 = 4; 2 + 4 = 6; 4 + 6 = 10.
MAPLE
with(numtheory): P:=proc(q, h) local a, b, k, n, t, v; v:=array(1..h);
for n from 2 to q do a:=n; b:=ilog10(a)+1; if b>1 then
for k from 1 to b do v[b-k+1]:=(a mod 10); a:=trunc(a/10); od; t:=b+1; v[t]:=add(v[k], k=1..b);
while v[t]<phi(n) do t:=t+1; v[t]:=add(v[k], k=t-b..t-1); od;
if v[t]=phi(n) then print(n); fi; fi; od; end: P(10^6, 1000);
MATHEMATICA
Select[Range[2, 10^5], (t = EulerPhi[#]; d = IntegerDigits[#]; While[Total[d] < t, d = Join[Rest[d], {Total[d]}]]; Total[d] == t) &] (* Robert Price, May 17 2019 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paolo P. Lava, Feb 24 2016
EXTENSIONS
a(38) from Lars Blomberg, Jan 18 2018
STATUS
approved