洛谷P1516 青蛙的约会（数论同余、扩欧）

之前遇到可以套用的模板直接就拿来用了，然鹅半年多了，水平还是这样差。所以争取对那些重要的只是都理解下来。今天看了一下午，研究这个题，终于还是理解的比较透彻了，虽然时间花的有点多，但是收获仿佛更大。

题意：emmmmmmmmmmmmmm，很简单求最小整数k使(n-m)*kx-y (mod l)；//还要满足n,m,k,x,y,l都是整数

分析：貌似这种同余都可以化成解方程的形式。=>(n-m)*k+z*l=x-y(zZ);令a*x+b*y=c,即a=n-m,x=k,b=l,y=z,c=x-y。

扩欧算法：

//扩展欧几里得算法 
int exgcd(int a，int b,int &x,int &y){
	if(b==0){
		x=1;y=0;
		return a;//到达递归边界开始向上一层返回 
	}
	int g=exgcd(b,a%b,x,y);
	int temp=y;//吧x,y变成上一层的
	y=(x-(a/b)*y);
	x=temp;
	return g;//得到a b的最大公因数 
}

解释扩欧算法原理：

一直递推用a=b,b=a%b直到求出a*x+b*y=gcd(a,b)的解x0=1,y0=0;然后再推回去（怎么退回去呢？就是上面的递归调用之后的操作，一会解释），得到x1的一个特解。然后可以根据这个特解求出一个x1。不过这个x为a*x+b*y=gcd(a,b)的解，所以a*x+b*y=c的特解实际上是x1*c/gcd(a,b)。求最小值的操作，emmmm，还有点点小不懂。

解释怎么用a*1+y*0=gcd(a,b)递推回去得到x1:

加入已经是状态（b,a%b）了，即b*x+(a%b)*y

又a%b=a-(a/b)*b(往回退，a,b增大)代入得到

=b*x+(a-(a/b)*b)*y

=a*y+b*(x-(a/b)*y)

......

一直这样递推回去（次数等于“正向”递推次数）（按照x=y,y=x-(a/b)*y）然后就可以得到x1

举个栗子：40*x+14*y

(a,b)=(40,14)=(14,12)=(12,2)=(2,0)得到g=gcd(a,b)=2。然后往回：(x,y)=(1,0)=(0,1)=(1,-1)=(-1,3)

可以验证，均满足a*b+b*y=gcd(a,b);如1*14+（-1）*12=2

如何用a*x+b*y=c特解x1得到最小值（非负数）:

if(c%g!=0) puts("Impossible");//c/g并不是整数，说明得不到整数 
else printf("%lld",(x1*c/g%(b/g)+b/g)%(b/g);//取最小非负整数 

第二行取最小非负的常规操作倒是挺好理解，比如mod=5

x1为负数（-7%5+5）

x1为正数  7%5

综合:（x1%mod+mod）%mod;//正负通吃

但是为什么mod=b/g呢？？？

如何由一个解得到其它解呢？有一个恒等式a(x+db)+b(y-da)=c

在保证db,da都是整数的情况下，我们让d最小，就可以得到所有的整数解，那么d=1/g​

如果解出的x>0，那么最小非负整数解等于x%(b/g)​；否则等于x%(b/g)+b/g;​

代码就可以直接写(x%(b/g)+b/g)%(b/g)

//其实这里还是有一点点不理解，先记住吧。

so,代码：

#include<bits/stdc++.h>
#define ll long long
#define pi acos(-1)
#define mst(a,i) memset(a,i,sizeof(a))
using namespace std;
template<class t>
inline void read(register t &x){register char c;while(c=getchar(),c>57||c<48);x=c^48;while(c=getchar(),c>47&&c<58)x=10*x+(c^48);}
const int inf=0x3f3f3f3f;
const int maxn=2e5+5;
const int mod=1e9+7;

ll x,y,n,m,L,a,xx,b,yy,c,g,tmp;

void exgcd(ll a,ll b){
	if(b==0){
		xx=1,yy=0;
		g=a;
		return ;
	}
	exgcd(b,a%b);
	ll tmp=yy;
	yy=xx-(a/b)*yy;
	xx=tmp;
	return ;//总要return ?可以不么 
}
int main(){
    ios::sync_with_stdio(false);
    int T=1;
    //read(T);
    while(T--){
    	read(x),read(y),read(m),read(n),read(L);
    	a=n-m,b=L,c=x-y;
    	if(a<0) a=-a,c=-c;
		//一定要有这一步，正数求gcd。两边同时取负，为啥b不变？因为b已经是正的了，解出来y的符号不一样罢了，无影响 
    	exgcd(a,b);
    	if(c%g) puts("Impossible");
    	else printf("%lld",(c/g*xx%(b/g)+b/g)%(b/g));
    }
    return 0;
}

终于还是AC了。不知道为神魔，竟然不能定义y0,y1(最后改为yy才正确)。