数论学习笔记

$\text{Part 1 gcd,lcm}$

$\gcd (x,y)$表示$x,y$的最大公约数，代码实现：

int gcd(int x,int y){return y?gcd(y,x%y):x;}

从代码中可以发现，每次$x,y$中有一个数至少减了一半，所以时间复杂度为$\Theta (\log \min (x,y))$

$\text{lcm}(x,y)$表示$x,y$的最小公倍数
$$\text{lcm}(x,y)=\frac{x\times y}{\gcd (x,y)}$$
即
$$\text{lcm}(x,y)\times \gcd (x,y)=x\times y$$

$\text{Part 2 exgcd}$

$(1)$贝祖定理：对于给定的正整数$a,b$，关于$x,y$的方程$ax+by=c$有整数解的充要条件是$c$是$\gcd (a,b)$的整数倍。

$(2)\text{exgcd}$正是贝祖定理中方程的解。

$(3)$我们发现每次递归的时候$x=y1,y=x1-\lfloor a/b\rfloor \times y1$，而当递归的$y=0$的时候易得该方程的一组解是$x=0,y=1$，然后逐层回溯得到其它解。

代码实现：

void exgcd(int a,int b)
{
	if (b==0){x=1,y=0;return;}
	exgcd(b,a%b);
	int t=x;
	x=y;
	y=t-(a/b)*y;
}

注意，此处的$x,y$是全局变量。

$\text{Part 3 }$莫比乌斯函数、反演

$(1)$莫比乌斯函数表示为$\mu (i)$，函数值：
$$\mu (n)=\{\begin{array}{r}\text{1 }(n=1)\\ (-1{)}^k\forall v_p(n)\le 2\\ 0\exists v_p(n)>2\end{array}$$
$(2)$性质：
$$\sum _{d|n}\mu (d)=[n=1]$$
$(3)$筛前缀和：欧拉筛$($时间复杂度$\Theta (n))$、杜教筛$($时间复杂度$\Theta (n^{\frac{2}{3}}))$

代码$($欧拉筛$)$：

for (int i=2;i<=n;++i)
{
	if (!bz[i]) prime[++tot]=i,mu[i]=-1;
	for (int j=1;j<=tot && i*prime[j]<=n;++j)
	{
		bz[i*prime[j]]=1;
		if (i%prime[j]==0{mu[i*prime[j]]=0;break;}
		mu[i*prime[j]]=-mu[i];	
	}
}
for (int i=1;i<=n;++i) sum[i]=sum[i-1]+mu[i];

代码$($杜教筛$)$：

int S(int x)
{
	if (x<N) return mu[x];
	if (sum[x]) return sum[x];
	int o=1;
	for (int i=2,lt;i<=x;i=lt+1) lt=x/(x/i),o-=(lt-i+1)*S(x/i);
	sum[x]=o;
	return o;
}

$(4)$莫比乌斯反演定理：若$F(n),f(n)$为算数函数且满足
$$F(n)=\sum _{d|n}f(n)$$
则有
$$f(n)=\sum _{d|n}\mu (d)F(\frac{n}{d})=\sum _{d|n}\mu (\frac{n}{d})F(d)$$
$(5)$较常见的莫比乌斯反演

$$\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)=k]$$
$$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }[\gcd (i,j)=1]$$
$$\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }\sum _{d|\gcd (i,j)}\mu (d)$$
$$\sum _{d=1}^{\min }\mu (d)\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }1$$
$$\sum _{d=1}^{\min }\mu (d)\lfloor \frac{n}{k}\rfloor \lfloor \frac{m}{k}\rfloor$$
$\mu (d)$的前缀和可以欧拉筛、杜教筛随便乱搞，$\lfloor \frac{n}{k}\rfloor \lfloor \frac{m}{k}\rfloor$随便数论分块一下就行了，时间复杂度$\Theta (n+\sqrt{n})$或$\Theta (n^{\frac{2}{3}}+\sqrt{n})$

$(6)$用$\mu$的前缀和求$\phi$的前缀和：
$$\sum _{i=1}^n\phi (i)$$
$$\sum _{i=1}^n\sum _{j=1}^{i-1}[\gcd (i,j)=1]$$
$$\frac{{\sum }_{i=1}^n{\sum }_{j=1}^n[\gcd (i,j)=1]}{2}$$
为了方便先单独挑出分子进行分解
$$\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)=1]$$
$$\sum _{i=1}^n\sum _{j=1}^n\sum _{d|\gcd (i,j)}\mu (d)$$
$$\sum _{d=1}^n\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{d}\rfloor }1$$
$$\sum _{d=1}^n\mu (d){\lfloor \frac{n}{d}\rfloor }^2$$
所以
$$\sum _{i=1}^n\phi (i)=\frac{{\sum }_{d=1}^n\mu (d){\lfloor \frac{n}{d}\rfloor }^2}{2}$$$\mu (d)$用欧拉筛、杜教筛照搞，后面照样数论分块，时间复杂度$\Theta (n+\sqrt{n})$或$\Theta (n^{\frac{2}{3}}+\sqrt{n})$

剩余待更$......$