约数

算数基本定理的推论

N的正约数集合：
{$p_1^{b1}$$p_2^{b2}$…$p_m^{bm}$}，其中0<=bi<=ci
N的正约术个数为：
$$(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 ...\ast (1+p_m+p_m^2+...p_m^{c_m})=\prod _{i=1}^m(\sum _{j=0}^{c_i}(p_i{)}^j)$$

求N的正约数集合

试除法

int factor[1600], m = 0;
void func(int n) {
	for (int i = 1; i * i <= n; i++) {
		if (n % i == 0) {
			factor[++m] = i;
			if (i != n / i)factor[++m] = n / i;
		}
	}
	for (int i = 1; i <= m; i++)
		cout << factor[i] << endl;
}

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

求1~N每个数的正约数集合

倍数法

vector<int>factor[500010];
void func(int n) {
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n / i; j++)
            factor[i * j].push_back(i);
    for (int i = 1; i <= n; i++) {
        for (int j = 0; j < factor[i].size(); j++)
            cout << factor[i][j] << " ";
        cout << endl;
    }
}

例题

AcWing198
AcWing199

最大公约数

欧几里得算法

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

互质与欧拉函数

互质的定义：

gcd(a,b)=1,a与b互质

欧拉函数

$$\phi (N)=N\ast \frac{p_1-1}{p_1}\ast \frac{p_2-1}{p_2}\ast ...\ast \frac{p_m-1}{p_m}=N\ast \prod _{质数p|N}(1-\frac{1}{p})$$

int phi(int n) {
    int ans = n;
    for (int i = 2; i <= sqrt(n); i++) {
        if (n % i == 0) {
            ans = ans / i * (i - 1);
            while (n % i == 0)n /= i;
        }
    }
    if (n > 1)ans = ans / n * (n - 1);
    return ans;
}

性质1~2

1.$\forall$n>1,1~n中与n互质的数的和为n*$\phi (n)$/2。
2.若a,b互质，则$$\phi (ab)=\phi (a)\phi (b)$$

积性函数

如果当a,b互质时，有f(a,b)=f(a)*f(b),那么f为积性函数

性质3~6

3.若f是积性函数，且在算术基本定理中 $$f(n)=\prod _{i=1}^m{p}_i^{c_i}$$
4.设p为质数，若p|n且$p^2|n$,则$$\phi (n)=\phi (n/p)\ast p$$
5.设p为质数，若p|n且$p^2\nmid n$,则$$\phi (n)=\phi (n/p)\ast (p-1)$$
6.$$\sum _{d|n}\phi (d)=n$$

例题：

AcWing201
可见的点

void euler(int n) {
    for (int i = 2; i <= n; i++) phi[i] = i;
    for (int i = 2; i <= n; i++)
        if (phi[i] == i)
            for (int j = i; j <= n; j += i)
                phi[i] = phi[i] / i * (i - 1);
}