【8.5】代码源 - 【GCD】【序列中位数】【最大连边数量】

#871. GCD

题意：给定 $1\le a<m\le 10^{10}$ ，计算 ${\sum }_{x=0}^{m-1}[\gcd (a,m)=\gcd (a+x,m)]$ 。

题解：(区间互质) 代码源每日一题 Div1 GCD

思路：推一下式子：

容易知道 $k=\gcd (a,m)$ 一定是 $m$ 的因子。我们先计算出 $a+x$ 的范围内是 $k$ 的倍数的左右边界，倍数范围即为 $[l,r]=[\lceil \frac{a}{k}\rceil ,\lfloor \frac{a+m-1}{k}\rfloor ]$ 。那么 ${\sum }_{x=0}^{m-1}[\gcd (a,m)=\gcd (a+x,m)]={\sum }_{i=l}^r[\gcd (i,\frac{m}{k})=1]$ ，那么问题转化为求 $l\sim r$ 中有多少个数与 $\frac{m}{k}$ 互质。

一个经典的问题就是求 $1\sim n$ 有多少个数与 $m$ 互质，方法为：对 $m$ 分解质因子，然后二进制枚举容斥计算，时间复杂度 $O(\sqrt{n}+2^{11})$ ，这里的 $11$ 是指 $10^{10}$ 内的数字最多有 $11$ 个质因子。

总结了一些关于 $\gcd$ 的性质：

$\gcd (a,b)=\gcd (a,a+b)=\gcd (a,a-b)$ ，两个数的和或差的 $\gcd$ 相等
$\gcd (k\times a,k\times b)=k\times \gcd (a,b)$
$\gcd (a,b)\Rightarrow p\times \gcd (\frac{a}{p},\frac{b}{p})$ ，要保证 $p$ 是 $a,b$ 的因子
$\gcd (a^k,b^k)={\gcd }^k(a,b)$
$c|a\times b\Rightarrow \frac{c}{\gcd (b,c)}|a$

模板：

LL cal(LL n, LL m){ // numbers of i in 1 ~ n : gcd(i, m) == 1
    vector<LL> vec;
    for(LL i = 2; i * i <= m; i ++ ){
        if(m % i == 0){
            while(m % i == 0) m /= i;
            vec.pb(i);
        }
    }
    if(m > 1) vec.pb(m);
    LL res = 0;
    for(int bit = 1; bit < 1 << vec.size(); bit ++ ){
        LL mul = 1;
        rep(i, 0, (int)vec.size() - 1) if(bit >> i & 1) mul *= vec[i];
        if(__builtin_popcount(bit) & 1) res += n / mul;
        else res -= n / mul;
    }
    return n - res;
}

AC代码：http://oj.daimayuan.top/submission/321288

---

#877. 序列中位数

题意：给定正整数 $N$ , 求出 $1\sim N-1$ 中所有与 $N$ 互质的数构成的序列 的中位数.

我们定义 : 一个长度为 $K$ 的序列的中位数是序列中第 $\lfloor \frac{K+1}{2}\rfloor$ 大的数字. 且两个正整数 $a$ 与 $b$ 互质当且仅当 $\gcd (a,b)=1$ 。

思路：打表可以发现规律。证明的话： $\gcd (i,x)=\gcd (x-i,x)$ ，大概就是这个式子。

AC代码：http://oj.daimayuan.top/submission/322050

---

#916. 最大连边数量

题意：$1\le n\le 10^5$

题解：(树状数组/DP) 代码源每日一题 最大连边数量

思路：树状数组优化 DP 。

我们按照下标枚举 $a$ 里的元素，去枚举 $b$ 里对应的元素。我们先定义 $dp(i,j)$ 表示 $a$ 前 $i$ 个元素和 $b$ 前 $j$ 个元素的最大匹配，转移方程为 $dp(i,j)={\max }_{1\le i^{\prime }<i,1\le j^{\prime }<j}dp(i^{\prime },j^{\prime })+1$ 。

可以发现这是个二维偏序，第一维能被优化。对于和当前枚举的 $i$ 对应的 $j$ ，转移方程为 $dp(j)=\max (d{p}^{\prime }(j),{\max }_{k=1}^{j-1}d{p}^{\prime }(k)+1)$ ，用树状数组维护第二维的前缀最值即可即可。

AC代码：http://oj.daimayuan.top/submission/322797