中国剩余定理(51nod)

先看古代大将的故事

韩信点兵:韩信让士兵先后以三人一排,五人一排,七人一排地变换队形,而他只要每次掠一眼队伍的排尾就知道总人数了

而中国剩余定理就是韩信的计算方法;


用数学表示:

对于一个数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;

}

 

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值