素数的相关问题

线性素数筛法

void judge(int n)
{
	int m=sqrt(n+0.5);
	for(int i=2;i<=m;i++)
	if(vis[i])
	{
		for(int j=i*i;j<=n;j+=i)
		vis[j]=1;
	}
}

上述代码看起来并不是线性的，但事实上它避免了重复筛查，从i*i开始筛查，因为2*i，3*i.....(i-1)*i之前已经筛查过了，所以上述算法基本可以在O（n）的时间复杂度内筛出素数。

唯一分解定理

for(int i=2;i*i<=A;)
{
	if(A%i==0)
	{
		p[++k]=i;
		n[k]=0;
		while(!(A%i))
		{
			n[k]++;
			A/=i;
		}
	}
	if(i==2) i++;
	else
	i+=2;
	if(A!=1)
	{
		p[++k]=A;//记录素数
		n[k]++;//记录次数
	}
}

威尔逊定理

(p-1)!=-1(mod p）

（p-1)!+1一定为p的倍数，再利用sin（（pi*（n-1）！+1）/n）

函数值为0的点就是素数所在的点

费马小定理

若p为素数，a为正整数，且a和p互质，则a^（p-1）=1(mod p)

首先，a,2*a,3*a,....(p-1)*a中没有一个是p的倍数

其次，a,2*a,3*a,.....(p-1)*a中没有一个是同余p

a*2*a*3*a*.....*(p-1)*a=1*2*3*4*...*(p-1)mod p;

a^p-1*(p-1)!=(p-1)!(mod p)

a^p=a(mod p);

Miller-Rabin测试

a^n%n=a;

int modular_exp(int a,int m,int n)//快速幂取模
{
	if(m==0) return 1;
	if(m==1) return a%n;
	long long w=modular_exp(a,m/2,n);
	w=w*w%n;
	if(m&1) w=w*a%n;
	return w;
}
bool Miller_Rabin(int n)//费马小定理
{
	if(n==2) return true;
	for(int i=0;i<5;i++)
	{
		int a=rand()%(n-2)+2;
		if(modular_exp(a,n,n)!=a) return false;
	}
	return true;
}
int main()
{
	srand(time(NULL));
	int n;
	scanf("%d",&n);
	if(Miller_Rabin(n)) printf("YES");
	else cout<<"NO";
	return 0;
}

欧拉函数的求法

void get_phi(int maxn)
{
	phi[1]=1;
	for(int i=2;i<=maxn;i++)
	{
		if(!vis[i])//质数 
		p[++p[0]]=i,phi[i]=i-1;//p用来存质数 
		for(int j=1;i*p[j]<=maxn&&j<=p[0];j++)
		{
			vis[i*p[j]]=1;
			if(i%p[j]==0)
			{
			phi[i*p[j]]=phi[i]*p[j];
		    break;//避免重复，提高效率，i应为质数 
			}
			else
			phi[i*p[j]]=phi[i]*(p[j]-1);
		}
	}
}

欧拉定理

a^phi[m]=1(mod m);

欧拉函数的线性素数筛法

void getphi()
{
	phi[1]=1;
	int i,j;
	for(int i=2;i<=n;i++)
	{
		if(!vis[i])
		{
			phi[i]=i-1;	
		    prime[++prime[0]]=i;
		}
		for(int j=1;j<=tot;j++)
	{
		if(i*prime[j]>N) break;
		mark[i*prime[j]]=1;
		if(i%prime[j]==0)
		{
		phi[i*prime[j]]=phi[i]*prime[j];
        break; 
		} 
		else
		phi[i*prime[j]]=phi[i]*(prime[j]-1);
	}
	}
}