质数的应用

1.试除法求质数 

bool is_prime(int x)
{
	if(x<2)
	return false;
	for(int i=2;i<=x/i;i++)
	if(x%i==0)
	return false;
	return true;
}

2.分解质因子

对于每个正整数 x，按照从小到大的顺序输出其分解质因数后，每个质因数的底数和指数，每个底数和指数占一行。 

void divide(int x)
{
	for(int i=2;i<=x/i;i++)
	{
		if(x%i==0)
		{
			int s=0;
			while(x%i==0)
			{
				x/=i,s++;
			}
			cout<<i<<' '<<s<<endl;
		}
	}
	if(x>1)
	cout<<x<<' '<<1<<endl;
	cout<<endl;
}

 最后的x>1的判断，是把没有除尽的数输出，这个数一定是质数。如21，在i=3时筛掉3，剩余7，跳出循环，但7也是21的质因子，所以也需要输出。

3.筛质数

（1）朴素筛法

//prime数组存储1~n中的质数
void get_prime(int n)
{
    for(int i=2;i<=n;i++)
	{
		if(!st[i])
		{
			prime[cnt++]=i;
		}
		for(int j=i+i;j<=n;j=j+i)//倍数一定不是质数
		{
			st[j]=true;
		}
    }	
} 

 （2）线性筛法

void get_prime(int n)
{
	for(int i=2;i<=n;i++)
	{
	    if(!st[i])
	    primes[cnt++]=i;
	    for(int j=0;primes[j]<=n/i;j++)
	    {
	    	st[primes[j]*i]=true;
	    	if(i%prime[j]==0)
	    	break;
		}
	}
} 

核心为：一个数只会且一定会被自身最小质因子筛掉

情况1：prime[j]是i的最小质因子，则prime是i*prime[j]的最小质因子

情况2：prime[j]不是i的最小质因子，则一定比i的最小质因子小，所以筛掉i*prime[j]也是正确的