POJ2142 The Balance

问题描述
Iyo姐姐是一个医学系大学生，她需要用一些不可名状之砝码+一个天枰测量一些不可名状之药。如图示她可以把砝码放在两边的秤盘（我怎么觉得这不符合实验室原则）进行称量。现在她需要知道，如何用最少的砝码称量指定重量的药品。

输入
多组数据。一行三个整数a、b、c，分别表示两种砝码的质量和指定的药品质量。输入0 0 0表示输入结束。

输出
对应每一行输入需要分别输出a砝码和b砝码的重量。且不应输出多余的格式符。

样例输入
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

代码

//By Chorolop
#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;

int abss(int a){
    return a>0?a:-a;
}

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

int gcd(int x,int y){
    return y==0?x:gcd(y,x%y);
}

int main()
{
    int a,b,c,cc = 0;
    scanf("%d%d%d",&a,&b,&c);
    while(a && b && c){
        int x,y,d;
        d = gcd(a,b);
//      if(c%d){printf("Error!\n");scanf("%d%d%d",&a,&b,&c);continue;}

        a /= d;
        b /= d;
        c /= d;

        ex_gcd(a,b,d,x,y);

        int vx = x*c; //此处的求的是x最接近0的特解
        vx = (vx%b + b) %b;
        int vy = (c - a*vx) /b;
        vy = vy<0?(-vy):vy;

        y = y*c; //此处求的是y最接近0的特解
        y = (y%a + a) %a;
        x = (c - b*y) /a;
        x = x<0?(-x):x;

//      if(cc++) printf("\n");

        if(x+y > vx+vy) x = vx,y = vy; //两相比较取其轻

        printf("%d %d\n",x,y);

        scanf("%d%d%d",&a,&b,&c);
    }
}

评价
一道标准的扩展欧几里德，我们可以直接把这个看作ax+by==c，那么这道题其实是求给定丢番图方程的一个特解。这个特解需要满足这样的条件：|x|+|y|达到最小。x和y都可以取负，因为显然他们取负时表示他们所对应的砝码放在天枰右边的托盘。值得注意的是这道题比较关注格式。