简单算法实现之《中国剩余定理》

中国剩余定理
这个定理之前没见过,是在RSA算法原理(一)中提到的,想深入了解的可以去读下。

故事背景
  • 韩信点兵

秦朝末年,楚汉相争。一次,韩信将1500名将士与楚王大将李锋交战。苦战一场,楚军不敌,败退回营,汉军也死伤四五百人,于是韩信整顿兵马也返回大本营。当行至一山坡,忽有后军来报,说有楚军骑兵追来。只见远方尘土飞扬,杀声震天。汉军本来已十分疲惫,这时队伍大哗。韩信兵马到坡顶,见来敌不足五百骑,便急速点兵迎敌。他命令士兵3人一排,结果多出2名;接着命令士兵5人一排,结果多出3名;他又命令士兵7人一排,结果又多出2名。韩信马上向将士们宣布:我军有1073名勇士,敌人不足五百,我们居高临下,以众击寡,一定能打败敌人。汉军本来就信服自己的统帅,这一来更相信韩信是“神仙下凡”、“神机妙算”。于是士气大振。一时间旌旗摇动,鼓声喧天,汉军步步进逼,楚军乱作一团。交战不久,楚军大败而逃。

题目就是:将军点兵,三三数余2,五五数余3,七七数余2。问兵几何?

  • 孙子算经

今有物,不知其数,三三数之,剩二,五五数之,剩三,七七数之,剩二,问物几何?
答曰:二十三  
术曰:三三数之剩二,置一百四十,五五数之剩三,置六十三,七七数之剩二,置三十,并之,得二百三十三,以二百一十减之,即得。凡三三数之剩一,则置七十,五五数之剩一,则置二十一,七七数之剩一,则置十五,即得。

《孙子算经》里给出了解决办法

步骤一
除3余2 -> 取140
除5余3 -> 取63
除7余2 -> 取30
步骤二:
对上面取到的数求和
140 + 63 + 30 = 233
步骤三
用上面求出的结果减去210,就是结果
233 - 210 = 23
结果也就是23了

不过那140,63,30,还有210是怎么来的呢?这就涉及到中国剩余定理的算法原理了。

中国剩余定理原理

#题目:
设要求的数为x,每个除法的结果分别为k1,k2,k3,那么就有:
[1]  x / 3 = k1 + 2
[2]  x / 5 = k2 + 3
[3]  x / 7 = k3 + 2
1.对步骤一中的式子1求解过程为:取式子2和式子3中除数(5和7)的最小公倍数LCM(35),
然后把这个LCM作为x带入式子1。如果LCM不能适应式子1,就对LCM再加上一个LCM,
然后带入式子1中......直到满足式子1为止。

按照这个求解过程,分别对3个式子求解,取到的结果为:35,63,30
2.把上面求到的结果求和

35 + 63 + 30 = 128
3.这里对3个式子的除数取最小公倍数

3,5,7的最小公倍数为105
然后用上面的结果减去这个最小公倍数,就是结果了
128 - 105 = 23

这里估计有人要问了:这个35也不是140,而且105也不等于210啊?! 这个不用纠结,这个跟古代数学的发展史有关嘛! 当然,这是我瞎说的,想要深挖的可以去研究下相关资料。

下面具体说代码实现

首先创建个model类,类里面有除数和余数2个属性,然后提供一个便利构造器,方便使用

@interface CRTModel : NSObject

@property (assign, nonatomic) NSInteger divider;

@property (assign, nonatomic) NSInteger remain;

+(instancetype)modelWithDivider:(NSInteger)divider
                         remain:(NSInteger)remain;

@end

@implementation CRTModel

+(instancetype)modelWithDivider:(NSInteger)divider remain:(NSInteger)remain
{
    CRTModel *model = [[self alloc] init];
    model.divider = divider;
    model.remain = remain;
    return model;
}

@end

还需要用到2个函数,最大公约数GCD和最小公倍数LCM

-(NSInteger)gcdOf:(NSInteger)a and:(NSInteger)b
{
    //这里使用欧几里德算法
    static const NSInteger(^GCDRecursionBlock)(NSInteger,NSInteger) 
    = ^(NSInteger ra, NSInteger rb){
        if (!ra || !rb) return MAX(ra, rb);
        return GCDRecursionBlock(rb,ra%rb);
    };
    return GCDRecursionBlock(a,b);
}

//最小公倍数
-(NSInteger)lcmOf:(NSInteger)a and:(NSInteger)b
{
    return a * b / [self gcdOf:a and:b]; //最小公倍数等于两数之积除以最大公约数
}

然后开始实现算法

0.先创建3个Model,把题目"抄"一遍

CRTModel *model1 = [CRTModel modelWithDivider:3 remain:2];
CRTModel *model2 = [CRTModel modelWithDivider:5 remain:3];
CRTModel *model3 = [CRTModel modelWithDivider:7 remain:2];
1.取到算法中步骤1的结果:

-(NSInteger)getSubMinNumberOfDivider1:(NSInteger)divider1
                              divider2:(NSInteger)divider2
                           andDivider3:(NSInteger)divider3
                             remainOf3:(NSInteger)remainOf3
{
    NSInteger lcm = [self lcmOf:divider1 and:divider2];
    NSInteger result = lcm;
    while ((result % divider3) != remainOf3) {
        result += lcm;
    }
    return result;
}

NSInteger a = [self getSubMinNumberOfDivider1:model1.divider divider2:model2.divider andDivider3:model3.divider remainOf3:model3.remain];
NSInteger b = [self getSubMinNumberOfDivider1:model2.divider divider2:model3.divider andDivider3:model1.divider remainOf3:model1.remain];
NSInteger c = [self getSubMinNumberOfDivider1:model3.divider divider2:model1.divider andDivider3:model2.divider remainOf3:model2.remain];
2.求和

NSInteger sum = a + b + c;
3.求全部除数的最小公倍数,然后求结果

//使用"更相减损术"
NSInteger lcmOfAll = 
    [self lcmOf:[self lcmOf:model1.divider and:model2.divider] and:model3.divider];
    
//求结果
NSInteger result = sum - lcmOfAll;
NSLog(@"计算结果为:%ld + %ld * k (k为自然数)",result,lcmOfAll);    

算法代码已经上传GitHub:
https://github.com/Bluelich/MyBlogCode/tree/master/CRTAlgorithm

就这么多了。这是我昨晚回顾RSA加密时候看到的,刚好就学习下了,其实想再复习下欧拉函数什么的,都忘光了! 有时间就搞!!

你可能感兴趣的:(简单算法实现之《中国剩余定理》)