一元线性同余方程
定义:a,b是整数,m是正整数,形如ax≡b(mod m),且x是未知整数的同余式称为一元线性同余方程。
定理:a,b,m是整数且m>0,gcd(a,m)=d,如果d|b,则方程恰有d个模m不同余的解,否则方程无解。
由同余方程式的定义可知:ax+my=b(y为整数),这个方程称为二元一次不定方程。
(1)解一元线性同余方程
设d=gcd(a,m),由定理可知若不满足d|b,那么方程无解,否则有:a=d×ao,m=d×mo;
那么方程变为:ao*x+mo*y=b/d;
由于此时gcd(ao,mo)=1,因此可以运用扩展欧几里得算法得出方程ao*x+mo*y=b/d的解x,虽然x不唯一,但是属于一个模m剩余系,由定理可知,共有d个模m剩余类满足方程,其代表元分别为:
x,x+mo,x+2*mo,...x+(d-1)*mo;
(2)一元线性同余方程求解代码实现
//扩展欧几里得算法
typedef long long LL;
void exgcd(LL a,LL b,LL &d,LL &x,LL &y)
{
if(!b)
{
x=1,y=0,d=a;
return ;
}
else
{
exgcd(b,a%b,d,x,y);
LL temp=x;
x=y;
y=temp-(a/b)*y;
}
}
//求方程的所有解
LL ans[maxn];
int solve(int a,int b,int m)
{
exgcd(a,m,d,x,y)
if(b%d) return -1; //代表无解
x=x*(b/d)%m;
for(int i=1;i<=d;i++)
ans[i]=(x+(i-1)*m/d)%m;
}
解一元线性同余方程组
由上一讲我们知道,任何一元同余方程都可以变成若干个形如x≡b(mod m)的方程。对于模线性方程组,可以进行方程组合并,求出合并后的方程的解,这样就可以很快推出方程的最终解,不管这样的方程有