题解 P3951 【小凯的疑惑】

infinityedge

2017-11-19 21:37:12

Solution

#### 说明 这个做法是我在考场上想出来的 ex\_gcd 做法,当时并没有想到打表,所以向发表一下让大家彻底摆脱小学奥数带来的阴影,我现在仍然认为,ex\_gcd 是这道题的标算。 #### 题解 首先,我们发现这两个数是**互质**的,并且有**无限**个。很容易想到不定方程 $ax + by = gcd(a, b)$ ,其中 $gcd(a, b) = 1$ 。 然后我们考虑,对于所有可行的能被 $a$ 和 $b$ 表示出来的数 $k$ ,都存在 $x \geq 0, y \geq 0,ax + by = k$。 现在我们要构造的是**最大的不合法的数**,显然,这个数 $+ 1$ 一定是一个合法的数,转化成了求**最大的减一后不合法的数**。 由于这个数 $k$ 一定是合法的,所以满足性质 $$\exists x \geq 0, \exists y \geq 0,ax + by = k$$ 那么如果 $k-1$ 合法,那么 $k - 1$ 可以表示成 $$a \left ( x - x' \right ) + b \left ( y - y' \right ) = k$$ 或 $$a \left ( x - x'' \right ) + b \left ( y - y'' \right ) = k$$ 其中 $x',y'$ 表示 $ax + by = 1$ 的 $x$ 最小且非负的整数解; $x'',y''$ 表示 $ax + by = 1$ 的 $y$ 最小且非负的整数解。 那么现在只需要让 $x - x'< 0$ 并且 $y - y'' < 0$ 即可 那么最后的**最大的减一后不合法的数**就是 $$a \left ( x' - 1 \right ) + b \left ( y'' - 1 \right ) $$ 那么最后的**答案**就是 $$a \left ( x' - 1 \right ) + b \left ( y'' - 1 \right ) - 1$$ #### 证明: 首先充分性成立。 然后证明必要性:若$a \left ( x' - 1 \right ) + b \left ( y'' - 1 \right ) $不是**最大的减一后不合法的数**,那么一定存在一个更大的数,显然该数的 $a$ 的系数大于 $x' - 1$ 或 $b$ 的系数大于 $y'' - 1$ (如果都小于等于,那么该数不会比当前数大)。显然,减一后仍然是合法的。所以必要性成立。 综上, $a \left ( x' - 1 \right ) + b \left ( y'' - 1 \right ) - 1$ 是最大的不合法的数。 #### 代码: ```cpp #include <iostream> #include <cstring> #include <cstdio> #include <algorithm> #include <cmath> using namespace std; typedef long long ll; ll gcd(ll a, ll b){ return b == 0 ? a : gcd(b, a % b); } void ex_gcd(ll a, ll b, ll &x, ll &y){ if(b == 0){ x = 1, y = 0; return; } ex_gcd(b, a % b, y, x); y -= (a / b) * x; } ll a, b; int main(){ cin >> a >> b; if(a > b) swap(a, b); ll x, y; ex_gcd(a, b, x, y); if(x > 0){ swap(a, b); swap(x, y); } ll tmp = (-x) / b; x = x + tmp * b; y = y - tmp * a; while(x < 0) x = x + b, y = y - a; while(x > 0) x = x - b, y = y + a; ll ans; ll xx2 = x + b; ans = a * (xx2 - 1) + b * (y - 1); cout << ans - 1 << endl; return 0; } ```