线性同余方程初步应用分析

一元线性同余方程

定义：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)的方程。对于模线性方程组，可以进行方程组合并，求出合并后的方程的解，这样就可以很快推出方程的最终解，不管这样的方程有