“韩信点兵”问题的数学真经

金创想

于 2024-04-06 16:38:56 发布

阅读量1.2k

点赞数 25

分类专栏：信息学奥赛 # 算法文章标签：算法韩信点兵同余方程组中国剩余定理

本文链接：https://blog.csdn.net/jjmhx/article/details/137431644

版权

信息学奥赛同时被 2 个专栏收录

114 篇文章 30 订阅

订阅专栏

算法

45 篇文章 0 订阅

订阅专栏

【题目描述】

相传韩信才智过人，从不直接清点自己军队的人数，只要让士兵先后以三人一排、五人一排、七人一排地变换队形，而他每次只掠一眼队伍的排尾就知道总人数了。输入包含多组数据，每组数据包含3个非负整数a，b，c，表示每种队形排尾的人数（a＜3，b＜5，c＜7），输出总人数的最小值（或报告无解）。已知总人数不小于10，不超过100。输入到文件结束为止。

【样例输入】

2 1 6

2 1 3

【样例输出】

Case 1: 41

Case 2: No answer

【题目来源】

刘汝佳《算法竞赛入门经典第2版》习题2-2　韩信点兵（hanxin）

【解析】

一、题目求解：穷举法

思路很简单：从10到100逐个检查是否存在对3、5、7取余分别为a、b、c的数。

代码如下：

#include<stdio.h>
int main(){
    int a, b, c, kase=0;
    while(scanf("%d%d%d", &a, &b, &c) == 3){
        int flag=0;
        kase ++;
        for(int i=10; i<=100; i++){
            if(i%3==a && i%5==b && i%7==c) {
                flag = 1;
                printf("Case %d: %d\n", kase, i);
                break;
            }
        }
        if(0==flag) printf("Case %d: No answer\n", kase);
    }
    return 0;
}

本题的关键只有两个：

(1)结束输入问题。题目要求“输入到文件结束为止”，但为了程序的健壮性，最好还是用代码中的形式，而不要用下面的代码：

while (scanf("%d %d %d", &a, &b, &c) != EOF)

这个代码也符合题目要求，但是一旦输入个字母或者小数，你将瞬间体验到自动刷屏的快感。

(2)多组数据变量重置问题。此问题在《“多组数据”题的注意事项，天杀的“鲁棒性”》一文中已详细阐述过。变量flag是参与到每组的计算之中的，所以应该放到while循环之内进行初始化。

题目到这里就解完了，但是老金不由感慨：咱们用二十一世纪最强大的计算工具求解的问题，早在两千多年前那一位仅凭口算就能搞定。

将军真乃神人也！

相传故事是这样的：在楚汉相争时期，韩信曾率1500名将士与楚王大将李锋交战，战后返回大本营时遭遇楚军骑兵追击。韩信开启巧妙点兵模式，迅速计算出己方兵力，并成功击退敌军。

二、韩信点兵问题的高级算法

用现代数学的话说，韩信点兵的问题，本质上是一个同余方程组问题。

如果整数a和b对正整数m取模，得到的结果相同，那么我们就说a和b对模m同余。

以余数为2 1 6为例，韩信点兵问题的同余方程组如下：

x≡2(mod3)

x≡1(mod5)

x≡6(mod7)

这里的“≡”表示同余关系，“mod m”则表示对m取模。比如第一个同余方程，就表示x和2对模3同余，解这个同余方程就是要求x。

方程组的解法源自一本真经，叫《孙子算经》。从这本真经我得得出一个定理，叫中国剩余定理（Chinese Remainder Theorem，简称CRT）。

因为源自《孙子算经》，此定理又称孙子定理。它的意思是，对于任意一组两两互质的正整数m1, m2, ..., mk（模数3、5、7），以及任意一组整数a1, a2, ..., ak（余数2、1、6），都存在一个整数x，满足对于所有的i（1 ≤ i ≤ k），x对mi取模的结果等于ai，即x%mi = ai。这就是中国剩余定理的核心内容。

中国剩余定理的公式可以表示为：x = a1 * M1 * y1 + a2 * M2 * y2 + ... + ak * Mk * yk，其中Mi = M/mi，M是m1, m2, ..., mk的乘积；yi是Mi关于mi的模逆元，即Mi * yi ≡ 1 (mod mi)。

这个公式理解起来也很简单，因为x对任何一个mi取余，等号右边第i项之外的其它单项式都含有因子mi，因而都能被mi整除，故只有ai*Mi*yi这一项会产生余数。而咱们预先设定Mi*yi对mi的余数为1，所以ai*Mi*yi%m=ai，即上面的公式永远满足x%mi=ai。

对应到韩信点兵问题，x便是要求的一个解，mi就是模数3、5、7，ai是就x除以mi的余数。

以题目中的样例输入2 1 6为例：

按3人一排排列时，队尾剩下2人，即 x≡2(mod3)；

按5人一排排列时，队尾剩下1人，即 x≡1(mod5)；

按7人一排排列时，队尾剩下6人，即 x≡6(mod7)。

现在，我们应用中国剩余定理来求解这个同余方程组。步骤如下：

(1)分别求出上面公式等号右侧中的每一项，即每个ai*Mi*yi的值；

先计算所有模数的最小公倍数（LCM）：M=3*5*7＝105。

①对于第一个方程 n≡2(mod3)，a1=2，M1= M/m1=M/3=35。

因为35×2≡1(mod3)，所以y1=2。