The Football Season（扩欧）

题意简化：

方程求解

 

x * w + y * d = p

x + y + z = n

其中方程解 x, y, z >= 0

 

Input

第一行输入4个数 n, p, w and d(1 <= n <= 10^{12}, 0 <= p <= 10^{17}, 1 <= d < w <= 10^{5}).

 

Output

如果方程无解, 输出-1

 

有解即按照顺序输出 x, y, z, (x,y,z均为非负数)如果有多个解, 输出任意解即可

 

x * w + y * d = p

x + y + z = n

Examples

Input

30 60 3 1

Output

17 9 4

Input

10 51 5 4

Output

-1

Input

20 0 15 5

Output

0 0 20

Note

第一个案例17 9 4 是方程一个解

 

第二个案例{10 51 5 4} , 方程无非负数解

 

题目大意：一场比赛，赢了得w分，平局得d分，如果输掉比赛，则不得分，现一共进行了n场比赛，得到了p分，w和d给出。问能否求出满足题意的解

 

 

1：赢了得w分，平局得d分，那么当赢d局和当平局w局的时候这两种情况得分是相同的。

2：d<w

 

由上面的1和2两个条件我们可以假设平局的数量为y，如果说y小于w的话，我们可以遍历y从0到w，因为w的范围是1e5，所以不会超时，当y大于等于w的时候，我们可以把平局的分数看做赢的分数，也就是赢得局面变成了w*(x+d)，这个时候x+y变成了x+d+y-w<x+y，但是总分数是不变的。

在有解的情况下，肯定存在一组解是y<w,如果y从0到w的情况下没有找到的话，就说明没有解。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

int main()
{
    ll n,p;
    int w,d,flag=1;
    cin>>n>>p>>w>>d;
    for(int i=0;i<w;++i)
    {
        int x=(p-d*i)/w;
        if(x>=0&&(p-d*i)%w==0&&x+i<=n)
        {
            flag=0;
            printf("%d %d %d\n",x,i,n-x-i);
            break;
        }
    }
    if(flag) printf("-1\n");
    return 0;
}

这道题直接可以使用扩展欧几里得来做，要求的是类似于ax+by=c的x、y的一组解，根据上面的方法可以想到x和y要尽量的小，那么问题就可以转化为求方程ax+by=c是否存在最小正整数，并且对于x和

y满足x+y<=n(题意中有给出)

 

提取问题：求是否存在最小正整数解

转化：先求一组解，通过这一组解求得最小正整数解

转化：根据扩欧公式里面的x和y求一组解

 

接下来分析过程

接触过扩展欧几里得的人就会知道扩展欧几里得函数带入四个变量ll(ll a,ll b,ll &x,ll &y)。这个式子是用来解ax+by=gcd(a,b)的，返回的值就是a和b的gcd，运算后得到的x，y就是我们要求的一组解，不过要注意这组解不是用来解ax+by=c的。

后者的解为 x0=x*c/gcd(a,b) y0=y*c/gcd(a,b)

x0和y0就是我们要求的一组解，接下来是求最小正整数解。

long long s=b/gcd(a,b)
long long xx=(x0%s+s)%s
long long yy=(c-a*xx)/b

long long s=a/gcd(a,b)
long long yy=(y0%s+s)%s
long long xx=(c-b*yy)/a

当这两个运行完后置xx和yy总存在小于0的数的话，就说明不存在最小正整数解，反之存在。

 

针对这一道题的代码：

#include <iostream>
#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
#define inf 0x3f3f3f3f
#define ll long long
#define mod 10000007
using namespace std;
ll gcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0)
    {
        x=1,y=0;
        return a;
    }
    ll r=gcd(b,a%b,x,y);
    ll temp=y;
    y=x-(a/b)*y;
    x=temp;
    return r;
}
int main()
{
    ll n,p,w,d,x,y;
    scanf("%lld %lld %lld %lld",&n,&p,&w,&d);
    ll r=gcd(w,d,x,y);
    if(p%r!=0)
        printf("-1\n");
    else
    {
        ll q=w/r;
        ll o=d/r;
        x=x%o*((p/r)%o);
        y=y%q*((p/r)%q);
//通过特解，求通解，在这里进行取余操作，防止爆long long

//        if(y<0)
//            y=y%q+q;
//        else
//            y%=q;
//直接合并
        y=(y%q+q)%q;
//直接根据wx+dy=p求x和y
        x=(p-y*d)/w;
        if(x<0||y<0||x+y>n)
            printf("-1\n");
        else
            printf("%lld %lld %lld\n",x,y,n-x-y);
    }
}