题面
【问题描述】
你有一个天平(天平左右两边都可以放砝码)与重量为 \(a,b\)(\(1\le a,b\le 10000\)) 的两种砝码。让你求出一种方案称出重为 \(c\) (\(1\le c\le 50000\))的物品,如有多种方案,请输出两种砝码需要数量的总和最小的方案。
【输入格式】
有若干行,每行三个数,\(a,b,c\)。
结束时用0 0 0表示。
【输出格式】
若干行,每行两个数,表示每个询问中 \(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;
}