C Looooops [exgcd]

最新推荐文章于 2021-11-18 13:52:17 发布

原创最新推荐文章于 2021-11-18 13:52:17 发布 · 151 阅读

0 ·

CC 4.0 BY-SA版权

扩展欧几里得专栏收录该内容

4 篇文章

订阅专栏

本文深入探讨了扩展欧几里得算法的原理及应用，通过实例代码展示了如何求解线性同余方程组，特别关注于解决特定条件下的最小非负解。文章详细解释了算法步骤，并提供了完整的C++实现，适用于解决数学和计算机科学中的多项式复杂度问题。

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

传送门

$t*C-y*2^k=B-A$

#include<bits/stdc++.h>
#define LL long long
using namespace std;
int A,B,C,k; LL x,y,gcd;
void exgcd(LL a,LL b){
	if(!b){x=1,y=0,gcd=a; return;}
	exgcd(b,a%b); LL x1=y,y1=x-a/b*y; x=x1,y=y1;
}
int main(){
	while(1){
		scanf("%d%d%d%d",&A,&B,&C,&k);
		if(A==0) break;
		if(A==B){printf("0\n"); continue;}
		LL b=(LL)1<<k,c=B-A,a=C;
		exgcd(a,b); if(c%gcd){printf("FOREVER\n"); continue;}
		x *= c/gcd; x %= (b/gcd); x = (x+b/gcd) % (b/gcd); printf("%lld\n",x);
	}return 0;
}