POJ 2142 && HDU 1356 The Balance（一元线性同余方程）

Description
给出两种质量的砝码，和一件物品，问要多少砝码能称出物品，砝码的个数满足两种砝码的数量和尽量小，如果还有多解满足两种砝码的总质量最少
Input
多组用例，每组用例包括三个整数a，b和d分别表示两种砝码的质量以及物品的质量，以0 0 0结束输入
Output
对于每组用例，输出能称出物品质量的砝码个数，砝码的个数满足两种砝码的数量和尽量小，如果还有多解满足两种砝码的总质量最少
Sample Input
700 300 200
500 200 300
500 200 500
275 110 330
275 110 385
648 375 4002
3 1 10000
0 0 0
Sample Output
1 3
1 1
1 0
0 3
1 1
49 74
3333 1
Solution
题意即为求解方程ax+by=d的解x，y，使其满足x,y均最小且x+y最小,我们可以先通过求解一元线性同余方程ax=d (mod b)得到满足条件的最小正整数x值，然后将其带入原方程解出y值，之后再通过解by=d (mod a)得到满足条件的最小正整数y值，然后将其带入方程解出x值，两个x+y值作比较即得到满足条件的解
Code

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector> 
using namespace std;
typedef long long ll;
ll extend_gcd(ll a,ll b,ll & x,ll & y)//扩展欧几里得,求ax+by=gcd(a,b) 
{ 
    ll d=a;
    if(b!=0)
    {
        d=extend_gcd(b,a%b,y,x);
        y-=(a/b)*x;
    }
    else 
    {
        x=1;
        y=0;
    }
    return d;
}
ll linear(ll a,ll b,ll c)//一元线性同余方程ax=b(mod c) 
{
    ll x,y;
    ll g=extend_gcd(a,c,x,y);
    if(b%g)
    {

        return -1;
    }
    x=x*(b/g);
    ll mod=c/g;
    x=(x%mod+mod)%mod;
    return x;
}
int main()
{
    ll a,b,d;
    while(scanf("%lld%lld%lld",&a,&b,&d),a||b||d)
    {
        ll x1=linear(a,d,b);
        ll y1=(d-a*x1)/b;
        if(y1<0)
            y1=-y1;
        ll y2=linear(b,d,a);
        ll x2=(d-b*y2)/a;
        if(x2<0)
            x2=-x2;
        if(x1+y1>x2+y2)
            x1=x2,y1=y2;
        printf("%lld %lld\n",x1,y1);
    }
    return 0;
}