斩杀线计算大师

链接：https://ac.nowcoder.com/acm/contest/4853/D
来源：牛客网

题目描述
算术能力是每个炉石玩家必不可少的，假设现在有三种伤害卡，伤害值分别是a,b,ca,b,ca,b,c。并且每种伤害卡的数量你可以认为是无限的。现在牛牛想知道是否存在一种方式可以刚好造成kkk点伤害，输出x,y,zx,y,zx,y,z分别表示三种伤害卡的使用个数。
数据保证一定存在解。如果存在多组解，输出任意一组。

输入描述:
一行四个整数分别表示a,b,c,ka,b,c,ka,b,c,k.

输出描述:
一行输出三个整数分别表示x,y,zx,y,zx,y,z.

示例1
输入

3 4 5 20

输出

4 2 0

备注:
1≤a,b,c≤1e51 \leq a, b, c \leq 1e51≤a,b,c≤1e5
0≤k≤1e120 \leq k \leq 1e120≤k≤1e12

ax+by+cz=k可以写成ax+by=k-cz,令ccc=k-z*c
所以就是扩展欧几里得定理了

#include<iostream>
#include<cstring>
using namespace std;
typedef long long LL;

LL exgcd(LL a, LL b, LL &x, LL &y)
{
    if (!b)
    {
        x = 1, y = 0;
        return a;
    }
    LL d = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return d;
}

int main()
{
    LL a,b,c,k;
    cin>>a>>b>>c>>k;
    for(LL z=0;z<k/c;z++)
    {
        LL ccc=k-z*c;
        LL x,y,d;
        d=exgcd(a,b,x,y);
        if(ccc%d) continue;
        x=x*ccc/d;y=y*ccc/d;  扩大 ccc/ gcd 倍
        x = (x % (b / d) + (b / d)) % (b / d);//计算最小正整数
        y = (ccc - x * a) / b;
        if(x >= 0 && y >= 0)
        {
            cout<<x<<" "<<y<<" "<<z<<endl;  
            break;
        }
    }
    return 0;
}