HDU 1356 The Balance && POJ2142

重新复习了一遍扩展欧几里德（一开始学的都忘了，丧

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