中国剩余定理模板

最新推荐文章于 2023-01-17 11:15:28 发布
最新推荐文章于 2023-01-17 11:15:28 发布
阅读量214 收藏 3
点赞数
分类专栏： 算法学习 数论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_41934068/article/details/98394330
版权

裸中国剩余定理

中国剩余定理的内容：

对于两两互质的一组整数m_{1},m_{2},.......,m_{n}，对于一组整数a_{1},a_{2},.......,a_{n}有x \equiv a_{i} ( mod m_{i})，则可以唯一地确定一个在m_{1}m_{2}......m_{n}内的解x使得此同余方程组成立。

方法：考虑两个方程的情况：即n%m1 = a，n%m2 = b，显然如果不在模意义下的话，那么a+(b-a)显然可以成为第二个方程的解，然而b-a的偏移量不满足模意义，我们可以求出m1关于m2的逆元w，则显然a+(b-a)m1w可以成为方程的解。

按照这种构造方法，同理，我们可以关于m1，m2，......,mn进行类似的构造，即可得出方程组的解。

POJ 1006

#include <cstdio>
#include <stack>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <map>
#include <vector>
#include <queue>
#include <set>
#define eps 1e-8
typedef long long ll;
const double PI = acos(-1.0);
const int maxn = 1e6;
const int INF = 0x3f3f3f;
const ll linf = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9+7;
using namespace std;
int a[maxn],m[maxn];
int n;
int exgcd(int a, int b, int &x, int &y)
{
    if(!b)
    {
        x = 1;
        y = 0;
        return a;
    }
    int d = exgcd(b,a%b,x,y);
    int z = x;
    x = y;
    y = z - (a/b)*y;
    return d;
}
int CRT()
{
    int M = 1,ans = 0;
    int x,y;
    n = 3;
    for(int i = 0; i<n; i++)
        M *= m[i];
    for(int i = 0; i<n; i++)
    {
        int Mi = M/m[i];
        exgcd(Mi,m[i],x,y);
        ans = (ans+a[i]*Mi*x)%M;
    }
    ans = (ans+M)%M;
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    m[0] = 23;
    m[1] = 28;
    m[2] = 33;
    int p,e,i,d;
    int cnt = 0;
    while(cin>>p>>e>>i>>d)
    {
        if(p+e+i+d == -4)
            break;
        a[0] = p;
        a[1] = e;
        a[2] = i;
        int ans = CRT();
        if(ans <= d) ans += 21252;
        cout<<"Case "<<++cnt<<": the next triple peak occurs in "<<ans-d<<" days."<<endl;
    }
    return 0;
}

扩展中国剩余定理

当mi不一定满足两两互质的条件时，则中国剩余定理不再适用。可以考虑使用数学归纳法，假设已急求出了前k-1个方程构成的方程组的一个解x。记m = \prod _{i = 1}^{k-1}m_{i}，则x+t*m是前k个方程的通解。

所以，只需考虑第k个方程，即求出一个整数t，使得x+tm\equiv a_{k} (mod m_{k})，该方程等价于tm\equiv a_{k}-x (mod m_{k})即可。若有解，则x' = x+t*m就是前k个方程构成的方程组的一个解。

HDU 3579

