每组数据有两个数n(0usingnamespacestd;voidexgcd(inta,intb,int&x,int&y)//扩展gcd，可以求出gcd(a,b)以及ax+by=gcd(a,b)中x

#includeusingnamespacestd;voidexgcd(inta,intb,int&x,int&y){if(b==0){x=1;y=0;}else{exgcd(b,a%b,y,x);y-

voidexgcd(inta1,intb,int&x,int&y){if(b==0){x=1;y=0;return;}exgcd(b,a1%b,x,y);intt=x;x=y;y=t-(a1/b)*y;

分析：非常板的容斥题。。。考场上时间多点应该还是写得出来的。。。转换一下题目，就是求满足x≡1(modai)且x≡0(modaj)x≡1(modai)且x≡0(modaj)的x的个数（x≤Nx≤N）。由于N非常大，无法判断求解，只能算贡献。就是对于一对数ai,ajai,aj求A∗ai=B∗aj+1A∗ai=B∗aj+1中A的解的个数。然而算贡献最大的弊病就是会重复，所以要去重。然后就可以用容斥原理

l（跳了y1圈）列出方程：x+m*x1-(y+n*x1)=y1*l;化简为（n-m）*x1+y1*l=x-y;代码：#include#includeusingnamespacestd;longlongexgcd

a*x1+b*y1=1;代码：#include#includeusingnamespacestd;longlongexgcd(longlonga,longlongb,longlong&x,longlong

如果c不整除最大公约数，则一定没有整数解~#include#include#includeusingnamespacestd;inta,b,c,x,y;intexgcd(inta,intb,int&x,

移项后有ax-km=b，转化到了EXGCD的解法。

intgcd(inta,intb){if(a%b==0)returnb;elsegcd(b,a%b);}intmain(){inta,b;cin>>a>>b;coutusingnamespacestd;intexGcd

includeusingnamespacestd;intmod=17,inv[30];voidinit(){inv[1]=1;for(inti=2;i#defineintlonglongusingnamespacestd;intb,a;voidexgcd

C.TheFootballSeason题意：给定四个数字n,p,w,dn,p,w,dn,p,w,d,要求求出x,yx,yx,y，使得x+y#definelllonglongusingnamespacestd;llexgcd

描述两只青蛙在网上相识了，它们聊得很开心，于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上，于是它们约定各自朝着对方那里跳，直到碰面为止。可是它们出发之前忘记了一件很重要的事情，既没有问清楚对方的特征，也没有约定见面的具体位置。不过青蛙们都是很乐观的，它们觉得只要一直朝着某个方向跳下去，总能碰到对方的。但是除非这两只青蛙在同一时间跳到同一点上，不然是永远都不可能碰面的。为了帮助这两只

1.逆元是个什么东西首先这里有个式子:(((ab\frac{a}{b}ba)))modmodmodppp,这个式子的答案怎么求呢？？A直接求不行吗？Q别想了不能。我们发现一个分数mod一个整数时是不能直接模运算的，但是可以进行乘法运算，我们就要用到逆元（大概可以看做一个数%p意义下的倒数吧）除以一个数模ppp，就等于乘以一个数在modp意义下的逆元也就是说，它的定义可以表示为ax≡1ax≡1ax≡

记或(定义了剩余系中的除法)求逆元的方法：1.扩展欧几里得ax≡1(modp)可以等价的转化为ax−py=1然后套用exgcd解方程，并检查gcd(a,p)是否等于1如果gcd(a,p)=1，把x调整到

是建立在扩展欧几里得上的算法用来求线性模方程的解即给出许多模数mi，然后给出余数ai，求原数x代码如下structChina{#defineN5LLA[N],p[N],T[N],m[N],T1[N],M;intn;voidexgcd

乘法逆元首先，假设你会扩展欧几里得（exgcd）的话，那么乘法逆元就很简单了。首先看一下定义：在(modp)意义下（p是素数），如果a*a'=1，那么我们就说a'是a的逆元。

每组数据有两个数n(0usingnamespacestd;typedeflonglongll;voidexgcd(lla,llb,ll&x,ll&y)//x,y有取地址符号，所以可以在主函数直接用//该函数还可以同时求出

举个例子，DeepinC大神说自己T1忘了exgcd，证了一个多小时，T2忘了CRT，证了一个小时，T3看出来个结论，证了半小时。然后一出分，80+40+10，我们姑且不算夸大的成分（。。。。。。）

说明：1.模板中maxn表示最大数据规模，可以用定义，其中数为数值2.对于含有模板的模板，用类似于STL中的map，bitset的方法定义数学模块扩展欧几里得算法说明：用于计算方程其中的一组解voidexgcd

目录引言欧几里得算法（gcd）代码模板：扩展欧几里得（exgcd）代码模板：通解形式为：最小正整数解：引言什么是扩展欧几里得，听起来好高深，别急先从欧几里得下手；欧几里得算法（gcd）欧几里得算法的用处

a:gcd(b,a%b));}intexgcd(inta,intb,int&x,int&y){intd=a;if(b!

gcd(b,a%b):a;}LLlcm(LLa,LLb){LLc=gcd(a,b);returna/c*b;}例题：hrbust1178hdu2028LowestCommonMultiplePlus二、exgcd

//exgcd,注意最后面求出最小解的公式...还有就是判断有没有解，还有就是把正负颠倒一下...

includeusingnamespacestd;intn,A,B;intans[1000005];intgcd(intx,inty){if(y==0)returnx;elsereturngcd(y,x%y);}longlongexgcd

https://ac.nowcoder.com/acm/contest/4853/D1#include2#include3usingnamespacestd;4usingll=longlong;5voidexgcd

文章目录数论排序高精度字符串数据结构组合数学多项式数论gcdCRTBSGSexgcdexCRTexBSGS整除分块线性筛欧拉函数线性筛莫比乌斯函数特殊数列的前缀最小公倍数排序堆排序归并排序快速排序高精度高精度加法高精度减法高精度乘法字符串

若ax+by=1,gcd(a,b)=1；则a的逆元为x%b；#include#includeusingnamespacestd;intl,m,n,s,p;voidexgcd(inta,intb,int&

拓扑排序O(n^3^)树的直径归并排序O(nlogn)LCS记录路径O(n*m)Dijkstra求单源最短路KMPtarjan+缩点网络流求最小割即最大流EKdinictarjan求lca扩展欧几里得llexgcd

有了前置结论，再用上扩展欧几里得求乘法逆元就ok代码#include#definelllonglongusingnamespacestd;constintMAXN=1e5+100;vectorvec;llexgcd

a:gcd(b,a%b);}扩展欧几里得核心代码llexgcd(lla,llb,ll&x,ll&y){if(b==0){x=1;y=0;retur

https://blog.csdn.net/sdutstudent/article/details/78795643附上一个超nice的exgcd的讲解，这里就是记录一下学习过程中卡住的一些地方我们根据某些追击问题的经验可以得出

EXGCD这个版块其实没有专门去打题，感觉对于题目来说就是要抽象出对应的ax+by=c这样的方程吧！

**~~EXGCD知识点总结**~~学任意模数NTT发现忘了一堆数论知识一下追溯到这里，复习知识点用令d=gcd（a,b）x最小整数姐x=(x%(b/d)+(b/d))%(b/d)求a对模b的逆元时，要求

看了T2，发觉好像没什么想法，可能用exgcd判一判吧。又看T3，完了，又没什么想法，莫非只会T1？！还好我心态比较稳健，很快打完T1过了大样例之后就去打T2。判了倍数

适合单个的，费马小定理，exgcd，都是不错的选择，利用积性函数的方法和欧拉筛的方法适合批量求，但是论时间和空间的话，还是积性函数的方法比较好用，线性的。

Code1Code1\\求解单个数在模意义下的逆元（exgcd）vo

负数取模还是负数，但范围+取模的数就是正数代码如下#include#include#include#includeusingnamespacestd;intexgcd(longlonga,longlongb

本来直接暴力枚举就可以搞完的事，当场脑子抽了想写EXGCD，浪费时间不说，还WA了一次（没看到均为非负数）。

假设方程为：aX+bY=c;首先判断是否有整数解：c%gcd(a,b)==0则有解代码：#include"iostream"#includeusingnamespacestd;voidexgcd(inta

>n)输出他们的最大公约数，同时输出s,t(m*s+n*t=1)#-*-coding:utf-8-*-#扩展欧几里得算法#输入mn#输出mn的最大公约数还有s,t##默认m>nimportsysdefexgcd

二元一次不定方程的整数解（扩展欧几里得算法）（不得不说这是一堂数学*信竞课）整数解解法c(modb)或ax+by=c有整数解当且仅当（a,b）|c一般意义下的解法：欧拉函数扩展欧几里得算法代码实现exgcd

：D1T1：最短路+并查集/kruskal重构树（今年复赛会不会考一道数据结构+图论综合在一起的一道题呢）D1T2：卡特兰数+DP（初赛也出现了卡特兰数感觉卡特兰数出现的概率有点大呀接下来）D2T1：exgcd

include#include#include#definelonglongLL;#defineINF0x3f3f3f3f#definePIacos(-1.0)usingnamespacestd;voidexgcd

给我们ab问能不能构造出前面两个分数，看到这个首先联想到exgcd拓欧至于构造，不妨利用b的约数在此我取的是b的最小质因子的全部幂次的数作为一个分母，另一个数直接用b除就好了，然后问题转为求解fx+dy

F.FractionConstructionProblem(exgcd&数论)题意：给定正整数a,ba,ba,b，构造出正整数c,d,e,fc,d,e,fc,d,e,f，满足cd−ef=ab,(d,f1gcd

F.F.F.题意：给定a,ba,ba,b，求出满足cd−ef=ab\frac{c}{d}-\frac{e}{f}=\frac{a}{b}dc−fe=ba的任意一组c,d,e,fc,d,e,fc,d,e,f。数据范围：a,b≤2×106a,b\leq2\times10^6a,b≤2×106，1≤d1g=gcd(a,b)>1g=gcd(a,b)>1，令b′=bg,a′=agb'=\frac{b}{g}

=B;variable+=C)statement;已知所有的数均要mod2^k(即0usingnamespacestd;typedeflonglongLL;LLexgcd(LLa,LLb,LL&x,LL

a:gcd(b,a%b);}LLexgcd(LLa,LLb,LL&x,LL&y){if(b==0){x=1;y=0

根据RSA一堆原理，写了个实现的代码，亲测可用，哈哈记录一下#include#include#includeusingnamespacestd;intexgcd(inta,intb,int&x,int&

2intexgcd(inta,intb,int&x,int&y){inti,j;if(b==0){x=1;y=0;ret