先看古代大将的故事
韩信点兵:韩信让士兵先后以三人一排,五人一排,七人一排地变换队形,而他只要每次掠一眼队伍的排尾就知道总人数了
而中国剩余定理就是韩信的计算方法;
用数学表示:
对于一个数n,有n%a=m1,n%b=m2,n%c=m3;求满足条件的最小值n
首先想想怎样使n同时满足n%a==m1,n%b=m2,n%c=m3这些条件
我们 令k1满足k1%(b*c)==0,k1%a==1,如果n=k1*m1,那么n%a==m1,这样的n就满足了第一个条件(这里为什么要使k1分别被b,c整除,看到后面就会理解了,下面的k2,k3一样) 将这里的n称为第一种n
同理,令k2满足k2%(a*c)==0,k1%b==1,如果n=k2*m2,那么n%b==m2,这样的n就满足了第二个条件
将这里的n称为第二种n
同理,令k3满足k3%(a*b)==0,k1%c==1,如果n=k3*m3,那么n%c==m2,这样的n就满足了第三个条件
将这里的n称为第三种n
然后我们将三种n都加到到一起成为新的n,n=k1*m1+k2*m2+k3*m3;
则因为k2与k3都能被(a*c)与(a*b)整除,故(k2*m2+k3*m3)%a==0,则n%a==m1;
同理:因为k1与k3都能被(b*c)与(a*b)整除,故(k1*m1+k3*m3)%b==0,则n%b==m2;
同理:因为k1与k2都能被(b*c) 与 (a*c)整除,故(k1*m1+k2*m2)%c==0, 则n%c==m3;
因此n=k1*m1+k2*m2+k3*m3就满足所有条件
现在只要n最小就行,由于n都是被a,b,c操作,故将n减去若干个a*b*c不会影响条件不成立,
因此最后答案ANSWER=(k1*m1+k2*m2+k3*m3)%(a*b*c)
中国剩余定理的一般化 链接
韩信点兵除以的数只有a,b,c三个数,如果将其扩展到t个数,我们可以以同样的方法求最小的n;(将除以的数称为因子)
韩信点兵,有三个条件,通过分别求出满足第一个条件的n,满足第二个条件的n,满足第三个条件的n,(由于每种n不能受其他的因子的影响,因此要求都能被其他所有因子整除); 然后将3种n加起来,就形成满足3个条件的n,这样的n并不一定是最小的,因此将其取余所有因子的乘积,就是最小的 | 对于一般的情况,有t个因子, 首先要求出满足第一个条件的n,一直求到满足第t个条件的n 然后将这t种n加起来形成新的的n(这样的n满足所有条件); 然后n%(所有因子的乘积)就最小的n |
对于多个因子的情况的代码:
#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
#define ll long long
int main()
{
int t;
int a[10];
ll k[10];//用来存满足第i个条件的n
int b[10];
while(scanf("%d",&t)==1)
{
ll sum=1;
for(int i=0;i<t;i++)
{
scanf("%d%d",&a[i],&b[i]);
sum*=a[i];//所有因子的乘积
}
for(int i=0;i<t;i++)
{
for(int j=1;;j++)//求出最小的满足第i个条件的n,存放在k[i]里
if((sum/a[i])*j%a[i]==1)//sum/a[i]就是除a[i]因子外其他所有因子的乘积
{
k[i]=sum/a[i]*j;
break;
}
}
ll res=0;
for(int i=0;i<t;i++)
{
res=(res+k[i]*b[i]);
}
//res就是满足所有条件的n
printf("%lld\n",res%sum);
}
return 0;
}