#624285 (C++14) No.186 中華風 (Easy)

• 問題: No.186 中華風 (Easy)
• ユーザー: nonpro3
• 提出日時: 2021-03-02 05:17:04
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 2 ms / 2,000 ms
• コード長: 2,236 bytes
• コンパイル時間: 720 ms
• コンパイル使用メモリ: 78,928 KB
• 実行使用メモリ: 5,248 KB
• 最終ジャッジ日時: 2024-10-03 01:32:02
• 合計ジャッジ時間: 1,464 ms
• ジャッジサーバーID （参考情報）: judge3 / judge1

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,248 KB
testcase_02	AC	2 ms; 5,248 KB
testcase_03	AC	2 ms; 5,248 KB
testcase_04	AC	2 ms; 5,248 KB
testcase_05	AC	2 ms; 5,248 KB
testcase_06	AC	2 ms; 5,248 KB
testcase_07	AC	2 ms; 5,248 KB
testcase_08	AC	2 ms; 5,248 KB
testcase_09	AC	2 ms; 5,248 KB
testcase_10	AC	2 ms; 5,248 KB
testcase_11	AC	2 ms; 5,248 KB
testcase_12	AC	2 ms; 5,248 KB
testcase_13	AC	2 ms; 5,248 KB
testcase_14	AC	2 ms; 5,248 KB
testcase_15	AC	2 ms; 5,248 KB
testcase_16	AC	2 ms; 5,248 KB
testcase_17	AC	2 ms; 5,248 KB
testcase_18	AC	2 ms; 5,248 KB
testcase_19	AC	2 ms; 5,248 KB
testcase_20	AC	2 ms; 5,248 KB
testcase_21	AC	2 ms; 5,248 KB
testcase_22	AC	2 ms; 5,248 KB

ソースコード

#include<iostream>
#include<algorithm>
#include<vector>
#include<cmath>

using namespace std;
typedef long long ll;

ll gcd(ll a, ll b) {
	if (a % b == 0)return b;
	return gcd(b, a % b);
}

ll lcm(ll a, ll b) {
	return a * b / gcd(a, b);
}

//e_gcd()自体はaとbの最大公約数を返す。
//xとyは、a*x+b*y=gcd(a,b)となる(x,y)の1つの組を返す。
//a*x+b*y=gcdの時、a=b*p+(a%b)を代入して、b*p*x+(a%b)*x+b*y=b*(p*x+y)+(a%b)*x
//次のx->p*x+y
//次のy->x
//次のa->b
//次のb->a%b
ll e_gcd(ll a, ll b, ll& x, ll& y)
{
	if (b == 0)
	{
		x = 1;
		y = 0;
		return a;
	}
	ll gcd_res = e_gcd(b, a % b, x, y);
	ll p = a / b;
	ll nx = x, ny = y;
	x = ny;
	y = nx - p * x;
	return gcd_res;
}

std::pair<ll, ll> CRT_Simple(const vector<ll>& val, const vector<ll>& mod) {
	//x == val[i] (mod[i])がn本ある感じ
	if (val.size() != mod.size())return make_pair(-1, -1);
	int n = val.size();
	ll x = 0;//最初は0 (mod 1)で始める。

	ll mlcm = 1;

	for (int i = 0; i < n; i++) {
		//m[i] % mgcdが等しくない場合は、解がない。

		ll p, q;
		ll mgcd = e_gcd(mlcm, mod[i], p, q);//pとqの符号は逆になる
		if ((val[i] - x) % mgcd != 0) {
			return make_pair(-1, -1);
		}
		
		//実際のpを求め(e_gcdのは右辺がgcd(mlcm, mod[i])でval[i]-xと定数倍の差がある。
		//その差を補正して、mod[i] / mgcdであまりを取ると正で一番小さいpを得られる。
		//pは常に正でqは常に負として考える。
		ll min_p = (p * (val[i] - x) / mgcd) % (mod[i] / mgcd);
		if (min_p < 0)min_p += mod[i] / mgcd;
		x += min_p * mlcm;//これはpでの計算

		//以下のは現在バグってる
		//以下のようにqで求めてもいい。
		//ll max_q = (q * (val[i] - x) / mgcd) % (mlcm / mgcd);
		//if(max_q > 0)max_q -= mlcm / mgcd;
		//x = val[i] - max_q * mod[i];
		
		mlcm = lcm(mlcm, mod[i]);
	}
	x %= mlcm;
	return make_pair(x, mlcm);
}

int main() {
	ll X1, X2, X3, Y1, Y2, Y3;
	cin >> X1 >> Y1 >> X2 >> Y2 >> X3 >> Y3;

	auto ret = CRT_Simple({ X1, X2, X3 }, { Y1, Y2, Y3 });

	if (ret.first == -1) {
		cout << -1 << endl;
		return 0;
	}
	else if (ret.first == 0) {
		cout << ret.second << endl;
		return 0;
	}
	cout << ret.first << endl;

	return 0;
}