POJ 2115 扩展欧几里得

MakingMaker

于 2012-10-14 00:03:55 发布

阅读量505

点赞数

文章标签：扩展 c

本文链接：https://blog.csdn.net/MakingMaker/article/details/8068678

版权

本文介绍了如何使用扩展欧几里得算法解决特定形式的线性方程，并提供了代码实现。重点在于理解算法原理及其在实际问题中的应用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

解线性方程：2^k * x + c*y=b-a;

#include<cstdio>
using namespace std;
typedef long long ll;
ll a,b,c,d,e,k,x,y;

ll Extended_Euclid(ll a,ll b)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	ll m=Extended_Euclid(b,a%b);
	ll temp=x;
	x=y;
	y=temp-a/b*y;
	return m;
}

void cal()
{
	d=1LL<<k;
	e=Extended_Euclid(c, d);
	if((b-a)%e!=0)
	{
		printf ("FOREVER\n");
		return ;
	}
	ll t=d/e;
	x=( x*(b-a)/e )%t+t;
	x%=t;
	printf("%lld\n",x);
}

int main()
{

	while(~scanf("%lld%lld%lld%lld",&a,&b,&c,&k))
	{
		if(a+b+c+k==0)
		{
			break;
		}
		cal();
	}
	return 0;
}