#814009 (C++17) No.2117 中国剰余定理入門

• 問題: No.2117 中国剰余定理入門
• ユーザー: kiliarin
• 提出日時: 2022-11-04 21:48:35
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 2 ms / 2,000 ms
• コード長: 593 bytes
• コンパイル時間: 172 ms
• コンパイル使用メモリ: 32,256 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-07-18 19:32:49
• 合計ジャッジ時間: 788 ms
• ジャッジサーバーID （参考情報）: judge3 / judge2

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

ソースコード

#include <stdio.h>
#include <string.h>
#include <stdlib.h>

#define LOG(FMT...) fprintf(stderr, FMT)

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

int main() {
	int b1, c1, b2, c2;
	scanf("%d%d%d%d", &b1, &c1, &b2, &c2);
	int x, y, d;
	d = exgcd(b1, b2, x, y);
	if ((c1 - c2) % d) {
		puts("NaN");
	} else {
		int t = (c1 - c2) / d;
		x *= t, y *= t;
		int g = b1 * b2 / d;
		int m = c2 + y * b2;
		int k = -m / g;
		if (k * g < -m) {
			k += 1;
		}
		printf("%d\n", k * g + m);
	}
	return 0;
}