数论算法学习笔记（一）：素数

整除

素数在整除的定义上衍生而来。于是就得先看一下整除的定义：
对于 $a,b\in \mathbb{Z},a\ne 0$，若 $\exists r\in \mathbb{Z}$，使得 $b=ar$，就可以说 $b$ 可以被 $a$ 整除，记作 $a\mid b$，否则，则称 $a$ 不能整除 $b$，记作 $a\nmid b$。
那么若 $a\mid b$，则称 $a$ 是 $b$ 的约数，$b$ 是 $a$ 的倍数。特殊地，$0$ 是任何非 $0$ 整数的倍数。
显而易见的是，对于 $b\ne 0$，$b$ 只有有限个约数。
对于 $b\ne 0$，我们称 $\pm 1,\pm b$ 为 $b$ 的平凡约数。特殊地，当 $b=\pm 1$ 时，$b$ 只有 $2$ 个平凡约数。$b$ 的其余约数被称为 $b$ 的真约数。
在具体问题中，如果没有特殊说明，我们所指的约数总是指正约数。
若对于 $a,b\in \mathbb{Z}$，$a$ 和 $b$ 的最大公约数 $\gcd (a,b)=1$，则称 $a,b$ 互素。

同余

素数往往涉及到同余，所以在写素数之前，先要写一点同余。
先给出余数的定义：
设 $a,b\in \mathbb{Z},a\ne 0$，设 $d\in \mathbb{Z}$，则一定存在唯一的一对 $q$ 和 $r$，满足 $b=qa+r,d\le r<|a|+d$，那么 $q$ 就是商，$r$ 就是余数。在一般情况下，$d=0$，则此时的 $r$ 被称为最小非负余数。如果没有特殊说明，余数总是指最小非负余数。
同余的定义：
设 $m\in \mathbb{Z}$，若 $m\mid (a-b)$，则可以记作 $a\equiv b(\mathrm{mod}m)$，这样的等式就被叫做同余式。其中，$m$ 为模（数），$a$ 同余于 $b$ 模 $m$，$b$ 是 $a$ 对模 $m$ 的剩余。如果没有特殊说明，模数总是正整数。

素数

设整数 $p\ne 0,\pm 1$，若 $p$ 除了平凡约数之外没有其他的约数，则称 $p$ 是一个素数。否则，$p$ 是一个合数。
$p$ 和 $-p$ 总是同为素数或同为合数。但若无特殊说明，素数总是指正的素数。
若一个整数的约数是一个素数，那这个素数被称作这个整数的素约数（素因数）。

算数基本定理（唯一分解定理）

对于整数 $a$，在不计次序的情况下，有唯一等式成立：
$$a=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}\cdots p_n^{{\alpha }_n}$$
其中 $p_1,p_2,\cdots ,p_n$ 是素数。这个等式被称为整数 $a$ 的标准素因数分解式。

素数的判定

知道了什么是素数，那么我们就要知道如何判断这个数是不是素数。
常用试除法判断：

