character*28c,d;1 read*,c;read(c,*,err=2)n;print'(i0" ="$)',n;do i=2,21;k=n**(1./i);if(n==k**i)print'(xi0"^"i0$)',k,i enddo;print*;goto1;2 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).