试除法分解质因数

试除法分解质因数

实现：

        {\\ \ \ \ \ \ \ \ }        枚举2~$\sqrt{\text{N}}$中的质数$\text{d}$，如果$\text{N}$能被$\text{d}$整除，则从$\text{N}$中将所有因子$\text{d}$除掉，并记录$\text{d}$的数量。

证明:只需枚举2~$\sqrt{\text{N}}$:

        {\\ \ \ \ \ \ \ \ }        $\text{N}$大于$\sqrt{\text{N}}$的质因子至多有一个，反证若有两个或以上则他们之积首先应是$\text{N}$的因子，但是却超过$\text{N}$。或者说，达不到分解的目的(因子之积等于$\text{N}$)。

证明:不需要筛素数:

        {\\ \ \ \ \ \ \ \ }        对于$\text{N}$的任意一个合数因子$\text{M}$，在枚举到它之前就已经将它从$\text{N}$中除去，因为$\text{M}$的质因子必是$\text{N}$的质因子，而这些质因子早已被枚举并除掉了。

代码:

#include<cstdio>
#include<cmath>
int p[100100],c[100100],len=0;
//p[]存入质因子,c[i]存入第i个质因子的个数 
void div(int n)
{
    for(int i=2;i<=sqrt(n);i++)/* sqrt(n)会在线更新， 
	按照证明,数a ( sqrt(n') < a <= sqrt(n) )不需要试除 */
	{
		if(n%i==0)
		{
			p[++len]=i;
			while(n%i==0) n/=i,c[len]++;
		}
	}
	if(n>1) p[++len]=n,c[len]++;
}
int main()
{
	int n;
	scanf("%d",&n);
	div(n); 
	printf("%d = %d^%d",n,p[1],c[1]);
	for(int i=2;i<=len;i++)
			printf(" * %d^%d",p[i],c[i]);
	return 0;
}

        {\\ \ \ \ \ \ \ \ }        分解质因数还有更高级的算法：Miller-Rabin和Pollard-Rho算法，笔者用时间会学的。