倒酒（拓展欧几里得）

// 对gcd的理解...

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int gcd(int x,int y){
return !y?x:gcd(y,x%y);
}
int exgcd(int a,int b,int& x,int& y){
if(b==0){
x=1; y=0;
return a;
}
int res=exgcd(b,a%b,x,y);
int t=x;
x=y;
y=t-(a/b)*y;
return res;
}
int main(){
int a,b,x,y;
scanf("%d%d",&a,&b);
int mod=gcd(a,b);
cout<<mod<<endl;
a/=mod; b/=mod;
exgcd(a,b,x,y);
while(x>0||y<0) x-=b,y+=a;
while(x+b<=0&&y-a>=0) x+=b,y-=a;
printf("%d %d\n",-x,y);
return 0;
}

倒酒（拓展欧几里得）相关推荐

poj1061-青蛙的约会(拓展欧几里得java)
题目: 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚 ...
拓展欧几里得模板/求逆元模板（java）
拓展欧几里得模板参考:哈尔滨理工大学ACM培训资料汇编/ACM-ICPC培训资料汇编* 基本原理 :设 a 和 b 不全为 0,则存在整数 x,y 使得 xa yb=gcd(a,b)=c 对于辗转相 ...
杭电2669拓展欧几里得
杭电2669 给a,b求Xa Yb = 1.如果没有则输出sorry. 可以通过拓展欧几里得指导Xa Yb = gcd(a,b). 不言而喻要判断gcd(a,b)是否等于1.如果不等于1,那么就是so ...
A/B HDU - 1576 （逆元或拓展欧几里得或数学公式）多解法求大数结果
题意:求(A/B)%9973,但由于A很大,我们只给出n(n=A%9973)(我们给定的A必能被B整除,且gcd(B,9973) = 1). 思维:(1)逆元+扩展欧几里得算法:满足a*k≡1 (mo ...
拓展欧几里得小结（初级理解）
什么是拓展欧几里得?简单的说,就是求关于x,y的方程 ax + by = gcd(a,b) 的所有整数解现在我们令g = gcd(a,b)则方程变成了ax + by = g 假如我们现在知道了关于这 ...
拓展欧几里得 [Noi2002]Savage
对于一个野人,他(她?)所在的位置,(C[i]+x*p[i])%ans,是的,暴力枚举每一个ans,用拓展欧几里得求出每两个wildpeople(wildrage?)相遇的年份,如果小于最小的寿限(就 ...
AcWing 878. 线性同余方程（拓展欧几里得）
题目链接 https://www.acwing.com/problem/content/880/ 思路 ai×xi≡bi(modmi)a_i\times x_i≡b_i(mod \ m_i)ai×x ...
数学--数论--欧几里得定理和拓展欧几里得定理
欧几里得定理: gcd(a, b) = gcd(b, a%b) 证明: 我们首先约定:m = gcd(a,b) , n = gcd(b, q) , a = b*p +q.(这里的gcd含义跟上面一样, ...
POJ1061青蛙的约会（拓展欧几里得）
青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 146847 Accepted: 34169 Description ...
从欧几里得到拓展欧几里得
欧几里得定理: 公式表述 gcd(a,b)=gcd(b,a mod b) 证明:a可以表示成a = kb + r,则r = a mod b 假设d是a,b的一个公约数,则有 d|a, d|b,而r ...

倒酒（拓展欧几里得）

倒酒（拓展欧几里得）相关推荐

最新文章

热门文章