扩展欧几里得算法

（1）扩展欧几里得算法
裴蜀定理：对于任意正整数a，b，一定存在非零整数x，y，使得ax+by=gcd（a，b）

/*
给定 n 对正整数 ai,bi，对于每对数，求出一组 xi,yi，使其满足 ai ×xi + bi ×yi = gcd(ai,bi)。
*/
#include<bits/stdc++.h>
using namespace std;
const int N = 1e5+10;
int exgcd(int a,int b,int &x,int &y)
{
	if(!b)
	{
		x=1,y=0;
		return a;
	}
	int d=exgcd(b,a%b,y,x);
	y=y-a/b*x;
	return d;
}
int main()
{
	int n;
	scanf("%d",&n);
	while(n--)
	{
		int a,b,x,y;
		scanf("%d%d",&a,&b);
		exgcd(a,b,x,y);
		printf("%d %d\n",x,y);
	}
	return 0;
} 

（2）线性同余方程
a × x ≡ b (mod m) 等价于a除以m的模为b，也等价于存在整数y，使得ax=my+b，即ax-my=b，令y’=-y，则变为ax+my’=b（变为了扩展欧几里得的问题，x有解的前提的b是gcd（a，m）的倍数）

/*
给定 n 组数据 ai,bi,mi，对于每组数求出一个 xi，使其满足 ai ×xi ≡ bi (mod mi)，如果无解则输出 impossible。
*/
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
int exgcd(int a,int b,int &x,int &y)
{
	if(!b)
	{
		x=1,y=0;
		return a;
	}
	int d=exgcd(b,a%b,y,x);
	y=y-a/b*x;
	return d;
}
int main()
{
	int n;
	scanf("%d",&n);
	while(n--)
	{
		int a,b,m,x,y;
		scanf("%d%d%d",&a,&b,&m);
		int d=exgcd(a,m,x,y);
		if(b%d) printf("impossible\n");//如果b不是gcd（a，m）的倍数，则无解 
		else printf("%lld\n",(LL)x*(b/d)%m);//x要扩大（b/d）倍再模m 
	}
	return 0;
}