中国剩余定理模板+习题

最新推荐文章于 2023-07-25 18:59:37 发布
最新推荐文章于 2023-07-25 18:59:37 发布
阅读量8k 收藏 9
点赞数 3
分类专栏： 数论 文章标签： ACM
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_42264485/article/details/88292114
版权

除数互质

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);           
    x = (x % n + n) % n;
    return x;
}
int china()  
{
    int p = 1,ans = 0;
    for(int i=0;i<n;i++)    p *= m[i];
    for(int i=0;i<n;i++){
        int w = p / m[i];
        ans = (ans + a[i] * w * inv(w,m[i])) % p;         
    }
    return ans;
}

除数不互质

void extgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(!b){d = a,x = 1,y = 0;}
	else{extgcd(b,a%b,d,y,x); y -= x*(a/b);}
}
ll China(){
	ll a1 = a[1],  n1 = n[1];   
	for(int i=2;i<=num;i++){
		ll a2 = a[i],	n2 = n[i];
		ll c = a2 - a1;
		ll k1,k2,d;
		extgcd(n1,n2,d,k1,k2);    //求解n1k1 - n2k2 = a2 - a1
		if(c % d)		return -1;
		c /= d,		n2 /= d;
		k1 *= c;
		k1 = (k1 % n2 + n2) % n2;  //找出最小的整数解
		a1 += n1 * k1;
		n1 *= n2;
	}
	return a1;
}
/*防止溢出*/
__int128 China(){
	__int128 a1 = a[1],  n1 = n[1];   
	for(int i=2;i<=num;i++){
		__int128 a2 = a[i],	n2 = n[i];
		__int128 c = a2 - a1;
		__int128 k1,k2,d;
		extgcd(n1,n2,d,k1,k2);
		if(c % d)		return -1;
		k1 = ((k1 * c / d) % (n2 / d) + (n2 / d)) % (n2 / d);
		a1 += n1 * k1;
		n1 = n1 * n2 / d;
	}
	return a1;
}

HDU3579

题解：输出最小正整数

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int const N = 20;
ll a[N],n[N];
int num;
void extgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(!b){d = a,x = 1,y = 0;}
	else{extgcd(b,a%b,d,y,x); y -= x*(a/b);}
}
ll China(){
	ll a1 = a[1],  n1 = n[1];   
	for(int i=2;i<=num;i++){
		ll a2 = a[i],	n2 = n[i];
		ll c = a2 - a1;
		ll k1,k2,d;
		extgcd(n1,n2,d,k1,k2);    //求解n1k1 - n2k2 = a2 - a1
		if(c % d)		return -1;
		c /= d,		n2 /= d;
		k1 *= c;
		k1 = (k1 % n2 + n2) % n2;  //找出最小的整数解
		a1 += n1 * k1;
		n1 *= n2;
	}
	return a1 > 0 ? a1 : n1;  //输出最小整数解，如果a1 = 0，那么a1 % n1还为0，所以a1 + n1.
}
int main(){
	int T,caser = 0;
	scanf("%d",&T);
	while(T--){
		scanf("%d",&num);
		for(int i=1;i<=num;i++)	scanf("%lld",&n[i]);
		for(int i=1;i<=num;i++)	scanf("%lld",&a[i]);   
		printf("Case %d: %lld\n",++caser,China());
	}
	return 0;
}

HDU1573

题解：输出≤n的个数

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int const N = 10 + 5;
ll a[N],n[N],high;
int num;
void extgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(!b){d = a,	x = 1,	y = 0;}
	else{extgcd(b,a%b,d,y,x);	y -= x*(a/b);}
}
ll solve(){
	ll a1 = a[1],	n1 = n[1];
	for(int i=2;i<=num;i++){
		ll a2 = a[i],	n2 = n[i],	c = a2 - a1;
		ll k1,k2,d;
		extgcd(n1,n2,d,k1,k2);
		if(c % d)	return 0;
		n2 /= d,	c /= d;
		k1 *= c;
		k1 = (k1 % n2 + n2) % n2;
		a1 += k1 * n1;
		n1 *= n2;
	}
	if(a1 > high)	return 0;
	else if(a1 == 0)	return high / n1;    //必须要正整数
	else return (high - a1) / n1 + 1;
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--){
		scanf("%lld%d",&high,&num);
		for(int i=1;i<=num;i++)	scanf("%lld",&n[i]);
		for(int i=1;i<=num;i++)	scanf("%lld",&a[i]);
		printf("%lld\n",solve());
	}
	return 0;
}

