51nod 1079 中国剩余定理模板

中国剩余定理就是同余方程组除数为质数的特殊情况

我直接用同余方程组解了。

记得exgcd后x要更新

还有先更新b1再更新m1,顺序不能错!!(不然会影响到b1的更新)

#include<cstdio>
#include<cctype>
#define REP(i, a, b) for(int i = (a); i < (b); i++)
#define _for(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;

typedef long long ll;
void exgcd(ll a, ll b, ll& d, ll& x, ll& y)
{
	if(!b) { d = a; x = 1; y = 0; }
	else { exgcd(b, a % b, d, y, x); y -= x * (a / b); }
}

void read(ll& x)
{
	ll f = 1; x = 0; char ch = getchar();
	while(!isdigit(ch)) { if(ch == '-1') f = -1; ch = getchar(); }
	while(isdigit(ch)) { x = x * 10 + ch - '0'; ch = getchar(); }
	x *= f;
}

int main()
{
	ll n, m1, b1, m2, b2;
	read(n); read(m1); read(b1);
	_for(i, 2, n)
	{
		read(m2); read(b2);
		ll A = m1, B = m2, K = b2 - b1, d, x, y;
		exgcd(A, B, d, x, y);
		x = (x * (K / d) % (B / d) + (B / d)) % (B / d);
		b1 = m1 * x + b1;
		m1 = m1 / d * m2; 
	}
	printf("%lld\n", b1);
	return 0;
}

 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值