不平等 - 2019年9月13日バランス

フェイス質問

[問題の説明]

あなたはバランス（左右のバランスが体重をかけることができます）の重量持つ\（A、B \） （\（1 \ルA、B \ 10000ル\））2つの重量。あなたが呼ばれるスキームの重量を見つけるしてみましょう\（C \） （\（1 \ルC \ル50000 \） ）の項目、様々なプログラムは、出力重量の2種類の数の最小和を必要としてください場合プログラム。

[入力形式]

3つの数の複数の行があり、\（A、B、C \） 。

最後に000で表されます。

[出力形式]

各質問を表す2つの数値のいくつかの行、\（A \）番号と\（B \）番号

解無し出力する場合 no solution

[サンプル入力】

700 300 200
500 200 300
500 200 500
275 110 330
275 110 385
648 375 4002
3 1 10000
0 0 0

[サンプル出力]

1 3
1 1
1 0
0 3
1 1
49 74
3333 1 

ソリューション

絶対値関数パラメータの最大値を含みます。

プログラム

#include <iostream>
#include <vector>
using namespace std;
typedef long long ll;
 
#define INF 0x3f3f3f3f
 
ll exgcd(ll a, ll b, ll& x, ll& y) {
    if (b) {
        ll xp, yp;
        ll d = exgcd(b, a % b, xp, yp);
        x = yp; y = xp - (a / b) * yp;
        return d;
    } else {
        x = 1; y = 0;
        return a;
    }
}

ll f(ll k, ll x, ll b, ll d, ll y, ll a) { 
    return abs(x + b * k / d) + abs(y - a * k / d);
}
 
int main() {
    ll a, b, c;
    while (cin >> a >> b >> c && a && b && c) {
        ll d, x, y;
        d = exgcd(a, b, x, y);
        if (c % d) {
            cout << "no solution" << endl;
            continue;
        }
        x *= c / d; y *= c / d;
        // assert(a * x + b * y == c);
        vector<ll> cand;
        cand.push_back(-d * x / b);
        cand.push_back(-d * x / b + 1);
        cand.push_back(-d * x / b - 1);
        cand.push_back(d * y / a);
        cand.push_back(d * y / a + 1);
        cand.push_back(d * y / a - 1);
        ll bval = INF, xf, yf;
        for (const int& k : cand) {
            if (f(k, x, b, d, y, a) < bval) {
                bval = f(k, x, b, d, y, a);
                xf = x + b * k / d;
                yf = y - a * k / d;
            }
        }
        cout << abs(xf) << " " << abs(yf) << endl;
    }
    return 0;
}