TJOI 2017 (Pollard_rho,MR素数测试)
题意是求分数 B A \frac{B}{A} AB在模M下的逆元,没有则输出-1。
本题的两大核心是:Pollard_rho算法,和MR素数测试算法。
B A \frac{B}{A} AB在模M下有逆元当且仅当 g c d ( A , M ) = 1 gcd(A,M)=1 gcd(A,M)=1,因为分数不是最简的,所以只要分数可以化简到 B ′ A ′ \frac{B'}{A'} A′B′,其中 g c d ( A ′ , M ) = 1 gcd(A',M)=1 gcd(A′,M)=1即有逆元。
我一开始的思路是先将 B A \frac{B}{A} AB化到最简分数 B ′ A ′ \frac{B'}{A'} A′B′,再判断A’是否有逆元。由于B’和A‘超出了264的范围,所以不能直接求 g c d gcd gcd以及约分。因为 A = a 1 ⋅ a 2 ⋯ a m , B = b 1 ⋅ b 2 ⋯ b m A=a_1\cdot a_2\cdots a_m,B=b_1\cdot b_2\cdots b_m A=a1⋅a2⋯am,B=b1⋅b2⋯bm,所以可以通过依次求ai,bi的质因子分解式,将A,B的质因子分解式求出,在通过A,B的质因子分解式去化简 B A \frac{B}{A} AB,再求M的质因子分解式,这样就可以判断是否 g c d ( A ′ , M ) = 1 gcd(A',M)=1 gcd(A′,M)=1了。
思路是简单明了的,答案能也保证正确,但是TLE。
在参考了别人的代码后,我发现并不用求A,B的质因子分解式,甚至也不用求M的质因子分解式,只需要求M的质因子集合即可。其实并不需要将 B A \frac{B}{A} AB化到最简,只需要将 B A \frac{B}{A}