#include <cstdio>
#include <stack>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <map>
#include <vector>
#include <queue>
#include <set>
#define eps 1e-8
typedef long long ll;
const double PI = acos(-1.0);
const int maxn = 1e6;
const int INF = 0x3f3f3f;
const ll linf = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9+7;
using namespace std;
ll a[maxn],m[maxn];
int n;
ll gcd(ll a, ll b)
{
    return b==0?a:gcd(b,a%b);
}
ll exgcd(ll a, ll b, ll &x, ll &y)
{
    if(!b)
    {
        x = 1;
        y = 0;
        return a;
    }
    ll d = exgcd(b,a%b,x,y);
    ll z = x;
    x = y;
    y = z-(a/b)*y;
    return d;
}
ll exCRT()
{
    ll lcm = 1;
    for(int i = 0; i<n; i++)
        lcm = m[i]/gcd(m[i],lcm)*lcm;
    for(int i = 1; i<n; i++)
    {
        ll A = m[0],B = m[i],C = a[i]-a[0];
        ll x,y;
        ll d = exgcd(A,B,x,y);
        if(C%d)
            return -1;
        ll Mod = m[i]/d;
        ll K = ((x*C/d)%Mod+Mod)%Mod;
        a[0] += m[0]*K;
        m[0] *= m[i]/d;
    }
    if(!a[0])
        return lcm;
    return a[0];

}
int main()
{
    //ios::sync_with_stdio(false);
    int T;
    cin>>T;
    int tt = 1;
    while(T--)
    {
        cin>>n;
        for(int i = 0; i<n; i++)
            cin>>m[i];
        for(int i = 0; i<n; i++)
            cin>>a[i];
        ll ans = exCRT();
        printf ("Case %d: ",tt++);
        cout<<ans<<endl;
    }

    return 0;
}

POJ 2891

/**

**/
#include <cstdio>
#include <stack>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <map>
#include <vector>
#include <queue>
#include <set>
#define eps 1e-8
typedef long long ll;
const double PI = acos(-1.0);
const int maxn = 1e6;
const int INF = 0x3f3f3f;
const ll linf = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9+7;
using namespace std;
ll a[maxn],m[maxn];
int n;
ll gcd(ll a, ll b)
{
    return b==0?a:gcd(b,a%b);
}
ll exgcd(ll a, ll b, ll &x, ll &y)
{
    if(!b)
    {
        x = 1;
        y = 0;
        return a;
    }
    ll d = exgcd(b,a%b,x,y);
    ll z = x;
    x = y;
    y = z-(a/b)*y;
    return d;
}
ll exCRT()
{
    ll lcm = 1;
    for(int i = 0; i<n; i++)
        lcm = m[i]/gcd(m[i],lcm)*lcm;
    for(int i = 1; i<n; i++)
    {
        ll A = m[0],B = m[i],C = a[i]-a[0];
        ll x,y;
        ll d = exgcd(A,B,x,y);
        if(C%d)
            return -1;
        ll Mod = m[i]/d;
        ll K = ((x*C/d)%Mod+Mod)%Mod;
        a[0] += m[0]*K;
        m[0] *= m[i]/d;
    }
    if(!a[0])
        return lcm;
    return a[0];

}

int main()
{
    //ios::sync_with_stdio(false);

        while(cin>>n)
        {
            for(int i = 0; i<n; i++)
                cin>>m[i]>>a[i];
            ll ans = exCRT();
            cout<<ans<<endl;
        }



    return 0;
}

HDU 1788

/**

**/
#include <cstdio>
#include <stack>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <map>
#include <vector>
#include <queue>
#include <set>
#define eps 1e-8
typedef long long ll;
const double PI = acos(-1.0);
const int maxn = 1e6;
const int INF = 0x3f3f3f;
const ll linf = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9+7;
using namespace std;
int n;
ll gcd(ll a, ll b)
{
    return b==0?a:gcd(b,a%b);
}
ll lcm(ll a, ll b)
{
    return a/gcd(a,b)*b;
}
int main()
{
    ios::sync_with_stdio(false);
    int n,a;
    while(cin>>n>>a && n+a)
    {
        ll ans = 1;
        for(int i = 0; i<n; i++)
        {
            ll m;
            cin>>m;
            ans = lcm(m,ans);
        }
        cout<<ans-a<<endl;
    }

    return 0;
}

求满足要求的解的个数

HDU 1573

