Avatar_small

乘法逆元:扩展gcd与欧拉phi函数

在计算一个很大的组合数modP时会用到乘法逆元。即把/a变成*(f(a)),其中f(a)为a在模P意义下的乘法逆元,即a*f(a) mod P=1 计算乘法逆元有两种方法,扩展gcd或基于欧拉公式的快速幂取模。 -----------------------------------------------------------------------------------------...

Avatar_small

Vim 的超级牛力--一次 debug 经历

Python 在数值计算方面,有一个非常好的功能,可以写一个 Fortran 程序,然后通过简单的 f2py 命令,便可变成一个 python 模块,以供在 python 里面调用。 在我用 f2py 把 Fortran 文件编译成 python 模块时,遇到一个奇怪的问题。当用其它的编辑器打开并修改  .f90 文件后,再使用 f2py 转到 python 模块,会发现模块里面的变量值没有...