初等数论——素数

基础

概念

整除：如果整数$b$可以表示成整数$a$的整数$k$倍，即$b=ka(a,b,c\in Z)$,则称$b$被$a$整除或$a$整除$b$，记作$a|b$

最大公约数：若$a,b\in N^{\ast }$ ,则$(a,b)$表示最大的正整数$d$，使得$d|a\wedge d|b$,成为$a,b$的最大公约数。特别的，规定任意正整数$a$,有$(a,0)=a$

素数

概念：除了$1$与其本身外，无其它约数的正整数。特别地，1既不是素数也不是合数。

判定

判定正整数$n(n>1)$是否为素数

• 枚举法

  从$2$枚举至$n-1$, 若其中存在数能整除$n$，则$n$为合数，否则$n$为素数。复杂度$\mathrm{O}(n)$

  小优化：若$d|n$,则必有$\frac{n}{d}|n$,因此$n$的约数成对存在，且乘积为$n$。任选一对约数$d1\cdot d2=n$,若$d1>\sqrt{n}\wedge d2>\sqrt{n}$，则$d1\cdot d2>n$与$d1\cdot d2=n$相矛盾，因此$n$的任意一对约数$d1\le d2$,必有$d1\le \sqrt{n}$. 又由于约数成对存在，因此我们仅需枚举$2$至$\sqrt{n}$判断是否存在$n$的约数即可。复杂度$O(\sqrt{n})$

• 素数测试

  费马素数测试

  费马小定理：对于任意素数$p$,若$p\nmid a$,则有$a^{p-1}\equiv 1(modp)$

容易想到可以通过选取随机数检验费马小定理来判定$n$是否是素数。但费马小定理逆定理并不成立，即使$a^{p-1}\equiv 1(modp)$,$p$也未必是素数，我们将满足该条件的合数称为伪素数或卡迈尔数

Miller-Rabin 素数测试

二次探测定理：若$p$为奇素数，则$x^2\equiv 1(modp)$解为$x\equiv 1$或$x\equiv p-1$.移项因式分解易证。

费马小定理结合二次探测定理，对于$n-1$,我们先分解成$n-1=n_0\cdot 2^t$,然后先求出$c=p^{n_0}$,然后对$c$进行$t$次自平方操作，若出现非平凡平方根(即$1$或$n-1$之外的数)，则不是素数，否则通过测试

参考代码：

inline bool Miller_Rabin(ll x) {
	if (x<3 || x%2==0) return (x==2);
	ll a=x-1, b=0;
	while (a%2==0) a>>=1, ++b;
	for (R ll i=1; i<=10; i++) {
		ll c=rand()%(x-1)+1;
		c=quick_pow(c, a, x);
		if (c==1) continue;
		ll j;
		for (j=0; j<b; j++) {
			if (c==x-1) break;
			c=(el)c*c%x;
		}
		if (j>=b) return false;
	}
	return true;
}

筛法

• Eratosthenes筛（埃筛）

  目的是筛出$2-n$内所有素数。素数即除$1$与其本身外没有其它约数的正整数。因此我们可以通过约数否决掉一个数是素数的可能性。每遍历到一个素数，枚举它的整数倍，将其倍数全部排除.复杂度$\mathrm{O}(nloglogn)$

  int prime[N], cnt;
  bool book[N];
  inline void get_prime(int n) {
      for (int i=2; i<=n; i++) {
          if (book[i]) continue;
          prime[++cnt]=i;
          for (int j=i; j<=n; j+=i) book[j]=true;
      }    
  }
  

• 欧拉筛

如果每个合数都只标记一次，那复杂度就是线性的了。详见$OI-wiki$代码解释.复杂度$\mathrm{O}(n)$

```cpp
inline void euler(){
    for(R ll i=2;i<=n;i++){
        if(!v[i]){
            v[i]=i;
            prime[++sum]=i;
        }
        for(R ll j=1;j<=sum&&prime[j]<=v[i]&&prime[j]*i<=n;j++){
            v[prime[j]*i]=prime[j];
            if (i%prime[j]==0) break;
        }
      }
}