数论

一、对于大整数$N$来说，不大于$N$的质数大约有$\frac{N}{\ln N}个$
二、 判断质数的方法

1. 试除法 复杂度$O(\sqrt{N})$ 绝对正确

2. Miller Robbin 非常高效，但有将合数判为质数的可能
   通过费马小定理，猜想当满足$a^{m-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$时$m$是质数（虽然这样不对，网上能查到不少反例）
   再辅以二次探测定理：

   $$若p是质数，x<p且x^2\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)则x=1或x=p-1$$

   进行验证。
   具体：取几个质数作为$a$，从$m-1$中除出尽可能多的2（）假设结果为$d$，个数为$s$算出$a^d$次方进行$s$平方，若某一刻模$p$的余数为$1$且平方前不为$1$或$m-1$则$m$不是质数
   代码

   long long a[8]={
        2,3,5,7,11,13,17,61};
   long long mul(long long x,long long y,long long p)
   {
        
   	long long ret=0;
   	x%=p;
   	y%=p;
   	while(y)
   	{
        
   		if(y&1)(ret+=x)%=p;
   		y>>=1;
   		(x*=2)%=p;
   	}
   	return ret%p;
   }
   
   long long qpow(long long x,long long y,long long p)
   {
        
   	long long ret=1;
   	while(y)
   	{
        
   		if(y&1)ret=mul(ret,x,p);
   		y>>=1;
   		x=mul(x,x,p);
   	}
   	return ret%p;
   }
   bool isprime(long long x)//miller robbin
   {
        
   	if(x==2)return true;
   	if(!(x&1))return false;
   	if(x<=1)return false;
   	long long d=x-1;
   	long long s=0;
   	while(!(d&1))
   	{
        
   		++s;
   		d>>=1;
   	}
   
   	for(int i=0;i<8;++i)
   	{
        
   		long long k=qpow(a[i],d,x);
   		for(int j=1;j<=s;++j)
   		{
        
   			long long y=mul(k,k,x);
   			if(y==1&&k!=1&&k!=x-1)return false;
   			k=y;
   		}
   	}
   	return true;
   }
   

三.质数筛法

1. Eratosthenes筛法 简单基础 复杂度大约$O(N{\log }^{2}N)$

   bool vis[100005];
   void Eratosthenes(int n)
   {
        
   	for(int i=2;i<=n;++i)
   	{
        
   		if(vis[i])continue;
   		printf("%d\n",i)l
   		for(int j=i;j<=n/i;++j)vis[i*j]=true;
   	}
   }	 
   

2. 线性筛法 复杂度$O(N)$

   int vis[100005];//v[i]:the smallest prime factor of i
   int cnt,prim[10005];
   void primes(int n)
   {
        
   	cnt=0;//the total number of prime numbers 
   	for(int i=2;i<=n;++i)
   	{
        
   		if(vis[i]==0)vis[i]=i,prim[++cnt]=i;
   		for(int j=1;j<=cnt;++j)
   		{
        
   			if(prim[j]>vis[i]||prim[j]>n/i)break;
   			vis[i*prim[j]]=prim[j];
   		}
   	}
   	for(int i=1;i<=cnt;++i)printf("%d\n",prim[i]);
   } 
   

   因为每个数只会被它的最小质因子筛一次 故是线性的

四、质因数分解

1. 试除法
2. Pollard’s Rho
   (To be continued…)

   ---

五、约数

1. 基本定理的推论

   若正整数N被唯一分解为$N=p_1^{c_1}{p}_2^{c_2}\dots p_m^{c_m}$，其中$c_i$都是整数，$p_i$都是质数，且满足$p_1<p_2<\cdots <p_m$ 那么
   $N$的正约数集合可写作 { p 1 b 1 p 2 b 2 … p m b m ， 0 ≤ b i ≤ c i } \{p_1^{b_1}p_2^{b_2}\dots p_m^{b_m}，0\leq b_i\leq c_i\} { p1b1​​p2b