character*28c,d;do1;read'(a)',c;read(c,*,err=3)n;print'(i0" ="$)',n do2 i=2,21;k=n**(1./i);2 if(n==k**i)print'(xi0"^"i0$)',k,i 1 print*;3 read*,d;print'(a,"= 10667507560"a)',c,"001^2",d,"159^2";end
Note that non-ascii characters in the above source code will be escaped (such as \x9f).