POJ2891

题解：模板题
代码：

#include <algorithm>
#include <cstdio>
#include <vector>
#include <iostream>
using namespace std;
typedef long long ll;
ll a1,a2,n1,n2;
int num;
void extgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(!b){d = a,x = 1,y = 0;}
	else{extgcd(b,a%b,d,y,x);	y -= a/b*x;}
}
int main(){
	while(~scanf("%d",&num)){
		bool flag = true;
		scanf("%lld%lld",&n1,&a1);
		for(int i=1;i<num;i++){
			scanf("%lld%lld",&n2,&a2);
			ll c = a2 - a1,k1,k2,d;
			extgcd(n1,n2,d,k1,k2);
			if(c % d)	flag = 0;
			n2 /= d,	c /= d;
			k1 *= c;
			k1 = (k1 % n2 + n2) % n2;
			a1 += n1 * k1;
			n1 *= n2;
		}
		printf("%lld\n",flag ? a1 : -1);	
	}
	return 0;
}

HDU1370

题解：还是模板题目，转换一下即可。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll const N =  21252;
ll a[5],n[5],m,d;
void extgcd(ll a,ll b,ll &x,ll &y){
	if(!b){x = 1,y = 0;}
	else{extgcd(b,a%b,y,x);	y -= x*(a/b);}
}
ll inv(ll a,ll n){
	ll x,y;
	extgcd(a,n,x,y);
	return (x % n + n) % n;
}
ll solve(){
	ll ans = 0;
	for(int i=1;i<=3;i++){
		ll w = m / n[i];
		ans = (ans + a[i] * w * inv(w,n[i])) % m;
	}
	return ans;
}
int main(){
	n[1] = 23,	n[2] = 28,	n[3] = 33;
	m = n[1] * n[2] * n[3];
	int caser = 0,T;
	scanf("%d",&T);
	while(cin>>a[1]>>a[2]>>a[3]>>d && d != -1){
		ll period = solve();
		while(period <= d)	period += m;   //第一次重合可能≤起始日期，下一个triple day，所以即使这一次是也不可以。
		printf("Case %d: the next triple peak occurs in %lld days.\n",++caser,period - d);
	}
	return 0;
}

HDU5668

题解：模拟一下。加入总共有n个人，一共数到K，现在统计的是第i个出局的人，从上一个出出局的人到这一个统计的个数为p。那么K % （n - i + 1） = p。每次出局一个做标记。这是一组除数为n - i + 1，余数为p的线性方程组。套模板。详情见代码。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int const N = 25;
int num[N],n,tmp,vis[N];
ll m[N],a[N];
void extgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(!b){d = a,x = 1, y = 0;}
	else{extgcd(b,a % b,d,y,x); y -= x * (a/b);}
}
ll solve(){
	ll a1 = a[1],	n1 = m[1];
	for(int i=2;i<=n;i++){
		ll a2 = a[i],	n2 = m[i],	c = a2 - a1;
		ll d,k1,k2;
		extgcd(n1,n2,d,k1,k2);
		if(c % d)	return -1;
		n2 /= d,	c /= d;
		k1 *= c;
		k1 = (k1 % n2 + n2) % n2;
		a1 += k1 * n1;
		n1 *= n2;
	}
	return a1;
}
int main(){
	int T;
	cin>>T;
	while(T--){
		cin>>n;
		memset(vis,true,sizeof(vis));   //表明在的人
		for(int i=1;i<=n;i++){
			cin>>tmp;   
			num[tmp] = i;   //第tmp个出局的是第i个人
		}
		int j = 1;    //现在是数到第j个人
		for(int i=1;i<=n;i++){	  //第i个出局的人
			int k = 1;  //统计余数
			while(num[i] != j){
				if(vis[j])	k++;
				j = j % n + 1;	 
			}
			vis[num[i]] = false;
 			m[i] = (ll)(n - i + 1);
			a[i] = (ll)k;
		}
		ll ans = solve();
		if(ans == -1)	printf("Creation August is a SB!\n");
		else 	printf("%lld\n",ans);
	}
}

HDU1930

