数论初步

$$数论初步$$

GDKOI 冬令营，算是正式开始面对数论吧。（我太菜了）

质数

1.1-N的质数个数（用于估计数组下标）：$$\pi (x)=\frac{x}{lnx}$$
2.试除法判断质数：复杂度$O(\sqrt{n})$
3 埃斯托尼克筛法：复杂度 $O(NloglogN)$
4.欧拉筛/线性筛：复杂度$O(N)$
附代码：

int prime[1000*1001],top;
bool pd[1000*10000+100];
void OL_find()
{
	for(register int i=2;i<=n;i++)
	{
		if(!pd[i]) prime[++top]=i;
		for(register int j=1;j<=top&&prime[j]*i<=n;j++)
		{
			pd[i*prime[j]]=true;
			if(i%prime[j]==0) break;//关键！！！
		}
	}
}

约数

1.x能被y整除：$y|x$
对于一个数N 有$$N=p_1^{c_1}\cdot p_2^{c_2}\cdot ......\cdot {p_m}^{c_m}$$
2.1~N的正约数个数：
$$(c_1+1)\ast (c_2+1)\ast ......\ast (c_m+1)=\prod _{i=1}^m(c_i+1)$$

3.1~N的正约数和：
$$(1+p_1+p_1^2+......p_1^{c_1})\ast ......\ast (1+p_m+p_m^2+\cdot \cdot \cdot \cdot \cdot \cdot p_m^{c_m})=\prod _{i=1}^m(\sum _{j=0}^{c_i}(p_i{)}^j)$$

最大公约数，最小公倍数

a,b的最大公约数： 记为 $gcd(a,b)$
a,b的最小公倍数： 记为 $lcm(a,b)$
性质：
$gcd(a,b)=gcd(b,a-b)=gcd(a,a-b)九章算术\cdot 更相减损术$
$gcd(a,b)\ast lcm(a,b)=a\ast b$
$gcd(m\cdot a,m\cdot b)=m\ast gcd(a,b)$
$gcd(a,b)=gcd(b,amodb)欧几里得算法——辗转相除法$

重头戏：欧拉函数

定义：$\phi (N)$表示1~N之间与N互质数的个数
$$\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 f(int n)
{
	int ans=n;
	for(int i=2;i<=sqrt(n);++i)
		if(n%i==0)
		{
			ans=ans*(i-1)/i;
			while(n%i==0) n/=i;
		}
	if(n>1) ans=ans*(n-1)/n;
	return ans;
}

此代码在对N试除时顺便求出了$\phi (n)$

·················································································································
未完待续
··················································································································