gcd

有序列表

1. LCM(a, b) = a / GCD(a, b) × b;//先除后乘 防止数字过大爆炸

2. 欧几里得cpp

void gcd(int  a,int  b)
{
    return b==0?a:gcd(b,a%b);
}

3. 扩展欧几里得

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

4. 线性模方程的解

 void mod_line(int a, int b, int n)
    {    
    //ax = b (mod n)
        int d, x, y, x0, i;
        exgcd(a, n, x, y,d);
        if(b%d == 0)
        {
            x0 = x * (b/d) % n + n;//x0为任意一个解
            for(i = 0;i < d;i ++)
               cout << (x0 + i*n/d) % n <<endl;//方程有d个解
        }
     }

一些结论（一下解都是基于整数）：

若：ax = b mod n;
转化后 ax-kn=b;
设：$d=gcd(a,n)$
上述方程的整数解可写成：
设一个解为$(x_0,y_0)$
$x=x_0+i*n/d，y=y_0+i*a/d$

1.求任意解

设：$s=\frac{b}{d}$
则任意解为：$(x_0*s$ mod $n+n)$ mod $n$

2.求通解

首先，明确一点通解共有d个
求通解前要先求得一个任意解：设：$x_0=(x_0*s$ mod $n+n)$ mod $n$
则通解为：$(x0 + i*n/d)$ mod n $0$≤i<$d$

3.求最小正整数解

设：$m=\frac{n}{d}$
同样要先求一个任意解$x_0$
则最小正整数解为：$x_0=(x_0*s$ mod $m+m)$ mod m m <script type="math/tex" id="MathJax-Element-1297">m</script>

解释一下，因为每个解都是以n/d间隔递增求模的，那么显然正解在（0，n/d）之间，所以mod n/d