BRCOCOLI的博客

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=3923 因为要取模，所以要乘法逆元 #include #include #include using namespace std; typedef long long LL; const int P=1000000007; int gcd(int a,int b) { if(b==0) retur

题目地址：http://poj.org/problem?id=1061 就是求（x-y)*t +（n-m)*k = -a 的式子中t的最小解，但注意k！=0 #include #include #include #include #include using namespace std; typedef long long LL; void gcd(LL a,LL b,LL &d,LL &x,

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=1576 p=9973 n=A%p; 求(A/B)%p; 因为( (A%p) /B )%p!=(A/B)%p 所以可以求出B的逆元NB 然后(n*NB)%p就是答案 因为B和p互质，所以可以用费马小定理 也即是B^(p-1)≡1 mod p，左右再同除B 也即是B^(p-2)≡1/B m

详细博客：逆元详解 ax ≡ 1 mod m 其中x就是a的逆元 x=a^-1 原式子也即是： a*a^-1 ≡ 1 mod m 求法： 1）根据扩展欧几里得算法算法求的 2）若a和m互素，根据费马小定理 a^(m-1)=1 mod m 所以 a^-1也即是a的逆元=a^(m-2)%m 有时候遇到一些题目： (a / b) % p != (a % p / b