求解单个同余方程的剩余系

//ax ≡ b (mod n)
#include <stdio.h>

void gcd(int a,int b,int& d,int& x,int& y){
    if(!b){
        d = a; x = 1; y = 0;
    }else{
        gcd(b, a % b, d, y, x);
        y -= x * (a / b);
    }
}

void linear_mod_equation(int a,int b,int n,int& d,int sol[]){
    // d = gcd(a,n)
    int x,y;
    gcd(a,n,d,x,y);
    if(b % d) d = 0;
    else{
        sol[0] = x * (b / d) % n;
        for(int i = 1; i < d; i++)
            sol[i] = (sol[i - 1] + n / d) % n;
    }
}

int sol[1000];
int main(void){
    int a,n,b,d;
    scanf("%d%d%d",&a,&n,&b);
    linear_mod_equation(a,b,n,d,sol);
    for(int i = 0; i < d; i++)
        printf("%d ",sol[i]);
    printf("
");
    return  0;
}