HDU 1356 The Balance

给出 两种砝码的重量 和要测量的重量
给定 a b c找到满足ax+by=c 的令|x|+|y|最小（等时令a|x|+b|y|最小）不妨a>b
先用扩展欧几里得算法求出 一组解 x0，y0
通解可以表示为x=x+k*b’ *t y=y-k*a’ 其中a’=a/d , b’=b/d;
|x|+|y|=|x+b/d *k |+|y-a/d *k| 这个函数的最小值
因为a b k x均为整数
故求|y-a/d *k|的最小值，分析可知会在y*d/a附件找到最值

我是在-x*d/b到y*d/a区间内找的最小值，水了点

#include <iostream> #include <cstdio> #include <cstdlib> using namespace std; //扩展欧几里得,返回a，b的gcd；求得x，y使ax+by=d int extend_gcd(int a,int b,int &x,int &y) { int t,m; if((b==0)&&(a==0)) return -1; if(b==0) { x=1;y=0; return a; } else { m=extend_gcd(b,a%b,x,y); t=x;x=y;y=t-a/b*y; } return m; } int main(void) { int a,b,x0,y0,c,d,x,y,t,s,sum,sum2,rx,ry,tx,ty; bool flag; while(cin>>a>>b>>c) { if(a==0&&b==0&&c==0) break; flag=false; if(a<b) { swap(a,b); flag=true; } d=extend_gcd(a,b,x0,y0); sum=sum2=1000000000; a/=d; b/=d; x=c/d*x0; y=c/d*y0; s=-x/b-1; t=y/a+1; for(int i=s;i<=t;++i) { if(abs(x+b*i)+abs(y-a*i)<sum||(abs(x+b*i)+abs(y-a*i)==sum&&a*abs(x+b*i)+b*abs(y-a*i)<sum2)) { sum=abs(x+b*i)+abs(y-a*i); sum2=a*abs(x+b*i)+b*abs(y-a*i); rx=abs(x+b*i); ry=abs(y-a*i); } } if(flag) swap(rx,ry); cout<<rx<<" "<<ry<<endl; } return 0; }