若ax≡1 mod f, 则称a关于模f的乘法逆元为x。也可表示为ax≡1(mod f)。 当a与f互素时，a关于模f的乘法逆元有唯一解。如果不互素，则无解。如果f为素数，则从1到f-1的任意数都与f互素，即在1到f-1之间都恰好有一个关于模f的乘法逆元。 （不会证明，想通了补） 首先a与f要互素， ...
分类：
其他好文 时间：
2016-04-03 20:16:24
阅读次数：
199
形如a*x+b*y=c 为不定方程，a,b>0其实无所谓，因为gcd(a,b)=gcd(|a|,|b|) //gcd为最大公约数 由数论的定理所知，当c%gcd==0，不定方程有解,现在我们来求这个解. gcd=gcd(a,b);a*b=gcd*lcm; //lcm为最小公倍数 a'=a/gcd;b ...
分类：
其他好文 时间：
2016-03-29 23:45:18
阅读次数：
278
题目大概说给一棵有点权的树，输出字典序最小的点对，使这两点间路径上点权的乘积模1000003的结果为k。 树的点分治搞了。因为是点权过根的两条路径的LCA会被重复统计，而注意到1000003是质数，所以这个用乘法逆元搞一下就OK了。还有要注意“治”的各个实现，把时间复杂度“控制”在O(nlogn)。 ...
分类：
其他好文 时间：
2016-03-27 18:06:33
阅读次数：
151
1. 题目描述求乘法逆元。2. 基本思路利用扩展gcd求逆元，模板题目。3. 代码 1 /* 3609 */ 2 #include <iostream> 3 #include <sstream> 4 #include <string> 5 #include <map> 6 #include <que
分类：
其他好文 时间：
2016-02-29 14:08:36
阅读次数：
158
定义：满足a*k≡1 (mod p)的k值就是a关于p的乘法逆元。 为什么要有乘法逆元呢？ 当我们要求（a/b） mod p的值，且a很大，无法直接求得a/b的值时，我们就要用到乘法逆元。 我们可以通过求b关于p的乘法逆元k，将a乘上k再模p，即(a*k) mod p。其结果与(a/b) mod p
分类：
其他好文 时间：
2016-02-26 21:53:38
阅读次数：
226
题目链接 求sigma(i : 1 to n)i^k。 为了做这个题这两天真是补了不少数论， 之前连乘法逆元都不知道... 关于拉格朗日插值法， 我是看的这里http://www.guokr.com/post/456777/， 还挺有趣... 根据题目给出的例子我们可以发现， k次方的通项公式的最高
分类：
其他好文 时间：
2016-02-25 21:10:24
阅读次数：
271
满足 a * k ≡ 1 (mod p) 的k 叫做 a关于p的乘法逆元。另一种表达方法是 k ≡ a-1 (mod p) 逆元在密码学中有广泛应用，AES密码体系的字节替代就是运用了逆元。（不知道说的smg） 应用： 我们知道(a+b)%p=(a%p+b%p)%p (a*b)%p=(a%p)*(b
分类：
其他好文 时间：
2016-02-16 20:36:02
阅读次数：
231
题目大概问小于等于m个的物品放到n个地方有几种方法。 即解这个n元一次方程的非负整数解的个数$x_1+x_2+x_3+\dots+x_n=y$，其中0<=y<=m。 这个方程的非负整数解个数是个经典问题，可以+1转化正整数解的个数用插板法解决：$C_{y+n-1}^{n-1}=C_{y+n-1}^y
分类：
其他好文 时间：
2016-02-07 17:22:41
阅读次数：
202
