中国剩余定理讲解及例题

凡有所学,皆成性格。
凡有所学,皆成性格。

目录🍒

  • 🍉==中国剩余定理==
    • 维基百科
    • [<font color=green>:cherries:形式描述:](#:cherries:形式描述:)
  • [<font color=red>:pear: 解法](#:pear: 解法)
    • [:strawberry:<font color=green>求解方法:](#:strawberry:求解方法:)
      • [:star:<font color=red>1.O(1)"大衍求一术"](#:star:1.O(1)“大衍求一术”)
      • [:star:<font color=lightgreen> 2.枚举O(n)](#:star: 2.枚举O(n))
      • [:star:<font color=lightred>3.逐级满足法](#:star:3.逐级满足法)
      • [:star:<font color=red>4.化为相同除数法](#:star:4.化为相同除数法)
  • [:dart:<font color=blude>参考例题](#:dart:参考例题)

🍉中国剩余定理

维基百科

🌽中国剩余定理,又称孙子定理或中国余数定理,是数论中的一个关于一元线性同余方程组的定理,说明了一元线性同余方程组有解的准则以及求解方法。


🍒形式描述:

一个整数除以三余二,除以五余三,除以七余二,求这个整数(x)。
≡ 就是全等于的意思   x m o   m 1 = = a 1   m o d   m 1 \equiv就是全等于的意思\,xmo\,m1==a1 \,mod\,m1 ≡就是全等于的意思xmom1==a1modm1
所求整数x和余数a就是同余的他们的余数都为a;

🍐 解法

2023-07-21 2023-07-23 2023-07-25 2023-07-27 2023-07-29 2023-07-31 2023-08-01 2023-08-03 已完成 进行中 解法一 解法二 中国剩余定理 解法

一阶线性同余方程组


S : = { a   m o d   d = e a   m o d   b = c a   m o d   f = g S:= \begin{cases}a\,mod\,d=e \\a\,mod\, b=c \\a\,mod\,f=g \end{cases} S:=⎩ ⎨ ⎧amodd=eamodb=camodf=g

类似这样求a的问题:一个整数除以三余二,除以五余三,除以七余二,求这个整数

🍓求解方法:

⭐️1.O(1)"大衍求一术"

mi除数,M除数相乘除了自己

|---|
| |

1 ) M i = M / m i , ∀   i ∈ { 1 , 2 , 3 , 4 , 5... , n } 是除了 m i 以外 n − 1 个数的乘积 1)Mi=M/mi,\forall\,i\in\lbrace1,2,3,4,5...,n\rbrace是除了mi以外n-1个数的乘积 1)Mi=M/mi,∀i∈{1,2,3,4,5...,n}是除了mi以外n−1个数的乘积

2 )计算出相应的数论倒数: t i = M i   m o d   m i ( 逆元) 2)计算出相应的数论倒数:ti=Mi\,mod\,mi(逆元) 2)计算出相应的数论倒数:ti=Mimodmi(逆元)

3 ) M i 与对应的 t i 相乘 3)Mi与对应的ti相乘 3)Mi与对应的ti相乘

4 )得到的结果与对应的余数相乘,结果相加为 a 4)得到的结果与对应的余数相乘,结果相加为a 4)得到的结果与对应的余数相乘,结果相加为a

5 )解的形式: x = a + k ∗ M , k ∈   Z 5)解的形式:x=a+k*M,k\in\,Z 5)解的形式:x=a+k∗M,k∈Z

⭐️ 2.枚举O(n)

枚举结果得出答案(笨办法):

比如一个整数除以三余二,除以五余三,除以七余二

就可以写成xmod2=3,xmod5=3,xmod7=2;

我们不知道的是x到底是2,5,7的几倍,因为都有余数。

所以设x是2,5,7的a,b,c倍

枚举a,b,c取得交集即为x

⭐️3.逐级满足法

先求出符合第一个式子得解,再向第二个式子和并一直到最后一个式子。

⭐️4.化为相同除数法

S : = { x   m o d   2 = 7 x   m o d   5 = 9 x   m o d   1 = 5 S:= \begin{cases}x\,mod\,2=7 \\x\,mod\, 5=9 \\x\,mod\,1=5 \end{cases} S:=⎩ ⎨ ⎧xmod2=7xmod5=9xmod1=5

这三个同余式,除数不同,分别为7、9、5,为了能利用同余式的和差特性,简化计算,先设法使它们的除数相同,为此:

X≡2 (mod 7 )两边都乘9x5,得Xx45≡2x45 (mod 7x45 ) →45 X≡90 (mod 315 ) ...(1)

X≡5 (mod 9 ) 两边都乘7x5,得Xx35≡5x35 (mod 9x35 ) →35 X≡175 (mod 315 ) ...(2)

X≡1 (mod 5 ) 两边都乘7x9,得Xx63≡1x63 (mod 9x63 ) →63 X≡ 63 (mod 315 ) ...(3)

​ (1)+(2)+(3)= (4) → 143 X≡328 (mod 315 ) ...4)

