real::a(51)=[0,12,20,40,60,100,120,140,143,5141,5144,5146,5156,5164,25564,61151,61163,61557,61564,& 62012,62412,62555,63412,64542,64544,66012,66145,66157,67151,67155,67165,67542,67546,67550,67557,70155,& 70157,71012,71141,71160,71163,71412,71557,71571,72012,72157,72412,73012,73145,73162,74165] do3 i=1,51;k=a(i);print'(i0h:$)',k;n=1;1 n=n+1;do2 while(k>1);if(k/n*n<k)goto1;print'(xi0$)',n;2 k=k/n;3 print*;end
Note that non-ascii characters in the above source code will be escaped (such as \x9f).