在《孙子算经》中有这样一个问题:“今有物不知其数,三三数之剩二(除以3余2),五五数之剩三(除以5余3),七七数之剩二(除以7余2),问物几何?”这个问题称为“孙子问题”,该问题的一般解法国际上称为“中国剩余定理”。具体解法分三步:
-
- 找出三个数:从3和5的公倍数中找出被7除余1的最小数15,从3和7的公倍数中找出被5除余1 的最小数21,最后从5和7的公倍数中找出除3余1的最小数70。
- 用15乘以2(2为最终结果除以7的余数),用21乘以3(3为最终结果除以5的余数),同理,用70乘以2(2为最终结果除以3的余数),然后把三个乘积相加(15*2+21*3+70*2)得到和233。
- 用233除以3,5,7三个数的最小公倍数105,得到余数23,即233%105=23。这个余数23就是符合条件的最小数。
在理解这个计算过程是非常简单的,但真正在代码中实现却不然,刚开始学知道模板对就套,现在有时间来详解一下中国剩余定理的模板
(感觉网上的大神们没有解释代码的(对他们来说太简单了!^^)),接下来我就详解一下中国剩余定理的模板,拿15*2来举例吧,代码中是这样呈现的
d=Extended_Euclid(w[i],m,x,y);
ret=(ret+y*m*a[i])%n;
ret=(ret+y*m*a[i])%n;
15是3和5的最小公倍数,现在我们需要找3和5的公倍数中除7余2的那个,即(15*x)%7==2,我们知道(15*x)%7,=gcd(15,7)中x为15的逆元
所以我们所求的x其实是15 的逆元,有因为所有数互质所以gcd(15,7)=1,所以15*x+7*y=gcd(15,7)所求除的值是15*x+7*y=1的解,
所以最后还需要乘以w[i](这里就是2了)
模板:
#include <iostream>
#include <cstdio>
using namespace std;
int Extended_Euclid(int a,int b,int &x,int &y) //扩展欧几里得算法
{
int d;
if(b==0)
{
x=1;
y=0;
return a;
}
d=Extended_Euclid(b,a%b,y,x);
y-=a/b*x;
return d;
}
int Chinese_Remainder(int a[],int w[],int len) //中国剩余定理 a[]存放余数 w[]存放两两互质的数
{
int i,d,x,y,m,n,ret;
ret=0;
n=1;
for (i=0; i<len; i++)
n*=w[i];
for (i=0; i<len; i++)
{
m=n/w[i];
d=Extended_Euclid(w[i],m,x,y);
ret=(ret+y*m*a[i])%n;
}
return (n+ret%n)%n;
}
int main()
{
int n,i;
int w[15],b[15];
while (scanf("%d",&n),n)
{
for (i=0; i<n; i++)
{
scanf("%d%d",&w[i],&b[i]);
}
printf("%d\n",Chinese_Remainder(b,w,n));
}
return 0;
}
#include <cstdio>
using namespace std;
int Extended_Euclid(int a,int b,int &x,int &y) //扩展欧几里得算法
{
int d;
if(b==0)
{
x=1;
y=0;
return a;
}
d=Extended_Euclid(b,a%b,y,x);
y-=a/b*x;
return d;
}
int Chinese_Remainder(int a[],int w[],int len) //中国剩余定理 a[]存放余数 w[]存放两两互质的数
{
int i,d,x,y,m,n,ret;
ret=0;
n=1;
for (i=0; i<len; i++)
n*=w[i];
for (i=0; i<len; i++)
{
m=n/w[i];
d=Extended_Euclid(w[i],m,x,y);
ret=(ret+y*m*a[i])%n;
}
return (n+ret%n)%n;
}
int main()
{
int n,i;
int w[15],b[15];
while (scanf("%d",&n),n)
{
for (i=0; i<n; i++)
{
scanf("%d%d",&w[i],&b[i]);
}
printf("%d\n",Chinese_Remainder(b,w,n));
}
return 0;
}