with(numtheory):T:=proc(w) local x, y, z; x:=w; y:=0;
for z from 1 to ilog10(x)+1 do y:=10*y+(x mod 10); x:=trunc(x/10); od; y; end:
P:=proc(q) local j, k, n; for n from 1 to q do
k:=0; j:=n; while j mod 2 <> 1 do k:=k+1; j:=j/2; od;
if sigma(2*n+1)+sigma(2*n-1)+sigma(n/2^k)*2^(k+1)-6*n-2=T(n);
then print(n); fi; od; end: P(10^9);
|