bool isprime(int a)
{
	if (a<2)
	{
		return false;
	}
	for (int i=2;i*i<=n;i++)
	{
		if (a%i==0)
		{
			return false;
		}
	}
	return true;

也可以使用概率算法 Fermat 和 Miller-Rabin：

int qpow(int a,int b,int p)
{
	int ans=1;
	while (b>0)
	{
		if (b&1)
		{
			ans=ans*a%p;
		}
		a=a*a%p;
		b>>=1;
	}
	return ans;
}
bool fermat(int n)
{
	if (n<3)
	{
		return n==2;
	}
	for (int i=1;i<=ttime;i++) //ttime>8, below is the same.
	{
		int a=rand()%(n-2)+2;
		if (qpow(a,n-1,n)!=1)
		{
			return false;
		}
	}
	return true;
}
bool millerrabin(int n)
{
	if (n<3||n%2==0)
	{
		return n==2;
	}
	int u=n-1,t=0;
	while (n%2==0)
	{
		u/=2;
		t++;
	}
	for (int i=1;i<=ttime;i++)
	{
		int a=rand()%(n-2)+2,v=qpow(a,u,n);
		if (v==1)
		{
			continue;
		}
		int flg=0;
		for (flg=0;s<t;s++)
		{
			if (v==n-1)
			{
				break;
			}
			v=(long long)v*v%n;
		}
		if (flg==t)
		{
			return false;
		}
	}
	return true;
}

筛法

既然我们知道了素数，那么我们就有必要知道怎么找到素数。筛法可以帮助我们找到自然数内小于等于 $n$ 的范围内的素数个数。

埃氏筛

考虑到任何素数的倍数（除了自己以外）都是合数，那么我们就可以标记出所有的合数，那么剩下的就是素数。

#define MAXN 
int isprime[MAXN];
int Enes(int n)
{
	int ans=0;
	for (int i=0;i<=n;i++)
	{
		isprime[i]=1;
	}
	isprime[0]=isprime[1]=0;
	for (int i=2;i<=n;i++)
	{
		if (isprime[i])
		{
			prime[ans++]=i;
			if ((long long)i*i<=n)
			{
				for (int j=i*i;;j<=n;j+=i)
				{
					isprime[j]=0;
				}
			}
		}
	}
	return ans;
}

可以证明，埃氏筛的时间复杂度为 $O(n\log \log n)$。

线性筛（欧拉筛）

如果只让每个合数被标记一次，那么时间复杂度就可以降到 $O(n)$ 了。

#define MAXN
bool vis[MAXN];
int pri[MAXN],cnt;
void euler(int n)
{
	for (int i=2;i<=n;i++)
	{
		if (!vis[i])
		{	
			pri[cnt++]=i;
		}
		for (int j=0;j<cnt;j++)
		{
			if (1ll*i*pri[i]>n)
			{
				break;
			}
			vis[i*pri[i]]=true;
			if (i%pri[i]==0)
			{
				break;
			}
		}
	}
}

筛法的应用

筛法有很多应用。

筛法求约数个数

按规矩，我们使用 $d[i]$ 表示 $i$ 的约数个数，$num[i]$ 表示 $i$ 的最小素约数出现的次数。于是就有：

约数个数定理
若
$$n=\prod _{i=1}^m{p}_i^{c_i}$$
则
$$d[i]=\prod _{i=1}^m(c_i+1)$$

代码：

void quantity()
{
	d[1]=1;
	for (int i=2;i<=n;i++)
	{
		if (!vis[i])
		{
			pri[++cnt]=i;
			d[i]=2;
			num[i]=1;
		}
		for (int j=1;j<=cnt&&i<=n/pri[j];j++)
		{
			vis[i*pri[j]]=true;
			if (i%pri[j]==0)
			{
				num[i*pri[j]]=num[i]+1;
				d[i*pri[j]]=d[i]/num[i*pri[j]]*(num[i*pri[j]]+1);
				break;
			}
			else
			{
				num[i*pri[j]]=1;
				d[i*pri[j]]=2*d[i];
			}
		}
	}
}			

筛法求约数和

令 $f[i]$ 是 $i$ 的约数和，$g[i]$ 表示最小素约数的 $p^0+p^1+\cdots +p^k$。

代码：

void sum()
{
	f[1]=g[1]=1;
	for (int i=2;i<=n;i++)
	{
		if (!vis[i])
		{
			vis[i]=true;
			pri[++cnt]=i;
			f[i]=i+1;
			g[i]=i+1;
		}
		for (int j=1;j<=cnt&&i<=n/pri[j];j++)
		{
			vis[i*pri[j]]=true;
			if (i%pri[j]==0)
			{
				g[i*pri[j]]=g[i]*pri[j]+1;
				f[i*pri[j]]=f[i]/g[i]*g[i*pri[j]];
				break;
			}
			else
			{
				f[i*pri[j]]=f[i]*pri[j];
				g[i*pri[j]]=pri[j]+1;
			}
		}
	}
}