重新复习了一遍扩展欧几里德(一开始学的都忘了,丧
http://acm.hdu.edu.cn/showproblem.php?pid=1356
这道题一眼就可以看出来是扩展欧几里德,但是将答案解出来了之后,还要找到符合他条件的答案,这样算出来的答案已经是符合条件1了,条件2是(x+y)最小,条件3是(ax+by)最小
算出来的x,y有正有负,正即放在天平左边,负即和物品放在天平右边
因为a|x|=b|y|+c,要求a|x|+b|y|最小,如果a|x|减小,(即|x|减小),因为天平两边一样重,所以|y|也要减小,当a|x|+b|y|最小时,就是|x|+|y|最小时(可能不是很严谨,这也是我脑糊的
于是题意就变成了求|x|+|y|的最小值了。x 和 y 是不定式 ax+by=c 的解。可由扩展欧几里德算出通式
x=x1+b/d*t ;
y=y1-a/d*t ;
|x|+|y|=|x0+b/d*t|+|y0-a/d*t|,我们规定 a>b(不行就交换),从这个式子中,我们可以发现:|x0+b/d*t|是单调递增的,|y0-a/d*t|是单调递减的,而由于我们规定了 a>b,那么减的速度是要大于增的速度,于是整个函数减少的要比增加的快,但是由于绝对值的符号的作用,最终函数还是递增的。也就是说,函数是凹的,先减小,再增大。那么什么时候最小呢?很显然是 y0-a/d*t==0 的时候,于是我们的最小值|x|+|y|也一定是在 t=y0*d/a附近了,只要在附近枚举几个值就能找到最优解了(我枚举了5个,应该判断一下2个就行了,但是我懒,不想想了。。。。。
//#include<bits/stdc++.h>
#include<iostream>
#include<cstdio>
#include<queue>
#include<stack>
#include<map>
#include<string>
#include<cstring>
#include<limits.h>
using namespace std;
typedef long long ll;
const ll INF=LONG_LONG_MAX;
const ll inf=LONG_LONG_MIN;
ll a1,b1,d1,x,y,r;
void exgcd(ll a,ll b){
if(b){
exgcd(b,a%b);
ll k=x;
x=y;
y=k-a/b*y;
}
else{
r=a;
x=1;
y=0;
}
}
bool flag=0;
int main(){
while(scanf("%lld%lld%lld",&a1,&b1,&d1)&&a1||b1||d1){
flag=0;
if(a1<b1){
swap(a1,b1);
flag=1;
}
exgcd(a1,b1);
if(d1%r!=0){
cout<<"no solution"<<endl;
}
else{
ll x1=x*(d1/r),y1=y*(d1/r);
//通解 x=x1+b1/r*t , y=y1-a1/r*t ,其中 t 为整数
ll resx=x1,resy=y1;
ll mid=y1*r/a1;
for(int i=mid-2;i<=mid+2;++i){
ll xx=x1+b1/r*i;
ll yy=y1-a1/r*i;
if(abs(xx)+abs(yy)<abs(resx)+abs(resy)){
resx=xx,resy=yy;
}
}
if(!flag){
cout<<abs(resx)<<' '<<abs(resy)<<endl;
}
else{
cout<<abs(resy)<<' '<<abs(resx)<<endl;
}
}
}
return 0;
}
附上我学习扩展欧几里德的网址https://www.cnblogs.com/hadilo/p/5914302.html