Modular Exponentiation by seichung

B=1e5;y=real(k<-200);y[1]=2;for(b in 1:2){for(n in 0:99){x=z=y;for(i in 1:b^2){r=0;for(m in 1:k){r=(p=y[1:m]%*%x[m:1]+r)%/%B;z[m]=p-r*B};y=z}}};cat(substr(B+y[k:1],2,6),sep="")

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

download

return to the top page