题解：耐心看完题目，发现很简单。输出注意一下，最后的空格不能输出。否则PE。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int const N = 50 + 10;
int n;
char mp[N];
ll m[5],a[5],ans[N];
ll num;
void extgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(!b){d = a,	x = 1,	y = 0;}
	else{extgcd(b,a%b,d,y,x);	y -= x*(a/b);}
}
ll China(){
	ll a1 = a[1],	n1 = m[1];
	for(int i=2;i<=4;i++){
		ll a2 = a[i],	n2 = m[i],	c = a2 - a1;
		ll k1,k2,d;
		extgcd(n1,n2,d,k1,k2);
		if(c % d)	return 0;
		n2 /= d,	c /= d;
		k1 *= c;
		k1 = (k1 % n2 + n2) % n2;
		a1 += k1 * n1;
		n1 *= n2;
	}
	return a1;
}
int main(){
	int T;
	for(int i=1;i<=26;i++)	mp[i] = i + 'A' - 1;
	mp[27] = ' ';
	cin>>T;
	while(T--){
		cin>>n;
		for(int i=1;i<=4;i++)	cin>>m[i];
		for(int i=1;i<=n;i++){
			cin>>num;
			for(int k=4;k>=1;k--){
				a[k] = num % 100;
				num /= 100;
			}
			ans[i] = China();
		}
		for(int i=1;i<=n;i++){	//080527
			string s;
			s[2] = mp[ans[i] % 100],	ans[i] /= 100;
			s[1] = mp[ans[i] % 100],	ans[i] /= 100;
			s[0] = mp[ans[i]];
			if(i != n)	cout<<s[0]<<s[1]<<s[2];
			else{
				if(s[0] != ' ')	cout<<s[0];
				if(s[1] != ' ')	cout<<s[1];
				if(s[2] != ' ')	cout<<s[2];
			}
		}
		cout<<endl;
	}
	return 0;
}

 

