题目来源:点击打开链接
判断一个线性方程组Ax+by+C=0是否有整数解,如果有,则输出整数解。
这个问题用到了欧几里得扩展求逆元的算法。
用这个算法可以求 X mod b=m 这种模线性方程(求的是X的通解,整数解)。变形一下,得 X+by=m,其中y为整数。通常写成 ax+by=m 的形式,然后解得x的通解,进而得到X=ax的通解。
其实就是在求GCD(a,b)的同时把解出满足 ax+by=m 的(x,y)通解。
然后判断是否是整数,这个好说,用最大公约数的性质就可以了。
#include <cmath>
#include <iostream>
using namespace std;
__int64 x,y;
__int64 gcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
/*if(a==0)
return b;*///这一句在标准算法中可以加上
if(b==0)
{
x=1;
y=0;
return a;
}
int result=gcd(b,a%b,x,y);
int temp=x-a/b*y;
x=y;
y=temp;
return result;
}
int main()
{
__int64 a,b,c;
cin>>a>>b>>c;
c=-c;
__int64 divgcd=gcd(a,b,x,y);
// cout<<divgcd<<" "<<x<<" "<<y<<" "<<endl;
if(c%divgcd!=0)
{
cout<<-1<<endl;
}
else
{
if(divgcd<0)
{
x=-x;
y=-y;
divgcd=-divgcd;
}
x=x*(c/divgcd);
y=y*(c/divgcd);
cout<<x<<" "<<y<<endl;
}
return 0;
}