当前位置: 代码迷 >> 综合 >> 模板--HihoCoder --1297--扩展欧几里得
  详细解决方案

模板--HihoCoder --1297--扩展欧几里得

热度:13   发布时间:2023-12-12 06:09:58.0

HihoCoder --1297--扩展欧几里得

求最大公倍数与不定方程的解;

LL gcd(LL a,LL b){return b?gcd(b,a%b):a;
}
void extgcd(LL a,LL b,LL &x, LL &y){if(b){extgcd(b,a%b,y,x);y-=(a/b)*x;}else  x=1,y=0;
}

完整代码:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
LL s1,s2,v1,v2,m;
LL A,B,C,D;
LL x,y;
LL gcd(LL a,LL b){return b?gcd(b,a%b):a;
}
void extgcd(LL a,LL b,LL &x, LL &y){if(b){extgcd(b,a%b,y,x);y-=(a/b)*x;}else  x=1,y=0;
}
int main(){scanf("%lld%lld%lld%lld%lld",&s1,&s2,&v1,&v2,&m);A=v1-v2;B=m;C=s2-s1;if(A<0)	A=A+m;D=gcd(A,B);if(C%D){printf("-1\n");return 0;} A=A/D;B=B/D;C=C/D;//缩小范围 extgcd(A,B,x,y);x=x*C%B;if(x<0)x=x+B;printf("%lld\n",x); return 0;
} 

 

  相关解决方案