月光下的魔术师 
关注
  • 3
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
秘密共享方案_中国剩余定理
11-15
#某电信安数学基础实验（3）基于中国剩余定理的秘密共享方案（更新）： #涉及位数的均为二进制位 #python #使用中国剩余定理 #使用（t,n）门限来控制秘密，即：n个子秘密中任取t个或以上即可计算得到秘密，而任意t-1...
几道中国剩余定理例题
m0_53524883的博客
10-15 8727
通过几道例题浅析解一次同余式组时遇到的常见问题
中国剩余定理（孙子定理）详解
热门推荐
彼岸花苏陌的博客
09-09 2万+
题目引出 在《孙子算经》中有这样一个问题：“今有物不知其数，三三数之剩二（除以3余2），五五数之剩三（除以5余3），七七数之剩二（除以7余2），问物几何？”这个问题称为“孙子问题”，该问题的一般解法国际上称为“中国剩余定理”。 解决方法 1.找出三个数：从3和5的公倍数中找出被7除余1的最小数15，从3和7的公倍数中找出被5除余1 的最小数21，最后从5和7的公倍数中找出除3余1的最小数70。 2.用15乘以2（2为最终结果除以7的余数），用21乘以3（3为最终结果除以5的余数），同理，用70乘以2（2为
中国剩余定理详解
weixin_44203780的博客
03-03 1940
title: 中国剩余定理详解 date: 2019-08-17 18:49:14 tags: - 数论 mathjax: true 中国剩余定理与扩展中国剩余定理 首先来看这个方程 {x≡a1(mod b1)x≡a2(mod b2)x≡a3(mod b3)⋯x≡an(mod bn)\begin{cases} x≡a_{1}(mod\ b_1)\\ x≡a...
中国剩余定理（孙子定理）
huatian5的博客
07-01 4145
中国剩余定理 以下内容摘自百度百科 例题解析 例一：一个数，除以5余1，除以3余2。问这个数最小是多少？ 采用通用的方法：逐步满足法 把除以5余1的数从小到大排列：1，6，11，16，21，26，…… 然后从小到大找除以3余2的，发现最小的是11. 所以11就是所求的数。 先满足一个条件，再满足另一个条件，所以称之为“逐步满足法”。 例二：一个数除以5余1，除以3也余
中国剩余定理总结+例题
weixin_43819578的博客
09-07 4781
解决类似n=b（mod a）的同余方程组的一种方法。 例题： 1.模板：P4777 【模板】扩展中国剩余定理（EXCRT）链接 给定 n组非负整数 a i,b i求解关于 xx 的方程组的最小非负整数解。 #include <cstdio> const int MAXN = 100010; typedef long long ll; int n; ll a[MAXN], b[MAXN...
C语言实现中国剩余定理
01-21
计算一次同余方程的解，并且调用miracl库函数，实现大数解的计算。
Python实现的中国剩余定理算法示例
09-21
主要介绍了Python实现的中国剩余定理算法,结合实例形式分析了中国剩余定理的概念、原理及具体算法实现技巧,需要的朋友可以参考下
2_中国剩余定理_信息安全数学基础_
09-30
中国剩余定理（Chinese Remainder Theorem，CRT）是数论中的一个重要概念，它在密码学，特别是信息安全领域，有着广泛的应用。这个定理解决了一类模线性同余方程组的问题，即给定一系列模数和相应的余数，找到一个数...
中国剩余定理 个人理解
weixin_53744560的博客
05-01 563
我们一定要养成随手保存文档的好习惯呀QAQ吐了好几百字的东西看到写错了一个字 按个 ctrl + z 结果全没了，心态给我搞崩了。 （csdn就是逊啦 bushi）
中国剩余定理讲解及例题
最新发布
A2105153335的博客
07-25 806
凡有所学，皆成性格。凡有所学，皆成性格。目录🍒。
【4.6 中国剩余定理详解】
LYS00Q的博客
07-16 508
概念 思想 例题 表达整数的奇怪方式
中国剩余定理+模板
qq_46667619的博客
05-01 285
中国剩余定理是一种用于求解诸如 形式的同余方程组的最小非负整数解x的定理，其中，a1,a2,…,an两两互质。 解法： 大佬博客 令 M=a1a2…an Mi= M / ai(1<=i<=n) Mi ti=1(mod ai)，ti为Mi在模ai下的逆元。 因为Mi与ai一定互质，但ai不一定为质数，可以用拓展欧几里得求逆元，不能用费马小定理求。 由此构造一个解 x0= b1 M1 t1+ b2 M2 t2+…+bn Mn tn；任意解即为x=x0+k*M,最小正整数解为x%M。 模板题：洛谷P
【数论系列】 中国剩余定理与拓展中国剩余定理
Hello World!
08-20 2426
在《孙子算经》中有这样一个问题：“今有物不知其数，三三数之剩二（除以3余2），五五数之剩三（除以5余3），七七数之剩二（除以7余2），问物几何？”这个问题称为“孙子问题”，该问题的一般解法国际上称为“中国剩余定理”。
中国剩余定理模板与证明+扩展（两种不同线性同余方程组的解法）
测试
09-12 9442
以x2为例，x2的组成是ai * M[2] * ti ，在这之中M[2]是 m/m2,所以M[2]是m1的倍数，因 此x2%m1的结果是0，x3~xn也是同理。5.最后得出来的x就是一个合法解，如果题目要求最小的非负整数解则可以让x对m取模让x落在 [0,m-1]的范围内。2.设ti是线性同余方程 M[i] * ti ≡ 1（mod mi）的一个解, 相当于求M[i]在mod（mi）下的逆元。3.由第2点以及模运算的性质可以结合上图推出以下公式。
密码学基础算法(二)中国剩余定理
洛奇看世界
09-25 3507
中国剩余定理推导和源码实现
【算法实践】1.4 编程-Biorhythms（中国剩余定理）
梁瑛平的博客
03-02 677
Description Some people believe that there are three cycles in a person’s life that start the day he or she is born. These three cycles are the physical, emotional, and intellectual cycles, and they h...
中国剩余定理——C++实现
m0_56473031的博客
03-04 1260
信息安全数学基础编程作业题2-2——实现中国剩余定理
中国剩余定理总结
zgdlxs的博客
02-14 2288
1.孙子问题 最早，在《孙子算经》中有这样一个问题：“今有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二，问物几何？用白话描述就是，现在有一个数不知道是多少，只知道这个数除以3余2，除以5余3，除以7余2， 问这个数是多少？ 上面的问题可以转换为以下这样一个方程组: 其中这个x就是我们要求的数。而中国剩余定理就是用来解决上述形式的方程组的x解的。 2.中国剩余定理公式 设正整数m1，m2,....,mk两两互素，则同余方程组: 3.解决例题 AcWing 204 表达整数的.

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值