解模线性方程ax=b(modn)

#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long LL;
// 求x和y使得ax+by=d并且|x|+|y|最小。其中d=gcd(a,b)
void exgcd(LL a,LL b,LL& d,LL& x,LL& y){
    if(!b) d = a,x = 1,y = 0;
    else{
        exgcd(b,a % b,d,y,x);
        y -= x * (a / b);
    }
}
// 用扩展欧几里得解模线性方程 ax=b(mod n)
bool modular_Linear_Equation(LL a,LL b,LL n){
    LL x,y,x0,i,d;
    exgcd(a,n,d,x,y);
    if(b % d) return false; //d不是b的约数时没有解
        x0 = x * (b / d) % n; // 得到一个特解
    for(i = 1 ; i <= d ; i++)
       printf("%d\n",(x0 + i * (n / d)) % n); // 恰好有d个解
    return true;
}
int main(){
    modular_Linear_Equation(3,6,6);
    return 0;
}