exgcd(扩展欧几里德算法)

题面描述

解不定方程Ax+By=K（得到的x和y只是其中一组解）
给出A、B、K，求出x和y，满足Ax+By=K。

【输入格式】
一行三个整数 A，B，K。 1 ≤ A，B，K ≤ 10^9

【输出格式】
一行两个整数 x,y。如果无解，输出"no solution!"

思路

$在欧几里得算法中，当b=0时，显然存在一对整数解x=1,y=0,使得a\ast 1+0\ast 0=gcd(a,0).$

$令ax+by=gcd(a,b)$
$\because gcd(b,amodb)=gcd(a,b);$
$\therefore a\ast x+b\ast y=gcd(a,b)=gcd(b,amodb)=b\ast tx+(amodb)\ast ty$
$\because amodb=a-\lfloor a/b\rfloor \ast b$
$\therefore a\ast x+b\ast y=b\ast tx+(a-\lfloor a/b\rfloor \ast b)\ast ty$
$\therefore a\ast x+b\ast y=a\ast ty+b\ast (tx-\lfloor a/b\rfloor \ast ty)$
$\therefore x=ty,y=tx-\lfloor a/b\rfloor \ast ty$
再用数学归纳法进行证明即可。
由于本题$K$不一定等同于$gcd(a,b)$，因此我们要先求到$gcd(a,b)$,后逆推得$x,y$，判断一下$K$是否被$gcd(a,b)$整除，再$x=x\ast k/d$，$y=(k-a\ast x)/b$

代码

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define ll long long 
using namespace std;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	if(!b)
	{
		x=1;y=0;
		return a;
	}
	else
	{
		ll tx,ty;
		ll d=exgcd(b,a%b,tx,ty);
		x=ty;y=tx-(a/b)*ty;
		return d;
	}
}
int main()
{
	ll a,b,k,x,y;scanf("%lld%lld%lld",&a,&b,&k);
	ll d=exgcd(a,b,x,y);
	if(k%d){puts("no solution!");return 0;}
	else
	{
		x=x*k/d;
		y=(k-a*x)/b;
	}
	printf("%lld %lld\n",x,y);
	return 0;
}