标签：style blog color strong for 问题 代码 div
【模运算的定义及概念】
模运算即求余运算。“模”是“Mod”的音译，模运算多应用于程序编写中。 Mod的含义为求余。模运算在数论和程序设计中都有着广泛的应用，从奇偶数的判别到素数的判别，从模幂运算到最大公约数的求法，从孙子问题到凯撒密码问题，无不充斥着模运算的身影。
【基本理论】
给定一个正整数p，任意一个整数n，一定存在等式 n = kp + r ；
其中k、r是整数，且 0 ≤ r < p，称呼k为n除以p的商，r为n除以p的余数。
对于正整数p和整数a,b，定义如下运算：
取模运算：a % p（或a mod p），表示a除以p的余数。
模p加法：(a + b) % p ，其结果是a+b算术和除以p的余数，也就是说，(a+b) = kp +r，则(a + b) % p = r。
模p减法：(a-b) % p ，其结果是a-b算术差除以p的余数。
模p乘法：(a * b) % p，其结果是 a * b算术乘法除以p的余数。
说明：
1. 同余式：正整数a，b对p取模，它们的余数相同，记做 a ≡ b % p或者a ≡ b (mod p)。
2. n % p得到结果的正负由被除数n决定,与p无关。例如：7%4 = 3， -7%4 = -3， 7%-4 = 3， -7%-4 = -3。
【基本性质】
（1）若p|(a-b)，则a≡b (% p)。例如 11 ≡ 4 (% 7)， 18 ≡ 4(% 7)
（2）(a % p)=(b % p)意味a≡b (% p)
（3）对称性：a≡b (% p)等价于b≡a (% p)
（4）传递性：若a≡b (% p)且b≡c (% p) ，则a≡c (% p)
【运算规则】
模运算与基本四则运算有些相似，但是除法例外。其规则如下：
（1）(a + b) % p = (a % p + b % p) % p
（2）(a - b) % p = (a % p - b % p) % p
（3）(a * b) % p = (a % p * b % p) % p
（4）(a^b) % p = ((a % p)^b) % p
（5）结合率： ((a+b) % p + c) % p = (a + (b+c) % p) % p
（6）((a*b) % p * c)% p = (a * (b*c) % p) % p
（7）交换率： (a + b) % p = (b+a) % p
（8）(a * b) % p = (b * a) % p
（9）分配率： ((a +b)% p * c) % p = ((a * c) % p + (b * c) % p) % p
（10）重要定理：若a≡b (% p)，则对于任意的c，都有(a + c) ≡ (b + c) (%p)；
（11）若a≡b (% p)，则对于任意的c，都有(a * c) ≡ (b * c) (%p)；
（12）若a≡b (% p)，c≡d (% p)，则 (a + c) ≡ (b + d) (%p)，(a - c) ≡ (b - d) (%p)，
(a * c) ≡ (b * d) (%p)，(a / c) ≡ (b / d) (%p)；
（13）若a≡b (% p)，则对于任意的c，都有ac≡ bc (%p)；
【基本应用】
1.判断素数
2.求最大公约数利用gcd(a,b)=gcd(b,a mod b);(欧几里得算法与扩展欧几里得算法）
3.模幂运算
利用模运算的运算规则，我们可以使某些计算得到简化。例如，我们想知道3333^5555的末位是什么。很明显不可能直接把3333^5555的结果计算出来，那样太大了。但我们想要确定的是3333^5555（%10），所以问题就简化了。
根据运算规则（4）a^b% p = ((a % p)^b) % p ，我们知道3333^5555（%10）= 3^5555（%10）。由于3^4 = 81，所以3^4（%10）= 1。
根据运算规则（3） (a * b) % p = (a % p * b % p) % p ，由于5555 = 4 * 1388 + 3，我们得到3^5555（%10）=（3^(4*1388) * 3^3）（%10）=（（3^(4*1388)（%10）* 3^3（%10））（%10）=（1 * 7）（%10）= 7。
下面的结论可以得到一种更好的算法：
如果N是偶数，那么X^N =（X*X）^[N/2]；
如果N是奇数，那么X^N = X*X^(N-1) = X *（X*X）^[N/2]；
其中[N]是指小于或等于N的最大整数。
下面我写出这一算法的代码：
1 int PowerMod( int x, int n, int p) 2 { 3 if(n==0) 4 return 1; 5 else 6 int cmp=PowerMod((x*x)%p,n/2,p);//递归调用 7 if(n&1!=0)//如果是奇数则多乘以一个x 8 cmp=(cmp*x)%p; 9 return cmp; 10 }
4.中国剩余定理
中国剩余定理可以描述为：
若某数x分别被d1、、…、dn除得的余数为r1、r2、…、rn，则可表示为下式：
x=R1r1+R2r2+…+Rnrn+RD 其中R1是d2、d3、…、dn的公倍数，而且被d1除后余数为1；（称为R1相对于d1的数论倒数）
R2 是d1,d3,d4,······dn的公倍数，而且被d2除后余数为1；
…
Rn是d1、d2、…、dn-1的公倍数，而且被dn除，余数为1；
D是d1、d2、…、dn的最小公倍数；R是任意整数（代表倍数），可根据实际需要决定；
且d1、、…、dn必须互质，以保证每个Ri(i=1,2，…，n)都能求得.
（注：因为R1对d1求余为1，所以R1r1对d1求余为r1，这就是为什么是R1对d1求余为1的目的，其次，R2r2,R3r3…Rnrn对d1求余都是0）
对于中国剩余定理有个简单理解并记忆的方法：
中国剩余定理的思想在于先找到一个满足条件的数，不管是不是最小的，如果不是最小的就不断减公倍数，中国剩余定理的巧妙在于把满足条件的数分成三个部分相加。
例如：
假设满足条件的数为：M=3*5*a+5*7*b+3*7*c
先让它满足条件1：除以3余1，我们看M的第一部分和第三部分能被3整除，只要第二部分除以3余1就行了，选择b=2就满足
再满足条件2：除以5余2，考虑第三部分就行了，选择c=2就满足
最后满足条件3：除以7余3，考虑第一部分就行了，选择a=3
这样M=3*5*3+5*7*2+3*7*2=157，比公倍数大，减去公倍数，157-105=52是满足条件最小数
以我个人理解写成下面这个形式（以3个数为例）
X被a,b,c处分别余r1,r2,r3。表示为：
X%a = r1 x%b = r2 x%c = r3
bc*k1 % a = 1 ac*k3 % b = 1 ab*k3 % c = 1
所以
x = bc * k1 * r1 + ac * k2 * r2 + ab * k3 * r3
下面我贴出代码：
1 int PowerMod( int d1, int d2, int d3，int s1,int s2,int s3)//d1,d2,d3是三个除数 2 { 3 int r1=d2*d3,r2=d1*d3,r3=d1*d2; 4 for(int i=1;1;i++) 5 if((r1*i)%d1==1) 6 { 7 r1=r1*i; 8 break; 9 } 10 for(int j=1;1;j++) 11 if((r2*j)%d2==1) 12 { 13 r2=r2*j; 14 break; 15 } 16 for(int z=1;1;z++) 17 if((r3*z)%d3==1) 18 { 19 r3=r3*z; 20 break; 21 } 22 cout<<r1<<‘ ‘<<r2<<‘ ‘<<r3<<endl; 23 int left=(r1*s1+r2*s2+r3*s3)%LCM(r1,r2,r3);//所求数字的最小值 24 return left; 25 }
标签：style blog color strong for 问题 代码 div
原文地址：http://www.cnblogs.com/khbcsu/p/3880053.html
