poj 2142

http://hi.baidu.com/5l2_/blog/item/a0381951b5738a1e367abef0.html 这个文章写的很好。。。我只是照办人家的。。。。TT 自己想不到的。。。。。

大概意思 给定 a b k找到满足ax+by=k 的令|x|+|y|最小（等时令a|x|+b|y|最小）不妨a 〉b

先用扩展欧几里得算法求出 一组解 x0，y0，通解可以表示为x=x0+b/d *t y=y0-a/d *t

|x|+|y|=|x0+b/d *t |+|y0-a/d *t| 这个关于t的函数的最小值在  t  =  y0*d/a  附近的两整点里取。故直接验证这两点即可。

因为   设a>b之后

|x0+b/d *t| 单调递增，|y0-a/d*t| 先递减再递增。因斜率a/d>b/d，所以总的|x0+b/d *t |+|y0-a/d *t| 先递减再递增，使y0-a/d*t0=0 的t0附近有最小值。

#include<iostream>
using namespace std;
int x, y, a1, b1;

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

int calx(int t) {
    return abs(x + a1 * t);
}

int caly(int t) {
    return abs(y - b1 * t);
}

int main() {
    int a, b, d, k;
    int tx1, tx2, ty1, ty2, r1, r2, r, ansx, ansy;
    bool flag = 1;
    while (scanf("%d%d%d", &a, &b, &k) && (a || b || k)) {
        flag = 1;
        if (a < b) flag = 0, swap(a, b);
        d = ex_gcd(a, b, x, y);
        x = x * (k / d);
        y = y * (k / d);
        a1 = b / d, b1 = a / d, r1 = y / b1;
        if (r1 * b1 - y >= 0) r1--;
        r2 = r1 + 1;
        tx1 = calx(r1),tx2 = caly(r1);
        ty1 = calx(r2),ty2 = caly(r2);
        if ((tx1 + tx2) < (ty1 + ty2)) r = r1;
        else if ((tx1 + tx2) > (ty1 + ty2)) r = r2;
        else {
            if ((tx1 * a + tx2 * b) <= (ty1 * a + ty2 * b)) r = r1;
            else r = r2;
        }
        ansx = calx(r);
        ansy = caly(r);
        if (!flag)
            printf("%d %d\n", ansy, ansx);
        else
            printf("%d %d\n", ansx, ansy);
    }
    return 0;
}