OFFSET
1,1
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..1000
EXAMPLE
91 - phi(91) = 91 - 72 = 19;
342 - phi(342) = 342 - 108 = 234.
MAPLE
with(numtheory): P:=proc(q) local a, b, c, d, k, n; for n from 1 to q do a:=n; b:=n-phi(n);
if ilog10(a)=ilog10(b) then c:=[]; d:=[]; for k from 1 to ilog10(n)+1 do
c:=[op(c), (a mod 10)]; a:=trunc(a/10); d:=[op(d), (b mod 10)]; b:=trunc(b/10); od;
c:=sort(c); d:=sort(d); if c=d then print(n); fi; fi; od; end: P(10^25);
MATHEMATICA
Select[Range[25000], Sort@ IntegerDigits@ # == Sort@ IntegerDigits[# - EulerPhi@ #] &] (* Michael De Vlieger, Jun 01 2016 *)
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Paolo P. Lava, May 31 2016
STATUS
approved