/**

**/
#include <cstdio>
#include <stack>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <iostream>
#include <map>
#include <vector>
#include <queue>
#include <set>
#define eps 1e-8
typedef long long ll;
const double PI = acos(-1.0);
const int maxn = 1e6;
const int INF = 0x3f3f3f;
const ll linf = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9+7;
using namespace std;
ll a[maxn],m[maxn];
int n;
ll M;
ll gcd(ll a, ll b)
{
    return b==0?a:gcd(b,a%b);
}
ll exgcd(ll a, ll b, ll &x, ll &y)
{
    if(!b)
    {
        x = 1;
        y = 0;
        return a;
    }
    ll d = exgcd(b,a%b,x,y);
    ll z = x;
    x = y;
    y = z-(a/b)*y;
    return d;
}
ll exCRT()
{
    ll lcm = 1;
    for(int i = 0; i<n; i++)
        lcm = m[i]/gcd(m[i],lcm)*lcm;
    for(int i = 1; i<n; i++)
    {
        ll A = m[0],B = m[i],C = a[i]-a[0];
        ll x,y;
        ll d = exgcd(A,B,x,y);
        if(C%d)
            return -1;
        ll Mod = m[i]/d;
        ll K = ((x*C/d)%Mod+Mod)%Mod;
        a[0] += m[0]*K;
        m[0] *= m[i]/d;
    }
    if(!a[0])
        return lcm;
    return a[0];
}
int main()
{
    ios::sync_with_stdio(false);
    int T;
    cin>>T;
    int tt = 1;
    while(T--)
    {
        cin>>M>>n;
        for(int i = 0; i<n; i++)
            cin>>m[i];
            ll lcm = 1;
        for(int i = 0; i<n; i++)
            lcm = m[i]/gcd(m[i],lcm)*lcm;
        for(int i = 0; i<n; i++)
            cin>>a[i];
        ll ans = exCRT();
        if(ans == -1 || ans>M)
        {
            cout<<0<<endl;
            continue;
        }
        cout<<(M-ans)/lcm+1<<endl;
    }

    return 0;
}

 

