OFFSET
1,1
COMMENTS
A051674 is a subsequence of this sequence.
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..5000
EXAMPLE
18036' = 36180, so 18036 is a term.
20852' = 22508, so 20852 is a term.
48872' = 74828, so 48872 is a term.
MAPLE
with(numtheory); A225902:= proc(i) local a, b, c, d, j, n, ok, p, pfs;
for n from 1 to i do b:=[]; c:=[];
a:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
while a>0 do b:=[op(b), a mod 10]; a:=trunc(a/10); od; a:=n;
while a>0 do c:=[op(c), a mod 10]; a:=trunc(a/10); od; d:=0;
if nops(b)=nops(c) then b:=sort(b); c:=sort(c); b:=b-c; ok:=1; for j from 1 to nops(b) do if b[j]<>0 then ok:=0; break; fi; od; if ok=1 then print(n); fi; fi;
od; end; A225902(10000);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paolo P. Lava, May 20 2013
STATUS
approved