模板——exCRT

#include<iostream>
using namespace std;

typedef long long ll;
const int N = 1e6 + 10;

int n;
ll b[N], p[N];

ll gcd(ll a, ll b) { return !b ? a : gcd(b, a % b); }

ll exgcd(ll a, ll b, ll& x, ll& y) {
	if (b == 0) {
		x = 1, y = 0;
		return a;
	}
	ll g = exgcd(b, a % b, y, x);
	y -= a / b * x;
	return g;
}

ll inv(ll a, ll p) {
	ll x, y; exgcd(a, p, x, y);
	return (x % p + p) % p;
}

ll exCRT() {
	for (int i = 2; i <= n; i++) {
		ll p1 = p[i - 1], p2 = p[i];
		ll b1 = b[i - 1], b2 = b[i];
		ll g = gcd(p1, p2), d = b2 - b1;
		if (d % g) return -1;
		p[i] = (p1 * p2) / g;
		b[i] = ((inv(p1 / g, p2 / g) * d / g) % (p2 / g)) * p1 + b1;
		b[i] = (b[i] % p[i] + p[i]) % p[i];
	}
	return b[n];
}
int main() {

	cin >> n;

	for (int i = 1; i <= n; i++) {
		cin >> p[i] >> b[i];
	}

	cout << exCRT() << endl;
	return 0;
}
  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值