浅析中国剩余定理(CRT)

中国剩余定理

用来求解同余方程组的最小非负整数解,其中 m_{1},m_{2},m_{3}......m_{n}都互质

       \left\{\begin{matrix} x\equiv a_{1}\, (\, mod\, m_{1})& & \\ x\equiv a_{2}\, (\, mod\, m_{2})& & \\ x\equiv a_{3}\, (\, mod\, m_{3})& & \\ \, \, \, ......& & \\ x\equiv a_{n}\, (\, mod\, m_{n})& & \end{matrix}\right.

首先让 M 等于所有 m_{i} 的最小公倍数,对于求解每一个的方程

先设一个 temp = M / m_{i},再求解其逆元 b_{i}

temp*b_{i}\equiv 1 (\, mod \, m_{i} )

则会有一组最小解 x=\sum_{i=1}^{n}a_{i}*b_{i}*temp_{i}

其通解就是x_{i}=x+i*M

如果没有看懂,可以看详细求解同余方程这一篇博客


代码

数论大法好,人间真善美。。。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
 
int a[1007] , m[1007] ;
int n , ans , M =1 ; 
 
int exgcd( int a , int b , int &x , int &y ) //扩展欧几里得
{
	if( b ==0 ) {
		x = 1;
		y = 0;
		return a;
	}
	else {
		int u = exgcd( b , a%b , y , x );
		y -= (a/b)*x ;
		return u;
	}
}
 
void China()
{
	for(int i = 1 ; i <= n ; ++ i ) {
		int temp = M / a[i] ;
		int x , y ;
		int gcd = exgcd( temp , a[i] , x , y );
		x = (x%a[i]+a[i])%a[i] ;
		ans = (ans + temp*x*m[i] ) % M ;
	}
}
 
int main()
{
	
	scanf("%d", &n );
	for(int i =1 ; i <= n ; ++ i ) {
		scanf("%d%d", &a[i] , &m[i] );
		M *= a[i] ;//最小公倍数
	}
	China();
	printf("%d", (ans%M+M) % M );
	return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值