整数与整除

整数与整除

扩展欧几里得算法

• 裴蜀定理：若 $a,b\in Z$，对于任意整数 $x,y$，$ax+by=k\gcd (a,b)$。特别的，存在$x,y$，使得 $ax+by=\gcd (a,b)$

  用来判断式子是否有解：$ax+by=m$ 有解 $⟺\gcd (a,b)|m$

  • 证明：

    设 $d=\gcd (a,b)$，则 $d|(ax+by)$

    不妨设 $s$ 为 $ax+by$ 的最小正值，则 $\because d|s\therefore d\le s$

    令 $q=a/s,r=a\%s$，即 $a=qs+r$.

    $\therefore r=a-q(ax+by)=a(1-qx)+b(-qy)$，即 $r$ 为 $a,b$ 的线性组合。

    $\because 0\le r<s且s为最小正值\therefore r=0,a=qs$，即 $s|a$.

    同理 $s|b$.

    故 $s$ 为 $a,b$ 的公约数，于是 $\because s|d\therefore s\le d$.

    于是 $s=d$ 得证。

  • 推论：$a,b$ 互质 $⟺$ 存在整数 $x,y$，使得 $ax+by=1$

• 扩展欧几里得：

  • 目的：计算 $ax+by=\gcd (a,b)$ 的解 $x,y$

  • 方法：根据普通辗转相除法，到达递归边界 b=0,a=gcd(a,b) 时，发现有解 $x=1,y=0$，即 $a\ast 1+b\ast 0=\gcd (a,b)$。

  • 注意：此时 $a,b$ 已不是开始的 $a,b$，于是根据递归算法和下一状态 $b,a\%b$ 的通解 $x_1,y_1$（即 $b\ast x_1+(a\%b)\ast y_1=\gcd$）求出当前状态的通解。
    $$\begin{gathered} \because a\%b=a-(a/b)\ast b \\ \therefore b\ast x_1+(a\%b)\ast y_1=b\ast x_1+(a-(a/b)\ast b)\ast y_1=a\ast y_1+b\ast (x_1-a/b\ast y_1)=\gcd \\ 即\{\begin{array}{l}x=y_1\\ y=x_1-a/b\ast y_1\end{array} \end{gathered}$$

    int exgcd(int a, int b, int &x, int &y){
        if(!b){
            x = 1, y = 0;
            return a;
        }
        int gcd = exgcd(b, a % b, x, y);
        int tmp = y;
        y = x - (a / b) * x;
        x = tmp;
        return gcd;
    }
    

例题：P2613 【模板】有理数取余

• 题目：计算有理数 $\frac{a}{b}(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的最小正数解

• 思路：即 $a\equiv bx(\mathrm{m}\mathrm{o}\mathrm{d}p)$，也就是 $bx\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的解的 $a$ 倍。处理大数据时只需要用 $a\%p,b\%p$ 代替 $a,b$ 即可。

  //计算 b * x (mod mod) == a 的最小正数解.即 b * x (mod mod) == 1 的解的 a 倍。
  #include<bits/stdc++.h>
  using namespace std;
  #define rep(i, a, b) for(int i = (a); i <= (b); i++)
  #define ll long long
  const int mod = 19260817;
  ll a = 0, b = 0, x = 0, y = 0;
  string A, B;
  ll exgcd(ll a, ll b, ll &x, ll &y){
      if(b == 0){
          x = 1, y = 0;
          return a;
      }
      ll gcd = exgcd(b, a % b, x, y);
      ll tmp = y;
      y = x - (a / b) * y;
      x = tmp;
      return gcd;
  }
  int main(){
      cin >> A >> B;
      rep(i, 0, A.length() - 1){
          (a *= 10) += (A[i] - '0');
          a %= mod;
      }
      rep(i, 0, B.length() - 1){
          (b *= 10) += (B[i] - '0');
          b %= mod;
      }
      if(b == 0 && a != 0){
          puts("Angry!");
          return 0;
      }
      exgcd(b, mod, x, y);
      x *= a;
      x = (x % mod + mod) % mod;
      cout << x << endl;
  }
  

例题：P1516 青蛙的约会

• 题目：计算 $xx+mx(\mathrm{m}\mathrm{o}\mathrm{d}L)=yy+nx(\mathrm{m}\mathrm{o}\mathrm{d}L)$ 的最小正数解 $x$

• 思路：
  $$\begin{array}{rl}(n-m)x& \equiv xx-yy(\mathrm{m}\mathrm{o}\mathrm{d}L)\\ ax& \equiv c(\mathrm{m}\mathrm{o}\mathrm{d}b)\\ ax+by& =c\end{array}$$

  • 有解 $⟺\gcd (a,b)|c$，其中 $a=n-m,b=L,c=xx-yy$

  • 若有解，可用扩展欧几里得求出满足 $a{x}_0+b{y}_0=\gcd (a,b)$ 的解。然后满足 $ax+by=c$ 的解即为 $x_1=\frac{c}{\gcd (a,b)}{x}_0$，则任意解为 $x=x_1+k\frac{b}{\gcd (a,b)}$

    要求最小正数的特解，只需要 $x=(x\%r+r)\%r,其中r=\frac{b}{\gcd (a,b)}$

  #include<bits/stdc++.h>
  using namespace std;
  #define rep(i, a, b) for(int i = (a); i <= (b); i++)
  #define per(i, a, b) for(int i = (a); i >= (b); i--)
  #define ll long long
  ll x, y, m, n, L, xx, yy, a, c, g;
  ll exgcd(ll a, ll b, ll &x, ll &y){  //求出 ax + by = gcd(a,b) 的特解
      if(b == 0){
          x = 1, y = 0;
          return a;
      }
      ll gcd = exgcd(b, a % b, x, y);
      ll tmp = y;
      y = x - (a / b) * y;
      x = tmp;
      return gcd;
  }
  int main(){
      cin >> xx >> yy >> m >> n >> L;
      m %= L, n %= L, xx %= L, yy %= L;
      a = n - m, c = xx - yy;
      if(a < 0) a *= -1, c *= -1;
      g = exgcd(a, L, x, y);
      if(c % g != 0){
          puts("Impossible");
          return 0;
      }
      x *= c / g;
      x = (x % (L / g) + (L / g)) % (L / g);
      cout << x << endl;
  }
  

逆元

• 单独求某个数逆元：扩展欧几里得 $O(\log n)$

• 统一求 $1\sim n$ 中每个数关于 $p$ 的逆元：线性求逆元 $O(n)$，单独求每个数花费约 $O(n^{\frac{1}{3}})$

  • 代码：

    inv[1] = 1;
    for (int i = 2; i <= n; ++i) {
      inv[i] = (long long)(p - p / i) * inv[p % i] % p;
    }
    

  • 证明：

    首先 $1^{-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$. 要求 $i^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

    令 $k=p/i,r=p\%i$，有 $p=ki+r$，则 $ki+r\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$

    两边同乘 $i^{-1}\times r^{-1}$：$i^{-1}\equiv -k{r}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$，即 $i^{-1}\equiv -(p/i)\ast (p\%i{)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$