题意简化:
方程求解
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);
}
}