斩杀线计算大师-----------------------数论(扩展欧几里得)

最新推荐文章于 2021-10-04 19:59:59 发布

AKone123456

最新推荐文章于 2021-10-04 19:59:59 发布

阅读量237

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/qq_43690454/article/details/105590642

版权

数论专栏收录该内容

50 篇文章 1 订阅

订阅专栏

在这里插入图片描述

解析:
ax+by+cz=k
ax+by=k-cz
这不就是形如ax+by=c
所以用扩展欧几里得求解就行了

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a,b,c,k;
ll x,y;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	if(!b)
	{
		x=1;y=0;
		return a;
	}
	ll d=exgcd(b,a%b,x,y);
	ll tmp=x;
	x=y;
	y=tmp-a/b*y;
	return d;
}
int main()
{
	 scanf("%lld%lld%lld%lld",&a,&b,&c,&k);
	for(int i=0;i<k/c;i++)
	{
		ll d=exgcd(a,b,x,y);
		ll p=k-c*i;
		if(p%d) continue;
		x=x*(p/d);//扩大若干倍
		y=y*(p/d);//扩大若干倍
		x=(x%(b/d)+(b/d))%(b/d); //最小正整数解
		y=(p-a*x)/b; //ax+by=p y=(p-a*x)/b
		if(x>=0&&y>=0)
		{
			printf("%lld %lld %lld\n",x,y,i);
				break;
		}
	
	}

 }