扩展欧几里得算法，解模线性方程，解ax+by=c的解集

typedef long long LL;
LL exgcd(LL a,LL b,LL &x,LL &y){
     if(a==0&&b==0) return -1;
     if(b==0) { x=1;y=0; return a;  }
     LL d=exgcd(b,a%b,y,x);
     y-=a/b*x;
     return d;
}
LL MLE(LL a,LL b,LL mod){  // 返回(a*x)%mod=b 的最小正整数解
    LL x,y;
    LL gcd=exgcd(a,mod,x,y);
    if(b%gcd!=0) return -1;  //若返回-1，则该方程无解
    x*=b/gcd;
    mod/=gcd;
    if(mod<0) mod=-mod;
    LL ans=x%mod;
    if(ans<=0) ans+=mod;  //其他解为： ans+i*mod （i为整数）
    return ans;
}
LL x0,y0,kx,ky;
bool LE(LL a,LL b,LL c){   //解线性方程ax+by=c
    LL x1,y1;
    LL gcd=exgcd(a,b,x1,y1);
    if(c%gcd)return false;  //无整数解
    x0=x1*c/gcd,y0=y1*c/gcd;
    kx=b/gcd,ky=-a/gcd;
    return true;   //有解，解集为：（x0+kx*t,y0+ky*t） t为整数
}

用扩展欧几里德算法求解模线性方程的方法&#x