中国剩余定理(非互质)

awing 204. 表达整数的奇怪方式

题意

给定2n个整数 a 1 , a 2 , … , a n a_1,a_2,…,a_n a1,a2,,an m 1 , m 2 , … , m n m_1,m_2,…,m_n m1,m2,,mn,
求一个最小的非负整x,满足 ∀i ∈ [ 1, n ], x ≡ m i ( % a i ) x ≡ m_i (\% \quad a_i) xmi(%ai)

题解

参考题解
先浏览参考题解,再看下面图片中的题解。

在这里插入图片描述

在这里插入图片描述
对于第一组方程的解,其实m[1]即为最小的解,lcm = a[1].
也即让x的初始值为m[1], lcm的初始值为a[1]即可。

我自己的模板如下:

#include<cstdio>
#include<queue>
#include<set>
#include<cstdlib>
#include<string.h>
#include<string>
#include<iostream>
#include<cmath>
#include<map>
#include<algorithm>
#define endl "\n"
#define IOS ios::sync_with_stdio(0), cin.tie(0),cout.tie(0)
#define ft first
#define sd second
#define pll pair<ll, ll>
#define pii pair<int, int>
#define ll long long int
#define mt(a,b) memset(a, b, sizeof a)
//#define int long long
const double PI = acos(-1.0);
const int inf = 0x3f3f3f3f;
const int INF = 0x7fffffff;
using namespace std;
const int N = 2e5 + 7, M = 1e6 + 10;
ll a[100], m[100];

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

int main()
{
	IOS;
	ll n; cin >> n;
	for (int i = 1; i <= n; i++)
		cin >> a[i] >> m[i];

	ll x = m[1], lcm = a[1];

	int f = 1;
	for (int i = 2; i <= n; i++)
	{
		ll c = m[i] - x;
		ll t, y;
		ll d = exgcd(lcm, a[i], t, y);
		if (c % d)
		{
			cout << "-1\n";
			return 0;
		}
		t *= c / d;
		ll mod = a[i] / d;
		t = (t % mod + mod) % mod;
		x = x + t * lcm;
		lcm = lcm / d * a[i];
	}
	cout << x  << endl;


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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

to cling

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值