lero1005
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
建议收藏算法基础课模板大全
08-12
基础算法 —— 代码模板链接 常用代码模板1——基础算法 ...中国剩余定理 高斯消元 组合计数 容斥原理 简单博弈论 动态规划 背包问题 线性DP 区间DP 计数类DP 数位统计DP 状态压缩DP 树形DP 记忆化搜索 贪心
常用算法模板_C++.zip
04-27
常用算法模板_C++.zip AC自动机，Dijkstra，Floyd，GCD，KMP，KMP扩展，Kruskal，LCM，LCS，LIS，Prim，...全排列，染色-二分图，树状数组，线段树-点更新，线段树-区间更新，向量，匈牙利算法，中国剩余定理，字典树
【中国剩余定理模板】中国剩余定理（crt）
blog
07-18 171
给出n个形同x ≡ ai(mod mi)的方程，保证mi互质的情况下，求出最小的x
中国剩余定理模板（互质与非互质）
稚始稚终的博客
03-19 148
互质 #include<bits/stdc++.h> using namespace std; int exgcd(int a,int b,int &x,int &y) { if(b==0) { x=1,y=0; return a; } int d=exgcd(b,a%b,y,x); y-=a/b*x; return d; } int Chinese_Remainder(int mod[],int p
中国剩余定理模板+习题
FlyingJack
03-07 8078
除数互质 void extgcd(int a,int b,int &x,int &y) { if(!b){x = 1,y = 0;} else{extgcd(b,a%b,y,x); y -= (a/b)*x;} } int inv(int a,int n) { int x,y; extgcd(a,n,x,y); ...
中国剩余定理模板poj1006
最新发布
s13166803785的博客
01-17 59
#include <cstdio> #include <iostream> #include <cstring> #include <cmath> #include <cstdlib> #include <algorithm> using namespace std; typedef long long ll; const int maxn = 20; ll exgcd(ll a, ll b, ll &x, ll &y
中国剩余定理模板题
king9666的博客
06-21 332
https://blog.csdn.net/wuswi0412/article/details/80112326
OI 信息奥赛 算法模板大全
02-16
- 数论：欧几里得算法求最大公约数，中国剩余定理，费马小定理等。 - 矩阵运算：快速幂、矩阵快速幂等，用于高次方运算的优化。 7. **搜索算法**： - A*搜索：一种启发式搜索算法，结合了Dijkstra算法和最佳优先...
Huantwofat之ACM模板库1
08-03
15. **中国剩余定理模板**： 中国剩余定理是数论中的一个重要结果，解决了同时解多个同余方程的问题，常用于大整数计算和密码学中。 16. **后缀数组（DC3算法）**： 后缀数组是字符串处理中的一种数据结构，用于...
几种博弈问题及算法
weixin_41934068的博客
08-14 9587
巴什博弈 巴什博弈是最基础的博弈游戏。 有一堆石子，共计n颗，规定二人每次拿1~m颗石头，先拿完者胜，求解先手是否能赢。 算法思想非常好考虑： 因为每次最多拿m颗石头，所以当n是m+1的倍数时，first无论拿多少颗，second只要在first之后拿能够凑够m+1颗石子就能胜利，所以first必败。 如果n不是m+1的倍数，假设n = (m+1)*s + r（1 <= r...
WARSHALL算法计算传递闭包
weixin_41934068的博客
05-10 2276
    考虑n+1个矩阵的序列M0，M1，……，Mn，将矩阵Mk的第i行记作Mk[i,j].    对于k=0,1,……，n，Mk[i,j]=1当且仅当在R的关系图中存在一条从xi到xj路径，并且这条路径除了端点外中间只经过{x1，x2，……，xk}中的结点。WARSHALL算法从M0开始，顺序计算M1，M2，……，直到Mn为止。    我们知道，对于关系R对应的关系矩阵M来讲，其传递闭包R（t）...
高斯消元总结
weixin_41934068的博客
09-10 795
什么是高斯消元 高斯消元是用于解形如的线性方程组的。 我们可以将系数存在系数和常数保存在增广矩阵中，通过行变换列变换求解。 高斯消元的思想是，将方程组中一方程中的未知数用含有另一未知数的方程表示，代入原方程，即可消去原方程的一个未知数。可以看出，高斯消元的时间复杂度是O（n^3）的。 我们来看一个例子： 求解线性方程组 首先，要将L1以下的等式中的x消除，然后再将L2以下的等式中的y...
数论基本定理典型例题
weixin_41934068的博客
09-03 654
1.HDU-2973（威尔逊定理） The math department has been having problems lately. Due to immense amount of unsolicited automated programs which were crawling across their pages, they decided to put Yet-Anoth...
prim算法和并查集算法解决简单最小生成树问题
weixin_41934068的博客
07-29 519
什么是最小生成树 在图{V,E}中，使用E中的边将V中所有的点联通，使其成为树，要求树中边权之和最小，这就是最小生成树。 prim算法 1.复杂度：O（n2） 2.适用范围：与dijkstra算法类似，不能处理负权边 3.思想：与dijkstra算法类似，松弛。找到距离源点（这里的源点任意）最近的点，如此循环。 4.代码实现： int mp[maxn][maxn],in[maxn]...
模意义下大整数乘法
weixin_41934068的博客
07-30 482
计算 与二进制取幂的思想一样，这次我们将其中的一个乘数表示为若干个 2 的整数次幂的和的形式。因为在对一个数做乘 2 并取模的运算的时侯，我们可以转化为加减操作防止溢出。 递归方法如下： ll p; ll qmul(ll a, ll b, ll p) { ll res = 0; a = ((a%p)+p)%p; b = ((b%p)+p)%p; if(a...
匈牙利算法和KM算法解决二分图匹配问题
weixin_41934068的博客
07-29 475
参考博文： https://www.cnblogs.com/logosG/p/logos.html https://www.cnblogs.com/wenruo/p/5264235.html 匈牙利算法 1.复杂度：O（n2） 2.适用范围：增广路径求最大匹配，无权图 3.思想：dfs进行匹配，如果匹配冲突，则检查能否改变冲突顶点的配对 4.代码实现： dfs代码模板： int ...
HDU-1847（裸SG）
weixin_41934068的博客
08-14 447
1、总共n张牌; 2、双方轮流抓牌； 3、每人每次抓牌的个数只能是2的幂次（即：1，2，4，8，16…） 4、抓完牌，胜负结果也出来了：最后抓完牌的人为胜者； 假设Kiki和Cici都是足够聪明（其实不用假设，哪有不聪明的学生~），并且每次都是Kiki先抓牌，请问谁能赢呢？ #include <cstdio> #include <algorithm> #...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值