根据同余式的加减性质,(1)+(2)+(3)得:

143 X≡328 (mod 315 ) 即 143 X≡13 (mod 315 ),化为带余除式为:

143 X÷315=K ...13 亦即 143X-13=315 K ,有(143X-13)÷315=K(整数)

(143X-13)÷315=K(整数)用通式表示为 (AX-B)÷S=K

解得 X=86、K=38 (实际上不用它,在此仅确认它是整数就行了),

通解为 X=86+315 N ( N= 1、2、3 ... )

或 X≡86 (mod 315 )

🎯参考例题

P1495 【模板】中国剩余定理(CRT)/ 曹冲养猪

问题描述

自从曹冲搞定了大象以后,曹操就开始捉摸让儿子干些事业,于是派他到中原养猪场养猪,可是曹冲满不高兴,于是在工作中马马虎虎,有一次曹操想知道母猪的数量,于是曹冲想狠狠耍曹操一把。举个例子,假如有

16 头母猪,如果建了 3个猪圈,剩下 1头猪就没有地方安家了。如果建造了 5 个猪圈,但是仍然有 1 头猪没有地方去,然后如果建造了 7 个猪圈,还有 2 头没有地方去。你作为曹总的私人秘书理所当然要将准确的猪数报给曹总,你该怎么办?


格式输入

第一行包含一个整数 n ------ 建立猪圈的次数,接下来 n 行,每行两个整数 ai

,bi ,表示建立了 ai个猪圈,有 bi头猪没有去处。你可以假定 a1∼a n互质。


格式输出

输出包含一个正整数,即为曹冲至少养母猪的数目。


样例输入

3

3 1

5 1

7 2


样例输出

16


评测用例规模与约定

1≤n≤10,0≤bi <ai ≤100000,1≤∏a i≤10^18


解析

就是中国剩余定理和拓展gcd的一个联合运用,洛谷有更详细的解释。


参考程序

c 复制代码
#include<iostream>
#include<cmath>
#include<vector>
using namespace std;
typedef long long ll;
void exgcd(int a,int b,int &x,int &y)
{
    if(b==0){x=1,y=0;return;}
    else{
        exgcd(b,a%b,x,y);
        int t=x-a/b*y;
        x=y,y=t;
    }
}
int crt(vector<int>a,vector<int> m,int n)
{
    int M=1,r=0;
    for(int i=0;i<n;i++)M*=m[i];
    for(int i=0,Mi,x,y;i<n;i++)
    {
        Mi=M/m[i];
        exgcd(Mi,m[i],x,y);
        r=(r+a[i]*Mi*x)%M;
    }
    if(r<0)r+=M;
    return r;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    ll n;cin>>n;
    vector<int>c,d;
    for(int i=0;i<n;i++){
       ll a,b;cin>>a>>b;c.push_back(a);
       d.push_back(b);
    }
  int rr=crt(d,c,n);
  cout<<rr;
}

参考文章:
中国剩余定理
中国剩余定理得五种解法

相关推荐
南东山人2 小时前
一文说清:C和C++混合编程
c语言·c++
LNTON羚通3 小时前
摄像机视频分析软件下载LiteAIServer视频智能分析平台玩手机打电话检测算法技术的实现
算法·目标检测·音视频·监控·视频监控
哭泣的眼泪4084 小时前
解析粗糙度仪在工业制造及材料科学和建筑工程领域的重要性
python·算法·django·virtualenv·pygame
Ysjt | 深5 小时前
C++多线程编程入门教程(优质版)
java·开发语言·jvm·c++
ephemerals__5 小时前
【c++丨STL】list模拟实现(附源码)
开发语言·c++·list
Microsoft Word5 小时前
c++基础语法
开发语言·c++·算法
天才在此5 小时前
汽车加油行驶问题-动态规划算法(已在洛谷AC)
算法·动态规划
一只小小汤圆5 小时前
opencascade源码学习之BRepOffsetAPI包 -BRepOffsetAPI_DraftAngle
c++·学习·opencascade
legend_jz6 小时前
【Linux】线程控制
linux·服务器·开发语言·c++·笔记·学习·学习方法