快乐地打牢基础（6）——约数

约数定义： 若整数$n$除以整数$d$的余数为0，即$d$能整除$n$，则称$d$是$n$的约数，$n$是$d$的倍数，记为 $d|n$ 。

一、整除的性质

---

设$a,b$是两个正整数，且$b!=0$，则存在唯一的整数$q$和$r$，使得:
$$a=qb+r，（0\le r<b）$$
这个式子叫做带余除法，并记余数$r=amodb$。
整除具有如下性质：

• 若$a|b$且$a|c$，则$\forall x,y,$有$a|(xb+yc)$。
• 若$a|b$且$b|c$，则$a|c$
• 设$m!=0$，则当且仅当$ma|mb$时，$a|b$
• 若$a|b$且$b|a$，$a=\pm b$

二、约数

---

1.算术基本定理（唯一分解定理）的推论
在算术基本定理中，若正整数$N$被唯一分解为$N=p_1^{c_1}{p}_2^{c_2}...p_m^{c_m}$，其中$c_i$都是正整数，$p_i$都是质数，且满足$p_1<p_2<...<p_m$，则$N$的正约数集合可写作：
{ p 1 b 1 p 2 b 2 . . . p m b m } ， 其 中 0 ≤ b i ≤ c i \{p_1^{b_1}p_2^{b_2}...p_m^{b_m}\}，其中0\leq b_i\leq c_i { p1b1​​p2b2​​...pmbm​​}，其中0≤bi​≤ci​
$N$的正约数个数为（$\prod$表示连乘）
$$(c_1+1)\ast (c_2+1)\ast ...\ast (c_m+1)=\prod _{i=1}^m(c_i+1)$$
$N$的所有正约数的和为
$$(1+p_1+p_1^2+...+p_1^{c_1})\ast (1+p_2+p_2^2+...+p_2^{c_2})\ast ...\ast (1+p_m+p_m^2+...+p_m^{c_m})=\prod _{i=1}^m(\sum _{j=0}^{c_i}(p_i{)}^j)$$

//使用算术基本定理求质因子代码
//prime数组是之前用筛法筛出的质数数组。
void divide(ll n)
{
   
    int len = prime.size();
    int cnt = 0;
    for(int i = 0; i < len; i++)
    {
   
        if(n % prime[i] == 0)
        {
   
            p[++cnt] = prime[i];
            c[cnt] = 0;
            while(n % prime[i] == 0 && n > 1)
            {
   
                n /= prime[i];
                c[cnt]++;
            }
        }
        printf("p = %d c = %d\n",p[cnt],c[cnt]);
        if(n == 1)
            break;
    }
    //如果遍历到了1e6的素数，n还没有变为1，那么剩下的n一定是一个素数。（前提，n<1e12）
    if(n > 1){
   
        p[++cnt] = n,c[cnt] = 1;
    }
}

2.求N的正约数集合——试除法

若$d\ge \sqrt{N}$是$N$的约数，则$N/d\le \sqrt{N}$也是$N$的约数。换言之，约数总是成对出现的（除了对于完全平方数，$\sqrt{N}$会单独出现）。
因此，只需要扫描$d=1～\sqrt{N}$，尝试$d$能否整除$N$，若能整除，则$N/d$也是$N$的约数。时间复杂度$O(\sqrt{N})$。

//试除法函数
int f[1600];
void divide(int n){
   
	m = 0;
	for(int i = 2; i * i <= n; i++){
   
		if(n % i == 0){
   
			f[++m] = i;
			if(i * i != n) f[++m] = n/i;
		}
	}	
	
}

试除法的推论
一个整数$N$的约数个数的上界为$2\sqrt{N}$。

3.求1~N每个数的正约数集合——倍数法

若用"试除法"分别求出$1～N$每个数的正约数集合，时间复杂度过高，为$O(N\sqrt{N})$。可以反过来考虑，对于每个数$d$，$1～N$中以$d$为约数的数就是$d$的倍数$d,2d,3d,...\lfloor N/d\rfloor \ast d$。以下程序采用"倍数法"求出$1～N$每个数的正约数集合：

//求1~N每个数的正约数集合——倍数法
vector<int> f2[500010];
void multiple(int n) {
   
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= n/i; j++)
			f2[i*j].push_back(i);
}

倍数法的推论
$1～N$每个数的约数个数的总和大约为$NlogN$

【例题1】BZOJ 1053 反素数

题意：
如果一个大于等于 1 的正整数 n，满足所有小于 n 且大于等于 1 的所有正整数的约数个数都小于 n 的约数个数，则 n 是一个反素数。譬如：1,2,4,6,12,24，它们都是反素数。
现给定一个数$N(1\le n\le 2\ast 10^9)$，请你计算不大于 N的最大反素数。
思路：
设$g(x)$为$x$的约数个数。

• 引理1： $1～N$中最大的反素数就是，其中约数个数最多的数中最小的一个。

证明： 设$m$是$1～N$中约数个数最多的数中最小的一个。根据$m$的定义，$m$显然满足： Ⅰ.$\forall x<m,g(x)<g(m)$ Ⅱ. ∀ x > m , g ( x ) ≤ g ( m ) \forall x > m,g(x)\leq g(m)