扩展欧几里得
本题要点:
1、这里需要解方程 ax + by = c 的所有正整数解。
不过,有个前提,a, b, c 都是正整数, 后面方便处理很多。
2、记 d = gcd(a, b), 如果 c % d != 0, 说明方程无整数解。
3、 用 扩展欧几里得 exgcd 算出 ax + by = gcd(a, b) 的 解 x0, y0;
那么原方程 ax + by = c 的所有 整数解可以表示为
x = (c/d)x0 + k(b/d)
y = (c/d)y0 - k(a/d), 其中,k取遍所有整数。
4、 记 m1 = b/d, m2 = a/d,
x 的最小正整数为: min_x = (x % m1 + m1) % m1,
y 的最小正整数为: min_y = (y % m2 + m2) % m2;
这里算出来的 min_x, min_y 可能为0, 只需加上 m1, m2 即可。
5、当 x 取得最小的正整数, 对应的 y 最大的正整数(如果存在的话)
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
typedef long long ll;
int T;
ll a, b, c;
ll exgcd(ll a, ll b, ll& x, ll& y)
{
if