Subfactorial FIXED by *yuko*

integer::a(99)=0,b(99)=0;a(99)=1;do i=1,51
do j=1,98;if(a(j)>0)exit;enddo;print'(99i0)',a(j:)
k=0;do j=99,1,-1;k=(a(j)+b(j))*i+k/10
a(j)=b(j);b(j)=mod(k,10);enddo;enddo;end

Note that non-ascii characters in the above source code will be escaped (such as \x9f).

download

return to the top page