poj2142

扩展欧几里得，注意细节

（“活着”的概念不是等着慢慢死去，而是要不断地奔跑，跑到很远的地方去看尽可能广大的世界，跑到筋疲力尽才不会后悔。——江南 《龙族》）

2015.10.30：

（错了就错了吧，一个做错的英雄，至少比什么都不做的笨蛋好(by昂热)     ——江南 《龙族III黑月之潮》）

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;

void exgcd(long long int *x,long long int *y,long long int a,long long int b,long long int *z){
    if(b==0){
        *x=1;
        *y=1;
        *z=a;

        return;
    }
    else{
        long long int x1,y1,z1;

        exgcd(&x1,&y1,b,a%b,&z1);
        *z=z1;
        *x=y1;
        *y=x1-a/b*y1;

        return;
    }
}

int main(){
    long long int a,b,d;
    long long int x,y,z;
    long long int ansx,ansy,ansxy,ansaxby;

    while(scanf("%lld%lld%lld",&a,&b,&d)!=EOF&&!(a==0&&b==0&&d==0)){
        /*if(a==0||b==0||d==0){
            printf("%d\n",1/0);
        }printf("%d\n",1/0);
        */

        exgcd(&x,&y,a,b,&z);

        x=x*d/z;
        y=y*d/z;

        long long int temp,tempx,tempy;

        tempx=0;
        if(x>0){//让x<=0
            tempx=x/(b/z);
            if(x%(b/z)){
                tempx++;
            }
        }
        if(y>0){
            tempy=-1;
        }
        else{//让y>=0;
            tempy=(-y)/(a/z);
            if((-y)%(b/z)){
                tempy++;
            }
        }
        temp=max(tempx,tempy);
        x=x-temp*(b/z);
        y=y+temp*(a/z);

        temp=min((-x)/(b/z),y/(a/z));
        x=x+temp*(b/z);
        y=y-temp*(a/z);
        ansx=-x;
        ansy=y;
        //printf("%lld %lldha1\n",ansx,ansy);
        ansxy=-x+y;
        ansaxby=(-x)*a+b*y;

        tempx=(-x)/(b/z);
        if((-x)%(b/z)){
            tempx++;
        }

        tempy=y/(a/z);
        if(y%(a/z)){
            tempy++;
        }

        temp=max(tempx,tempy);
        x=x+temp*(b/z);
        y=y-temp*(a/z);

        if(x-y<ansxy){
            ansx=x;
            ansy=-y;
            //printf("%lld %lldha2\n",ansx,ansy);
            ansxy=x-y;
            ansaxby=x*a+(-y)*b;
        }
        else if(x-y==ansxy&&x*a+(-y)*b<ansaxby){
            ansx=x;
            ansy=-y;
            //printf("%lld %lldha3\n",ansx,ansy);
            ansxy=x-y;
            ansaxby=x*a+(-y)*b;
        }

        tempy=(-y)/(a/z);//还有x，y都大于0时的情况
        if((-y)%(a/z)){
            tempy++;
        }

        if(x-tempy*(b/z)>=0){
            x=x-tempy*(b/z);
            y=y+tempy*(a/z);

            if(a<b){
                tempx=x/(b/z);
                x=x-tempx*(b/z);
                y=y+tempx*(a/z);//这里把a写成了b
            }

            if(x+y<ansxy){
                ansx=x;
                ansy=y;
                ansxy=x+y;
                ansaxby=x*a+y*b;
            }
            else if(x+y==ansxy&&x*a+y*b<ansaxby){
                ansx=x;
                ansy=y;
                ansxy=x+y;
                ansaxby=x*a+y*b;
            }
        }

        printf("%lld %lld\n",ansx,ansy);
    }

    return 0;
}