我们要求1到p-1在%p意义下的逆元:
ny[1]=1; for(int i=2;i<=n;i++) { ny[i]=(p-p/i)*ny[p%i]%p; }
其中nv为逆元的意思。
时间复杂度O(p)