logmod:计算 y 使得 mod(a^y, p^N) == x-matlab开发

认为p 是奇素数， ord(a) = p-1 GCD(x, p) = 1 N 是一个整数 >= 2 然后y = logmod(x, a, p, N) 返回 y 使得 a^y = x (mod p^N) 例子设 x= vpi(154), a = vpi(7), p = vpi(17), N = vpi(37) y = logmod(x, a, p, N) 返回y = 2088349219044680767324467844670001776975183904